SNP calling for Haploid Genomes

mikawinmikawin Member
edited October 2012 in Ask the GATK team


I used both haploid version (-ploidy flag set to 1) and the default version of UG to call SNPs between a reference assembly and sequence reads that come from the same haploid genome.

In the former case I get about 125 alt. hom calls and in the diploid mode of calling about 1,140 het and 120 alt. hom SNPs. Since the comparison is being done within haplotypes, I don’t expect any variation at all. What could be the basis for noise that is observed?

Thank you for your input!



Sign In or Register to comment.