This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Different total of SNPs and Indels after selection

Hi, I used HaplotypeCaller to call Indels and SNPs on my 30 samples and I separated the SNPs and Indels into 2 different files using SelectVariant tool.
I got 1954 SNPs and Indels in combination (in raw files), but after selection, I got only 277 indels and 1666 SNPs (in total only 1943, different from its raw file).
Why is it happening?
Many thanks before.

haplotypecaller gatk

Could you provide the commands that you used to split the variants and to count them?

Hi, Pawel.

I used the SelectVariants tool (-T SelectVariants) and -selectType INDEL for INDEL calling or -selectType SNP for SNP calling. But, fortunately I have found the answer, the 11 missing samples were mixed Indel/SNP sites, so the tool could not determine and split them. Thanks anyway.

0 answers

No answers yet.

Log in to answer this question.