fork/minimap2
1
0
Fork 0
mirror of https://github.com/lh3/minimap2.git synced 2026-10-06 15:48:12 +08:00
Code Issues Packages Projects Releases Wiki Activity

Merge remote-tracking branch 'remotes/origin/master'

Browse Source
This commit is contained in:
Heng Li
2025-12-12 11:39:33 -05:00
parent 9bb4d2bed4 6d49eb690f
commit e2542e6425
1 changed files with 13 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
misc/paftools.js
+13 -3
View File
@@ -1740,10 +1740,11 @@ function paf_gff2bed(args)
function paf_sam2paf(args) function paf_sam2paf(args)
{ {
var c, pri_only = false, long_cs = false, pri_pri_only = false; var c, pri_only = false, long_cs = false, pri_pri_only = false, allow_unmapped = false;
while ((c = getopt(args, "pPL")) != null) { while ((c = getopt(args, "pPUL")) != null) {
if (c == 'p') pri_only = true; if (c == 'p') pri_only = true;
else if (c == 'P') pri_pri_only = pri_only = true; else if (c == 'P') pri_pri_only = pri_only = true;
else if (c == 'U') allow_unmapped = true;
else if (c == 'L') long_cs = true; else if (c == 'L') long_cs = true;
} }
if (args.length == getopt.ind) { if (args.length == getopt.ind) {
@@ -1751,6 +1752,7 @@ function paf_sam2paf(args)
print("Options:"); print("Options:");
print(" -p convert primary or supplementary alignments only"); print(" -p convert primary or supplementary alignments only");
print(" -P convert primary alignments only"); print(" -P convert primary alignments only");
print(" -U convert unmapped reads as well");
print(" -L output the cs tag in the long form"); print(" -L output the cs tag in the long form");
exit(1); exit(1);
} }
@@ -1775,7 +1777,15 @@ function paf_sam2paf(args)
var flag = parseInt(t[1]); var flag = parseInt(t[1]);
if (t[9] != '*' && t[10] != '*' && t[9].length != t[10].length) if (t[9] != '*' && t[10] != '*' && t[9].length != t[10].length)
throw Error("at line " + lineno + ": inconsistent SEQ and QUAL lengths - " + t[9].length + " != " + t[10].length); throw Error("at line " + lineno + ": inconsistent SEQ and QUAL lengths - " + t[9].length + " != " + t[10].length);
if (t[2] == '*' || (flag&4) || t[5] == '*') continue; if (t[2] == '*' || (flag&4) || t[5] == '*') {
if (allow_unmapped) {
// emit an unmapped PAF line instead of skipping
// fields: qname, qlen, qstart, qend, strand, tname, tlen, tstart, tend, n_match, aln_len, mapq
var qlen_val = (t[9] == '*' ? 0 : t[9].length);
print([t[0], qlen_val, 0, 0, '*', '*', 0, 0, 0, 0, 0, 0].join("\t"));
}
continue;
}
if (pri_only && (flag&0x100)) continue; if (pri_only && (flag&0x100)) continue;
if (pri_pri_only && (flag&0x900)) continue; if (pri_pri_only && (flag&0x900)) continue;
var tlen = ctg_len[t[2]]; var tlen = ctg_len[t[2]];