How to calculate my own phred quality score
1
0
Entering edit mode
7.0 years ago
clear.choi ▴ 30

I am planning to calculate my own phred quality score. I need to share some idea how to calculate my own Q Score. I have sequences like below:

AAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAA
                    ....
TAAAAAAAAAAAAAAAAAAAAA

So in the above case, There is A: 100 and T: 1 on the first base. So that position has 99% accuracy. And from 99% accuracy how I can get Q = 20?

What if there is A: 151 and 1 T on the first position? Anyone knows how to calculate that score from the percentage?

Thank you

phred-quality-score fastq Qscore • 10k views
ADD COMMENT
2
Entering edit mode
7.0 years ago
Tej Sowpati ▴ 250

From Wikipedia page:

Q=-10\ \log _{10}P

Keep in mind that the original Phred score is intrinsic to the sequencing machine, which determines the score based on peak shape and resolution at every base. Assuming you still want to represent your percentages as Phred scores, you have to use the probability of the wrong base in the above equation. For example, if a base is 99% of the times 'A', then the probability of the base not being 'A' is 0.01. And Q score can then be derived as

Q = -10 * log10(0.01) = 20

Cheers

ADD COMMENT

Login before adding your answer.

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