bca->e = errmod_init(1. - theta);
bca->min_frac = 0.002;
bca->min_support = 1;
+ bca->per_sample_flt = 0;
return bca;
}
bca->bases = (uint16_t*)realloc(bca->bases, 2 * bca->max_bases);
}
// fill the bases array
- memset(r, 0, sizeof(bcf_callret1_t));
for (i = n = r->n_supp = 0; i < _n; ++i) {
const bam_pileup1_t *p = pl + i;
int q, b, mapQ, baseQ, is_diff, min_dist, seqQ;