If you have all regions smaller than 78 (1260 >> kPSS), with coverage grather or equal ot 4, you don't have any valid "pile". When you try to calculate the median, a segmentation fault is raised. This solution avoids this, but I think it is better to inform the user about the lack of overlap because none contigs will be generated.
If you have all regions smaller than 78 (1260 >> kPSS), with coverage grather or equal ot 4, you don't have any valid "pile". When you try to calculate the median, a segmentation fault is raised. This solution avoids this, but I think it is better to inform the user about the lack of overlap because none contigs will be generated.
Datasets: ERR4082025 and ERR4082026