r747: fixed a minor issue in the last (mis)commit

This commit is contained in:
Heng Li 2014-05-02 16:17:50 -04:00
parent 6db761e269
commit e8f28cb529
2 changed files with 3 additions and 3 deletions

View File

@ -248,11 +248,11 @@ mem_chain_v mem_chain(const mem_opt_t *opt, const bwt_t *bwt, const bntseq_t *bn
l_rep += e - b;
for (i = 0; i < aux->mem.n; ++i) {
bwtintv_t *p = &aux->mem.a[i];
int step, slen = (uint32_t)p->info - (p->info>>32); // seed length
int step, count, slen = (uint32_t)p->info - (p->info>>32); // seed length
int64_t k;
if (slen < opt->min_seed_len) continue; // ignore if too short or too repetitive
step = p->x[2] > opt->max_occ? p->x[2] / opt->max_occ : 1;
for (k = 0; k < p->x[2]; k += step) {
for (k = count = 0; k < p->x[2] && count < opt->max_occ; k += step, ++count) {
mem_chain_t tmp, *lower, *upper;
mem_seed_t s;
int rid, to_add = 0;

2
main.c
View File

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