This is a test version of Biostars. For the public version, visit https://www.biostars.org.
does bcftools has anything like vcf-check?

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
• 4,948 views • link
written by roalva1 roalva1

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

• 1 views • link
written by roalva1 roalva1

1 answer

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
• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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

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

and I get an empty file

• 1 views • link
written by roalva1 roalva1

Log in to answer this question.

More posts like this