This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Necessary to filter out nanopore (or PacBio) reads before mapping?

I have some nanopore reads I want to map to hg38 human genome using minimap2. In the past I would filter out mean QV<7 reads but lately I've been thinking: if the reads are long enough, there's little chance that (despite errors) they map to the wrong regions. Is it really necessary to filter these out? If the coverage is high enough, the consensus should cover up low quality bases. Anyone have any opinions?

nanopore minimap2 alignment mapping
• 2,386 views • link
updated by Biostar Biostar • written by BioinformaticsLad BioinformaticsLad

1 answer

I think there is indeed a growing consensus that you should just map them and if they map OK then use them. In order words, let the mapping do the 'filtering' rather than on basecall quality scores.

• 0 views • link
written by lieven.sterck lieven.sterck

In addition, PacBio Sequel don't use quality scores, and all raw reads have a Phred score of zero.

• 0 views • link
written by h.mon h.mon

Log in to answer this question.

More posts like this