Schema for gnomAD v3.1.1 - Genome Aggregation Database (gnomAD) Genome Variants v3.1.1
  Database: hg38    Primary Table: gnomadGenomesVariantsV3_1_1 Data last updated: 2021-06-11
Big Bed File: /gbdb/hg38/gnomAD/v3.1.1/genomes.bb
Item Count: 491,375,940
Format description: gnomAD v3.1.1 variant data and gnomAD v3.1 chrM data
fieldexampledescription
chromchr1Chromosome (or contig, scaffold, etc.)
chromStart166022215Start position in chromosome
chromEnd166022216End position in chromosome
name27cebde3972ed646c4441287735857cfmd5 of data used as a key into external data file
score0Score from 0-1000
strand.+ or -
thickStart166022215Start of where display should be thick (start codon)
thickEnd166022216End of where display should be thick (stop codon)
reserved95,95,95Used as itemRgb as of 2004-11-22
refTReference Sequence
altCAlternate Sequence
FILTERPASSFILTER tags from VCF
AC1Allele Count
AN152174Allele Number
AF6.57142e-06Allele Frequency
faf950.00000Filtering allele frequency (using Poisson 95% CI) for samples
nhomalt0Count of homozygous individuals in samples
rsIddbSnp rsID
genesRPS3AP10List of genes affected by variant
annototherAnnotation type: pLoF, missense, synonymous, or other
variation_typenon_coding_transcript_exon_variantVariant type(s)
_startPos166022216Unshifted chromStart position from VCF for link outs
_displayNamechr1-166022216-T-CgnomAD display name
_dataOffset383549341255Offset into gnomad.v3.1.1.details.tab.gz for line with more info
_dataLen935Length of the line in gnomad.v3.1.1.details.tab.gz

Sample Rows
 
chromchromStartchromEndnamescorestrandthickStartthickEndreservedrefaltFILTERACANAFfaf95nhomaltrsIdgenesannotvariation_type_startPos_displayName_dataOffset_dataLen
chr116602221516602221627cebde3972ed646c4441287735857cf0.16602221516602221695,95,95TCPASS11521746.57142e-060.000000RPS3AP10othernon_coding_transcript_exon_variant166022216chr1-166022216-T-C383549341255935
chr1166022223166022224bcddfcb19e539c2be7dd9dc10a783a5c0.16602222316602222495,95,95GAPASS11521426.57281e-060.000000rs751452093RPS3AP10othernon_coding_transcript_exon_variant166022224chr1-166022224-G-A383549342191936
chr1166022230166022231cb9a570e0f665bf8f7e647baf889077d0.16602223016602223195,95,95AGPASS11522146.56970e-060.000000RPS3AP10othernon_coding_transcript_exon_variant166022231chr1-166022231-A-G383549343128936
chr1166022231166022232631bde6ef7cad226888a87cb4d3fe7440.16602223116602223295,95,95CGPASS11520966.57479e-060.000000rs1007767227RPS3AP10othernon_coding_transcript_exon_variant166022232chr1-166022232-C-G383549345002936
chr116602223116602223273463e89c9196211da436ac5af28bb9c0.16602223116602223295,95,95ACAPASS21522141.31394e-052.18000e-060rs1390260903RPS3AP10othernon_coding_transcript_exon_variant166022231chr1-166022231-AC-A383549344065936
chr1166022231166022232dbbbd22848b77fd8d331597facebd76c0.16602223116602223295,95,95CTPASS11520966.57479e-060.000000rs1007767227RPS3AP10othernon_coding_transcript_exon_variant166022232chr1-166022232-C-T383549345939936
chr116602224216602224397618629defca8ca6be14af1ea00bf760.16602224216602224395,95,95CTPASS31521541.97169e-055.24000e-060rs1304269434RPS3AP10othernon_coding_transcript_exon_variant166022243chr1-166022243-C-T383549346876936
chr1166022248166022249add95b7f236652b4e13b7b8f0e338b250.16602224816602224995,95,95AGPASS11521846.57099e-060.000000rs868005923RPS3AP10othernon_coding_transcript_exon_variant166022249chr1-166022249-A-G383549347813936
chr11660222491660222503e4caac9c3ba81b8317ae343ca5692bb0.16602224916602225095,95,95TCPASS11521326.57324e-060.000000rs1325472800RPS3AP10othernon_coding_transcript_exon_variant166022250chr1-166022250-T-C383549348750936
chr116602225116602225229f37c0a22a48f7cfe7159d0f46cfa5b0.16602225116602225295,95,95GAPASS31520941.97246e-055.24000e-060rs1323076827RPS3AP10othernon_coding_transcript_exon_variant166022252chr1-166022252-G-A383549349687940

gnomAD v3.1.1 (gnomadGenomesVariantsV3_1_1) Track Description
 

Description

The gnomAD v3.1 track shows variants from 76,156 whole genomes (and no exomes), all mapped to the GRCh38/hg38 reference sequence. 4,454 genomes were added to the number of genomes in the previous v3 release. For more detailed information on gnomAD v3.1, see the related blog post.

On the gnomAD side, the gnomAD v3.1.1 track contains the same underlying data as v3.1, but with minor corrections to the VEP annotations and dbSNP rsIDs. On the UCSC side, we have now included the mitochondrial chromosome data that was released as part of gnomAD v3.1 (but after the UCSC version of the track was released). For more information about gnomAD v3.1.1, please see the related changelog.

The gnomAD v2 tracks show variants from 125,748 exomes and 15,708 whole genomes, all mapped to the GRCh37/hg19 reference sequence and lifted to the GRCh38/hg38 assembly. The data originate from 141,456 unrelated individuals sequenced as part of various population-genetic and disease-specific studies collected by the Genome Aggregation Database (gnomAD), release 2.1.1. Raw data from all studies have been reprocessed through a unified pipeline and jointly variant-called to increase consistency across projects. For more information on the processing pipeline and population annotations, see the following blog post and the 2.1.1 README.

gnomAD v2 data are based on the GRCh37/hg19 assembly. These tracks display the GRCh38/hg38 lift-over provided by gnomAD on their downloads site.

For questions on the gnomAD data, also see the gnomAD FAQ.

More details on the Variant type(s) can be found on the Sequence Ontology page.

Display Conventions and Configuration

gnomAD v3.1.1

The gnomAD v3.1.1 track version follows the same conventions and configuration as the v3.1 track, except as noted below.

  1. There are additional FILTER field filters: AS_VQSR, indel_stack (chrM only), and npg (chrM only).
  2. Where possible, variants overlapping multiple transcripts/genes have been collapsed into one variant, with additional information available on the details page, which has roughly halved the number of items in the bigBed.
  3. The bigBed has been split into two files, one with the information necessary for the track display, and one with the information necessary for the details page. For more information on this data format, please see the Data Access section below.
  4. The VEP annotation is shown as a table instead of spread across multiple fields.

gnomAD v3.1

By default, a maximum of 50,000 variants can be displayed at a time (before applying the filters described below), before the track switches to dense display mode.

Mouse hover on an item will display many details about each variant, including the affected gene(s), the variant type, and annotation (missense, synonymous, etc).

Clicking on an item will display additional details on the variant, including a population frequency table showing allele count in each sub-population.

Following the conventions on the gnomAD browser, items are shaded according to their Annotation type:

pLoF
Missense
Synonymous
Other

Label Options

By default, variants are labeled according to their chromosomal position followed by the reference and alternate alleles, for example "chr1:1234-1235 T/CAG".

dbSNP rsID's are also available as an additional label, if the variant is present in dbSnp.

Filtering Options

Three filters are available for these tracks:

  • FILTER: Used to exclude/include variants that failed Random Forest (RF), Inbreeding Coefficient (Inbreeding Coeff), or Allele Count (AC0) filters. The PASS option is used to include/exclude variants that pass all of the RF, InbreedingCoeff, and AC0 filters, as denoted in the original VCF.
  • Annotation type: Used to exclude/include variants that are annotated as Probability Loss of Function (pLoF), Missense, Synonymous, or Other, as annotated by VEP version 85 (GENCODE v19).
  • Variant Type: Used to exclude/include variants according to the type of variation, as annotated by VEP v85.
There is one additional configurable filter on the minimum minor allele frequency.

gnomAD v2.1.1

The gnomAD v2.1.1 track follows the standard display and configuration options available for VCF tracks, briefly explained below.
  • In mode, a vertical line is drawn at the position of each variant.
  • In mode, "ref" and "alt" alleles are displayed to the left of a vertical line with colored portions corresponding to allele counts. Hovering the mouse pointer over a variant pops up a display of alleles and counts.
  • Filtering Options

    Four filters are available for these tracks, the same as the underlying VCF:

    • AC0: Allele Count 0 after filtering out low confidence genotypes (GQ < 20; DP < 10; and AB < 0.2 for het calls))
    • InbreedingCoeff: Inbreeding Coefficient < -0.3
    • RF: Used to exclude/include variants that failed Random Forest filtering thresholds of 0.055272738028512555, 0.20641025579497013 (probabilities of being a true positive variant) for SNPs, indels)
    • Pass: Variant passes all 3 filters

    There are two additional filters available, one for the minimum minor allele frequency, and a configurable filter on the QUAL score.

    UCSC Methods

    For the v3.1 update, In order to cut down on the amount of displayed data, the following variant types have been filtered out, but are still viewable in the gnomAD browser:

    • Regulatory Region Variants
    • Downstream/Upstream Gene Variants
    • Transcription Factor Binding Site Variants

    For the full steps used to create the track at UCSC, please see the section denoted "gnomAD v3.1 update" in the hg38 makedoc.

    Data Access

    The raw data can be explored interactively with the Table Browser, or the Data Integrator. For automated analysis, the data may be queried from our REST API, and the genome annotations are stored in files that can be downloaded from our download server, subject to the conditions set forth by the gnomAD consortium (see below). Variant VCFs can be found in the vcf/ subdirectory. The v3.1 and v3.1.1 variants can be found in a special directory as they have been transformed from the underlying VCF.

    For the v3.1.1 variants in particular, the underlying bigBed only contains enough information necessary to use the track in the browser. The extra data like VEP annotations and CADD scores are available in the same directory as the bigBed but in the files gnomad.v3.1.1.details.tab.gz and gnomad.v3.1.1.details.tab.gz.gzi. The gnomad.v3.1.1.details.tab.gz contains the gzip compressed extra data in JSON format, and the .gzi file is available to speed searching of this data. Each variant has an associated md5sum in the name field of the bigBed which can be used along with the _dataOffset and _dataLen fields to get the associated external data, as show below:

    # find item of interest:
    bigBedToBed genomes.bb stdout | head -4 | tail -1
    chr1    12416    12417    854246d79dc5d02dcdbd5f5438542b6e    [..omitted for brevity..]    chr1-12417-G-A    67293    902
    
    # use the final two fields, _dataOffset and _dataLen (add one to _dataLen to include a newline), to get the extra data:
    bgzip -b 67293 -s 903 gnomad.v3.1.1.details.tab.gz
    854246d79dc5d02dcdbd5f5438542b6e    {"DDX11L1": {"cons": ["non_coding_transcript_variant",  [..omitted for brevity..]
    

    The data can also be found directly from the gnomAD downloads page. Please refer to our mailing list archives for questions, or our Data Access FAQ for more information.

    Credits

    Thanks to the Genome Aggregation Database Consortium for making these data available. The data are released under the ODC Open Database License (OBdL) as described here.

    References

    Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alfoldi J, Wang Q, Collins RL, Laricchia KM, Ganna A, Birnbaum DP et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. doi: https://doi.org/10.1101/531210.

    Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O'Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016 Aug 17;536(7616):285-91. PMID: 27535533; PMC: PMC5018207