- if (conf->vcf) {
- double f0, f, fpost, pref = -1.0; // the reference allele frequency
- int j, _ref, _alt, _ref0, depth, rms_q, _ref0b, is_var = 0, qref = 0;
+ if (conf->flag & MPLP_GLF) {
+ int _ref0, ref16;
+ bcf1_t *b = calloc(1, sizeof(bcf1_t));
+ _ref0 = (ref && pos < ref_len)? ref[pos] : 'N';
+ ref16 = bam_nt16_table[_ref0];
+ for (i = 0; i < n; ++i)
+ bcf_call_glfgen(n_plp[i], plp[i], ref16, bca, bcr + i);
+ bcf_call_combine(n, bcr, ref16, &bc);
+ bcf_call2bcf(tid, pos, &bc, b);
+ bcf_write(bp, bh, b);
+ //fprintf(stderr, "%d,%d,%d\n", b->tid, b->pos, b->l_str);
+ bcf_destroy(b);
+ } else if (conf->flag & MPLP_VCF) {
+ mc_rst_t r;
+ int j, _ref0, depth, rms_q, _ref0b, is_var = 0, qref = 0, level = 2, tot;