From: Heng Li Date: Mon, 8 Nov 2010 19:28:27 +0000 (+0000) Subject: fixed a silly bug in the indel caller X-Git-Url: https://git.donarmstrong.com/?p=samtools.git;a=commitdiff_plain;h=0917f66462dc2a0e7be6a33cb4f3e969e3bb1c9b fixed a silly bug in the indel caller --- diff --git a/bam2bcf.c b/bam2bcf.c index be8cf1f..c2d9944 100644 --- a/bam2bcf.c +++ b/bam2bcf.c @@ -183,7 +183,7 @@ int bcf_call2bcf(int tid, int pos, bcf_call_t *bc, bcf1_t *b, bcf_callret1_t *bc if (bc->a[i] < 0) break; if (i > 1) kputc(',', &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];