GNOMAD exome vcf field missing? number of hemizygous individuals
Entering edit mode
2.9 years ago
brismiller ▴ 50

Hey everyone,

I am searching the GNOMAD 2.1.1 exome vcf and I am having trouble finding a field that I am interested in. I am looking for the number of individuals that are hemizygous for that variant.

For genes that the variant is on a non sex chromosome, this information is easy to find, but not for variants on the sex chromosomes.

e.g. The number of individuals can be found in the field "nhomalt". For example look at the gene , you can see on the right that the number of homozygotes is shown.

Sadly though I cannot find the field that corresponds to the number of hemizygotes.

Any suggestions?

Thanks in advance, this is my first time using GNOMAD

gnomad exome hemizygous • 1.0k views
Entering edit mode
2.9 years ago
rbagnall ★ 1.8k

Lets look at a gene on the X chromosome, say F8

There is a variant at X-154065843-G-A, which is found in 103 males in the exome data and 92 males in the genome data, giving a total of 195 males who have this variant and are hemizygous.

Now lets get the male allele count (thus number of hemizygous males) from the gnomad exome VCF file. Specifically, we want the annotation called AC_male, which is the count of male with the variant.

tabix X:154065843-154065843 | tr ';' '\n' | grep '^AC_male'

This returns AC_male=103, which is the number of hemizygous males who have this variant in the exome VCF.

Entering edit mode

Thank you! That was exactly what I was looking for.


Login before adding your answer.

Traffic: 2241 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6