does bcftools has anything like vcf-check?
1
1
Entering edit mode
6.5 years ago
roalva1 ▴ 90

Hello,

I am merging 14k vcfs and I get this error:

[W::vcf_parse] INFO 'D' is not defined in the header, assuming Type=String
Error: The INFO field is not defined in the header: D

So, I think I have some vcf that is not well formatted for bcftools, would you know if bcftools has something like vcf-check? I haven't seen it. Then, I could remove this vcf.

Thanks

bcftools vcf • 3.9k views
ADD COMMENT
0
Entering edit mode

I've just discovered a truncated file. Sorry!!!!

ADD REPLY
2
Entering edit mode
6.5 years ago

here is a command to find the VCF missing the INFO/D in the header

$ find dir -name "*.vcf.gz" | while read F; do gunzip -c "$F" | grep "^##INFO=<ID=D," > /dev/null || echo $F ; done
ADD COMMENT
0
Entering edit mode

I think there is a problem with bcftools, becuase I do this:

zgrep   ";D=" ./*.vcf.gz > D.txt

and I get an empty file

ADD REPLY

Login before adding your answer.

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