<div dir="ltr">Hi,<div><br></div><div>This variant and its MAF is not reported by the VEP REST endpoint because it has failed our QC process, as you can see from the message in the grey box at the top of the variant page you linked to.</div><div><br></div><div>Currently there's no way to force inclusion of such failed variants using the REST service, though it is a simple fix to add the flag necessary to do this. We'll look into getting this incorporated soon.</div><div><br></div><div>Regards</div><div><br></div><div>Will McLaren</div><div>Ensembl Variation</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 31 July 2016 at 15:50, Wolf Beat <span dir="ltr"><<a href="mailto:Beat.Wolf@hefr.ch" target="_blank">Beat.Wolf@hefr.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Excellent question.<br>
<br>
I essentially expect the REST interface to return the same information as the website and the biomart API.<br>
Now, it could be that biomart and the website are both wrong in displaying the MAF number they display, but it should be coherent between all three ways to access.<br>
<br>
What i would expect personnaly is:<br>
As the 3 alleles are regrouped under the same rs number, i need to have the MAF of all 3 alleles combined. That what i would expect, as the rs number indicates that all 3 alleles map to it. I think thats a reasonably assumption?<br>
<br>
Kind regards<br>
<span class="HOEnZb"><font color="#888888"><br>
Beat Wolf<br>
________________________________________<br>
From: <a href="mailto:dev-bounces@ensembl.org">dev-bounces@ensembl.org</a> [<a href="mailto:dev-bounces@ensembl.org">dev-bounces@ensembl.org</a>] on behalf of Paul Flicek [<a href="mailto:flicek@ebi.ac.uk">flicek@ebi.ac.uk</a>]<br>
Sent: Sunday, July 31, 2016 3:38 PM<br>
To: Ensembl developers list<br>
Subject: Re: [ensembl-dev] Missing MAF numbers on VEP REST service<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
Hi Beat,<br>
<br>
Could you give me some insight about what you are looking for in the<br>
case of SNPs with more than two alleles.  For example, what do you<br>
expect the MAF to return in these cases?<br>
<br>
More specifically, MAF is not a property of a variant site itself.<br>
Frequency is a property of the various alleles.  We are currently<br>
loading and returning the dbSNP MAF for each SNP, which is defined as<br>
the frequency of the second most common allele, but in reality, the SNP<br>
has several alleles and each of them has a frequency.  As more<br>
individuals are assayed, more SNPs will be found to have 3 or 4 alleles,<br>
so this is an issue that will grow.<br>
<br>
Thanks,<br>
Paul<br>
<br>
<br>
<br>
On 31 Jul 2016, at 8:39, Wolf Beat wrote:<br>
<br>
> Hello,<br>
><br>
> i recently migrated a lot of my code from the biomart interface to the<br>
> ensembl rest service. I'm mainly using grch37 as my source.<br>
> I made the following observation:<br>
> When i query a multi allelic variant, the REST VEP service returns not<br>
> MAF numbers.<br>
> But when i query the biomart service or check the variant on the<br>
> website manually, the MAF numbers are present.<br>
> Just to give one example:<br>
> Variant rs6687605 has an MAF number of 0.43 when looked at on the<br>
> website and biomart:<br>
> <a href="http://grch37.ensembl.org/Homo_sapiens/Variation/Explore?r=1:25889132-25890132;v=rs6687605;vdb=variation;vf=106422589" rel="noreferrer" target="_blank">http://grch37.ensembl.org/Homo_sapiens/Variation/Explore?r=1:25889132-25890132;v=rs6687605;vdb=variation;vf=106422589</a><br>
><br>
> But the REST service returns no MAF number:<br>
> <a href="http://grch37.rest.ensembl.org/vep/hsapiens/id/rs6687605?content-type=application/json" rel="noreferrer" target="_blank">http://grch37.rest.ensembl.org/vep/hsapiens/id/rs6687605?content-type=application/json</a><br>
><br>
> I observed the same behaviour for:<br>
> rs2229546, rs6022, rs1058111<br>
><br>
> I also observed that the VEP rest service specifies:<br>
> <a href="http://grch37.rest.ensembl.org/documentation/info/vep_id_get" rel="noreferrer" target="_blank">http://grch37.rest.ensembl.org/documentation/info/vep_id_get</a><br>
> That the ExAC information is by default off. But i do recieve it<br>
> without specifying the option.<br>
> In my case i actually want it, but its still interesting to note.<br>
><br>
> The first problem is kind of a big issue for me, any help in solving<br>
> it would be greatly appreciated.<br>
><br>
> Kind regards<br>
><br>
> Beat Wolf<br>
><br>
> _______________________________________________<br>
> Dev mailing list    <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
> Posting guidelines and subscribe/unsubscribe info:<br>
> <a href="http://lists.ensembl.org/mailman/listinfo/dev" rel="noreferrer" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
> Ensembl Blog: <a href="http://www.ensembl.info/" rel="noreferrer" target="_blank">http://www.ensembl.info/</a><br>
<br>
_______________________________________________<br>
Dev mailing list    <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" rel="noreferrer" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" rel="noreferrer" target="_blank">http://www.ensembl.info/</a><br>
<br>
_______________________________________________<br>
Dev mailing list    <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" rel="noreferrer" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" rel="noreferrer" target="_blank">http://www.ensembl.info/</a><br>
</div></div></blockquote></div><br></div>