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);
}
}