How to calculate my own phred quality score
1
0
Entering edit mode
6.1 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 a sequences like below,

 

        AAAAAAAAAAAAAAAAAAAAAA

        AAAAAAAAAAAAAAAAAAAAAA

        AAAAAAAAAAAAAAAAAAAAAA

        AAAAAAAAAAAAAAAAAAAAAA

        AAAAAAAAAAAAAAAAAAAAAA

                            ....

 

        TAAAAAAAAAAAAAAAAAAAAA

 

 

So 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.

 

 

 

Qscore phred quality score quality score fastq • 8.9k views
ADD COMMENT
2
Entering edit mode
6.1 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: 1898 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