Make a 3-state allele HMM - allow transitions to netural state
get_allele_hmm_s3(pAD, DP, R, p_s, t, theta_min, gamma = 20, r = 0.015)
HMM object
integer vector Paternal allele counts
integer vector Total alelle counts
numeric vector Phase switch probabilities
numeric Transition probability between haplotype states
numeric Haplotype imbalance
numeric Overdispersion in the allele-specific expression
numeric Variant mapping bias