It looks like you're new here. If you want to get involved, click one of these buttons!
kumar35885
Posts: 2Member ✭
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
ebanks
Posts: 483 mod
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.
