zzh
/
bwa
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

dev-473: added a few assertions 

to make sure the new change works as is expected
This commit is contained in:
Heng Li 2014-04-10 21:03:13 -04:00
parent 8638cfadc8
commit f02cd42679
3 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
bwamem.c
View File

@ -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); a.rid = bns_pos2rid(bns, pos);
assert(a.rid == ar->rid);
a.pos = pos - bns->anns[a.rid].offset; a.pos = pos - bns->anns[a.rid].offset;
a.score = ar->score; a.sub = ar->sub > ar->csub? ar->sub : ar->csub; a.score = ar->score; a.sub = ar->sub > ar->csub? ar->sub : ar->csub;
free(query); free(query);

1
bwamem_extra.c
View File

@ -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; int64_t pos, rb = p->rb, re = p->re;
pos = bns_depos(bns, rb < bns->l_pac? rb : re - 1, &is_rev); pos = bns_depos(bns, rb < bns->l_pac? rb : re - 1, &is_rev);
rid = bns_pos2rid(bns, pos); rid = bns_pos2rid(bns, pos);
assert(rid == a->rid);
pos -= bns->anns[rid].offset; pos -= bns->anns[rid].offset;
kputs(s->name, &str); kputc('\t', &str); kputs(s->name, &str); kputc('\t', &str);
kputw(s->l_seq, &str); kputc('\t', &str); kputw(s->l_seq, &str); kputc('\t', &str);

2
main.c
View File

@ -4,7 +4,7 @@
#include "utils.h" #include "utils.h"
#ifndef PACKAGE_VERSION #ifndef PACKAGE_VERSION
#define PACKAGE_VERSION "0.7.8+dev-r472" #define PACKAGE_VERSION "0.7.8+dev-r473"
#endif #endif
int bwa_fa2pac(int argc, char *argv[]); int bwa_fa2pac(int argc, char *argv[]);