Updating equations for the error rates
Consider a site and a read
. The joint probability of the true nucleotide in the read,
, at the site being
and the observed nucleotide at the site
is:
Using Bayes Theorem, the probability of the true nucleotide in the read, , at the site being
, given that we observe
is:
Inserting 19.25 in 19.26 we get:
The equation 19.27 gives us the probabilities for a given read, , and site,
, given the observed nucleotide
, that the true nucleotide is
,
, given our current values for the frequency
(inserted for
) and error rates. Since we know the sequenced nucleotide in each read at each site, we can get new updated values for the error rate of producing an
nucleotide when the true nucleotide is
,
, for
by summing the probabilities of the true nucleotide being
for all reads across all sites for which the sequenced nucleotide is
, and dividing by the sum of all probabilities of the true nucleotide being a
across all reads and all sites:
![$\displaystyle e^*_q(N \rightarrow M) = \frac{\sum_h \sum_{i=1,...,{k_h}: n_i^h ...
...(r_i^k = N \vert n_i^h)}{\sum_h \sum_{i=1,...,{k_h}} P(r_i^k = N \vert n_i^h)} $](img168.gif)