struct __mc_aux_t;
typedef struct __mc_aux_t mc_aux_t;
+#define MC_PTYPE_FULL 1
+#define MC_PTYPE_COND2 2
+
#ifdef __cplusplus
extern "C" {
#endif
mc_aux_t *mc_init(int n);
+ void mc_init_prior(mc_aux_t *ma, int type, double theta);
void mc_destroy(mc_aux_t *ma);
double mc_freq0(int ref, int *n, const bam_pileup1_t **plp, mc_aux_t *ma, int *_ref, int *alt);
double mc_freq_iter(double f0, mc_aux_t *ma);