De novo fragment assembly with short mate-paired reads: Does the read length matter?
Increasing read length is currently viewed as the crucial condition for fragment assembly with next-generation sequencing technologies. However, introducing mate-paired reads (separated by a gap of length, GapLength) opens a possibility to transform short mate-pairs into long mate-reads of length approximately GapLength, and thus raises the question as to whether the read length (as opposed to GapLength) even matters. We describe a new tool, EULER-USR, for assembling mate-paired short reads and use it to analyze the question of whether the read length matters. We further complement the ongoing experimental efforts to maximize read length by a new computational approach for increasing the effective read length. While the common practice is to trim the error-prone tails of the reads, we present an approach that substitutes trimming with error correction using repeat graphs. An important and counterintuitive implication of this result is that one may extend sequencing reactions that degrade with length "past their prime" to where the error rate grows above what is normally acceptable for fragment assembly.
Chaisson
Mark J
MJ
Bioinformatics Program, University of California San Diego, La Jolla, California 92093, USA. mchaisso@bioinf.ucsd.edu
Brinza
Dumitru
D
Pevzner
Pavel A
PA
D000465:Algorithms; D020029:Base Pairing / Q000502:physiology; D001483:Base Sequence / Q000502:physiology*; D015982:Bias; D002874:Chromosome Mapping / Q000379:methods*; D019295:Computational Biology / Q000379:methods; D004926:Escherichia coli / Q000235:genetics; D005819:Genetic Markers; D016680:Genome, Bacterial; D008954:Models, Biological; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*
D005819:Genetic Markers