zzh
/
fast-bwa
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Emergent bugfix: wrong TLEN sign 

It is interesting that Picard did not find the issue.
This commit is contained in:
Heng Li 2013-03-09 18:03:15 -05:00
parent 1d132a546d
commit 9ea7f83974
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
bwamem.c
View File

@ -671,7 +671,7 @@ void bwa_hit2sam(kstring_t *str, const int8_t mat[25], int q, int r, int w, cons
if (mid == rid) { if (mid == rid) {
int64_t p0 = p->rb < bns->l_pac? p->rb : (bns->l_pac<<1) - 1 - p->rb; int64_t p0 = p->rb < bns->l_pac? p->rb : (bns->l_pac<<1) - 1 - p->rb;
int64_t p1 = m->rb < bns->l_pac? m->rb : (bns->l_pac<<1) - 1 - m->rb; int64_t p1 = m->rb < bns->l_pac? m->rb : (bns->l_pac<<1) - 1 - m->rb;
kputw(p0 - p1 + (p0 > p1? 1 : p0 < p1? -1 : 0), str); kputw(-(p0 - p1 + (p0 > p1? 1 : p0 < p1? -1 : 0)), str);
} else kputw(0, str); } else kputw(0, str);
kputc('\t', str); kputc('\t', str);
} else if (m && is_mapped(p)) { // then copy the position } else if (m && is_mapped(p)) { // then copy the position