Necessary to filter out nanopore (or PacBio) reads before mapping?
1
1
Entering edit mode
2.3 years ago

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 • 808 views
ADD COMMENT
3
Entering edit mode
2.3 years ago

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.

ADD COMMENT
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 1237 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