US Holiday notice: this Thursday and Friday (Nov 25-26) the forum will be unattended. Normal service will resume Monday Nov 29. Happy Thanksgiving!

GATK UGT -heterozygosity / -hets

kumar35885kumar35885 Posts: 2Member
edited October 2012 in Ask the GATK team

I have an inbred mouse strain that I am sequencing and there should be little to NO heterozygosity. Yet with the default settings of UGT -heterozygosity (which is 0.001) many homs are being called as hets. When 230/250 reads are alternate and 20/250 are reference, it calls a het, even though it should be homozygous alternate.

What do you recommendations for this setting for inbred animals?

thanks, GATK is great!


Post edited by Geraldine_VdAuwera on

Best Answer

  • ebanksebanks Posts: 684GATK Developer mod
    Answer ✓

    The heterozygosity is the probability that the allele for a single chromosome is non-reference (which in humans is 1/1000 for SNPs). So even for an "inbred human" you would still expect a variant on a given chromosome once every 1000 bases. It is not the probability that a diploid sample is heterozygous at a given position.

Sign In or Register to comment.