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

r870: separate <50% and <10%

This commit is contained in:
Heng Li 2018-11-04 19:31:25 -05:00
parent 9e27575387
commit 57ec73ec6c
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/paftools.js
View File

@ -1,6 +1,6 @@
#!/usr/bin/env k8 #!/usr/bin/env k8
var paftools_version = '2.13-r869-dirty'; var paftools_version = '2.13-r870-dirty';
/***************************** /*****************************
***** Library functions ***** ***** Library functions *****
@ -831,7 +831,7 @@ function paf_asmgene(args)
file.close(); file.close();
} }
var col1 = ["full_sgl", "full_dup", "frag", "part50+", "part50-"]; var col1 = ["full_sgl", "full_dup", "frag", "part50+", "part10+", "part10-"];
var rst = []; var rst = [];
for (var k = 0; k < col1.length; ++k) { for (var k = 0; k < col1.length; ++k) {
rst[k] = []; rst[k] = [];
@ -846,7 +846,8 @@ function paf_asmgene(args)
else if (gene[g][i][0] > 1) rst[1][i]++; else if (gene[g][i][0] > 1) rst[1][i]++;
else if (gene[g][i][1] >= opt.min_cov) rst[2][i]++; else if (gene[g][i][1] >= opt.min_cov) rst[2][i]++;
else if (gene[g][i][1] >= 0.5) rst[3][i]++; else if (gene[g][i][1] >= 0.5) rst[3][i]++;
else rst[4][i]++; else if (gene[g][i][1] >= 0.1) rst[4][i]++;
else rst[5][i]++;
} }
} }
print(header.join("\t")); print(header.join("\t"));