Merge VCF INFO tags
0
0
Entering edit mode
21 months ago
FGV ▴ 130

Dear all,

I have a VCF file that, after annotation from VEP and extraction of tags with bcftools +split-vep ends up like:

[...]|amr|4.05357e-01||0|-0.143||||||;gnomADg_AF=0.341347,0.341347;[...]

Is there any way to merge the gnomADg_AF tag? For example, choose the MIN() or MAX(), so that it ends up with a single number on that tag:

[...]|amr|4.05357e-01||0|-0.143||||||;gnomADg_AF=0.341347;[...]

thanks,

vcf tags info • 457 views
ADD COMMENT
0
Entering edit mode

My recommendation would be to do this in R or Python. VCFs are tabular, so you can easily manipulate the INFO field with strsplit(..., split = "|") or .split('|'), then apply whatever logic you want for the gnomAD annotation.

ADD REPLY

Login before adding your answer.

Traffic: 2616 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