From 081df6ac7df97c43f2d45a1d62840a62bdf03683 Mon Sep 17 00:00:00 2001 From: xdudiagnoa <47234985+xdudiagnoa@users.noreply.github.com> Date: Mon, 16 Sep 2019 17:21:40 +0800 Subject: [PATCH] Fix example.c seq read logic for every idx should map all input seqs --- example.c | 2 ++ 1 file changed, 2 insertions(+) diff --git a/example.c b/example.c index 80e493ce..ca1fdbe2 100644 --- a/example.c +++ b/example.c @@ -35,6 +35,8 @@ int main(int argc, char *argv[]) while ((mi = mm_idx_reader_read(r, n_threads)) != 0) { // traverse each part of the index mm_mapopt_update(&mopt, mi); // this sets the maximum minimizer occurrence; TODO: set a better default in mm_mapopt_init()! mm_tbuf_t *tbuf = mm_tbuf_init(); // thread buffer; for multi-threading, allocate one tbuf for each thread + gzrewind(f); + kseq_rewind(ks); while (kseq_read(ks) >= 0) { // each kseq_read() call reads one query sequence mm_reg1_t *reg; int j, i, n_reg;