Extract variants from genomic range from gnomAD GRCh38
0
0
Entering edit mode
2.1 years ago
j.lunger18 ▴ 30

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


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.

gnomad variants genome exome • 822 views
1
Entering edit mode

It may be hard to help you without knowing Y:XXXX-XXXX

0
Entering edit mode

Y = chromosome, X = genomic positions --> for p53 this would be something like 17:7661779-7687550

0
Entering edit mode

I know how coordinates work, thank you very much. It doesn't hurt to be specific, as such we can replicate your problem.

1
Entering edit mode

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

0
Entering edit mode

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