what does "mean base quality" mean?
1
0
Entering edit mode
@christiangriffioen-11642
Last seen 5.8 years ago

We did this quality control with Ioniser in R, with our paired-end nanopore FastQ data

it's about the Y-axis "mean base quality per read" in the quality graph

What does the "Mean base quality per read mean?? it's to low to be PHRED for instance

You can see it here:

https://gyazo.com/df7cd33fddd1b41a7fab138720fda634

(hope the link works)

base quality quality mean base quality • 2.5k views
ADD COMMENT
0
Entering edit mode
Mike Smith ★ 6.6k
@mike-smith
Last seen 9 hours ago
EMBL Heidelberg

This does indeed relate the phred quality score for the base-called data.  

If base calling was to some extent successful each fast5 file contains the called sequence in FASTQ format.  Ideally both strands of a molecule were read and are used together to give a more accurate call.  These reads end up in the '2D' category.  However a fair number of reads don't fall into this category and you have separate information for the 'template' and (potentially) 'complement' strands.  

For every read in each of those bins IONiseR calculates the mean phred score, and the plot you've generated is the boxplot of those mean values.  This of course discards any information about the quality of sequencing along the length of reads, it's a very high level look at the data quality.

The reliability of single base calls is really not that great with nanopore sequencing, and that plot doesn't look particularly worse than other data that I've seen, but you can clearly see that reported confidence in the 2D reads is far better, as one would expect when you information from both strands.


Also, as a bit of advice, it's a good idea to use the name of the package your question relates to as one of the tags when you create a post.  That way the package author is alerted to your comment, rather than just stumbling across it.

ADD COMMENT
0
Entering edit mode

Oh thank you very much!

Weird that they don't mention it on the Y-axis

 

ADD REPLY
0
Entering edit mode

It's my package, so thanks for the feedback.  I'll update the axis label to explicitly state the values are phred scores, and clarify the details in the manual page.

ADD REPLY

Login before adding your answer.

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