fork/minimap2
1
0
Fork 0
mirror of https://github.com/lh3/minimap2.git synced 2026-09-22 04:18:11 +08:00
Code Issues Packages Projects Releases Wiki Activity

exoneval may skip terminal exons in aa mode

Browse Source
This commit is contained in:
Heng Li
2025-05-17 23:06:41 -04:00
parent 1fd85be6e2
commit 005c9a1f6b
1 changed files with 8 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
misc/paftools.js
+8 -1
View File
@@ -2613,7 +2613,8 @@ function paf_junceval(args)
function paf_exoneval(args) // adapted from paf_junceval()
{
var c, l_fuzzy = 0, print_ovlp = false, print_err_only = false, first_only = false, chr_only = false, aa = false, is_bed = false, use_cds = false, eval_base = false;
while ((c = getopt(args, "l:epcab1ds")) != null) {
var skip_start = false, skip_last = false;
while ((c = getopt(args, "l:epcab1dsft")) != null) {
if (c == 'l') l_fuzzy = parseInt(getopt.arg);
else if (c == 'e') print_err_only = print_ovlp = true;
else if (c == 'p') print_ovlp = true;
@@ -2623,6 +2624,8 @@ function paf_exoneval(args) // adapted from paf_junceval()
else if (c == '1') first_only = true;
else if (c == 'd') use_cds = true;
else if (c == 's') eval_base = true;
else if (c == 'f') skip_start = true;
else if (c == 't') skip_last = skip_start = true;
}
if (args.length - getopt.ind < 1) {
@@ -2635,6 +2638,8 @@ function paf_exoneval(args) // adapted from paf_junceval()
print(" -e print erroreous overlapping exons");
print(" -c only consider alignments to /^(chr)?([0-9]+|X|Y)$/");
print(" -1 only process the first alignment of each query");
print(" -f skip the first exon in the miniprot mode");
print(" -t skip the first and the last exons");
print(" -b BED as input");
print(" -s compute base Sn and Sp (more memory)");
exit(1);
@@ -2764,6 +2769,8 @@ function paf_exoneval(args) // adapted from paf_junceval()
for (var i = tmp_exon.length - 1; i >= 0; --i)
exon.push([pos + (glen - tmp_exon[i][1]), pos + (glen - tmp_exon[i][0])]);
}
if (skip_start) exon.shift();
if (skip_last) exon.pop();
} else {
var tmp_st = pos;
while ((m = re_cigar.exec(cigar)) != null) {