Hi all,
I have Illumina sequencing results of a bacterial genome and a quality score of 35.89 is associated with these data. I know that a quality score of 30 is 99.99% of base calling accuracy based on this but what about the meaning of 35.89? Is there a scale used for quality score?
Thank you very much for your much appreciated help on this!
1 answer
The phred score describes the probability of a base call being incorrect. Wikipedia has a great explanation including the formula, which is -10xlog10(probability of being wrong).
phred = 35.89 # = -10*log10(probability)
35.89/-10 = -3.589
probability = 10^-3.589 = 0.0002576321
Therefore, a base call with a Phred score of 35.89 has a 0.02% chance of being wrong. Or, in other words, 0.02% of bases with a Phred score of 35.89 will have been wrongly assigned.
Log in to answer this question.
-
Why does including BQ<10 bases change Mutect2 calls in NovaSeq X Plus ctDNA data?
written by Anthony •
-
What is the relationship between the per-base sequence quality of a FASTQ file and the QUAL of a VC…
written by ManuelDB
-
Python program that reads a fastq file and calculate Phred score
written by sherrygeek •
-
Meaning of TOP/BOTTOM strand and QC implications
written by ognjen011
-
Recommended Depth of Coverage for Bacterial Genotyping
written by michael •
-
Understanding Trimmomatic Sliding Window Approach
written by nikelle.petrillo
-
Generating Read Length?
written by nikelle.petrillo
-
How to calculate my own phred quality score
written by clear.choi •
-
consistent problem of strand-specific information between TopHat/RSeQC/picard
written by pengchy
-
Estimating Assembly Error Rate For Given Sequencing Depth
written by Leszek
Can you elaborate a bit on how you obtained that quality score? Is that the average quality across all Illumina reads? Was it determined by a specific tool? If yes, how?
Thank you for your help. This quality score that I mentionned is the phred quality score. I would like to interpret this data properly.