Bug Bulletin: The recent 3.2 release fixes many issues. If you run into a problem, please try the latest version before posting a bug report, as your problem may already have been solved.

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!

Vivek

Post edited by Geraldine_VdAuwera on

Best Answer

  • ebanksebanks Posts: 679 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.