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

r79: drop bad hits

This commit is contained in:
Heng Li 2017-06-26 15:28:04 -04:00
parent 5b614ae828
commit 99c57b86c5
2 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
main.c
View File

@ -10,7 +10,7 @@
#include "minimap.h" #include "minimap.h"
#include "mmpriv.h" #include "mmpriv.h"
#define MM_VERSION "2.0-r78-pre" #define MM_VERSION "2.0-r79-pre"
void liftrlimit() void liftrlimit()
{ {

4
map.c
View File

@ -55,7 +55,7 @@ mm_tbuf_t *mm_tbuf_init(void)
{ {
mm_tbuf_t *b; mm_tbuf_t *b;
b = (mm_tbuf_t*)calloc(1, sizeof(mm_tbuf_t)); b = (mm_tbuf_t*)calloc(1, sizeof(mm_tbuf_t));
// b->km = km_init(); b->km = km_init();
b->sdb = sdust_buf_init(b->km); b->sdb = sdust_buf_init(b->km);
return b; return b;
} }
@ -369,6 +369,8 @@ static void *worker_pipeline(void *shared, int step, void *in)
bseq1_t *t = &s->seq[i]; bseq1_t *t = &s->seq[i];
for (j = 0; j < s->n_reg[i]; ++j) { for (j = 0; j < s->n_reg[i]; ++j) {
mm_reg1_t *r = &s->reg[i][j]; mm_reg1_t *r = &s->reg[i][j];
if (r->p && r->p->blen - r->p->n_ambi - r->p->n_diff < p->opt->min_score)
continue;
if (p->opt->flag & MM_F_OUT_SAM) mm_write_sam(&p->str, mi, t, j, r); if (p->opt->flag & MM_F_OUT_SAM) mm_write_sam(&p->str, mi, t, j, r);
else mm_write_paf(&p->str, mi, t, j, r); else mm_write_paf(&p->str, mi, t, j, r);
puts(p->str.s); puts(p->str.s);