Repository 'dr_disco'
hg clone https://toolshed.g2.bx.psu.edu/repos/erasmus-medical-center/dr_disco

Changeset 9:c3ce472a315f (2017-09-26)
Previous changeset 8:922b52cf9b35 (2017-09-14) Next changeset 10:f2c70d589596 (2017-11-02)
Commit message:
planemo upload for repository https://github.com/ErasmusMC-Bioinformatics/galaxytools-emc/tree/master/tools/dr-disco commit 78dd573f8aa888ee261580852bf049b4b3500b52
modified:
dr-disco_integrate.xml
macros.xml
test-data/classify_test_16.in.dbed
test-data/classify_test_16.test-01.out.dbed
test-data/classify_test_16.test-02.out.dbed
test-data/classify_test_16.test-03.out.dbed
test-data/classify_test_16.test-04.out.dbed
test-data/classify_test_16.test-05.out.dbed
test-data/classify_test_16.test-06.out.dbed
test-data/classify_test_16.test-07.out.dbed
test-data/classify_test_16.test-08.out.dbed
test-data/detect.txt
test-data/fix_test-01.fixed.bam
test-data/integrate_test_terg_s041.in.dbed
test-data/integrate_test_terg_s041.out.no-gtf.txt
test-data/integrate_test_terg_s041.out.txt
b
diff -r 922b52cf9b35 -r c3ce472a315f dr-disco_integrate.xml
--- a/dr-disco_integrate.xml Thu Sep 14 08:35:12 2017 -0400
+++ b/dr-disco_integrate.xml Tue Sep 26 08:27:47 2017 -0400
b
@@ -20,7 +20,7 @@
 
     <inputs>
         <param name="dr_disco_classify_output" type="data" format="tabular" label="Discordant alignment file of STAR (processed with dr-disco fix)" />
-        <param argument="--gtf" type="data" format="gtf" optional="True" label="GTF file (Optional: for predicting frame shifts)" help="This GTF file requires the following attributes annotated: gene_name, transcript_id and transcript_version"/>
+        <param argument="--gtf" type="data" format="gtf" optional="True" label="GTF file (Optional: for predicting frame shifts)" help="This GTF file requires the following attributes: gene_name, transcript_id and transcript_version"/>
     </inputs>
     
     <outputs>
b
diff -r 922b52cf9b35 -r c3ce472a315f macros.xml
--- a/macros.xml Thu Sep 14 08:35:12 2017 -0400
+++ b/macros.xml Tue Sep 26 08:27:47 2017 -0400
b
@@ -1,12 +1,12 @@
 <macros>
     <xml name="requirements">
         <requirements>
-            <requirement type="package" version="0.10.0">dr-disco</requirement>
+            <requirement type="package" version="0.11.0=py27_0">dr-disco</requirement>
             <yield/>
         </requirements>
     </xml>
 
-    <token name="@TOOL_VERSION@">0.10.0-g3</token>
+    <token name="@TOOL_VERSION@">0.11.0</token>
 
     <xml name="version_command">
         <version_command>dr-disco --version | head -n 1</version_command>
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.in.dbed
--- a/test-data/classify_test_16.in.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.in.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tunclassified\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tunclassified\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tunclassified\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tunclassified\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tunclassified\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n-chr21\t39851814\t+\t0\t6\tchr21\t39852667\t-\t6\t0\t853\tunclassified\tcircula'..b'ed_1_t:1,spanning_paired_2_t:1)\n+chr21\t39838727\t+\t0\t2\tchr21\t39863701\t-\t2\t0\t24974\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t56\t0\t1\t0\t0\t0.3333\t2\tchr21:39838727/39838728(+)->chr21:39863701/39863702(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39849621\t+\t0\t2\tchr21\t39850505\t-\t2\t0\t884\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tunclassified\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tunclassified\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tunclassified\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tunclassified\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tunclassified\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tunclassified\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tunclassified\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tunclassified\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tunclassified\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tunclassified\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tunclassified\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-01.out.dbed
--- a/test-data/classify_test_16.test-01.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-01.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-02.out.dbed
--- a/test-data/classify_test_16.test-02.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-02.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
@@ -1,3 +1,3 @@
-chr-A pos-A direction-A pos-A-acceptor pos-A-donor chr-B pos-B direction-B pos-A-acceptor pos-A-donor genomic-distance filter-status circRNA intronic/exonic score soft+hardclips n-split-reads n-discordant-reads n-edges n-nodes-A n-nodes-B n-splice-junc-A n-splice-junc-B entropy-bp-edge entropy-all-edges bp-pos-stddev entropy-disco-bps lr-A-slope lr-A-intercept lr-A-rvalue lr-A-pvalue lr-A-stderr lr-B-slope lr-B-intercept lr-B-rvalue lr-B-pvalue lr-B-stderr disco/split clips/score nodes/edge data-structure
-chr21 39817544 - 140 0 chr21 42880007 + 0 140 3062463 valid linear exonic 253 170 85 6 5 2 4 1 0 0.8111 0.8497 0.0000 0.9535 0.7247 63.7843 0.8440 0.0000 0.0559 0.3651 17.3030 0.8378 0.0000 0.0289 0.0706 0.3360 1.2000 chr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)
-chr21 39877811 - 98 0 chr21 42873374 + 0 98 2995563 valid linear intronic 91 42 21 28 1 1 1 0 0 0.9528 0.9528 0.0000 0.8200 2.9207 17.3314 0.9125 0.0000 0.1910 1.5880 68.1747 0.8769 0.0000 0.1269 1.3333 0.2308 2.0000 chr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)
+chr-A pos-A direction-A pos-A-acceptor pos-A-donor chr-B pos-B direction-B pos-A-acceptor pos-A-donor genomic-distance filter-status circRNA intronic/exonic score soft+hardclips n-split-reads n-discordant-reads alignment-score mismatches n-edges n-nodes-A n-nodes-B n-splice-junc-A n-splice-junc-B entropy-bp-edge entropy-all-edges bp-pos-stddev entropy-disco-bps lr-A-slope lr-A-intercept lr-A-rvalue lr-A-pvalue lr-A-stderr lr-B-slope lr-B-intercept lr-B-rvalue lr-B-pvalue lr-B-stderr disco/split clips/score nodes/edge data-structure
+chr21 39817544 - 140 0 chr21 42880007 + 0 140 3062463 valid linear exonic 253 170 85 6 10000 0 5 2 4 1 0 0.8111 0.8497 0 0.9535 0.7247 63.7843 0.844 0 0.0559 0.3651 17.303 0.8378 0 0.0289 0.0706 0.336 1.2 chr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)
+chr21 39877811 - 98 0 chr21 42873374 + 0 98 2995563 valid linear intronic 91 42 21 28 10000 0 1 1 1 0 0 0.9528 0.9528 0 0.82 2.9207 17.3314 0.9125 0 0.191 1.588 68.1747 0.8769 0 0.1269 1.3333 0.2308 2 chr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-03.out.dbed
--- a/test-data/classify_test_16.test-03.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-03.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-04.out.dbed
--- a/test-data/classify_test_16.test-04.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-04.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-05.out.dbed
--- a/test-data/classify_test_16.test-05.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-05.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-06.out.dbed
--- a/test-data/classify_test_16.test-06.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-06.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-07.out.dbed
--- a/test-data/classify_test_16.test-07.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-07.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/classify_test_16.test-08.out.dbed
--- a/test-data/classify_test_16.test-08.out.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/classify_test_16.test-08.out.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,37 +1,37 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t39817544\t-\t140\t0\tchr21\t42880007\t+\t0\t140\t3062463\tvalid\tlinear\texonic\t253\t170\t85\t6\t5\t2\t4\t1\t0\t0.8111\t0.8497\t0.0000\t0.9535\t0.7247\t63.7843\t0.8440\t0.0000\t0.0559\t0.3651\t17.3030\t0.8378\t0.0000\t0.0289\t0.0706\t0.3360\t1.2000\tchr21:39817544/39817545(-)->chr21:42880007/42880008(+):(spanning_paired_1:65,spanning_paired_2:65,spanning_singleton_1:2,spanning_singleton_1_r:3,spanning_singleton_2:2,spanning_singleton_2_r:3)&chr21:39817544/39817545(-)->chr21:42879876/42879877(+):(discordant_mates:10,spanning_paired_1:6,spanning_paired_1_t:1,spanning_paired_2:6,spanning_paired_2_t:1,spanning_singleton_1_r:2,spanning_singleton_2_r:2)&chr21:39817544/39817545(-)->chr21:42878371/42878372(+):(discordant_mates:2,spanning_paired_1:1,spanning_paired_1_t:1,spanning_paired_2:1,spanning_paired_2_t:1,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:39846044/39846045(-)->chr21:42879876/42879877(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:39817544/39817545(-)->chr21:42876293/42876294(+):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t39877811\t-\t98\t0\tchr21\t42873374\t+\t0\t98\t2995563\tvalid\tlinear\tintronic\t91\t42\t21\t28\t1\t1\t1\t0\t0\t0.9528\t0.9528\t0.0000\t0.8200\t2.9207\t17.3314\t0.9125\t0.0000\t0.1910\t1.5880\t68.1747\t0.8769\t0.0000\t0.1269\t1.3333\t0.2308\t2.0000\tchr21:39877811/39877812(-)->chr21:42873374/42873375(+):(discordant_mates:56,spanning_paired_1:12,spanning_paired_1_t:9,spanning_paired_2:12,spanning_paired_2_t:9)\n-chr21\t42861433\t+\t6\t18\tchr21\t42866505\t-\t18\t6\t5072\tn_discordant_reads=2<3\tcircular\texonic\t41\t30\t14\t2\t7\t5\t6\t1\t2\t1.0000\t0.9375\t18.7474\t0.7500\t26.9000\t3.0000\t0.9425\t0.0164\t5.5097\t20.1000\t42.8000\t0.8833\t0.0470\t6.1587\t0.1429\t0.3659\t1.5714\tchr21:42861433/42861434(+)->chr21:42866505/42866506(-):(discordant_mates:2,spanning_paired_1:2,spanning_paired_1_t:1,spanning_paired_2:2,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42861433/42861434(+)->chr21:42870116/42870117(-):(spanning_paired_1:3,spanning_paired_2:3,spanning_singleton_1_r:1,spanning_singleton_2_r:1)&chr21:42860568/42860569(+)->chr21:42867210/42867211(-):(spanning_paired_1_t:2,spanning_paired_2_t:2)&chr21:42860320/42860321(+)->chr21:42861520/42861521(-):(spanning_paired_1:1,spanning_paired_2:1,spanning_singleton_1:1,spanning_singleton_2:1)&chr21:42859986/42859987(+)->chr21:42863878/42863879(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:42860320/42860321(+)->chr21:42866505/42866506(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:42866446/42866447(+)->chr21:42870110/42870111(-):(discordant_mates:2)\n-chr21\t39836011\t+\t0\t8\tchr21\t39836545\t-\t8\t0\t534\tn_support=6<8,bp_pos_stddev=25.3065>15.0,log_ratio_slope=9.62>1.8,log_ratio_rvalue=8.96>0.4\tcircular\tintronic\t17\t12\t6\t0\t2\t2\t2\t0\t0\t0.4056\t0.6934\t25.3065\t0.0000\t1.5000\t29.5000\t0.7746\t0.2254\t0.8660\t0.0000\t93.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3529\t2.0000\tchr21:39836011/39836012(+)->chr21:39836545/39836546(-):(spanning_paired_1_t:4,spanning_paired_2_t:4)&chr21:39835707/39835708(+)->chr21:39839850/39839851(-):(spanning_paired_1_t:1,spanning_paired_2_t:1,spanning_singleton_1:1,spanning_singleton_2:1)\n-chr21\t39875880\t+\t0\t6\tchr21\t39876476\t-\t6\t0\t596\tentropy=0.0<0.6839,n_support=4<8\tcircular\tintronic\t12\t8\t4\t0\t2\t2\t2\t0\t0\t0.0000\t0.4056\t0.0000\t0.0000\t0.0000\t68.0000\t0.0000\t1.0000\t0.0000\t0.0000\t58.0000\t0.0000\t1.0000\t0.0000\t0.0000\t0.3333\t2.0000\tchr21:39875880/39875881(+)->chr21:39876476/39876477(-):(spanning_paired_1:3,spanning_paired_2:3)&chr21:39874280/39874281(+)->chr21:39876277/39876278(-):(spanning_paired_1'..b'rt=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t52\t0\t1\t0\t0\t74\t0\t1\t0\t0\t0.3333\t2\tchr21:39849621/39849622(+)->chr21:39850505/39850506(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39858845\t+\t0\t2\tchr21\t39859592\t-\t2\t0\t747\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t41\t0\t1\t0\t0\t84\t0\t1\t0\t0\t0.3333\t2\tchr21:39858845/39858846(+)->chr21:39859592/39859593(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t39863020\t+\t0\t2\tchr21\t39863629\t-\t2\t0\t609\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t101\t0\t1\t0\t0\t25\t0\t1\t0\t0\t0.3333\t2\tchr21:39863020/39863021(+)->chr21:39863629/39863630(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t39871034\t+\t0\t2\tchr21\t39877568\t-\t2\t0\t6534\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t31\t0\t1\t0\t0\t66\t0\t1\t0\t0\t0.3333\t2\tchr21:39871034/39871035(+)->chr21:39877568/39877569(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42836581\t+\t0\t2\tchr21\t42839672\t-\t2\t0\t3091\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t19\t0\t1\t0\t0\t107\t0\t1\t0\t0\t0.3333\t2\tchr21:42836581/42836582(+)->chr21:42839672/42839673(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)\n+chr21\t42853701\t+\t0\t2\tchr21\t42860029\t-\t2\t0\t6328\tentropy=0.0<0.6828,n_support=1<7\tcircular\tintronic\t3\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t22\t0\t1\t0\t0\t31\t0\t1\t0\t0\t0.3333\t2\tchr21:42853701/42853702(+)->chr21:42860029/42860030(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42866483\t-\t0\t2\tchr21\t42870116\t-\t2\t0\t3633\tentropy=0.0<0.6828,n_support=1<7\tlinear\tintronic\t3\t9\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t38\t0\t1\t0\t0\t88\t0\t1\t0\t0\t1.5\t2\tchr21:42866483/42866484(-)->chr21:42870116/42870117(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t42839824\t+\t0\t4\tchr21\t42840323\t-\t4\t0\t499\tn_support=2<7,log_ratio_slope=12.99>1.8,log_ratio_rvalue=9.21>0.4\tcircular\tintronic\t2\t2\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t1\t1\t0\t0.4056\t44\t82\t1\t0\t0\t0\t120\t0\t1\t0\t40\t0.5\t2\tchr21:42839824/42839825(+)->chr21:42840323/42840324(-):(discordant_mates:4)\n+chr21\t39770910\t+\t2\t0\tchr21\t39771416\t-\t0\t2\t506\tentropy=0.0<0.6827,n_support=1<7\tcircular\tintronic\t2\t2\t1\t0\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t75\t0\t1\t0\t0\t48\t0\t1\t0\t0\t0.5\t2\tchr21:39770910/39770911(+)->chr21:39771416/39771417(-):(spanning_singleton_1_r:1,spanning_singleton_2_r:1)\n+chr21\t39860702\t-\t4\t0\tchr21\t42843890\t+\t0\t4\t2983188\tentropy=0.0<0.6827,n_support=2<7\tlinear\tintronic\t2\t0\t0\t2\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t40\t0\t2\tchr21:39860702/39860703(-)->chr21:42843890/42843891(+):(discordant_mates:4)\n+chr21\t39781087\t+\t0\t2\tchr21\t39782101\t-\t2\t0\t1014\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39781087/39781088(+)->chr21:39782101/39782102(-):(discordant_mates:2)\n+chr21\t39846131\t+\t0\t2\tchr21\t39846834\t-\t2\t0\t703\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t39\t0\t1\t0\t20\t0\t2\tchr21:39846131/39846132(+)->chr21:39846834/39846835(-):(discordant_mates:2)\n+chr21\t39853321\t-\t2\t0\tchr21\t42851610\t+\t0\t2\t2998289\tentropy=0.0<0.6826,n_support=1<7\tlinear\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t126\t0\t1\t0\t20\t0\t2\tchr21:39853321/39853322(-)->chr21:42851610/42851611(+):(discordant_mates:2)\n+chr21\t42840341\t+\t0\t2\tchr21\t42842586\t-\t2\t0\t2245\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t124\t0\t1\t0\t0\t124\t0\t1\t0\t20\t0\t2\tchr21:42840341/42840342(+)->chr21:42842586/42842587(-):(discordant_mates:2)\n+chr21\t42843861\t+\t2\t0\tchr21\t42845272\t-\t0\t2\t1411\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t65\t0\t1\t0\t0\t60\t0\t1\t0\t20\t0\t2\tchr21:42843861/42843862(+)->chr21:42845272/42845273(-):(discordant_mates:2)\n+chr21\t42863757\t+\t0\t2\tchr21\t42864220\t-\t2\t0\t463\tentropy=0.0<0.6826,n_support=1<7\tcircular\tintronic\t1\t0\t0\t1\t10000\t0\t1\t1\t1\t0\t0\t0\t0\t0\t0\t0\t126\t0\t1\t0\t0\t123\t0\t1\t0\t20\t0\t2\tchr21:42863757/42863758(+)->chr21:42864220/42864221(-):(discordant_mates:2)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/detect.txt
--- a/test-data/detect.txt Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/detect.txt Tue Sep 26 08:27:47 2017 -0400
b
@@ -1,2 +1,2 @@
-chr-A pos-A direction-A pos-A-acceptor pos-A-donor chr-B pos-B direction-B pos-B-acceptor pos-B-donor genomic-distance filter-status circRNA intronic/exonic score soft+hardclips n-split-reads n-discordant-reads n-edges n-nodes-A n-nodes-B n-splice-junc-A n-splice-junc-B entropy-bp-edge entropy-all-edges bp-pos-stddev entropy-disco-bps lr-A-slope lr-A-intercept lr-A-rvalue lr-A-pvalue lr-A-stderr lr-B-slope lr-B-intercept lr-B-rvalue lr-B-pvalue lr-B-stderr disco/split clips/score nodes/edge data-structure
-chr21 39877811 - 6 0 chr21 42873374 + 0 6 2995563 unclassified linear intronic 9 6 3 0 1 1 1 0 0 1.0000 1.0000 0.0000 0.0000 17.0000 37.3333 0.9994 0.0216 0.5774 24.0000 42.3333 0.9954 0.0611 2.3094 0.0000 0.3333 2.0000 chr21:39877811/39877812(-)->chr21:42873374/42873375(+):(spanning_paired_1_t:3,spanning_paired_2_t:3)
+chr-A pos-A direction-A pos-A-acceptor pos-A-donor chr-B pos-B direction-B pos-B-acceptor pos-B-donor genomic-distance filter-status circRNA intronic/exonic score soft+hardclips n-split-reads n-discordant-reads alignment-score mismatches n-edges n-nodes-A n-nodes-B n-splice-junc-A n-splice-junc-B entropy-bp-edge entropy-all-edges bp-pos-stddev entropy-disco-bps lr-A-slope lr-A-intercept lr-A-rvalue lr-A-pvalue lr-A-stderr lr-B-slope lr-B-intercept lr-B-rvalue lr-B-pvalue lr-B-stderr disco/split clips/score nodes/edge data-structure
+chr21 39877811 - 6 0 chr21 42873374 + 0 6 2995563 unclassified linear intronic 9 6 3 0 642 0 1 1 1 0 0 1.0000 1.0000 0.0000 0.0000 17.0000 37.3333 0.9994 0.0216 0.5774 24.0000 42.3333 0.9954 0.0611 2.3094 0.0000 0.3333 2.0000 chr21:39877811/39877812(-)->chr21:42873374/42873375(+):(spanning_paired_1_t:3,spanning_paired_2_t:3)
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/fix_test-01.fixed.bam
b
Binary file test-data/fix_test-01.fixed.bam has changed
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/integrate_test_terg_s041.in.dbed
--- a/test-data/integrate_test_terg_s041.in.dbed Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/integrate_test_terg_s041.in.dbed Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,6 +1,6 @@\n-chr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n-chr21\t38445621\t-\t200\t0\tchr21\t41498118\t+\t0\t200\t3052497\tvalid\tlinear\texonic\t588\t378\t189\t21\t14\t7\t9\t2\t4\t0.7421\t0.8197\t0\t0.7162\t0.7344\t19.9882\t0.9657\t0\t0.0178\t0.7194\t21.976\t0.9648\t0\t0.0177\t0.1111\t0.3214\t1.1429\tchr21:38445621/38445622(-)->chr21:41498118/41498119(+):(discordant_mates:14,spanning_paired_1:114,spanning_paired_1_t:4,spanning_paired_2:114,spanning_paired_2_t:4)&chr21:38445621/38445622(-)->chr21:41508080/41508081(+):(spanning_paired_1:33,spanning_paired_2:33)&chr21:38445621/38445622(-)->chr21:41507949/41507950(+):(discordant_mates:8,spanning_paired_1:15,spanning_paired_2:15)&chr21:38474121/38474122(-)->chr21:41498118/41498119(+):(spanning_paired_1:14,spanning_paired_2:14)&chr21:38445621/38445622(-)->chr21:41480475/41480476(+):(discordant_mates:6,spanning_paired_1:2,spanning_paired_2:2)&chr21:38423561/38423562(-)->chr21:41498118/41498119(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:38445621/38445622(-)->chr21:41506444/41506445(+):(discordant_mates:6,spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:38445621/38445622(-)->chr21:41504366/41504367(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38446208/38446209(-)->chr21:41498118/41498119(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41480475/41480476(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41508080/41508081(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38423567/38423568(-)->chr21:41485824/41485825(+):(discordant_mates:4)&chr21:38445482/38445483(-)->chr21:41500346/41500347(+):(discordant_mates:2)&chr21:38474093/38474094(-)->chr21:41507966/41507967(+):(discordant_mates:2)\n-chr21\t38487350\t-\t200\t0\tchr21\t41479719\t+\t0\t200\t2992369\tvalid\tlinear\tintronic\t151\t96\t48\t7\t1\t1\t1\t0\t0\t0.869\t0.869\t0.2887\t0.8124\t2.1465\t5.6435\t0.9817\t0\t0.0571\t1.9162\t26.2818\t0.99\t0\t0.0376\t0.1458\t0.3179\t2\tchr21:38487350/38487351(-)->chr21:41479719/41479720(+):(discordant_mates:14,spanning_paired_1:41,spanning_paired_1_t:7,spanning_paired_2:41,spanning_paired_2_t:7)\n-chr21\t36338771\t+\t0\t10\tchr21\t36344707\t-\t10\t0\t5936\tvalid\tcircular\texonic\t117\t58\t38\t3\t2\t2\t2\t0\t0\t0.8536\t0.8591\t2.436\t0.871\t2.1714\t19.1829\t0.9605\t0\t0.1021\t2.0941\t32.6902\t0.9725\t0\t0.0814\t0.0789\t0.2479\t2\tchr21:36338771/36338772(+)->chr21:36344707/36344708(-):(discordant_mates:6,spanning_paired_1:14,spanning_paired_1_t:23,spanning_paired_2:14,spanning_paired_2_t:23)&chr21:36338841/36338842(+)->chr21:36341534/36341535(-):(spanning_paired_1:1,spanning_paired_2:1)\n-chr21\t29321220\t+\t10\t0\tchr21\t29329693\t-\t0\t10\t8473\tvalid\tcircular\texonic\t49\t35\t15\t4\t4\t4\t3\t0\t0\t0.8605\t0.8018\t0\t0.6667\t8.3531\t34.3077\t0.9863\t0\t0.4412\t9.4685\t14.5897\t0.9701\t0\t0.749\t0.2667\t0.3571\t1.75\tchr21:29321220/29321221(+)->chr21:29329693/29329694(-):(discordant_mates:4,spanning_paired_1:9,spanning_paired_1_t:1,spanning_paired_2:9,spanning_paired_2_t:1)&chr21:29326058/29326059(+)->chr21:29329693/29329694(-):(spanning_paired_1:4,spanning_paired_2:4)&chr21:29327324/29327325(+)->chr21:29329704/29329705(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:29321216/29321217(+)->chr21:29329647/29329648(-):(discordant_mates:4)\n-chr21\t33414887\t+\t5\t5\tchr21\t33432871\t-\t5\t5\t17984\tvalid\tcircular\texonic\t20\t12\t6\t2\t3\t3\t2\t1\t0\t0.8277\t0.9167\t0\t0.75\t12.6\t50.8\t0.7995\t0.1045\t5.4663\t12.8\t24.2\t0.8072\t0.0987\t5.4049\t0.3333\t0.3\t1.6667\tchr21:33414887/33414888(+)->chr21:33432871/33432872(-):(spanning_paired_1:2,spanning_paired_1_t:3,spanning_paired_2:2,spanning_paired_2_t:3)&chr21:33426883/33426884(+)->chr21:33432871/3343287'..b'lips\tn-split-reads\tn-discordant-reads\talignment-score\tmismatches\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tdata-structure\n+chr21\t38445621\t-\t200\t0\tchr21\t41498118\t+\t0\t200\t3052497\tvalid\tlinear\texonic\t588\t378\t189\t21\t10000\t0\t14\t7\t9\t2\t4\t0.7421\t0.8197\t0\t0.7162\t0.7344\t19.9882\t0.9657\t0\t0.0178\t0.7194\t21.976\t0.9648\t0\t0.0177\t0.1111\t0.3214\t1.1429\tchr21:38445621/38445622(-)->chr21:41498118/41498119(+):(discordant_mates:14,spanning_paired_1:114,spanning_paired_1_t:4,spanning_paired_2:114,spanning_paired_2_t:4)&chr21:38445621/38445622(-)->chr21:41508080/41508081(+):(spanning_paired_1:33,spanning_paired_2:33)&chr21:38445621/38445622(-)->chr21:41507949/41507950(+):(discordant_mates:8,spanning_paired_1:15,spanning_paired_2:15)&chr21:38474121/38474122(-)->chr21:41498118/41498119(+):(spanning_paired_1:14,spanning_paired_2:14)&chr21:38445621/38445622(-)->chr21:41480475/41480476(+):(discordant_mates:6,spanning_paired_1:2,spanning_paired_2:2)&chr21:38423561/38423562(-)->chr21:41498118/41498119(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:38445621/38445622(-)->chr21:41506444/41506445(+):(discordant_mates:6,spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:38445621/38445622(-)->chr21:41504366/41504367(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38446208/38446209(-)->chr21:41498118/41498119(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41480475/41480476(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41508080/41508081(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38423567/38423568(-)->chr21:41485824/41485825(+):(discordant_mates:4)&chr21:38445482/38445483(-)->chr21:41500346/41500347(+):(discordant_mates:2)&chr21:38474093/38474094(-)->chr21:41507966/41507967(+):(discordant_mates:2)\n+chr21\t38487350\t-\t200\t0\tchr21\t41479719\t+\t0\t200\t2992369\tvalid\tlinear\tintronic\t151\t96\t48\t7\t10000\t0\t1\t1\t1\t0\t0\t0.869\t0.869\t0.2887\t0.8124\t2.1465\t5.6435\t0.9817\t0\t0.0571\t1.9162\t26.2818\t0.99\t0\t0.0376\t0.1458\t0.3179\t2\tchr21:38487350/38487351(-)->chr21:41479719/41479720(+):(discordant_mates:14,spanning_paired_1:41,spanning_paired_1_t:7,spanning_paired_2:41,spanning_paired_2_t:7)\n+chr21\t36338771\t+\t0\t10\tchr21\t36344707\t-\t10\t0\t5936\tvalid\tcircular\texonic\t117\t58\t38\t3\t10000\t0\t2\t2\t2\t0\t0\t0.8536\t0.8591\t2.436\t0.871\t2.1714\t19.1829\t0.9605\t0\t0.1021\t2.0941\t32.6902\t0.9725\t0\t0.0814\t0.0789\t0.2479\t2\tchr21:36338771/36338772(+)->chr21:36344707/36344708(-):(discordant_mates:6,spanning_paired_1:14,spanning_paired_1_t:23,spanning_paired_2:14,spanning_paired_2_t:23)&chr21:36338841/36338842(+)->chr21:36341534/36341535(-):(spanning_paired_1:1,spanning_paired_2:1)\n+chr21\t29321220\t+\t10\t0\tchr21\t29329693\t-\t0\t10\t8473\tvalid\tcircular\texonic\t49\t35\t15\t4\t10000\t0\t4\t4\t3\t0\t0\t0.8605\t0.8018\t0\t0.6667\t8.3531\t34.3077\t0.9863\t0\t0.4412\t9.4685\t14.5897\t0.9701\t0\t0.749\t0.2667\t0.3571\t1.75\tchr21:29321220/29321221(+)->chr21:29329693/29329694(-):(discordant_mates:4,spanning_paired_1:9,spanning_paired_1_t:1,spanning_paired_2:9,spanning_paired_2_t:1)&chr21:29326058/29326059(+)->chr21:29329693/29329694(-):(spanning_paired_1:4,spanning_paired_2:4)&chr21:29327324/29327325(+)->chr21:29329704/29329705(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:29321216/29321217(+)->chr21:29329647/29329648(-):(discordant_mates:4)\n+chr21\t33414887\t+\t5\t5\tchr21\t33432871\t-\t5\t5\t17984\tvalid\tcircular\texonic\t20\t12\t6\t2\t10000\t0\t3\t3\t2\t1\t0\t0.8277\t0.9167\t0\t0.75\t12.6\t50.8\t0.7995\t0.1045\t5.4663\t12.8\t24.2\t0.8072\t0.0987\t5.4049\t0.3333\t0.3\t1.6667\tchr21:33414887/33414888(+)->chr21:33432871/33432872(-):(spanning_paired_1:2,spanning_paired_1_t:3,spanning_paired_2:2,spanning_paired_2_t:3)&chr21:33426883/33426884(+)->chr21:33432871/33432872(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:33414915/33414916(+)->chr21:33421509/33421510(-):(discordant_mates:4)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/integrate_test_terg_s041.out.no-gtf.txt
--- a/test-data/integrate_test_terg_s041.out.no-gtf.txt Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/integrate_test_terg_s041.out.no-gtf.txt Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,6 +1,6 @@\n-shared-id\tfusion\tchr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tfull-gene-disregulation\tframeshift=0\tframeshift=+1\tframeshift=+2\tdata-structure\n-1\tTMPRSS2->ERG\tchr21\t38487350\t-\t200\t0\tchr21\t41479719\t+\t0\t200\t2992369\tvalid\tlinear\tintronic\t151\t96\t48\t7\t1\t1\t1\t0\t0\t0.869\t0.869\t0.2887\t0.8124\t2.1465\t5.6435\t0.9817\t0\t0.0571\t1.9162\t26.2818\t0.99\t0\t0.0376\t0.1458\t0.3179\t2\t\t\t\t\tchr21:38487350/38487351(-)->chr21:41479719/41479720(+):(discordant_mates:14,spanning_paired_1:41,spanning_paired_1_t:7,spanning_paired_2:41,spanning_paired_2_t:7)\n-1\tTMPRSS2->ERG\tchr21\t38445621\t-\t200\t0\tchr21\t41498118\t+\t0\t200\t3052497\tvalid\tlinear\texonic\t588\t378\t189\t21\t14\t7\t9\t2\t4\t0.7421\t0.8197\t0\t0.7162\t0.7344\t19.9882\t0.9657\t0\t0.0178\t0.7194\t21.976\t0.9648\t0\t0.0177\t0.1111\t0.3214\t1.1429\tTMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000398897.5)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000453032.6)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000398897.5)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000453032.6)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000492833.5)-ensembl_havana\tTMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000288319.11)-ensembl_havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000398905.5)-havana,TMPRSS2(ENST00000332149.9)-'..b'ENST00000454499.5)-havana(+2)->ERG(ENST00000398919.6)-ensembl_havana(+0),TMPRSS2(ENST00000454499.5)-havana(+2)->ERG(ENST00000417133.6)-ensembl_havana(+0),TMPRSS2(ENST00000454499.5)-havana(+2)->ERG(ENST00000442448.5)-ensembl(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000288319.11)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398905.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398907.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398910.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398911.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398919.6)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000417133.6)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000442448.5)-ensembl(+0)\tchr21:38445621/38445622(-)->chr21:41498118/41498119(+):(discordant_mates:14,spanning_paired_1:114,spanning_paired_1_t:4,spanning_paired_2:114,spanning_paired_2_t:4)&chr21:38445621/38445622(-)->chr21:41508080/41508081(+):(spanning_paired_1:33,spanning_paired_2:33)&chr21:38445621/38445622(-)->chr21:41507949/41507950(+):(discordant_mates:8,spanning_paired_1:15,spanning_paired_2:15)&chr21:38474121/38474122(-)->chr21:41498118/41498119(+):(spanning_paired_1:14,spanning_paired_2:14)&chr21:38445621/38445622(-)->chr21:41480475/41480476(+):(discordant_mates:6,spanning_paired_1:2,spanning_paired_2:2)&chr21:38423561/38423562(-)->chr21:41498118/41498119(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:38445621/38445622(-)->chr21:41506444/41506445(+):(discordant_mates:6,spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:38445621/38445622(-)->chr21:41504366/41504367(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38446208/38446209(-)->chr21:41498118/41498119(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41480475/41480476(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41508080/41508081(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38423567/38423568(-)->chr21:41485824/41485825(+):(discordant_mates:4)&chr21:38445482/38445483(-)->chr21:41500346/41500347(+):(discordant_mates:2)&chr21:38474093/38474094(-)->chr21:41507966/41507967(+):(discordant_mates:2)\n+2\tchr21:36338771->chr21:36344707\tchr21\t36338771\t+\t0\t10\tchr21\t36344707\t-\t10\t0\t5936\tvalid\tcircular\texonic\t117\t58\t38\t3\t10000\t0\t2\t2\t2\t0\t0\t0.8536\t0.8591\t2.436\t0.871\t2.1714\t19.1829\t0.9605\t0\t0.1021\t2.0941\t32.6902\t0.9725\t0\t0.0814\t0.0789\t0.2479\t2\t\t\t\t\tchr21:36338771/36338772(+)->chr21:36344707/36344708(-):(discordant_mates:6,spanning_paired_1:14,spanning_paired_1_t:23,spanning_paired_2:14,spanning_paired_2_t:23)&chr21:36338841/36338842(+)->chr21:36341534/36341535(-):(spanning_paired_1:1,spanning_paired_2:1)\n+3\tchr21:29329693->chr21:29321220\tchr21\t29321220\t+\t10\t0\tchr21\t29329693\t-\t0\t10\t8473\tvalid\tcircular\texonic\t49\t35\t15\t4\t10000\t0\t4\t4\t3\t0\t0\t0.8605\t0.8018\t0\t0.6667\t8.3531\t34.3077\t0.9863\t0\t0.4412\t9.4685\t14.5897\t0.9701\t0\t0.749\t0.2667\t0.3571\t1.75\t\t\t\t\tchr21:29321220/29321221(+)->chr21:29329693/29329694(-):(discordant_mates:4,spanning_paired_1:9,spanning_paired_1_t:1,spanning_paired_2:9,spanning_paired_2_t:1)&chr21:29326058/29326059(+)->chr21:29329693/29329694(-):(spanning_paired_1:4,spanning_paired_2:4)&chr21:29327324/29327325(+)->chr21:29329704/29329705(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:29321216/29321217(+)->chr21:29329647/29329648(-):(discordant_mates:4)\n+4\tchr21:33432871<->chr21:33414887\tchr21\t33414887\t+\t5\t5\tchr21\t33432871\t-\t5\t5\t17984\tvalid\tcircular\texonic\t20\t12\t6\t2\t10000\t0\t3\t3\t2\t1\t0\t0.8277\t0.9167\t0\t0.75\t12.6\t50.8\t0.7995\t0.1045\t5.4663\t12.8\t24.2\t0.8072\t0.0987\t5.4049\t0.3333\t0.3\t1.6667\t\t\t\t\tchr21:33414887/33414888(+)->chr21:33432871/33432872(-):(spanning_paired_1:2,spanning_paired_1_t:3,spanning_paired_2:2,spanning_paired_2_t:3)&chr21:33426883/33426884(+)->chr21:33432871/33432872(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:33414915/33414916(+)->chr21:33421509/33421510(-):(discordant_mates:4)\n'
b
diff -r 922b52cf9b35 -r c3ce472a315f test-data/integrate_test_terg_s041.out.txt
--- a/test-data/integrate_test_terg_s041.out.txt Thu Sep 14 08:35:12 2017 -0400
+++ b/test-data/integrate_test_terg_s041.out.txt Tue Sep 26 08:27:47 2017 -0400
b
b'@@ -1,6 +1,6 @@\n-shared-id\tfusion\tchr-A\tpos-A\tdirection-A\tpos-A-acceptor\tpos-A-donor\tchr-B\tpos-B\tdirection-B\tpos-A-acceptor\tpos-A-donor\tgenomic-distance\tfilter-status\tcircRNA\tintronic/exonic\tscore\tsoft+hardclips\tn-split-reads\tn-discordant-reads\tn-edges\tn-nodes-A\tn-nodes-B\tn-splice-junc-A\tn-splice-junc-B\tentropy-bp-edge\tentropy-all-edges\tbp-pos-stddev\tentropy-disco-bps\tlr-A-slope\tlr-A-intercept\tlr-A-rvalue\tlr-A-pvalue\tlr-A-stderr\tlr-B-slope\tlr-B-intercept\tlr-B-rvalue\tlr-B-pvalue\tlr-B-stderr\tdisco/split\tclips/score\tnodes/edge\tfull-gene-disregulation\tframeshift=0\tframeshift=+1\tframeshift=+2\tdata-structure\n-1\tTMPRSS2->ERG\tchr21\t38487350\t-\t200\t0\tchr21\t41479719\t+\t0\t200\t2992369\tvalid\tlinear\tintronic\t151\t96\t48\t7\t1\t1\t1\t0\t0\t0.869\t0.869\t0.2887\t0.8124\t2.1465\t5.6435\t0.9817\t0\t0.0571\t1.9162\t26.2818\t0.99\t0\t0.0376\t0.1458\t0.3179\t2\t\t\t\t\tchr21:38487350/38487351(-)->chr21:41479719/41479720(+):(discordant_mates:14,spanning_paired_1:41,spanning_paired_1_t:7,spanning_paired_2:41,spanning_paired_2_t:7)\n-1\tTMPRSS2->ERG\tchr21\t38445621\t-\t200\t0\tchr21\t41498118\t+\t0\t200\t3052497\tvalid\tlinear\texonic\t588\t378\t189\t21\t14\t7\t9\t2\t4\t0.7421\t0.8197\t0\t0.7162\t0.7344\t19.9882\t0.9657\t0\t0.0178\t0.7194\t21.976\t0.9648\t0\t0.0177\t0.1111\t0.3214\t1.1429\tTMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000424093.5)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000455813.1)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000458356.5)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000398897.5)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000453032.6)-ensembl_havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000463138.1)-havana->ERG(ENST00000492833.5)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000398897.5)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000429727.6)-ensembl,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000453032.6)-ensembl_havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000468474.5)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000473107.1)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000481609.5)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000485493.1)-havana,TMPRSS2(ENST00000497881.5)-havana->ERG(ENST00000492833.5)-ensembl_havana\tTMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000288319.11)-ensembl_havana,TMPRSS2(ENST00000332149.9)-ensembl_havana->ERG(ENST00000398905.5)-havana,TMPRSS2(ENST00000332149.9)-'..b'ENST00000454499.5)-havana(+2)->ERG(ENST00000398919.6)-ensembl_havana(+0),TMPRSS2(ENST00000454499.5)-havana(+2)->ERG(ENST00000417133.6)-ensembl_havana(+0),TMPRSS2(ENST00000454499.5)-havana(+2)->ERG(ENST00000442448.5)-ensembl(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000288319.11)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398905.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398907.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398910.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398911.5)-havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000398919.6)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000417133.6)-ensembl_havana(+0),TMPRSS2(ENST00000458356.5)-havana(+2)->ERG(ENST00000442448.5)-ensembl(+0)\tchr21:38445621/38445622(-)->chr21:41498118/41498119(+):(discordant_mates:14,spanning_paired_1:114,spanning_paired_1_t:4,spanning_paired_2:114,spanning_paired_2_t:4)&chr21:38445621/38445622(-)->chr21:41508080/41508081(+):(spanning_paired_1:33,spanning_paired_2:33)&chr21:38445621/38445622(-)->chr21:41507949/41507950(+):(discordant_mates:8,spanning_paired_1:15,spanning_paired_2:15)&chr21:38474121/38474122(-)->chr21:41498118/41498119(+):(spanning_paired_1:14,spanning_paired_2:14)&chr21:38445621/38445622(-)->chr21:41480475/41480476(+):(discordant_mates:6,spanning_paired_1:2,spanning_paired_2:2)&chr21:38423561/38423562(-)->chr21:41498118/41498119(+):(spanning_paired_1:2,spanning_paired_2:2)&chr21:38445621/38445622(-)->chr21:41506444/41506445(+):(discordant_mates:6,spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:38445621/38445622(-)->chr21:41504366/41504367(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38446208/38446209(-)->chr21:41498118/41498119(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41480475/41480476(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38474121/38474122(-)->chr21:41508080/41508081(+):(spanning_paired_1:1,spanning_paired_2:1)&chr21:38423567/38423568(-)->chr21:41485824/41485825(+):(discordant_mates:4)&chr21:38445482/38445483(-)->chr21:41500346/41500347(+):(discordant_mates:2)&chr21:38474093/38474094(-)->chr21:41507966/41507967(+):(discordant_mates:2)\n+2\tchr21:36338771->chr21:36344707\tchr21\t36338771\t+\t0\t10\tchr21\t36344707\t-\t10\t0\t5936\tvalid\tcircular\texonic\t117\t58\t38\t3\t10000\t0\t2\t2\t2\t0\t0\t0.8536\t0.8591\t2.436\t0.871\t2.1714\t19.1829\t0.9605\t0\t0.1021\t2.0941\t32.6902\t0.9725\t0\t0.0814\t0.0789\t0.2479\t2\t\t\t\t\tchr21:36338771/36338772(+)->chr21:36344707/36344708(-):(discordant_mates:6,spanning_paired_1:14,spanning_paired_1_t:23,spanning_paired_2:14,spanning_paired_2_t:23)&chr21:36338841/36338842(+)->chr21:36341534/36341535(-):(spanning_paired_1:1,spanning_paired_2:1)\n+3\tchr21:29329693->chr21:29321220\tchr21\t29321220\t+\t10\t0\tchr21\t29329693\t-\t0\t10\t8473\tvalid\tcircular\texonic\t49\t35\t15\t4\t10000\t0\t4\t4\t3\t0\t0\t0.8605\t0.8018\t0\t0.6667\t8.3531\t34.3077\t0.9863\t0\t0.4412\t9.4685\t14.5897\t0.9701\t0\t0.749\t0.2667\t0.3571\t1.75\t\t\t\t\tchr21:29321220/29321221(+)->chr21:29329693/29329694(-):(discordant_mates:4,spanning_paired_1:9,spanning_paired_1_t:1,spanning_paired_2:9,spanning_paired_2_t:1)&chr21:29326058/29326059(+)->chr21:29329693/29329694(-):(spanning_paired_1:4,spanning_paired_2:4)&chr21:29327324/29327325(+)->chr21:29329704/29329705(-):(spanning_paired_1_t:1,spanning_paired_2_t:1)&chr21:29321216/29321217(+)->chr21:29329647/29329648(-):(discordant_mates:4)\n+4\tchr21:33432871<->chr21:33414887\tchr21\t33414887\t+\t5\t5\tchr21\t33432871\t-\t5\t5\t17984\tvalid\tcircular\texonic\t20\t12\t6\t2\t10000\t0\t3\t3\t2\t1\t0\t0.8277\t0.9167\t0\t0.75\t12.6\t50.8\t0.7995\t0.1045\t5.4663\t12.8\t24.2\t0.8072\t0.0987\t5.4049\t0.3333\t0.3\t1.6667\t\t\t\t\tchr21:33414887/33414888(+)->chr21:33432871/33432872(-):(spanning_paired_1:2,spanning_paired_1_t:3,spanning_paired_2:2,spanning_paired_2_t:3)&chr21:33426883/33426884(+)->chr21:33432871/33432872(-):(spanning_paired_1:1,spanning_paired_2:1)&chr21:33414915/33414916(+)->chr21:33421509/33421510(-):(discordant_mates:4)\n'