Greetings, thank you for the tool. I have been trying to compare results to other sweep detectors and have run into an issue with my reference genome. The reference fasta has several gap regions in each chromosome that are represented by "NNNNNNNNNNN" sequences that range from 1000-1,000,000 bp. These are the regions that RASID identifies as standing out with high muVar. I have tried to exclude the regions using the -X command, which seems to help reduce the signal but it still identifies these regions. As an example, the highest hit in the output is below, despite the exclude file listing "chr5 58097111 59097111" as a region.
58599927 58096412 59103441 366.60 1.89 2.00 0.00
Is there a better way to handle this issue? Thank you
Greetings, thank you for the tool. I have been trying to compare results to other sweep detectors and have run into an issue with my reference genome. The reference fasta has several gap regions in each chromosome that are represented by "NNNNNNNNNNN" sequences that range from 1000-1,000,000 bp. These are the regions that RASID identifies as standing out with high muVar. I have tried to exclude the regions using the -X command, which seems to help reduce the signal but it still identifies these regions. As an example, the highest hit in the output is below, despite the exclude file listing "chr5 58097111 59097111" as a region.
58599927 58096412 59103441 366.60 1.89 2.00 0.00
Is there a better way to handle this issue? Thank you