* the BQ tag is now 64 shifted, not 33 shifted
bq = (uint8_t*)bam_aux_get(p->b, "BQ");
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;
}
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);
}
#endif
#ifndef PACKAGE_VERSION
-#define PACKAGE_VERSION "0.1.9-16 (r822)"
+#define PACKAGE_VERSION "0.1.9-16 (r823)"
#endif
int bam_taf2baf(int argc, char *argv[]);