title: Mapping single molecule sequencing reads using basic local alignment with successive refinement (BLASR): application and theory.
abstract: BACKGROUND
  Recent methods have been developed to perform high-throughput sequencing of DNA by Single Molecule Sequencing (SMS). While Next-Generation sequencing methods may produce reads up to several hundred bases long, SMS sequencing produces reads up to tens of kilobases long. Existing alignment methods are either too inefficient for high-throughput datasets, or not sensitive enough to align SMS reads, which have a higher error rate than Next-Generation sequencing.
  
  
  RESULTS
  We describe the method BLASR (Basic Local Alignment with Successive Refinement) for mapping Single Molecule Sequencing (SMS) reads that are thousands of bases long, with divergence between the read and genome dominated by insertion and deletion error. The method is benchmarked using both simulated reads and reads from a bacterial sequencing project. We also present a combinatorial model of sequencing error that motivates why our approach is effective.
  
  
  CONCLUSIONS
  The results indicate that it is possible to map SMS reads with high accuracy and speed. Furthermore, the inferences made on the mapability of SMS reads using our combinatorial model of sequencing error are in agreement with the mapping accuracy demonstrated on simulated reads.
authors.affiliation: Department of Mathematics, University of California, San Diego, 9500 Gilman Dr, CA, La Jolla, USA.
authors.forename: Mark J
authors.initials: MJ
authors.lastname: Chaisson
authors.forename: Glenn
authors.initials: G
authors.lastname: Tesler
mesh_terms: D000465:Algorithms; D001483:Base Sequence; D004247:DNA / Q000235:genetics*; D016678:Genome; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D012984:Software
chemical_list: D004247:DNA