Error and Error Mitigation in Low-Coverage Genome Assemblies
Name
Hubisz-2011-Error and Error Miti.pdf
Size
565.76 KB
Format
Adobe PDF
Checksum (MD5)
44ce39ad44317c0c3dc7de171eee290f
Author(s) • • •
Hubisz, Melissa J.
Lin, Michael F.
Kellis, Manolis
Siepel, Adam
Date Issued
February 2011
Journal
PLoS ONE
Publisher
Public Library of Science
Citation
Hubisz, Melissa J. et al. “Error and Error Mitigation in Low-Coverage Genome Assemblies.” Ed. Thomas Mailund. PLoS ONE 6.2 (2011) : e17034.
Version
Final published version
Abstract
The recent release of twenty-two new genome sequences has dramatically increased the data available for mammalian comparative genomics, but twenty of these new sequences are currently limited to ~2× coverage. Here we examine the extent of sequencing error in these 2× assemblies, and its potential impact in downstream analyses. By comparing 2× assemblies with high-quality sequences from the ENCODE regions, we estimate the rate of sequencing error to be 1–4 errors per kilobase. While this error rate is fairly modest, sequencing error can still have surprising effects. For example, an apparent lineage-specific insertion in a coding region is more likely to reflect sequencing error than a true biological event, and the length distribution of coding indels is strongly distorted by error. We find that most errors are contributed by a small fraction of bases with low quality scores, in particular, by the ends of reads in regions of single-read coverage in the assembly. We explore several approaches for automatic sequencing error mitigation (SEM), making use of the localized nature of sequencing error, the fact that it is well predicted by quality scores, and information about errors that comes from comparisons across species. Our automatic methods for error mitigation cannot replace the need for additional sequencing, but they do allow substantial fractions of errors to be masked or eliminated at the cost of modest amounts of over-correction, and they can reduce the impact of error in downstream phylogenomic analyses. Our error-mitigated alignments are available for download.
MIT Department
Massachusetts Institute of Technology. Computer Science and Artificial Intelligence Laboratory
Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science
Terms of Use
Creative Commons Attribution
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1371/journal.pone.0017034