- 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);
- double mc_freq_post(mc_aux_t *ma);
- double mc_ref_prob(mc_aux_t *ma);
+ int mc_cal(int ref, int *n, const bam_pileup1_t **plp, mc_aux_t *ma, mc_rst_t *rst, int level);