BRCA1 E23Vfs*17 is a highly recurrent germline variant which is present in gnomAD (https://gnomad.broadinstitute.org/variant/17-41276044-ACT-A) but has no gnomAD frequency in this query:
This is strange. Maybe after all, myvariantinfo won't work? @leexgh please investigate.
@inodb i can't find it on myVariantInfo, could you check please?
@tmazor yikes - excellent catch.
It is in VEP:
http://grch37.rest.ensembl.org/vep/human/hgvs/chr17:g.41276045_41276046delCT?content-type=application/json
but not in myvariant.info:
http://myvariant.info/v1/variant/chr17:g.41276045_41276046delCT
@leexgh could you report to the myvariant.info issue tracker that this variant is missing from their gnomad data? Maybe it's an easy fix for them
https://github.com/biothings/myvariant.info/issues?q=is%3Aissue+is%3Aopen+sort%3Aupdated-desc
I am wondering if this is not just one of the many. Can we check all brca1 and brca2 variants in gnomad and see how many are missing i myvariant.info?
I report the issue to myVariantInfo (https://github.com/biothings/myvariant.info/issues/75)
I second @jjgao 's concern that this could be indicative of a large problem of missing data. Might be worth checking a few other genes also (in addition to BRCA1/2).
@tmazor @jjgao @inodb the myVariantInfo team replays me with this:
We have this variant in MyVariant.info. You could access it by: http://myvariant.info/v1/variant/chr17:g.41276045_41276046del. We normally remove the deleted nucleotides when representing a variant in HGVS id.
Thanks, @leexgh. Let's use change it on our end when calling myvariant then.
@leexgh @inodb : please make sure we code them correctly for other mutation types.
@jjgao Sure I'm working on it
@leexgh is this solved?
@jjgao yes it's done, maybe we can close this issue?