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

sequence retrieval working

Browse Source
This commit is contained in:
Heng Li
2017-06-23 14:11:56 -04:00
parent 326d91deb0
commit 9cd313eae1
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
align.c
+1 -1
View File
@@ -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);
}