|
From: Rachel H. <ha...@so...> - 2011-01-25 00:19:30
|
I am using the vcf-annotate tool and I am using a command similar to the
one suggested on the VCF tools project website:
cat in.vcf | vcf-annotate -a annotations.gz \
-d key=INFO,ID=ANN,Number=1,Type=Integer,Description='My custom
annotation' > out.vcf
I get this error:
No ID in description?
[key=INFO,ID=ANN,Number=1,Type=String,Description=My custom annotation]
at ./vcf-annotate line 20
main::error('No ID in description?
[key=INFO,ID=ANN,Number=1,Type=String,D...') called at ./vcf-annotate
line 66
main::annotate('HASH(0x76e4220)') called at ./vcf-annotate line 11
When I looked at the vcf-annotate perl script, it looks like it is
parsing the description string given to the -d option splitting on ";"
but the items in this string are comma-separated. If I change the code
to split on commas, I still get errors. How can I get this to work?
Thanks.
Rachel
|