An MCMC algorithm for haplotype assembly from whole-genome sequence data

(Downloading may take up to 30 seconds. If the slide opens in your browser, select File -> Save As to save it.)

Click on image to view larger version.

Figure 6.
Figure 6.

Results of running the MCMC algorithm with different Γ on a fragment matrix with n = 200 columns (from chromosome 22 of HuRef genome). (A) A comparison of the HASH algorithm against two other MCMC algorithms: (1) ℳ(Γ1) and (2) ℳ(Γ) where Γ was computed using the recursive graph-partitioning algorithm G(X). All algorithms were initialized with a random haplotype pair. (B) Comparison of HASH algorithm initialized with a random haplotype against ℳ(Γ) (graph-partitioning) initialized with a good haplotype. Note that we are zooming in on the first 10,000 steps in the iteration.

This Article

  1. Genome Res. 18: 1336-1346

Preprint Server