This is a test version of Biostars. For the public version, visit https://www.biostars.org.
ANNOVAR annotation

Hi,

I'm trying to annotate VCF file:

enter image description here

I run the following command

table_annovar.pl \
  --vcfinput eczema.vcf \
  humandb/ \
  -buildver hg19 \
  -out ecze \
  -remove \
  -protocol refGene,cytoBand,exac03,avsnp147,dbnsfp30a \
  -operation g,r,f,f,f \
  -nastring . \
  -polish

but I get the following error :

NOTICE: Running with system command <convert2annovar.pl  -includeinfo -allsample -withfreq -format vcf4 eczema.vcf > ecze.avinput >
Error: invalid record in VCF file: the GT specifier is not present in the FORMAT string: < 1       754182       .       A       G   .> 
Error running system command: <convert2annovar.pl  -includeinfo -allsample -withfreq -format vcf4 eczema.vcf > ecze.avinput>

I can't seem to understand what is the error would be happy for help

vcf annovar
• 1,140 views • link
updated by Ram Ram • written by Eliza • Eliza

1 answer

: the GT specifier is not present in the FORMAT string

The is a FORMAT column. When this column exists, a VCF is expected to have some other columns for the genotypes.

The following column could fix the problem.

cut -f1-8 < eczema.vcf > new.vcf

and test with 'new.vcf'

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this