Question: filtering of mapped Nanopore reads
gravatar for m.eitel
6 months ago by
m.eitel0 wrote:


I am trying to filter a bam file. I mapped raw nanopore reads to a reference assembly and want to keep only reads that map (almost) full length. bwa mem and minimap2 softmask bases at both ends that do not match to the reference and I want to keep only those reads that map full-length and allowing for 40bp unmatched/softmasked bp on both ends. Tried Samclip, but it does not work (outputting text file instead of true sam file). Any ideas on how to filter otherwise?

Thanks Michael

sequencing next-gen assembly • 216 views
ADD COMMENTlink written 6 months ago by m.eitel0

please, post a SAM example with a few reads + header on

ADD REPLYlink written 6 months ago by Pierre Lindenbaum116k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1353 users visited in the last hour