* samtools-0.1.9-17 (r824)

[samtools.git] / bam_md.c

diff --git a/bam_md.c b/bam_md.c
index 47ec97c124b3fae015bb126b3d99c936308daa00..29b315f4e9235f4e83516d6fb544821185d19a82 100644 (file)
--- a/bam_md.c
+++ b/bam_md.c
@@ -220,7 +220,7 @@ int bam_prob_realn_core(bam1_t *b, const char *ref, int write_bq)
                        else if (op == BAM_CDEL) x += l;
                }
                if (write_bq) {
-                       for (i = 0; i < c->l_qseq; ++i) bq[i] = bq[i] - qual[i] + 33;
+                       for (i = 0; i < c->l_qseq; ++i) bq[i] = bq[i] - qual[i] + 64;
                        bam_aux_append(b, "BQ", 'Z', c->l_qseq + 1, bq);
                        free(bq);
                }