From a832a42f6f1e20a8824e3cb85a8e807a27776787 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Fri, 11 Apr 2025 01:08:16 -0400 Subject: [PATCH] r1270: fixed ts tag for PE reads this also fixed the wrong --write-junc --- map.c | 4 ++++ minimap.h | 2 +- 2 files changed, 5 insertions(+), 1 deletion(-) diff --git a/map.c b/map.c index e2cb1c5..85b812c 100644 --- a/map.c +++ b/map.c @@ -463,6 +463,10 @@ static void worker_for(void *_data, long i, int tid) // kt_for() callback r->qs = qlens[j] - r->qe; r->qe = qlens[j] - t; r->rev = !r->rev; + if (r->p) { + if (r->p->trans_strand == 1) r->p->trans_strand = 2; + else if (r->p->trans_strand == 2) r->p->trans_strand = 1; + } } } if (mm_dbg_flag & MM_DBG_PRINT_QNAME) diff --git a/minimap.h b/minimap.h index b53a0ed..1e9a720 100644 --- a/minimap.h +++ b/minimap.h @@ -5,7 +5,7 @@ #include #include -#define MM_VERSION "2.28-r1269-dirty" +#define MM_VERSION "2.28-r1270-dirty" #define MM_F_NO_DIAG (0x001LL) // no exact diagonal hit #define MM_F_NO_DUAL (0x002LL) // skip pairs where query name is lexicographically larger than target name