Alignment of 1000 Genomes Project reads to reference assembly GRCh38.

Journal Article (Journal Article)

The 1000 Genomes Project produced more than 100 trillion basepairs of short read sequence from more than 2600 samples in 26 populations over a period of five years. In its final phase, the project released over 85 million genotyped and phased variants on human reference genome assembly GRCh37. An updated reference assembly, GRCh38, was released in late 2013, but there was insufficient time for the final phase of the project analysis to change to the new assembly. Although it is possible to lift the coordinates of the 1000 Genomes Project variants to the new assembly, this is a potentially error-prone process as coordinate remapping is most appropriate only for non-repetitive regions of the genome and those that did not see significant change between the two assemblies. It will also miss variants in any region that was newly added to GRCh38. Thus, to produce the highest quality variants and genotypes on GRCh38, the best strategy is to realign the reads and recall the variants based on the new alignment. As the first step of variant calling for the 1000 Genomes Project data, we have finished remapping all of the 1000 Genomes sequence reads to GRCh38 with alternative scaffold-aware BWA-MEM. The resulting alignments are available as CRAM, a reference-based sequence compression format. The data have been released on our FTP site and are also available from European Nucleotide Archive to facilitate researchers discovering variants on the primary sequences and alternative contigs of GRCh38.

Full Text

Duke Authors

Cited Authors

  • Zheng-Bradley, X; Streeter, I; Fairley, S; Richardson, D; Clarke, L; Flicek, P; 1000 Genomes Project Consortium,

Published Date

  • July 1, 2017

Published In

Volume / Issue

  • 6 / 7

Start / End Page

  • 1 - 8

PubMed ID

  • 28531267

Pubmed Central ID

  • PMC5522380

Electronic International Standard Serial Number (EISSN)

  • 2047-217X

Digital Object Identifier (DOI)

  • 10.1093/gigascience/gix038


  • eng

Conference Location

  • United States