Question: VCF to ANNOVAR format
0
gravatar for reza.jabal
3.5 years ago by
reza.jabal320
United Kingdom
reza.jabal320 wrote:

Here, I've got a VCF file I want to convert to ANNOVAR format through running the following command, but apparently it is not working!

table_annovar.pl 1_23_intersect.vcf /scratch/GENE_DATABASES/humandb/HG19_APR13/ -protocol refGene,phastConsElements44way,genomicSuperDups,esp6500si_all,1000g2012apr_all,snp135,ljb_all -operation g,r,r,f,f,f,f -vcfinput .

I was wondering if anyone has any suggestion on how to sort this out.

next-gen software error • 3.1k views
ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by reza.jabal320

Could you solve the issue? Even I am facing the same problem. Everytime I use this table_annovar.pl, I get an error, saying its a syntax issue. 

ADD REPLYlink written 3.5 years ago by Samar0
1
./convert2annovar.pl -format vcf4 input.vcf > output.vcf
ADD REPLYlink written 19 months ago by Karma200
1
gravatar for ebrown1955
3.5 years ago by
ebrown1955290
United States
ebrown1955290 wrote:

The -vcfinput option is relatively new, and if it's not working then you might have an older version of Annovar. Try updating to the most current release.

If all else fails, use convert2annovar.pl, i.e.

convert2annovar.pl -format vcf4 1_23_intersect.vcf > 1_23_intersect.annovar

Then you can use table annovar in some variant of this:

table_annovar.pl 1_23_intersect.annovar /scratch/GENE_DATABASES/humandb/HG19_APR13/ -protocol refGene,phastConsElements44way,genomicSuperDups,esp6500si_all,1000g2012apr_all,snp135,ljb_all -operation g,r,r,f,f,f,f -csvout .
ADD COMMENTlink written 3.5 years ago by ebrown1955290
1
gravatar for cpad0112
3.5 years ago by
cpad011211k
India
cpad011211k wrote:

Down load the converter from here: http://www.openbioinformatics.org/annovar/download/convert2annovar.pl

Try running simple command like this on your vcf file:

perl convert2annovar.pl -format vcf4 ex2.vcf

You can download ex2.vcf from here: https://raw.githubusercontent.com/Bioconductor-mirror/VariantAnnotation/master/inst/extdata/ex2.vcf

This perl script seems to be working and output i got was:

WARNING to old ANNOVAR users: this program no longer does line-to-line conversion for multi-sample VCF files. If you want to include all variants in output, use '-format vcf4old' or use '-format vcf4 -allsample -withfreq' instead.
20    1110696    1110696    A    G    het    67    6
20    1110696    1110696    A    T    het    67    6
20    1234568    1234569    TC    -    het    50    4
NOTICE: Finished reading 27 lines from VCF file
NOTICE: A total of 5 locus in VCF file passed QC threshold, representing 5 SNPs (2 transitions and 2 transversions) and 2 indels/substitutions
NOTICE: Finished writing 2 SNP genotypes (1 transitions and 1 transversions) and 1 indels/substitutions for 1 sample (but input contains 3 samples)
WARNING: 1 invalid alternative alleles found in input file

Run your VCF file through VCF validating tools. Upgrade/update perl script tool as mentioned above. First try running with simple commands.

 

ADD COMMENTlink modified 3.5 years ago • written 3.5 years ago by cpad011211k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1162 users visited in the last hour