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

r415: bug in the new (optional) mapQ computation 

I may use the new method as the default. Testing needed.
This commit is contained in:
Heng Li 2013-09-09 17:51:05 -04:00
parent ebb7b02e9b
commit 7144a0cefc
2 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
bwamem.c
View File

@ -807,7 +807,7 @@ int mem_approx_mapq_se(const mem_opt_t *opt, const mem_alnreg_t *a)
} else if (opt->mapQ_coef_len > 0) {
double tmp;
tmp = l < opt->mapQ_coef_len? 1. : opt->mapQ_coef_fac / log(l);
tmp *= a->seedcov < l? (double)a->seedcov / identity : 1;
tmp *= identity * identity;
mapq = (int)(6.02 * (a->score - sub) / opt->a * tmp * tmp + .499);
} else {
mapq = (int)(MEM_MAPQ_COEF * (1. - (double)sub / a->score) * log(a->seedcov) + .499);

2
main.c
View File

@ -3,7 +3,7 @@
#include "utils.h"
#ifndef PACKAGE_VERSION
#define PACKAGE_VERSION "0.7.5a-r414"
#define PACKAGE_VERSION "0.7.5a-r415"
#endif
int bwa_fa2pac(int argc, char *argv[]);