Calculate the chance that an allele is not sequenced?
1
0
Entering edit mode
4 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 • 171 views
0
Entering edit mode

Or is the probability just 1 in 30?

3
Entering edit mode
4 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

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.