Schema for GWAS Variants - GWAS Variants 2022-03-08 (lifted) |
|---|
| |
Database: hub_567047_hs1 Primary Table: hub_567047_gwasSNPs2022-03-08
VCF File Download: /gbdb/hs1/gwasSNPs2022-03-08/chm13v2.0_GWASv1.0rsids_e100_r2022-03-08.vcf.gz
Format description: The fields of a Variant Call Format data line
| field | description |
|---|
| chrom | An identifier from the reference genome | | pos | The reference position, with the 1st base having position 1 | | id | Semi-colon separated list of unique identifiers where available | | ref | Reference base(s) | | alt | Comma separated list of alternate non-reference alleles called on at least one of the samples | | qual | Phred-scaled quality score for the assertion made in ALT. i.e. give -10log_10 prob(call in ALT is wrong) | | filter | PASS if this position has passed all filters. Otherwise, a semicolon-separated list of codes for filters that fail | | info | Additional information encoded as a semicolon-separated series of short keys with optional comma-separated values | | format | If genotype columns are specified in header, a semicolon-separated list of of short keys starting with GT | | genotypes | If genotype columns are specified in header, a tab-separated set of genotype column values; each value is a colon-separated list of values corresponding to keys in the format column |
|
| |
|
|
|
Sample Rows
|
|---|
| |
| chrom | pos | id | ref | alt | qual | filter | info |
|---|
| chr1 | 261839 | rs2977608 | A | C,T | . | PASS | COMMON;FREQ=1000Genomes:0.5106,0.4894,.|GENOME_DK:0.15,0.85,.|GoNL:0.2565,0.7435,.|KOREAN:0.5648,0.4352,.|Korea1K:0.5885,0.4115, ... |
| chr1 | 273986 | rs141175086 | C | T | . | PASS | FREQ=1000Genomes:0.989,0.01098|GnomAD:0.9934,0.006578|KOREAN:0.9747,0.02527|Korea1K:0.9738,0.0262|SGDP_PRJ:0.5,0.5|TOMMO:0.9425, ... |
| chr1 | 375395 | rs2272756 | G | A,C | . | PASS | COMMON;FREQ=1000Genomes:0.8632,0.1368,.|ALSPAC:0.7532,0.2468,.|Estonian:0.7415,0.2585,.|GENOME_DK:0.625,0.375,.|GoNL:0.7255,0.27 ... |
| chr1 | 384464 | rs13303065 | T | A,C | . | MismatchedRefAllele | AttemptedAlleles=C*->A,T;AttemptedLocus=chr1:384464-384464;COMMON;FREQ=1000Genomes:0.5467,.,0.4533|ALSPAC:0.343,.,0.657|Estonian ... |
| chr1 | 387924 | rs115438739 | G | A | . | PASS | COMMON;FREQ=1000Genomes:0.9922,0.007788|ALSPAC:0.9779,0.02206|GnomAD:0.9868,0.01318|GoESP:0.982,0.01801|NorthernSweden:0.955,0.0 ... |
| chr1 | 387978 | rs13303010 | A | G | . | MismatchedRefAllele | AttemptedAlleles=G*->A;AttemptedLocus=chr1:387978-387978;COMMON;FREQ=1000Genomes:0.3656,0.6344|ALSPAC:0.08381,0.9162|ExAC:0.1828 ... |
| chr1 | 389115 | rs13303327 | A | G | . | MismatchedRefAllele | AttemptedAlleles=G*->A;AttemptedLocus=chr1:389115-389115;COMMON;FREQ=1000Genomes:0.3636,0.6364|GENOME_DK:0.025,0.975|GnomAD:0.35 ... |
| chr1 | 389681 | rs112905931 | T | A,C | . | MismatchedRefAllele | AttemptedAlleles=C*->A,T;AttemptedLocus=chr1:389681-389681;COMMON;FREQ=1000Genomes:0.6222,.,0.3778|ALSPAC:0.3581,.,0.6419|Estoni ... |
| chr1 | 391276 | rs201385366 | C | G,T | . | PASS | FREQ=1000Genomes:0.9994,.,0.000599|Estonian:0.9987,.,0.001339|ExAC:0.9988,.,0.001164|FINRISK:0.9897,.,0.01034|GnomAD:0.9985,.,0. ... |
| chr1 | 393854 | rs140019196 | A | G | . | PASS | FREQ=ExAC:0.9998,0.0001597|GnomAD:0.9996,0.0004434|GnomAD_exomes:0.9999,0.0001284|GoESP:0.9996,0.0003847|TOPMED:0.9996,0.0003967 ... |
|
| |
|
|
|
GWAS Variants (hub_567047_gwasSNPs2022-03-08) Track Description
|
|---|
| |
Description
A GWAS Catalog VCF was generated by intersecting RefSeq IDs in the GWAS Catalog associations v1.0 file (ebi.ac.uk/gwas/api/search/downloads/full, accessed 2022-03-08) with all dbSNP build 155 variants on the primary contigs for Chromosomes 1-22, Chromsome X and Chromosome Y. These variants were lifted over from GRCh38 to the T2T-CHM13 assembly. This track contains GWAS Catalog variants that lifted over from GRCh38 to the T2T-CHM13v2.0 assembly. This includes variants that lifted over perfectly, as well as variants that failed initial liftover due to reference/alternative allele swaps but were recovered on subsequent liftover, with reference and alternative alleles swapped appropriately.
These two sets of variants are included together in this track. If you are interested in downloading these sets separately (i.e., variants that lifted over perfectly vs. recovered variants with ref/alt allele swaps) they can be accessed here: https://s3-us-west-2.amazonaws.com/human-pangenomics/index.html?prefix=T2T/CHM13/assemblies/annotation/liftover/.
Methods
We performed liftover using the GATK release 4.1.9 LiftoverVcf (Picard Version 2.23.3) tool with the default parameters. This successfully lifts over variants that map exactly from GRCh38 to T2T-CHM13v2.0 but does not recover variants with swapped reference and alternative alleles. To recover variants with swapped reference/alternative alleles, we ran LiftoverVCF again, with the RECOVER_SWAPPED_REF_ALT flag. Notably, this feature does not recover multiallelic variants, so to recover these variants, we first separated them into multiple biallelic variants, performed liftover using the RECOVER_SWAPPED_REF_ALT tag, and converted them back to their multiallelic representations.
Contacts
References
Van der Auwera GA & O'Connor BD. (2020). Genomics in the Cloud: Using Docker, GATK, and WDL in Terra (1st Edition). O'Reilly Media.
| |
|
|
|
US Server