fork/minimap2
1
0
Fork 0
mirror of https://github.com/lh3/minimap2.git synced 2026-09-24 15:08:11 +08:00
Code Issues Packages Projects Releases Wiki Activity

support inserting RG lines

Browse Source
This commit is contained in:
Heng Li
2017-08-17 23:34:09 +08:00
parent b5f5929bf9
commit bbb37d95f2
5 changed files with 83 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files
map.c
+1 -5
View File
@@ -386,11 +386,7 @@ int mm_map_file(const mm_idx_t *idx, const char *fn, const mm_mapopt_t *opt, int
if (pl.fp == 0) return -1;
pl.opt = opt, pl.mi = idx;
pl.n_threads = n_threads, pl.mini_batch_size = mini_batch_size;
if (opt->flag & MM_F_OUT_SAM) {
uint32_t i;
for (i = 0; i < idx->n_seq; ++i)
printf("@SQ\tSN:%s\tLN:%d\n", idx->seq[i].name, idx->seq[i].len);
}
if (opt->flag & MM_F_OUT_SAM) sam_write_sam_header(idx);
kt_pipeline(n_threads == 1? 1 : 2, worker_pipeline, &pl, 3);
free(pl.str.s);
mm_bseq_close(pl.fp);