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

correctly parse ##PAF lines in GFF3 (for miniprot)

This commit is contained in:
Heng Li 2022-10-07 19:33:27 -04:00
parent 6c2cbf7903
commit 85a5260408
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
misc/paftools.js
View File

@ -2422,9 +2422,11 @@ function paf_junceval(args)
var re_cigar = /(\d+)([MIDNSHP=XFGUV])/g; var re_cigar = /(\d+)([MIDNSHP=XFGUV])/g;
while (file.readline(buf) >= 0) { while (file.readline(buf) >= 0) {
var m, t = buf.toString().split("\t"); var m, t = buf.toString().split("\t");
var ctg_name = null, cigar = null, pos = null, qname = t[0]; var ctg_name = null, cigar = null, pos = null, qname;
if (t[0].charAt(0) == '@') continue; if (t[0].charAt(0) == '@') continue;
if (t[0] == "##PAF") t.shift();
qname = t[0];
if (is_bed) { if (is_bed) {
ctg_name = t[0], pos = parseInt(t[1]), cigar == null; ctg_name = t[0], pos = parseInt(t[1]), cigar == null;
} else if (t[4] == '+' || t[4] == '-' || t[4] == '*') { // PAF } else if (t[4] == '+' || t[4] == '-' || t[4] == '*') { // PAF