Any Idea why I'm not getting any gnomAD results for NM_001029.3:c.-22C>G?
Here it is in gnomAD:
my VEP results:
You could try wAnnovar which seems to work for that variant.
This variant is listed as suspect by dbSNP:
So by default, the VEP only reports information about the co-located variant from the HGMD database:
To retrieve information from the gnomAD database, you can select the 'Include flagged variants' option in the VEP submission form.
1) Does this mean it is suspected by dbSNP for being a false positive variant in the gnomad?
2) How would I 'Include flagged variants' via the command line tool?
1) The variant was flagged as suspect by dbSNP version 151, which is the dbSNP version available when using the GRCh37 Ensembl browser and VEP. The variant is not flagged as suspect in the latest version of dbSNP (dbSNP 153), which is available when using the GRCh38 Ensembl browser and VEP tool:
So you can retrieve the GnomAD data if you use the GRCh38 Ensembl browser and VEP:
I'm not sure why dbSNP flagged this variant a suspect in version 151. It may be for a number of different reasons.
2) You can include flagged variants in the command line VEP tool by using the --failed 1 flag:
You can find the number of failed variants using MySQL from the failed_variation table:
or using the Perl API:
thanks again :)
I have a very large variant set and I must use GRCh37.
I'm using the command line tool to annotate.
I'm not sure I want to include failed variants for the entire set because that would include the up-to-date failed variants as well.
Before I dig in, perhaps you can save me the effort if you happen to have a very rough estimate of how many variants I'm dealing with? (failure status changed from 151 to 153?)
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy