bcftools info and filter error
1
0
Entering edit mode
7 months ago

Hi. I am trying to do a comparative analysis of my vcf file against the vcf files of ExAc. I'm using the bcftool isec here. I am getting an error that says:

[W::vcf_parse_info] INFO 'HOM_CONSANGUINEOUS' is not defined in the header, assuming Type=String
[W::vcf_parse_filter] FILTER 'AC_Adj0_Filter' is not defined in the header

I'm not sure how to annotate the header and resolve this. [I'm an amateur in this with very rudimentary knowledge of programming]

ExAc bcftools vcf • 383 views
ADD COMMENT
0
Entering edit mode
7 months ago

use awk to insert a filter defintion in the VCF

awk '/^#CHROM/ {printf("##FILTER=<ID=HOM_CONSANGUINEOUS,Description=\"blabla\">\n##FILTER=<ID=AC_Adj0_Filter,Description=\"blabla\">\n");} {print}'  in.vcf
ADD COMMENT

Login before adding your answer.

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