* samtools-0.1.10-4 (r840)

[samtools.git] / bam2bcf_indel.c

diff --git a/bam2bcf_indel.c b/bam2bcf_indel.c
index 7ce6cd4294d8f0e2a83b78c3e0f7ac0d86e5602a..819c1e7b1371296b7da5536939fab691f53c1722 100644 (file)
--- a/bam2bcf_indel.c
+++ b/bam2bcf_indel.c
@@ -286,10 +286,10 @@ int bcf_call_gap_prep(int n, int *n_plp, bam_pileup1_t **plp, int pos, bcf_calla
                                        const uint8_t *qual = bam1_qual(p->b), *bq;
                                        uint8_t *qq;
                                        qq = calloc(qend - qbeg, 1);
-                                       bq = (uint8_t*)bam_aux_get(p->b, "BQ");
+                                       bq = (uint8_t*)bam_aux_get(p->b, "ZQ");
                                        if (bq) ++bq; // skip type
                                        for (l = qbeg; l < qend; ++l) {
-                                               qq[l - qbeg] = bq? qual[l] + (bq[l] - 33) : qual[l];
+                                               qq[l - qbeg] = bq? qual[l] + (bq[l] - 64) : qual[l];
                                                if (qq[l - qbeg] > 30) qq[l - qbeg] = 30;
                                                if (qq[l - qbeg] < 7) qq[l - qbeg] = 7;
                                        }