I am using Variant effect predictor to enrich my VarScan results which in turn are a results of tumor-normal variant calling in real human cancer patient data. One goal is to distinguish mutations, that overlap known protein domains and therefore we use --domains argument in VEP. However, so far we haven't been able to get any protein domain information using VEP, even though we have a test data with more than 25 000 missense mutations in known genes. These mutations are correctly annotated with PolyPhen and SIFT (using VEP with corresponding arguments). I think it is not possible to have such a large amount of mutations without a single one of them overlapping protein domains. I'm using VEP in offline mode and I have downloaded the large data sets including PolyPhen and SIFT data. Can you think of any reasons why VEP is unable to find overlapping protein domains?
0 answers
No answers yet.
Log in to answer this question.
-
Variant effect prediction
written by wrab425 •
-
how to use VEP and DrGAP tool
written by huynguyen96.dnu •
-
VEP output and stop gain mutation
written by seta
-
Large number of overlapped genes/transcripts reported by VEP
written by newbio17
-
How can I get other parameters from Ensembl VEP for vcf annotation
written by mittu1602
-
How to filter mutant allele frequency in Variant Effect Predictor (VEP)?
written by kin182 •
-
VEP vs webserver Polyphen scores
written by regrant •
-
Ensembl VEP filtering most_severe keeping Symbol, Feature,Consequence etc.
written by t2
-
Analysis Of Low-Complexity Regions
written by Khader Shameer
-
Missense Mutations Predictions Using Sift And Polyphen-2
written by Kamila001
You are right -- this is a bit strange. I would check few cases and confirm that the variants are on the domains and not on unassigned regions and then contact VEP.
Thanks for your reply. Could you please suggest a procedure for finding variants that occur on protein domains manually. We have obviously large amounts of data and we suppose that most of the variation occurs outside domains.