filtering of mapped Nanopore reads
Hi.
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
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Minimap2 segfaults when mapping Nanopore reads to a very large reference database
written by firefox91 •
-
PacBio amplicon reads partially aligned using minimap2 – library or analysis issue?
written by louisflower1999 •
-
Restrictive alignment with minimap2
written by njornet •
-
filter reads with >90% match
written by anuradharavi10 •
-
demultiplex paired end fastq file with mixed samples with barcode
written by hellocita •
-
BWA mem minimum read length
written by yolek64754 •
-
How to properly assemble the virus genome using map to reference strategy
written by xatabadich •
-
Mapping to custom reference minimap2/IGV
written by sendhelp •
-
Necessary to filter out nanopore (or PacBio) reads before mapping?
written by BioinformaticsLad
-
3 peaks in per base coverage plot
written by m.eitel •
please, post a SAM example with a few reads + header on gist.github.com