Question: Necessary to filter out nanopore (or PacBio) reads before mapping?
1
gravatar for BioinformaticsLad
14 months ago by
BioinformaticsLad160 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 14 months ago by Biostar ♦♦ 20 • written 14 months ago by BioinformaticsLad160
3
gravatar for lieven.sterck
14 months ago by
lieven.sterck8.3k
VIB, Ghent, Belgium
lieven.sterck8.3k 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 14 months ago by lieven.sterck8.3k
1

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

ADD REPLYlink written 14 months ago by h.mon30k
Please log in to add an answer.

Help
Access

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