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

r422: matesw hits not sorted

This commit is contained in:
Heng Li 2013-11-21 14:43:50 -05:00
parent ff4762f3c7
commit 29aa855432
2 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
bwamem_pair.c
View File

@ -108,6 +108,7 @@ void mem_pestat(const mem_opt_t *opt, int64_t l_pac, int n, const mem_alnreg_v *
int mem_matesw(const mem_opt_t *opt, int64_t l_pac, const uint8_t *pac, const mem_pestat_t pes[4], const mem_alnreg_t *a, int l_ms, const uint8_t *ms, mem_alnreg_v *ma) int mem_matesw(const mem_opt_t *opt, int64_t l_pac, const uint8_t *pac, const mem_pestat_t pes[4], const mem_alnreg_t *a, int l_ms, const uint8_t *ms, mem_alnreg_v *ma)
{ {
extern int mem_sort_and_dedup(int n, mem_alnreg_t *a, float mask_level_redun);
int i, r, skip[4], n = 0; int i, r, skip[4], n = 0;
for (r = 0; r < 4; ++r) for (r = 0; r < 4; ++r)
skip[r] = pes[r].failed? 1 : 0; skip[r] = pes[r].failed? 1 : 0;
@ -166,6 +167,7 @@ int mem_matesw(const mem_opt_t *opt, int64_t l_pac, const uint8_t *pac, const me
} }
++n; ++n;
} }
if (n) ma->n = mem_sort_and_dedup(ma->n, ma->a, opt->mask_level_redun);
if (rev) free(rev); if (rev) free(rev);
free(ref); free(ref);
} }

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.5a-r421" #define PACKAGE_VERSION "0.7.5a-r422"
#endif #endif
int bwa_fa2pac(int argc, char *argv[]); int bwa_fa2pac(int argc, char *argv[]);