This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is it possible to add variant type "VT" to each variant in vcf?

I tried the following line:

bcftools +fill-tags test.vcf -Ov -o testVT.vcf -- -t VT

But get an error

Error parsing "--tags VT": the tag "VT" is not supported

Is it possible to add variant type annotation to vcf? How to do that? My INFO field contains things like VQSRTrancheSNP99.80to100.00 or VQSRTrancheINDEL99.00to100.00 which indicates the variant type (SNP and INDEL), maybe it's possible to extract it?

Thanks guys.

bcftools
• 1,939 views • link
updated by Ram Ram • written by optimistsso4co3 • optimistsso4co3

You should be able to use bcftools query to extract the VQSRTranche... field (along with CHR POS REF and ALT) and then use bcftools annotate to add a modified version of the extracted field as a new annotation.

bcftools might also have an intrinsic SNP/INDEL annotation protocol, plus using REF and ALT lengths might be an option too.

• 0 views • link
written by Ram Ram

0 answers

No answers yet.

Log in to answer this question.

More posts like this