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

tunable overlap ratio

This commit is contained in:
Heng Li 2017-07-26 18:30:29 -04:00
parent 323293fbda
commit 9f728dd96a
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
misc/pbsim-eval.js
View File

@ -36,9 +36,10 @@ var getopt = function(args, ostr) {
return optopt; return optopt;
} }
var c, max_mapq = 60, mode = 0, err_out_q = 256, print_err = false; var c, max_mapq = 60, mode = 0, err_out_q = 256, print_err = false, ovlp_ratio = 0.333;
while ((c = getopt(arguments, "Q:")) != null) { while ((c = getopt(arguments, "Q:r:")) != null) {
if (c == 'Q') err_out_q = parseInt(getopt.arg), print_err = true; if (c == 'Q') err_out_q = parseInt(getopt.arg), print_err = true;
else if (c == 'r') ovlp_ratio = parseFloat(getopt.arg);
} }
var file = arguments.length == getopt.ind? new File() : new File(arguments[getopt.ind]); var file = arguments.length == getopt.ind? new File() : new File(arguments[getopt.ind]);
@ -61,7 +62,7 @@ function is_correct(s, b)
o = (s[2] < b[2]? s[2] : b[2]) - s[1]; o = (s[2] < b[2]? s[2] : b[2]) - s[1];
l = (s[2] > b[2]? s[2] : b[2]) - b[1]; l = (s[2] > b[2]? s[2] : b[2]) - b[1];
} }
return o/l > .333? true : false; return o/l > ovlp_ratio? true : false;
} }
function count_err(qname, a, tot, err, mode) function count_err(qname, a, tot, err, mode)