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] - 64) : qual[l];