This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Duplicate variants in same position: bcftools

I create a vcf file by bcftools mpileup command. But I got the huge number of duplicate variants. Duplicate means it shows different variants in the same position. How can I handle this situation?

Example:

chr1 3993567 T A
chr1 3993567 tc t
bcftools vcf
• 5,003 views • link
updated by sami • sami • written by Mbillah Mbillah

Could you please provide the exact command you have used and some complete vcf lines?

• 1 views • link
written by finswimmer finswimmer

POS-----VARIANTS---------REF---------------------Allele Variations

99162-----DEL/MNP---TCGGTGTGCGCGG-----TCGG

99166-----SNP-------------------T-------------------------------C

I have on my VCF: DELETION in position 99162 and SNP in position 99166.

It is impossible to have deletion in 99162 and snp in 99166 in the same time. How can i explain that ?? Should i keep snp or deletion in my vcf? Please help me.

Thank you for your answer.

• 0 views • link
written by sami • sami

You've already asked a question. Please do not spam other threads - this behavior will lead to your account being suspended.

Plus, this is not an answer. I'm moving it to a comment for now, but will be deleting any similar posts in the future.

• 0 views • link
written by Ram Ram

@RamRS don't let me deactivate my account and leave this group forever !

• 0 views • link
written by sami • sami

Please stop adding answers. Use Add Comment or Add Reply as appropriate instead.

• 0 views • link
written by Ram Ram

1 answer

How can I handle this situation?

try to use bcftools norm

https://samtools.github.io/bcftools/bcftools.html#norm (...) recover multiallelics from multiple rows.

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this