Add utility function to count heterozygotes