hap_obs_prob.RdThis will ultimately not be exported. It gives the probability that H2 is observed from a haplotype that is truly H1 given a super simple model of independent sequencing errors. This is not intended to be vectorized, at the moment.
hap_obs_prob(H1, H2, e)the DNA bases at the SNPs in the first haplotype
the DNA bases at the SNPs in the second haplotype
per snp error rates