Illumina Q score
1
1
Entering edit mode
6 weeks ago
A_heath ▴ 60

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!

Illumina quality score • 362 views
ADD COMMENT
1
Entering edit mode

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?

ADD REPLY
1
Entering edit mode

Thank you for your help. This quality score that I mentionned is the phred quality score. I would like to interpret this data properly.

ADD REPLY
5
Entering edit mode
6 weeks ago

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.

ADD COMMENT
1
Entering edit mode

Thank you very much for this clear explanation!!

ADD REPLY

Login before adding your answer.

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