Figure 3

Genotype call error rate as a function of call confidence score for bins of 10,000 calls in simulated Ciona intestinalis genome data. The stippled blue region shows 95% confidence intervals of the Bayesian posterior probability distribution of the underlying error rate computed from the Beta distribution Beta(n _e  + 1, n _c  − n _e  + 1) conjugate the assumed binomial distribution of observed errors, where n _e and n _c are the number of errors and number of calls in each bin respectively. The red curve shows the best fit error model ${\mathit{a}}_1{\mathit{e}}^{-\frac{{\mathit{s}}^{{\mathit{c}}_1}}{{\mathit{b}}_1}}+{\mathit{a}}_2{\mathit{e}}^{-\frac{{\mathit{s}}^{{\mathit{c}}_2}}{{\mathit{b}}_2}}$, with parameter values a ₁ = 0.49, b ₁ = 2.08, c ₁ = 1.26, a ₂ = 5.47, b ₂ = 0.17, c ₂ = 0.16. χ ² reduced = 0.82.