fixed a silly bug in the indel caller

[samtools.git] / bam2bcf.c

diff --git a/bam2bcf.c b/bam2bcf.c
index be8cf1fe529a735800947b80b2a52c7fbaadb2d4..c2d994479a9094ea330291db34290e15faca8850 100644 (file)
--- 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];