With the new release of gnomAD, I want to pull variants from several genes for which I have their genomic coordinates. I was successfully able to do this for gnomad GRCh37, so I'm trying to do the same with the following terminal command:
module load google-cloud-sdk
tabix -h gs://gnomad-public/release/3.0/vcf/genomes/gnomad.genomes.r3.0.sites.chr1.vcf.bgz Y:XXXX-XXXX > OUTPUTfile.vcf
While my output file has all of the gnomAD meta-data in the correct vcf format, there are no listed variants... I couldn't find any documentation on how to pull specific locations, this code was just given to me by someone else, so any help would be great.
•
written
by
j.lunger18 •
1 answer
I created a python package based on SQLite databases, where you can easily query all gnomAD variants for GRCh37/38. https://github.com/KalinNonchev/gnomAD_DB I have precomputed SQLite databases for gnomAD WGS for GRCh37/38 in the description of the package. Please take a look there.
Log in to answer this question.
-
The allele with index 2 is not defined in the REF/ALT columns" error after AF filtering with gnomAD…
written by erikatatianacs •
-
How to extract and analyze variant sequences from gnomAD v4 exome VCF
written by Rohan •
-
gnomAD v3 link not working for download
written by peter •
-
How to calculate popmax values for gnomAD variants?
written by Brady •
-
How to extract variants falling in a genomic window from a VCF?
written by francois •
-
Gnomad genome vs Gnomad exom (extremely rare Ref case)
written by Eugene A
-
Pipe output from curl to tabix to slice .vcf file
written by j.lunger18 •
-
vcfanno: issues to annotate VCF with gnomAD r3.0
written by Nicolas Rosewick
-
Import VCF from GnomAD
written by j.lunger18 •
-
Missing SNP in gnomAD hg38 liftover vcfs
written by Nicolas Rosewick
It may be hard to help you without knowing
Y:XXXX-XXXXY = chromosome, X = genomic positions --> for p53 this would be something like 17:7661779-7687550
I know how coordinates work, thank you very much. It doesn't hurt to be specific, as such we can replicate your problem.
Are you using intervals in the format chrY:XXXX-XXXX ?
The older versions of gnomAD were mapped to hg19, with chromosome names like 1,2,3,X,Y The new version 3 of gnomAD is mapped to GRCh38, with chromosome names like chr1, chr2, chrX, chrY
check the chromosome notation for query and indexed vcf (chrY vs Y) ? If they match, expand the coordinates, to check if there are any variants at all at said range. @ j.lunger18