Given the estimated error rate E, the read depth N and the count of alternative allele X, a binomial model was used to compute the P-value for the event that more than X alternative alleles were observed when the null hypothesis is true. When a position is covered by both forward and reverse reads, we used the Stouffer's Z-score method (weighted by the read depths) to combine P-values of individual tests into a single P-value.
