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

r282: set min split_len to read length

This commit is contained in:
Heng Li 2013-02-25 17:29:35 -05:00
parent 30cc8a95d1
commit 77b5b586ad
2 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
bwamem.c
View File

@ -199,6 +199,7 @@ static void mem_insert_seed(const mem_opt_t *opt, kbtree_t(chn) *tree, smem_i *i
{ {
const bwtintv_v *a; const bwtintv_v *a;
int split_len = (int)(opt->min_seed_len * opt->split_factor + .499); int split_len = (int)(opt->min_seed_len * opt->split_factor + .499);
split_len = split_len < itr->len? split_len : itr->len;
while ((a = smem_next(itr, split_len, opt->split_width)) != 0) { // to find all SMEM and some internal MEM while ((a = smem_next(itr, split_len, opt->split_width)) != 0) { // to find all SMEM and some internal MEM
int i; int i;
for (i = 0; i < a->n; ++i) { // go through each SMEM/MEM up to itr->start for (i = 0; i < a->n; ++i) { // go through each SMEM/MEM up to itr->start

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.6.2-r281-beta" #define PACKAGE_VERSION "0.6.2-r282-beta"
#endif #endif
static int usage() static int usage()