zzh
/
minimap2
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

r238: bugfix to cs - rev sequence not complemented

This commit is contained in:
Heng Li 2017-08-01 10:33:21 -04:00
parent cd105b47f2
commit 12cea727b8
2 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
format.c
View File

@ -72,8 +72,10 @@ static void write_cs(void *km, kstring_t *s, const mm_idx_t *mi, const mm_bseq1_
for (i = r->qs; i < r->qe; ++i) for (i = r->qs; i < r->qe; ++i)
qseq[i - r->qs] = seq_nt4_table[(uint8_t)t->seq[i]]; qseq[i - r->qs] = seq_nt4_table[(uint8_t)t->seq[i]];
} else { } else {
for (i = r->qs; i < r->qe; ++i) for (i = r->qs; i < r->qe; ++i) {
qseq[r->qe - i - 1] = seq_nt4_table[(uint8_t)t->seq[i]]; uint8_t c = seq_nt4_table[(uint8_t)t->seq[i]];
qseq[r->qe - i - 1] = c >= 4? 4 : 3 - c;
}
} }
for (i = q_off = t_off = 0; i < r->p->n_cigar; ++i) { for (i = q_off = t_off = 0; i < r->p->n_cigar; ++i) {
int j, op = r->p->cigar[i]&0xf, len = r->p->cigar[i]>>4; int j, op = r->p->cigar[i]&0xf, len = r->p->cigar[i]>>4;

2
main.c
View File

@ -8,7 +8,7 @@
#include "minimap.h" #include "minimap.h"
#include "mmpriv.h" #include "mmpriv.h"
#define MM_VERSION "2.0rc1-r237-dirty" #define MM_VERSION "2.0rc1-r238-dirty"
void liftrlimit() void liftrlimit()
{ {