How To Prepare Gatk .Vcf Files Output For Annovar ?
2
1
Entering edit mode
8.9 years ago
newDNASeqer ▴ 730

I have managed to generate vcf files by using bwa + PiCard + GATK pipeline, and would like to use ANNOVAR to do the variant calling. I checked the vcf content output from GATK, and think it is not ready for feeding into ANNOVAR. I was wondering if there's any handy tool to prepare the GATK vcf files for ANNOVAR? thanks for reply

variant calling gatk annovar vcf • 5.4k views
ADD COMMENT
3
Entering edit mode
8.9 years ago

From the ANNOVAR manual:

perl convert2annovar.pl infile.vcf -format vcf4 > outfile

As an aside, I think you are slightly confused when you say use ANNOVAR to do "variant calling". What you have in your vcf file are variants that were called using GATK. ANNOVAR does ANNOtation of VARiants by either locating interrupted or closest genomic features (genes, transcription factors, etc.) or filtering the variants based on another set of variant genotypes you supply.

ADD COMMENT
0
Entering edit mode

Thank you for your answer, appreciate it.

ADD REPLY
0
Entering edit mode

Also, is there any way to convert the Annovar file (containing INDELs info) to vcf file (using hg19 as ref genome).

ADD REPLY
0
Entering edit mode
5 months ago

For those who find this 8 years later... table_annovar.pl takes vcf input directly with the --vcfinput argument.

ADD COMMENT

Login before adding your answer.

Traffic: 1379 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6