We present a computational method for analyzing deep sequencing data obtained from a genetically diverse sample. The set of reads obtained from a deep sequencing experiment represents a statistical sample of the underlying population. We develop a generative probabilistic model for assigning observed reads to unobserved haplotypes in the presence of sequencing errors. This clustering problem is solved in a Bayesian fashion using the Dirichlet process mixture to define a prior distribution on the unknown number of haplotypes in the mixture. We devise a Gibbs sampler for sampling from the joint posterior distribution of haplotype sequences, assignment of reads to haplotypes, and error rate of the sequencing process to obtain estimates of the local haplotype structure of the population. The method is evaluated on simulated data and on experimental deep sequencing data obtained from HIV samples. © Springer-Verlag Berlin Heidelberg 2009.
CITATION STYLE
Zagordi, O., Geyrhofer, L., Roth, V., & Beerenwinkel, N. (2009). Deep sequencing of a genetically heterogeneous sample: Local haplotype reconstruction and read error correction. In Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics) (Vol. 5541 LNBI, pp. 271–284). https://doi.org/10.1007/978-3-642-02008-7_21
Mendeley helps you to discover research relevant for your work.