Make a 2-state allele HMM - no transitions to netural state
get_allele_hmm_s2(pAD, DP, R, p_s, theta, gamma = 20, r = 0.015)
HMM object
integer vector Paternal allele counts
integer vector Total alelle counts
numeric vector Phase switch probabilities
numeric Haplotype imbalance
numeric Overdispersion in the allele-specific expression