how to capture vcf's INFO field with bcftools
1
0
Entering edit mode
5.0 years ago
genya35 ▴ 40

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

software error • 3.2k views
ADD COMMENT
1
Entering edit mode
5.0 years ago

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

grep INFO your.vcf | grep TAG

##INFO=<ID=TAG,Number=...
ADD COMMENT
0
Entering edit mode

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;

ADD REPLY

Login before adding your answer.

Traffic: 1814 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6