I am quite confused about to extract the information I want from the Genomad annotation database (downloaded .vcf ). I've download the .vcf from Genomad and annotated using SnPEff, which resulted in maaany informations, for instance, the allele frequencies from each population, from controls, non_cancer, topmed, etc.
However, what interests me is the total allele frequency, the same as it appears in the browser version. The only annotation I could see that approximates from this is the "AF_raw", but when I compared to the browser version, I saw it was not the same for some variants. Could anyone tell me which INFO could I extract from my vcf to obtain this information?
Do you mean gnomAD? Please give us a couple of example loci where you see a difference between the browser and your VCF. For those loci, also paste the entire
INFOvalues here so we can see what's going on.
Can I ask you why did you have to annotate the vcf file? I have the vcf file from the genomes of gnomAD and they already contain the information you mentioned ("the allele frequencies from each population, from controls, non_cancer, topmed, etc.").
Maybe it was an older version and when I downloaded the new one? (I downloaded the data the first week of March)
On top of the INFO rows I can see for example:
Is the last one the one you are interested in?
I think it's a gap in understanding - the OP meant that they downloaded the gnomAD VCF and used it to annotate their own VCF.
Also, this post should not be an answer. It's not answering the top level question; it's a request for clarification and hence should be a comment. I'm moving it to a comment now, but please be more careful in the future.
Hi Lecob, it is exactly what RamRS has said. I donwloaded the gnomAD VCF in order to use it to annotate MY vcf, with my data. And if you pay attention, this AF doesn't match with the AF described in the browser. The one would match with the browser it usually the adjusted one, i. e. AF_adj, but I couldn't find it. So, not sure I'm right.
You should really just post your SnpEff command if you're having trouble annotating your vcf file, since we know the value you want shoud be in the gnomAD vcf under
AF. But you're saying your gnomAD file has the raw value under AF instead?
Did you get your gnomAD file directly from gnomAD or some other source? I've never seen a gnomAD field named
AF_adj, rather it should be
AFso to me this suggests you're not using an official gnomAD file. Try downloading the file directly from gnomAD and using that instead, or post a few sample lines of your gnomAD file so we can see how the INFO is denoted.
Thank you for trying to help. yes, I downloaded my gnomAD file from gnomAD website. I haven't seen AF_adj for gnomAD either, it was just what I was expecting to find, as there is a field alike for ExAC for example.
Thanks your comments I found the AF field and managed to annotate it. However, comparing both AF_raw and AF fields to what it says in the browser version, it is still different. There is no match. Of course the values are very similar, yet not the same.
I suppose I download the last version: gnomad.exomes.r2.1.sites.vcf.bgz