q = (uint32_t)bam1_qual(p->b)[p->qpos];
if (q < bm->min_baseQ) continue;
x |= (uint32_t)bam1_strand(p->b) << 18 | q << 8 | p->b->core.qual;
q = (uint32_t)bam1_qual(p->b)[p->qpos];
if (q < bm->min_baseQ) continue;
x |= (uint32_t)bam1_strand(p->b) << 18 | q << 8 | p->b->core.qual;