Accuracy and quality of massively parallel DNA pyrosequencing
Data file 2: These sequences correspond to those included in Additional data files 3-5 but include only the final sequence information in fasta format. (3.305Mb)
Data file 3: Text translation of the original GS20 binary output (sff) file for the first of three PicoTiterPlate™ regions used in the analysis. (36.58Mb)
Data file 4: Text translation of the original GS20 binary output (sff) file for the second of three PicoTiterPlate™ regions used in the analysis. (37.93Mb)
Huse, Susan M.
Huber, Julie A.
Morrison, Hilary G.
Sogin, Mitchell L.
Mark Welch, David B.
MetadataShow full item record
Massively parallel pyrosequencing systems have increased the efficiency of DNA sequencing, although the published per-base accuracy of a Roche GS20 is only 96%. In genome projects, highly redundant consensus assemblies can compensate for sequencing errors. In contrast, studies of microbial diversity that catalogue differences between PCR amplicons of ribosomal RNA genes (rDNA) or other conserved gene families cannot take advantage of consensus assemblies to detect and minimize incorrect base calls. We performed an empirical study of the per-base error rate for the Roche GS20 system using sequences of the V6 hypervariable region from cloned microbial ribosomal DNA (tag sequencing). We calculated a 99.5% accuracy rate in unassembled sequences, and identified several factors that can be used to remove a small percentage of low-quality reads, improving the accuracy to 99.75% or better. By using objective criteria to eliminate low quality data, the quality of individual GS20 sequence reads in molecular ecological applications can surpass the accuracy of traditional capillary methods.
© 2007 Huse et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The definitive version was published in Genome Biology 8 (2007): R143, doi:10.1186/gb-2007-8-7-r143.Additional data file 1 is a fasta file of the 43 known sequences used. Additional data file 2 is a gzip-compressed fasta file of the sequences output by the GS20. These sequences correspond to those included in Additional data files 3, 4, 5 but include only the final sequence information. Additional data files 3, 4, 5 are three compressed text files representing the text translations of the original GS20 binary output (sff) files for all of the sequencing used in the analysis, including sequence, flowgram and other run information. GS20 data are reported by region of the PicoTiterPlate™; we sequenced three plate regions.
The following license files are associated with this item: