Restrictive alignment with minimap2
0
0
Entering edit mode
3 months 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!

nanopore alignment minimap2 • 314 views
ADD COMMENT
0
Entering edit mode

I want to align only the reads that have a very high identity

Filtering after the fact based on MAPQ or CIGAR's an option?

ADD REPLY
0
Entering edit mode

If I have the resulting SAM it is right?

ADD REPLY

Login before adding your answer.

Traffic: 1992 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6