Reads length Distribution of ONT reads and Error rate
2
0
Entering edit mode
4.6 years ago
midox ▴ 270

Hello,

I want to know the distribution of the reads and the error rate (mismatch, insertion, deletion) of the Nanopore Sequencing Technology.

Thank you

nanopore NGS error rate reads Distribution • 2.0k views
0
Entering edit mode
4.6 years ago

BBMap can report those statistics, though reads longer than 6kbp will be broken into 6kbp segments. For this purpose, that shouldn't matter, though it's worth noting that the shorter the fragments, the lower the apparent error rate is.

mapPacBio.sh in=reads.fastq maxlen=6000 out=mapped.sam ref=reference.fasta


That will report the rates of mismatches, insertions, and deletions.

0
Entering edit mode

in the case of PacBio reads or i can also use Nanopore?

0
Entering edit mode

The error rates of PacBio and Nanopore are similar (both extremely high) so I use the same error profile. Though in practice, Nanopore error rates seem to be much higher than PacBio, so to map as many reads as possible, you can reduce the kmer length and increase sensitivity. But... nobody really cares about reads with under 75% identity, so I'd just classify those as junk, and ignore them.

0
Entering edit mode

I am looking for a distribution that already exists of number of nanopore reads according to their lengths and the different error rates. If I can do it with BBmap I try to map the raw reads on the reference genome and extract this informations.

0
Entering edit mode

You can't directly get a distribution of error rates by read length with BBMap with reads over 6kbp, since it chops reads. It's easy to post-process, though, since the chopped reads still contain their name.

0
Entering edit mode

Please brian. When I do mapping sequences on the reference genomes using BBmap I consider the overall mapping rate in "pct bases" ?? Thanks

1
Entering edit mode

@Brian referred to this recently here.

0
Entering edit mode

Are you talking to Brian or your brain?

0
Entering edit mode

oups sorry. To Brian. and everyone

0
Entering edit mode
4.6 years ago

It's impossible to make an accurate statement about the distribution of read lengths, since these are dependent on the length of the input DNA and how careful you treat the library. With the most common library prep you have reads of 6-10kb, with a long tail up to ~200kb.
However, with an adapted protocol, very careful handling and old-school DNA extraction you can get reads of up to 970kb, as reported by Nick Loman and Josh Quick.