From 9cd313eae12e859c6ada77fa562b328b7b5054d8 Mon Sep 17 00:00:00 2001 From: Heng Li Date: Fri, 23 Jun 2017 14:11:56 -0400 Subject: [PATCH] sequence retrieval working --- align.c | 2 +- minimap.h | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/align.c b/align.c index e56fc99..7f9e138 100644 --- a/align.c +++ b/align.c @@ -30,7 +30,7 @@ void mm_align_skeleton(void *km, const mm_mapopt_t *opt, const mm_idx_t *mi, int qseq = &qseq0[rev][qs]; ret = mm_idx_getseq(mi, rid, rs, re, tseq); assert(ret > 0); - fprintf(stderr, "===> %d,%d (%s:%d-%d) <===\n", reg, i, mi->seq[rid].name, rs, re); + fprintf(stderr, "===> [%d,%d] %d-%d %c (%s:%d-%d) <===\n", reg, i, qs, qe, "+-"[rev], mi->seq[rid].name, rs, re); for (k = 0; k < re - rs; ++k) fputc("ACGTN"[tseq[k]], stderr); fputc('\n', stderr); for (k = 0; k < qe - qs; ++k) fputc("ACGTN"[qseq[k]], stderr); fputc('\n', stderr); } diff --git a/minimap.h b/minimap.h index 5ab9c6a..53273d3 100644 --- a/minimap.h +++ b/minimap.h @@ -74,8 +74,8 @@ typedef struct mm_tbuf_s mm_tbuf_t; struct bseq_file_s; -#define mm_seq4_set(s, i, c) ((s)[(i)>>4] |= (uint32_t)(c) << (((i)&7)<<2)) -#define mm_seq4_get(s, i) ((s)[(i)>>4] >> (((i)&7)<<2) & 0xf) +#define mm_seq4_set(s, i, c) ((s)[(i)>>3] |= (uint32_t)(c) << (((i)&7)<<2)) +#define mm_seq4_get(s, i) ((s)[(i)>>3] >> (((i)&7)<<2) & 0xf) #ifdef __cplusplus extern "C" {