Open Access Highly Accessed Open Badges Methodology article

Estimation of sequencing error rates in short reads

Xin Victoria Wang12, Natalie Blades3, Jie Ding12, Razvan Sultana14 and Giovanni Parmigiani12*

Author affiliations

1 Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA 02215, USA

2 Department of Biostatistics, Harvard School of Public Health, Boston, MA 02115, USA

3 Department of Statistics, Brigham Young University, Provo, UT 84602, USA

4 Department of Bioinformatics, Boston University, Boston, MA 02215, USA

For all author emails, please log on.

Citation and License

BMC Bioinformatics 2012, 13:185  doi:10.1186/1471-2105-13-185

Published: 30 July 2012



Short-read data from next-generation sequencing technologies are now being generated across a range of research projects. The fidelity of this data can be affected by several factors and it is important to have simple and reliable approaches for monitoring it at the level of individual experiments.


We developed a fast, scalable and accurate approach to estimating error rates in short reads, which has the added advantage of not requiring a reference genome. We build on the fundamental observation that there is a linear relationship between the copy number for a given read and the number of erroneous reads that differ from the read of interest by one or two bases. The slope of this relationship can be transformed to give an estimate of the error rate, both by read and by position. We present simulation studies as well as analyses of real data sets illustrating the precision and accuracy of this method, and we show that it is more accurate than alternatives that count the difference between the sample of interest and a reference genome. We show how this methodology led to the detection of mutations in the genome of the PhiX strain used for calibration of Illumina data. The proposed method is implemented in an R package, which can be downloaded from∼vwang/shadowRegression.html webcite.


The proposed method can be used to monitor the quality of sequencing pipelines at the level of individual experiments without the use of reference genomes. Furthermore, having an estimate of the error rates gives one the opportunity to improve analyses and inferences in many applications of next-generation sequencing data.