Table 3.

Human and mouse assemblies

Assemblies:HumanMouse
Assembly no.:123456
Sequence data:IlluminaIlluminaABI3730IlluminaIlluminaABI3730
Program:ALLPATHS-LGSOAPCeleraALLPATHS-LGSOAPARACHNE
Completeness
 Covered, %91.174.396.288.786.294.2
 Captured, %6.618.61.38.68.03.8
 Uncaptured, %2.37.02.52.75.72.0
 Segmental duplication coverage, %41.112.162.242.327.965.7
 Exon bases covered, %95.181.296.296.792.497.3
Continuity
 Contig N50, kb245.5109161625
 Scaffold N50, kb11,54339917,6467,15634016,871
Contig accuracy
 Ambiguous bases, %0.08000.0400
 1-kb chunks vs. referenceNA12878GRCGRCGRCB6B6B6
 (I) perfect77.188.676.878.0
 (II) ≤0.1% error rate8.72.52.97.0
 (III) ≤1%10.25.76.111.7
 (IV) ≤10%3.13.65.53.62.811.82.4
 (V) >10%0.40.40.70.50.22.40.3
 Base quality, from I–IIIQ33Q36Q35Q33
 Misassembly % of 1-kb chunks, from IV–V3.54.06.24.13.014.22.7
Scaffold accuracy
 Validity at 100 kb, %99.199.599.799.098.899.1
  • An evaluation of human and mouse assemblies is shown. Contigs of size <1 kb were excluded from the analysis. Reference sequences are described in SI Materials and Methods. Assembly no.: Assemblies 1, 4, and 5 are from the data of this paper and are deposited in DDBJ/EMBL/GenBank under accession nos. AEKP00000000, AEKQ00000000, and AEKR00000000, respectively. The versions described in this paper are the first versions, AEKP01000000, AEKQ01000000, and AEKR01000000. For each ambiguity {x1, … , xn}, we inserted x1 into the fasta sequence and referred to x2, … , xn in a note at the locus. Assemblies 2, 3, and 6 are from refs. 3, 12, and 19). Completeness: Contigs were aligned to the reference sequence, with each contig assigned to at most one location. The covered fraction of a genome consists of the fraction of total bases in the reference (exclusive of gaps) that lie under a contig. The captured fraction consists of those bases that lie within a gap in a scaffold. All other bases are uncaptured. Exon coverage was computed from RefSeq gene annotations (http://genome.ucsc.edu/cgi-bin/hgTables). Segmental duplication coverage was computed from http://humanparalogy.gs.washington.edu/build36/oo.weild10kb.join.all.cull.xwparse and http://mouseparalogy.gs.washington.edu/She2008_download/WGAC.tab.gz. Continuity: We report the N50 sizes of contigs and scaffolds, excluding gaps in the latter case. Contig accuracy: We first report the fraction of bases labeled as ambiguous (SI Materials and Methods). We then divide the contigs into 1-kb chunks (as in ref. 9, which, however, used a chunk size of 10 kb). Each chunk is then aligned to the reference sequence using the Smith–Waterman algorithm, seeded on perfect 100-mer matches, to find the optimal placement, and the number of errors (mismatch plus indel bases) is computed. (Contigs having no 100-mer match were treated as novel sequence and ignored for purposes of this analysis. There was <1% of novel sequence in all cases.) The contig is then assigned to one of five mutually exclusive classes on the basis of its error rate. The percentages of chunks landing in each class are listed. Note that for assembly 1, contig accuracy was calculated with respect to two reference sequences. Base quality: Inferred Phred quality (20) of bases in chunk classes I–III. Misassembly %: Total fraction of bases in chunk classes IV–V. Scaffold accuracy: Validity at 100 kb (9): We report the probability that two 100-base sequences in the assembly, separated by 100 kb, and also present in the reference, have the same orientation and are separated by 100 kb ± 10%.