Allow p to be varied in error model used to calculate phread likelihoods in genotypes.
Allow beta bionomial to be used in place of bionomial for error model.
Initial set of errmod tests for my sanity