Schema for gnomAD v3.1 - Genome Aggregation Database (gnomAD) Genome Variants v3.1
  Database: hg38    Primary Table: gnomadGenomesVariantsV3_1 Data last updated: 2021-01-13
Big Bed File: /gbdb/hg38/gnomAD/v3.1/variants/genomes.bb
Item Count: 798,965,229
Format description: Browser extensible data (9 fields), plus gnomAD related fields.
fieldexampledescription
chromchr1Chromosome (or contig, scaffold, etc.)
chromStart165970949Start position in chromosome
chromEnd165970950End position in chromosome
namechr1:165970949-165970950 (G/C)Name of item
score0Score from 0-1000
strand.+ or -
thickStart165970949Start of where display should be thick (start codon)
thickEnd165970950End of where display should be thick (stop codon)
reserved95,95,95Used as itemRgb as of 2004-11-22
refGReference Sequence
altCAlternate Sequence
FILTERAC0,AS_VQSRFILTER tags from VCF
AC0Allele Count
AN133970Allele Number
AF0.00000Allele Frequency
faf950.00000Filtering allele frequency (using Poisson 95% CI) for samples
nhomalt0Count of homozygous individuals in samples
rsIddbSnp rsID
genesList of genes affected by variant
annototherAnnotation type: pLoF, missense, synonymous, or other
variation_typeintergenic_variantVariant type(s)
hgvscHGVS c. terms
hgvspHGVS p. terms
popmaxN/APopulation with maximum AF
AC_popmaxN/AAllele count in the population with the maximum AF
AN_popmaxN/ATotal number of alleles in the population with the maximum AF
AF_popmaxN/AMaximum allele frequency across populations (excluding samples of Ashkenazi, Finnish, and indeterminate ancestry)
AC_afr0Alternate allele count for samples of African-American/African ancestry
AN_afr35480Total number of alleles in samples of African-American/African ancestry
AF_afr0.00000Alternate allele frequency in samples of African-American/African ancestry
nhomalt_afr0Count of homozygous individuals in male samples of African-American/African ancestry
AC_ami0Alternate allele count for samples of Amish ancestry
AN_ami862Total number of alleles in samples of Amish ancestry
AF_ami0.00000Alternate allele frequency in samples of Amish ancestry
nhomalt_ami0Count of homozygous individuals in samples of Amish ancestry
AC_amr0Alternate allele count for samples of Latino/Admixed American ancestry
AN_amr13226Total number of alleles in samples of Latino/Admixed American ancestry
AF_amr0.00000Alternate allele frequency in samples of Latino/Admixed American ancestry
nhomalt_amr0Count of homozygous individuals in samples of Latino/Admixed American ancestry
AC_asj0Alternate allele count for samples of Ashkenazi Jewish ancestry
AN_asj3262Total number of alleles in samples of Ashkenazi Jewish ancestry
AF_asj0.00000Alternate allele frequency in samples of Ashkenazi Jewish ancestry
nhomalt_asj0Count of homozygous individuals in samples of Ashkenazi Jewish ancestry
AC_eas0Alternate allele count for samples of East Asian ancestry
AN_eas4626Total number of alleles in samples of East Asian ancestry
AF_eas0.00000Alternate allele frequency in samples of East Asian ancestry
nhomalt_eas0Count of homozygous individuals in samples of East Asian ancestry
AC_fin0Alternate allele count for samples of Finnish ancestry
AN_fin7526Total number of alleles in samples of Finnish ancestry
AF_fin0.00000Alternate allele frequency in samples of Finnish ancestry
nhomalt_fin0Count of homozygous individuals in samples of Finnish ancestry
AC_mid0Alternate allele count for samples of Middle Eastern ancestry
AN_mid292Total number of alleles in samples of Middle Eastern ancestry
AF_mid0.00000Alternate allele frequency in samples of Middle Eastern ancestry
nhomalt_mid0Count of homozygous individuals in samples of Middle Eastern ancestry
AC_nfe0Alternate allele count for samples of Non-Finnish European ancestry
AN_nfe63000Total number of alleles in samples of Non-Finnish European ancestry
AF_nfe0.00000Alternate allele frequency in samples of Non-Finnish European ancestry
nhomalt_nfe0Count of homozygous individuals in samples of Non-Finnish European ancestry
AC_sas0Alternate allele count for samples of South Asian ancestry
AN_sas3896Total number of alleles in samples of South Asian ancestry
AF_sas0.00000Alternate allele frequency in samples of South Asian ancestry
nhomalt_sas0Count of homozygous individuals in samples of South Asian ancestry
AC_oth0Alternate allele count for samples of Other ancestry
AN_oth1800Total number of alleles in samples of Other ancestry
AF_oth0.00000Alternate allele frequency in samples of Other ancestry
nhomalt_oth0Count of homozygous individuals in samples of Other ancestry
cadd_phred2.99500Cadd Phred-like scores ('scaled C-scores') ranging from 1 to 99, based on the rank of each variant relative to all possible 8.6 billion substitutions in the human reference genome. Larger values are more deleterious
revel_scoreN/AdbNSFP's Revel score from 0 to 1. Variants with higher scores are predicted to be more likely to be deleterious
splice_ai_max_dsN/AIllumina's SpliceAI max delta score. Interpreted as the probability of the variant being splice-altering
splice_ai_consequenceN/AThe consequence term associated with the max delta score in 'splice_ai_max_ds'
primate_ai_scoreN/APrimateAI's deleteriousness score from 0 (less deleterious) to 1 (more deleterious)
_startPos165970950Unshifted chromStart position from VCF for link outs

Sample Rows
 
chromchromStartchromEndnamescorestrandthickStartthickEndreservedrefaltFILTERACANAFfaf95nhomaltrsIdgenesannotvariation_typehgvschgvsppopmaxAC_popmaxAN_popmaxAF_popmaxAC_afrAN_afrAF_afrnhomalt_afrAC_amiAN_amiAF_aminhomalt_amiAC_amrAN_amrAF_amrnhomalt_amrAC_asjAN_asjAF_asjnhomalt_asjAC_easAN_easAF_easnhomalt_easAC_finAN_finAF_finnhomalt_finAC_midAN_midAF_midnhomalt_midAC_nfeAN_nfeAF_nfenhomalt_nfeAC_sasAN_sasAF_sasnhomalt_sasAC_othAN_othAF_othnhomalt_othcadd_phredrevel_scoresplice_ai_max_dssplice_ai_consequenceprimate_ai_score_startPos
chr1165970949165970950chr1:165970949-165970950 (G/C)0.16597094916597095095,95,95GCAC0,AS_VQSR01339700.000000.000000otherintergenic_variantN/AN/AN/AN/A0354800.00000008620.0000000132260.000000032620.000000046260.000000075260.00000002920.0000000630000.000000038960.000000018000.0000002.99500N/AN/AN/AN/A165970950
chr1165970949165970950chr1:165970949-165970950 (G/T)0.16597094916597095095,95,95GTPASS21339761.49280e-052.48000e-060rs1571083851otherintergenic_variantNon-Finnish European2630023.17450e-050354820.00000008620.0000000132280.000000032620.000000046260.000000075260.00000002920.0000002630023.17450e-050038960.000000018000.0000002.86000N/AN/AN/AN/A165970950
chr1165970950165970951chr1:165970950-165970951 (G/C)0.16597095016597095195,95,95GCAC001348320.000000.000000otherintergenic_variantN/AN/AN/AN/A0356400.00000008760.0000000132560.000000032760.000000046840.000000076860.00000002880.0000000633200.000000039900.000000018160.0000003.00500N/AN/AN/AN/A165970951
chr1165970950165970951chr1:165970950-165970951 (G/T)0.16597095016597095195,95,95GTPASS21348341.48331e-052.46000e-060otherintergenic_variantAfrican/African American2356425.61136e-052356425.61136e-05008760.0000000132560.000000032760.000000046840.000000076860.00000002880.0000000633200.000000039900.000000018160.0000002.86900N/AN/AN/AN/A165970951
chr1165970952165970953chr1:165970952-165970953 (G/C)0.16597095216597095395,95,95GCPASS11340787.45835e-060.000000otherintergenic_variantLatino/Admixed American1132667.53807e-050353660.00000008600.0000001132667.53807e-050032680.000000046460.000000075960.00000002940.0000000630660.000000039240.000000017920.0000003.01700N/AN/AN/AN/A165970953
chr1165970953165970954chr1:165970953-165970954 (C/G)0.16597095316597095495,95,95CGPASS11341087.45668e-060.000000otherintergenic_variantAfrican/African American1353822.82630e-051353822.82630e-05008720.0000000132660.000000032580.000000046720.000000075900.00000002960.0000000630660.000000039120.000000017940.0000003.32500N/AN/AN/AN/A165970954
chr1165970953165970954chr1:165970953-165970954 (C/T)0.16597095316597095495,95,95CTAS_VQSR11341007.45712e-060.000000rs1486985498otherintergenic_variantNon-Finnish European1630581.58584e-050353840.00000008720.0000000132640.000000032580.000000046720.000000075900.00000002960.0000001630581.58584e-050039120.000000017940.0000003.81900N/AN/AN/AN/A165970954
chr1165970954165970955chr1:165970954-165970955 (G/A)0.16597095416597095595,95,95GAPASS639381172300.5454060.54186317196rs35625097otherintergenic_variantAfrican/African American18287290280.62997818287290280.62997855844047800.517949975340117900.4529261219165030220.545996432164041300.397094292302656880.5319978291742620.6641225630976575680.5380778100154633960.45524133189515660.5715202563.45700N/AN/AN/AN/A165970955
chr1165970955165970956chr1:165970955-165970956 (C/G)0.16597095516597095695,95,95CGPASS11304387.66648e-060.000000otherintergenic_variantSouth Asian137360.0002676660340720.00000008600.0000000129660.000000031960.000000045840.000000075400.00000002940.0000000615080.000000137360.0002676660016820.0000003.47000N/AN/AN/AN/A165970956
chr1165970956165970957chr1:165970956-165970957 (T/C)0.16597095616597095795,95,95TCAC001317980.000000.000000otherintergenic_variantN/AN/AN/AN/A0345640.00000008620.0000000131260.000000032060.000000045960.000000080040.00000002960.0000000616980.000000037300.000000017160.0000005.45900N/AN/AN/AN/A165970957

gnomAD v3.1 (gnomadGenomesVariantsV3_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.

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.

GnomAD Genome Mutational Constraint is based on v3.1.2 and is available only on hg38. It shows the reduced variation caused by purifying natural selection. This is similar to negative selection on loss-of-function (LoF) for genes, but can be calculated for non-coding regions too. Positive values are red and reflect stronger mutation constraint (and less variation), indicating higher natural selection pressure in a region. Negative values are green and reflect lower mutation constraint (and more variation), indicating less selection pressure and less functional effect. Briefly, for any 1kbp window in the genome, a model based on trinucleotide sequence context, base-level methylation, and regional genomic features predicts expected number of mutations, and compares this number to the observed number of mutations using a Z-score (see preprint in the Reference section for details). The chrX scores were added as received from the authors, as there are no de novo mutation data available on chrX (for estimating the effects of regional genomic features on mutation rates), they are more speculative than the ones on the autosomes.

The gnomAD Predicted Constraint Metrics track contains metrics of pathogenicity per-gene as predicted for gnomAD v2.1.1 and identifies genes subject to strong selection against various classes of mutation. This includes data on both the gene and transcript level.

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.
  5. Intergenic variants have not been pre-filtered.

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

To maintain consistency with the gnomAD website, variants are by default labeled according to their chromosomal start position followed by the reference and alternate alleles, for example "chr1-1234-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

The gnomAD v3.1.1 data is unfiltered.

For the v3.1 update only, 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.

The mutational constraints score was updated in October 2022 from a previous, now deprecated, pre-publication version. The old version can be found in our archive directory on the download server. It can be loaded by copying the URL into our "Custom tracks" input box.

Credits

Thanks to the Genome Aggregation Database Consortium for making these data available. The data are released under the ODC Open Database License (ODbL) 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

Chen S, Francioli L, Goodrich J, Collins R, Wang Q, Alfoldi J, Watts N, Vittal C, Gauthier L, Poterba T, Wilson M A genome-wide mutational constraint map quantified from variation in 76,156 human genomes. Biorxiv 2022