Calculate the chance that an allele is not sequenced?
1
0
Entering edit mode
5 weeks ago
vpsev3 ▴ 10

Hello,

With 30X depth NGS sequencing, what is the chance at one position that an allele is not sequenced at all with 30 reads ?

Should we do 2^30 = 1 073 741 824 ?

Thank you.

sequencing • 174 views
ADD COMMENT
0
Entering edit mode

Or is the probability just 1 in 30?

ADD REPLY
3
Entering edit mode
5 weeks ago
Fabio Marroni ★ 2.8k

The probability is a number between 0 and 1, so by definition it cannot be 2^30. However, I think you were close.

First of all: I assume you want to know the probability that, in a heterozygous individual, sequencing 30 reads, you always sequence allele 1 and never allele 2.

Each read has a 0.5 probability of sequencing allele 1, and 0.5 of sequencing allele 2. The probability that with 2 reads you always sequence allele 1 is 0.5 * 0.5 (0.5^2), with 3 is 0.5^3, and so on... thus the probability that with 30 reads you always sequence allele 1 is 0.5^30, i.e. 9^-10.

I hope this helps

ADD COMMENT
1
Entering edit mode

While the above is valid in theory, in practice, and depending on the situation, I think the values ought to be corrected for multiple comparisons i.e. of all the variants we found, what is the likelihood of observing the imbalance that we see in this case.

ADD REPLY

Login before adding your answer.

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