New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
ExAC plugin returning only ExAC_AF, not population-specific AF #8
Comments
I can't recreate this error:
|
I have the same Issues using VEP 81. I have seen that the problem occurs if the --vcf VEP option is chosen. perl ~/vep/variant_effect_predictor.pl \
X 1544525 . G A . . CSQ=A|missense_variant|MODERATE|ASMTL|ENSG00000169093|Transcript|ENST00000381317|protein_coding|8/13||ENST00000381317.3:c.953C>T|ENSP00000370718.3:p.Pro318Leu|986/2027|953/1866|318/621|P/L|cCc/cTc||1||-1|SNV|HGNC|751|YES||CCDS43917.1|ENSP00000370718|ASML_HUMAN||UPI0000141AF9|NM_004192.3|tolerated(0.05)|benign(0.212)|PROSITE_profiles:PS51683&hmmpanther:PTHR11746&hmmpanther:PTHR11746:SF69&Gene3D:1.10.10.10&Superfamily_domains:SSF46785||||||||||||||||||||||0.000712||||,A|missense_variant|MODERATE|ASMTL|ENSG00000169093|Transcript|ENST00000381333|protein_coding|7/12||ENST00000381333.4:c.905C>T|ENSP00000370734.4:p.Pro302Leu|994/2035|905/1818|302/605|P/L|cCc/cTc||1||-1|SNV|HGNC|751|||CCDS55363.1|ENSP00000370734|ASML_HUMAN||UPI0000071671|NM_001173474.1|tolerated(0.06)|possibly_damaging(0.844)|PROSITE_profiles:PS51683&hmmpanther:PTHR11746&hmmpanther:PTHR11746:SF69&Gene3D:1.10.10.10&Superfamily_domains:SSF46785||||||||||||||||||||||0.000712||||,A|missense_variant|MODERATE|ASMTL|ENSG00000169093|Transcript|ENST00000534940|protein_coding|8/13||ENST00000534940.1:c.779C>T|ENSP00000446410.1:p.Pro260Leu|1005/2048|779/1692|260/563|P/L|cCc/cTc||1||-1|SNV|HGNC|751|||CCDS55362.1|ENSP00000446410|ASML_HUMAN||UPI0001CA7E5F|NM_001173473.1|tolerated(0.17)|benign(0.212)|Superfamily_domains:SSF46785&Gene3D:1.10.10.10&hmmpanther:PTHR11746:SF69&hmmpanther:PTHR11746&PROSITE_profiles:PS51683||||||||||||||||||||||0.000712||||,A|missense_variant|MODERATE|ASMTL|ENSG00000169093|Transcript|ENST00000416733|protein_coding|7/12||ENST00000416733.2:c.725C>T|ENSP00000410578.2:p.Pro242Leu|923/1966|725/1638|242/545|P/L|cCc/cTc||1||-1|SNV|HGNC|751||||ENSP00000410578||E7ER97_HUMAN|UPI0001AE6F52|||benign(0.15)|PROSITE_profiles:PS51683&hmmpanther:PTHR11746&hmmpanther:PTHR11746:SF69&Gene3D:1.10.10.10&Superfamily_domains:SSF46785||||||||||||||||||||||0.000712|||| GT:AD:DP 0/1:116,89:. 0/0:2654,7:. perl ~/vep/variant_effect_predictor.pl -plugin ExAC,/Users/alessandro/.vep/homo_sapiens/ExAC/ExAC.r0.3.sites.vep.vcf.gz --offline --gencode_basic --everything --total_length --allele_number --no_escape --check_existing --xref_refseq --shift_hgvs 1 --check_alleles --quiet --no_stats --assembly GRCh37 --dir $VEP_DATA --fasta $VEP_DATA/homo_sapiens/81_GRCh37/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa --input_file s_IM_pcnsl_011_T_vs_s_Normal_Pooled_FFPE_1.vcf --output_file test.vep.vcp X_1544525_G/A X:1544525 A ENSG00000169093 ENST00000381317 Transcript missense_variant 986/2027 953/1866 318/621 P/L cCc/cTc - ExAC_AF=0.000712;ExAC_AF_AMR=8.75e-05;ExAC_AF_OTH=0.00114;ExAC_AF_SAS=6.14e-05;ExAC_AF_FIN=0.000152;ExAC_AF_NFE=0.00124;ALLELE_NUM=1;ExAC_AF_AFR=0;ExAC_AF_EAS=0;IMPACT=MODERATE;STRAND=-1;VARIANT_CLASS=SNV;SYMBOL=ASMTL;SYMBOL_SOURCE=HGNC;HGNC_ID=751;BIOTYPE=protein_coding;CANONICAL=YES;CCDS=CCDS43917.1;ENSP=ENSP00000370718;SWISSPROT=ASML_HUMAN;UNIPARC=UPI0000141AF9;RefSeq=NM_004192.3;SIFT=tolerated(0.05);PolyPhen=benign(0.212);EXON=8/13;DOMAINS=PROSITE_profiles:PS51683,hmmpanther:PTHR11746,hmmpanther:PTHR11746:SF69,Gene3D:1.10.10.10,Superfamily_domains:SSF46785;HGVSc=ENST00000381317.3:c.953C>T;HGVSp=ENSP00000370718.3:p.Pro318Leu |
Hi Will I have update my ExAC.pm plugin, I unfortunately still have the same issues. This is my veep run output with and without the --vcf parameter. perl variant_effect_predictor.pl --plugin ExAC,/Users/alessandro/.vep/homo_sapiens/ExAC/ExAC.r0.3.sites.vep.vcf.gz --quiet --offline --no_stats --everything --shift_hgvs 1 --check_existing --check_alleles --total_length --allele_number --no_escape --xref_refseq --assembly GRCh37 --fasta /Users/alessandro/.vep/homo_sapiens/81_GRCh37/Homo_sapiens.GRCh37.75.dna.primary_assembly.fa --dir /Users/alessandro/.vep/ --input_file example_GRCh37.vcf --output_file example_GRCh37.vep.vcf --vcf |
Thanks for the report, c295b79 should fix it. You may also need to update your ensembl-tools folder (just re-download, or run "perl INSTALL.pl --UPDATE") |
I need to dig through the source, but I am seeing empty columns for everything except ExAC_AF. The source ExAC VCF is directly from the ExAC download page.
The text was updated successfully, but these errors were encountered: