Question: Necessary to filter out nanopore (or PacBio) reads before mapping?
gravatar for BioinformaticsLad
4 months ago by
BioinformaticsLad140 wrote:

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?

ADD COMMENTlink modified 4 months ago by Biostar ♦♦ 20 • written 4 months ago by BioinformaticsLad140
gravatar for lieven.sterck
4 months ago by
VIB, Ghent, Belgium
lieven.sterck6.1k wrote:

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 COMMENTlink written 4 months ago by lieven.sterck6.1k

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

ADD REPLYlink written 4 months ago by h.mon27k
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: 2547 users visited in the last hour