This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to Turn Multiple ALT alleles to Genotype Calls

I have a VCF file and many lines have multiple ALT alleles such as:

1       104159  .       CA      C,TA,TTT

I have an algorithm that takes genotype data in a 0,1,2 matrix form. I am wondering if there are any standards or best practices in how to turn VCF lines with multiple ALT alleles (0/0,0/1,0/2,1/2,2/3) into 0,1,2 genotype form.

genotype
• 3,421 views • link
written by jimkozubek • jimkozubek

did not work xd Can you share the solution if you have?

If I use the command below to divide them into biallelic, it does not work correctly.

Another solution is maybe to get rid of them but I do not lean towards.. Would be appreciated if you share the solution. Thanks!

• 0 views • link
written by obidobi • obidobi

1 answer

Yes, you can split multi-allelic calls with 

bcftools norm -Ov -m-any MyVariants.vcf > MyVariantsSplit.vcf ;

A useful addition is to also set the variants such that the REF allele matches that of a chosen reference genome, for example:

bcftools norm -Ov -m-any MyVariants.vcf | bcftools norm -Ov -f human_g1k_v37.fasta > MyVariantsSplitRefChecked.vcf ;

Kevin

• 1 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this