From: peterjc Date: Fri, 23 Mar 2012 12:51:32 +0000 (+0000) Subject: Cope with CIGAR D or P to mark gaps in an embedded padded reference X-Git-Url: https://git.donarmstrong.com/?p=samtools.git;a=commitdiff_plain;h=fe8fbc9eae570d7afc3e0e91df735e6b6fdc329d Cope with CIGAR D or P to mark gaps in an embedded padded reference --- diff --git a/padding.c b/padding.c index ad68728..b330326 100644 --- a/padding.c +++ b/padding.c @@ -42,10 +42,10 @@ static void unpad_seq(bam1_t *b, kstring_t *s) ++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); } }