This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Annotating a VCF [FORMAT not INFO field]

Hi,

I'm annotating a vcf and have done this successfully in the past but I'm curious whether I can add a custom annotation in the FORMAT field rather than the INFO field.

cat file.vcf | vcf-annotate -a annotation.gz -d key=FORMAT,ID=anno,Number=.,Type=String,Description='custom annotation'  -c CHROM,FROM,TO,FORMAT/anno >  new.vcf

That is, where the GT:PL:GL information is usually.

Or can't this be done? The above code does not work. Though does work when you swap FORMAT with INFO.

Thanks!

vcf vcftools annotation

The FORMAT field describes the tags for each value in each genotype column. If your annotation file is going to add more tags, how are you assigning values per sample in the genotype columns? So whatever you're trying to do, can't be done with vcf-annotate. You're most likely going to need vcf-merge. Can you describe your overall goal, so we can help troubleshoot?

0 answers

No answers yet.

Log in to answer this question.