bcftools icon indicating copy to clipboard operation
bcftools copied to clipboard

Bcftools isec GT aware

Open Rashesh7 opened this issue 2 years ago • 1 comments

Hello,

Is there a way to run bcftools isec with taking genotypes into consideration? I want to compare multiple VCF files and check even in the same variant call, how many have a shift from heterozygous to homozygous.

I am running "bcftools isec -p Unique -n-1 -c none $FILE1 $FILE2 $FILE3 ...." currently. But it just checks for the exact Alt allele. I know vcf-compare has an option for comparing genotypes "--cmp-genotypes ". Is there something similar that I can use?

Many Thanks, Rashesh

Rashesh7 avatar Aug 28 '22 12:08 Rashesh7

I wonder if the GTisec or GTsubset plugins would help. I did not write these so can't really tell.

pd3 avatar Sep 08 '22 14:09 pd3