Annovar using VCF file
1
0
Entering edit mode
6.9 years ago
Samar • 0

Hi all, I have a VCF file that I want to annotate using Annovar. I tried using convert2annovar.pl command but since I have a multi-sample file, it generates separate avinput files. I want one single file with all the sample information.

So I now used the command:

table_annovar.pl \
  /archive/samarpana/2_Unified_Genotyper/Family2/Family2.vcf humandb/ \
  -buildver hg19 \
  myanno \
  -remove \
  -protocol refGene,phastConsElements46way,esp6500siv2_all,1000g2015aug_all,avsnp144,exac03,ljb26_all,avsift \
  -operation g,r,f,f,f,f,f,f \
  -nastring .\
   -csvout

and also

table_annovar.pl \
  /archive/samarpana/2_Unified_Genotyper/Family2/Family2.vcf humandb/ \
  -buildver hg19 \
  myanno \
  -remove \
  -protocol refGene,phastConsElements46way,esp6500siv2_all,1000g2015aug_all,avsnp144,exac03,ljb26_all,avsift \
  -operation g,r,f,f,f,f,f,f \
  -nastring NA \
  -vcfinput

But none of the commands are working. I am getting error message saying Syntax error. Can anyone suggest what I am doing wrong?

Thanks

software-error next-gen ANNOVAR • 2.6k views
ADD COMMENT
0
Entering edit mode

Can you try

table_annovar.pl /archive/samarpana/2_Unified_Genotyper/Family2/Family2.vcf humandb/ -buildver hg19 -out myanno -remove -protocol refGene,phastConsElements46way,esp6500siv2_all,1000g2015aug_all,avsnp144,exac03,ljb26_all,avsift -operation g,r,f,f,f,f,f,f -nastring NA -vcfinput
#___________________________________________________________________________________________________^^^^
ADD REPLY
1
Entering edit mode
6.7 years ago

As @Sam pointed out, you're missing the -out before myanno. The -vcfinput argument (which you had in the second command but not your first) is also necessary, and omitting it wouldn't cause the "Syntax error" message, but instead all your variants would get discarded (and written to myanno.invalid_input)

I think your main problem is that Annovar complains when -nastring is not '.' when using -vcfinput. This will give you an Syntax Error message like you described.

table_annovar.pl /archive/samarpana/2_Unified_Genotyper/Family2/Family2.vcf humandb/ \
    -buildver hg19 \
    -out myanno \
    -remove \
    -protocol refGene,phastConsElements46way,esp6500siv2_all,1000g2015aug_all,avsnp144,exac03,ljb26_all,avsift \
    -operation g,r,f,f,f,f,f,f \
    -nastring . \
    -vcfinput

The above line should work.

ADD COMMENT

Login before adding your answer.

Traffic: 1361 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