<div dir="ltr">Hi Alexey,<div><br></div><div>I'm afraid it's not easily possible to do as you request, and you're right, in cases where the 1000 genomes alt allele is actually the most common allele this can be misleading.</div>
<div><br></div><div>It would, however, be possible to implement a plugin to do the switching process - if you are reasonably confident in Perl it shouldn't be too big a task.</div><div><br></div><div>It's something we'll look at incorporating into the core VEP code in future.</div>
<div><br></div><div>Regards</div><div><br></div><div>Will McLaren</div><div>Ensembl Variation</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 28 July 2014 20:45, Alexey Larionov <span dir="ltr"><<a href="mailto:alexey_larionov@hotmail.com" target="_blank">alexey_larionov@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-GB" link="blue" vlink="purple"><div><p class="MsoNormal">VEP provides “co-localised MAF” for the detected variants.  This may not be sufficient to estimate the 1k frequency of the detected variant, for instance if there are multiple alleles at this location.  <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Furthermore, VEP provides “--filter_common” option to exclude variants co-located with existing common allele (MAF > 1%).  This is commonly interpreted as excluding the common variants.  Unfortunately the above interpretation is completely wrong for many hundreds of variants, where the minor (very rare) allele is incorporated in reference genome (e.g. rs6672356, rs4274008, rs7545802, rs198400, rs1763642, rs6429745 etc).  On my experience, this this “co-localised MAF” clause is very confusing for many biologists, who merely want excluding the common variants.  <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal">Is it possible to use VEP for retrieving 1k allele frequency for the DETECTED variant, in addition/instead to the frequency co-localised minor allele?  <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><u></u> <u></u></p></div></div><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" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
<br></blockquote></div><br></div>