We use a new statistical test, which I find beautiful in its simplicity. The idea is to exploit the length of the reads to look at multiple loci at the same time. The p-value of a *group of* variants being due only to sequencing errors is 0.1^(ab)*C(m,b)*C(n,a), C denoting binomial coefficients.