diff --git a/bwamem.c b/bwamem.c index fe83f93..346715c 100644 --- a/bwamem.c +++ b/bwamem.c @@ -1036,6 +1036,7 @@ mem_aln_t mem_reg2aln2(const mem_opt_t *opt, const bntseq_t *bns, const uint8_t } } a.rid = bns_pos2rid(bns, pos); + assert(a.rid == ar->rid); a.pos = pos - bns->anns[a.rid].offset; a.score = ar->score; a.sub = ar->sub > ar->csub? ar->sub : ar->csub; free(query); diff --git a/bwamem_extra.c b/bwamem_extra.c index aee1eb4..2eb87d6 100644 --- a/bwamem_extra.c +++ b/bwamem_extra.c @@ -90,6 +90,7 @@ void mem_reg2ovlp(const mem_opt_t *opt, const bntseq_t *bns, const uint8_t *pac, int64_t pos, rb = p->rb, re = p->re; pos = bns_depos(bns, rb < bns->l_pac? rb : re - 1, &is_rev); rid = bns_pos2rid(bns, pos); + assert(rid == a->rid); pos -= bns->anns[rid].offset; kputs(s->name, &str); kputc('\t', &str); kputw(s->l_seq, &str); kputc('\t', &str); diff --git a/main.c b/main.c index 49cc9ae..d78832b 100644 --- a/main.c +++ b/main.c @@ -4,7 +4,7 @@ #include "utils.h" #ifndef PACKAGE_VERSION -#define PACKAGE_VERSION "0.7.8+dev-r472" +#define PACKAGE_VERSION "0.7.8+dev-r473" #endif int bwa_fa2pac(int argc, char *argv[]);