Record actual sequence in unpad_seq (fixes apparently harmless bug)

diff --git a/padding.c b/padding.c
index 6c8b31ff3de5d01afa79bb94428bd4e711078116..2b432756d29805b4fe102d376d8ac1810eb38c20 100644 (file)
--- a/padding.c
+++ b/padding.c
@@ -38,8 +38,7 @@ static void unpad_seq(bam1_t *b, kstring_t *s)
                op = bam_cigar_op(cigar[k]);
                ol = bam_cigar_oplen(cigar[k]);
                if (op == BAM_CMATCH || op == BAM_CEQUAL || op == BAM_CDIFF) {
-                       for (i = 0; i < ol; ++i) s->s[s->l++] = bam1_seqi(seq, j);
-                       ++j;
+                       for (i = 0; i < ol; ++i, ++j) s->s[s->l++] = bam1_seqi(seq, j);
                } else if (op == BAM_CSOFT_CLIP) {
                        j += ol;
                } else if (op == BAM_CDEL || op == BAM_CPAD) {