This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Better splitting of mutliallelic sites then bcftools norm --m-any

Hey everyone,

I am trying to split the multiallelic sites of my VCF. I used bcftools norm --m-any. However, the result is not really reasonable to me. Here's an example.

Let's say, I have this multiallelic site:

REF     ALT     GT1     GT2     GT3
A       C,G     1/2     0/2     0/1

After splitting I get these two:

REF     ALT     GT1     GT2     GT3
A       C       1/0     0/0     0/1
A       G       0/1     0/1     0/0

So, the results for the "unused" ALT allele for a specific row is just set to REF. Is there a way to change this behavior, since I don't think it's reasonable to do it this way, at least for my analysis. I would like my result to be more like this:

REF     ALT     GT1     GT2     GT3          GT1     GT2     GT3
A       C       1/.     0/.     0/1    or    ./.     ./.     0/1
A       G       ./1     0/1     0/.          ./.     0/1     ./.

Or similar. At least I don't want to have REF where there was an ALT before.

bcftools
• 2,243 views • link
written by gernophil • gernophil

Because this question just received an upvote: I used vt decompose -s back in the days, but I think bcftools now also makes you choose the behavior you need.

• 0 views • link
written by gernophil • gernophil

1 answer

There isn’t a standard way to handle half-missing genotype calls. Often, conversion of half-missing to fully-missing genotypes is the best that can be done.

As a consequence, your suggested split isn’t really information-preserving. The bcftools split operation is analytically problematic, but at least it’s reversible.

• 0 views • link
written by chrchang523 chrchang523

Log in to answer this question.

More posts like this