how to capture vcf's INFO field with bcftools

Question

[genya35]

Hello,

I'm trying to print the content of a vcf file but having problem with capturing the 'INFO' column.

Here is the code.

my $cmd2 = `bcftools query $vcf_indexed -f '%CHROM\t%POS\t%REF\t%ALT\t%INFO/TAG[\t%SAMPLE=%AD]\n' --samples '$sample'`;print $cmd2;

I get "Error: no such tag defined in the VCF header: INFO/TAG"

Thanks

Answer 1

[Pierre Lindenbaum]

it only works if TAG is declared in the VCF header:

grep INFO your.vcf | grep TAG

##INFO=<ID=TAG,Number=...

[genya35]

Thanks, sorry for a simple question but how do I redirect the output to a text file of my choice?

Sorry, i got it. my $cmd2 = bcftools query $vcf_indexed -f '%CHROM\t%POS\t%REF\t%ALT\t%INFO/ANN][\t%SAMPLE=%AD]\n' --samples '$sample' > $output_file;print $cmd2;