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.