[Vcftools-help] vcf-annotate example annotates wrong column

[Powell, E. <EPo...@me...>]

I am trying to replicate the vcf-annotate example.

My command line is:
-----------------------------------------------------------------------------
/mihg/users/epowell1/gatk/vcftools_0.1.4a/bin/vcf-annotate -a annotate.txt.gz -c FROM,TO,CHROM,INFO/HM2,INFO/GN,INFO/DP -d key=INFO,ID=HM2,Number=0,Type=Flag,Description=HapMap2 membership -d key=INFO,ID=GN,Number=1,Type=String,Description=Gene Name -d key=INFO,ID=DP,Number=0,Type=Integer,Description=Depth,etc /mihg/users/epowell1/gatk/vcftools_0.1.4a/examples/concat-a.vcf

The format of annotate.txt.gz is:
-----------------------------------------------------------------------------
[epowell1@mendel1 test_annotate_example]$ zcat annotate.txt.gz
100     100     1       HM2     gene1   5
110     110     1       0       .       6
130     130     1       HM2     .       7
140     140     1       0       .       8
150     150     1       0       .       9
110     150     2       HM2     gene2   .
160     160     2       0       gene3   11

Unfortunately, while I do not get any errors, the output is not correct:
-----------------------------------------------------------------------------
[epowell1@dfa  test_annotate_example]$ perl $vcfannotate -a $annot.gz     -c FROM,TO,CHROM,INFO/HM2,INFO/GN,INFO/DP     -d key=INFO,ID=HM2,Number=0,Type=Flag,Description="HapMap2 membership"     -d key=INFO,ID=GN,Number=1,Type=String,Description="Gene Name"     -d key=INFO,ID=DP,Number=0,Type=Integer,Description="Depth,etc" $invcf
##fileformat=VCFv4.0
##INFO=<ID=DP,Number=1,Type=Integer,Description="Total Depth">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">
##FORMAT=<ID=DP,Number=1,Type=Integer,Description="Read Depth">
##FILTER=<ID=q10,Description="Quality below 10">
##INFO=<ID=HM2,Number=0,Type=Flag,Description="HapMap2 membership">
##INFO=<ID=GN,Number=1,Type=String,Description="Gene Name">
##source_20110305.1=/mihg/users/epowell1/gatk/vcftools_0.1.4a/bin/vcf-annotate -a annotate.txt.gz -c FROM,TO,CHROM,INFO/HM2,INFO/GN,INFO/DP -d key=INFO,ID=HM2,Number=0,Type=Flag,Description=HapMap2 membership -d key=INFO,ID=GN,Number=1,Type=String,Description=Gene Name -d key=INFO,ID=DP,Number=0,Type=Integer,Description=Depth,etc /mihg/users/epowell1/gatk/vcftools_0.1.4a/examples/concat-a.vcf
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  A
;GN=gene1;HM2   GT:GQ:DPGTTT    0/1:409:3506    q10     DP=5
1       GT:GQ:DP.       0/1:245:32      1792    PASS    DP=6
1       120     .       GA      G       628     q10     DP=21   GT:GQ:DP        1/1:21:21
;HM2    GT:GQ:DP.       0/1:212:22      1016    PASS    DP=7
1       GT:GQ:DP.       0/1:150:30      727     PASS    DP=8
1       GT:GQ:DP.       1/2:12:10A,T    246     PASS    DP=9
1       160     .       TAAAA   TA,T    246     PASS    DP=10   GT:GQ:DP        1/2:12:10
2       100     .       GTTT    G       1806    q10     DP=35   GT:GQ:DP        0/1:409:35
;GN=gene2;HM2   GT:GQ:DPCAAA    0/1:245:3292    PASS    DP=.
;GN=gene2;HM2   GT:GQ:DPGA      1/1:21:2128     q10     DP=.
;GN=gene2;HM2   GT:GQ:DPGAA     0/1:212:2216    PASS    DP=.
;GN=gene2;HM2   GT:GQ:DPGT      0/1:150:307     PASS    DP=.
;GN=gene2;HM2   GT:GQ:DPTAAAA   1/2:12:1046     PASS    DP=.
;GN=gene360     GT:GQ:DPTAAAA   1/2:12:1046     PASS    DP=11