[ensembl-dev] VEP: How to get frequency of detected variant instead of the "co-localised MAF"?

Will McLaren wm2 at ebi.ac.uk
Tue Jul 29 13:11:21 BST 2014


Hi Alexey,

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.

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.

It's something we'll look at incorporating into the core VEP code in future.

Regards

Will McLaren
Ensembl Variation


On 28 July 2014 20:45, Alexey Larionov <alexey_larionov at hotmail.com> wrote:

> 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.
>
>
>
> 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.
>
>
>
> 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?
>
>
>
>
>
> _______________________________________________
> Dev mailing list    Dev at ensembl.org
> Posting guidelines and subscribe/unsubscribe info:
> http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20140729/c69522c2/attachment.html>


More information about the Dev mailing list