Is it possible to add variant type "VT" to each variant in vcf?
0
0
Entering edit mode
23 months ago

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 • 604 views
ADD COMMENT
1
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

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