Entering edit mode
12 weeks ago
njornet
▴
20
I need to align nanopore reads from human to a custom reference genome and I want to align only the reads that have a very high identity. I tried changing the missmatch penalty but if I increase it too much minimap raises a Segmentation fault error. Should I change something else in the configuration of minimap2 or should I filter the aligned reads by score or something and how would I do that?
Thank you!
Filtering after the fact based on MAPQ or CIGAR's an option?
If I have the resulting SAM it is right?