++j;
} else if (op == BAM_CSOFT_CLIP) {
j += ol;
- } else if (op == BAM_CDEL) {
+ } else if (op == BAM_CDEL || op == BAM_CPAD) {
for (i = 0; i < ol; ++i) s->s[s->l++] = 0;
} else {
- fprintf(stderr, "[depad] ERROR: Didn't expect CIGAR op %c in %s\n", BAM_CIGAR_STR[op], bam1_qname(b));
+ fprintf(stderr, "[depad] ERROR: Didn't expect CIGAR op %c in embedded reference %s\n", BAM_CIGAR_STR[op], bam1_qname(b));
assert(-1);
}
}
int n2 = 0, m2 = 0, *posmap = 0;
h = bam_header_read(in);
+ /* TODO - The reference sequence lengths in the BAM + SAM headers should be updated */
bam_header_write(out, h);
b = bam_init1();
r.l = r.m = q.l = q.m = 0; r.s = q.s = 0;