From 9c0e2c67f8aa499c104ea09b42ab56e8eacd8fe1 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Tue, 21 Apr 2020 13:21:03 -0400 Subject: [PATCH] r986: don't estimate dv with --qstrand --- main.c | 2 +- map.c | 3 ++- 2 files changed, 3 insertions(+), 2 deletions(-) diff --git a/main.c b/main.c index e84a7f6..15074bd 100644 --- a/main.c +++ b/main.c @@ -7,7 +7,7 @@ #include "mmpriv.h" #include "ketopt.h" -#define MM_VERSION "2.17-r985-qs-dirty" +#define MM_VERSION "2.17-r986-qs-dirty" #ifdef __linux__ #include diff --git a/map.c b/map.c index 1cc4a05..1def5fd 100644 --- a/map.c +++ b/map.c @@ -358,7 +358,8 @@ void mm_map_frag(const mm_idx_t *mi, int n_segs, const int *qlens, const char ** i == regs0[j].as? 0 : ((int32_t)a[i].y - (int32_t)a[i-1].y) - ((int32_t)a[i].x - (int32_t)a[i-1].x)); chain_post(opt, max_chain_gap_ref, mi, b->km, qlen_sum, n_segs, qlens, &n_regs0, regs0, a); - if (!is_sr) mm_est_err(mi, qlen_sum, n_regs0, regs0, a, n_mini_pos, mini_pos); + if (!is_sr && !(opt->flag&MM_F_QSTRAND)) + mm_est_err(mi, qlen_sum, n_regs0, regs0, a, n_mini_pos, mini_pos); if (n_segs == 1) { // uni-segment regs0 = align_regs(opt, mi, b->km, qlens[0], seqs[0], &n_regs0, regs0, a);