if (theta <= 0.) theta = CALL_DEFTHETA;
bca = calloc(1, sizeof(bcf_callaux_t));
bca->capQ = 60;
- bca->openQ = 40; bca->extQ = 20; bca->tandemQ = 80;
+ bca->openQ = 40; bca->extQ = 20; bca->tandemQ = 100;
bca->min_baseQ = min_baseQ;
bca->e = errmod_init(1. - theta);
return bca;
kputc(ref[pos], &s);
for (i = 1; i < 4; ++i) {
if (bc->a[i] < 0) break;
- if (i > 1) kputc(',', &s);
+ if (i > 1) {
+ kputc(',', &s); kputc(ref[pos], &s);
+ }
if (bca->indel_types[bc->a[i]] < 0) { // deletion
- for (j = -bca->indel_types[bc->a[i]]; j < bca->indelreg; ++i)
+ for (j = -bca->indel_types[bc->a[i]]; j < bca->indelreg; ++j)
kputc(ref[pos+1+j], &s);
} else { // insertion; cannot be a reference unless a bug
char *inscns = &bca->inscns[bc->a[i] * bca->maxins];