Hi, Everybody,
I am trying to add an annotation column to my vcf file, after calling variants with the Samtools 1.x pipeline. I am using bcftools annotate, but I keep getting the same error regarding one of the headers: Could not parse the header line: "##FEATURE=<web_id=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100">" Could not parse bcftools_annots.hdr: ##FEATURE=<web_id=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100"> I am not sure what is wrong with that line; this is the content of the file I am calling with the -h parameter: ##FEATURE=<ID=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100"> ##FEATURE=<Name=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100"> ##FEATURE=<description=STRING_TAG,Number=1,Type=STRING,Description="erythrocyte+membrane+protein+1%2C+PfEMP1+%28VAR%29"> ##FEATURE=<size=NUMERIC_TAG,Number=2,Type=Integer,Description="7617"> ##FEATURE=<web_id=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100"> ##FEATURE=<locus_tag=STRING_TAG,Number=1,Type=STRING,Description="PF3D7_0100100"> ##FEATURE=<Alias=STRING_TAG,Number=1,Type=STRING,Description="VAR-UPSB1,124505645"> and this is an example of the 'FEATURE' field in my annotation tab-delimited file: ID=PF3D7_0100100;Name=PF3D7_0100100;description=erythrocyte+membrane+protein+1%2C+PfEMP1+%28VAR%29;size=7617;web_id=PF3D7_0100100;locus_tag=PF3D7_0100100;size=7617;Alias=VAR-UPSB1,124505645,MAL1P4.01,VAR,PF3D7_0100100,7670005,PFA0005w Have you got any idea on what I am doing wrong? Thanks for your kind attention, Massimiliano S. Tagliamonte Graduate Student University of Florida College of Veterinary Medicine Department of Infectious Diseases and Pathology
------------------------------------------------------------------------------
_______________________________________________ Samtools-help mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/samtools-help
