Fix example.c seq read logic

for every idx should map all input seqs
This commit is contained in:
xdudiagnoa 2019-09-16 17:21:40 +08:00 committed by Heng Li
parent a3e7a575fb
commit 081df6ac7d
1 changed files with 2 additions and 0 deletions

View File

@ -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;