This function is a work in progress. As such, it is not exported yet; but it allows for some leniency on the sequencing error parameter, hence "SESlack"

pmf_betabinomMix_SEslack(
  ref_counts,
  var_counts,
  probshift,
  piRR,
  piVV,
  pr,
  pv,
  prv,
  theta_RR,
  theta_VV,
  theta_het
)