论文信息 - 2SNV: Quasispecies reconstruction from PacBio reads

2SNV: Quasispecies reconstruction from PacBio reads

Pacific Biosciences (PacBio) sequencing is providing thousands of reads with the length up to 10,000 bases. In most cases this length is enough to cover entire region of interest however this technology has high (≈15%) error rate. We propose a method for viral haplotype reconstruction using long SMRT reads. The proposed method based on correlation between SNVs and maximum likelihood for frequency estimation. When applied to PacBio reads from an Influenza A Virus (IAV) sample with ten variants, our method was able to reconstruct the nine most frequent.

[1] Nicholas Eriksson,et al. ShoRAH: estimating the genetic diversity of a mixed sample from next-generation sequencing data , 2011, BMC Bioinformatics.

[2] Richard Durbin,et al. Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[3] Alice McCarthy. Third generation DNA sequencing: pacific biosciences' single molecule real time technology. , 2010, Chemistry & biology.

[4] Ion I. Mandoiu,et al. kGEM: An EM-based algorithm for local reconstruction of viral quasispecies , 2013, 2013 IEEE 3rd International Conference on Computational Advances in Bio and medical Sciences (ICCABS).