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

wrong logic: paired mapQ should <q_pe

This commit is contained in:
Heng Li 2013-02-22 12:34:46 -05:00
parent c5ce72f593
commit c0093264de
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_pair.c
View File

@ -269,7 +269,7 @@ int mem_sam_pe(const mem_opt_t *opt, const bntseq_t *bns, const uint8_t *pac, co
if (q_pe > 60) q_pe = 60; if (q_pe > 60) q_pe = 60;
// the following assumes no split hits // the following assumes no split hits
if (z[0] == 0 && z[1] == 0) { // the best hit if (z[0] == 0 && z[1] == 0) { // the best hit
q_pe = q_pe > q_se[0] + q_se[1]? q_pe : q_se[0] + q_se[1]; q_pe = q_pe < q_se[0] + q_se[1]? q_pe : q_se[0] + q_se[1];
if (q_pe > 60) q_pe = 60; if (q_pe > 60) q_pe = 60;
q_se[0] = is_tandem[0]? q_se[0] : q_pe; q_se[0] = is_tandem[0]? q_se[0] : q_pe;
q_se[1] = is_tandem[1]? q_se[1] : q_pe; q_se[1] = is_tandem[1]? q_se[1] : q_pe;