Repository 'control_freec'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/control_freec

Changeset 2:f923c54a17ee (2021-11-13)
Previous changeset 1:2c6349fb175c (2020-08-18)
Commit message:
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/freec commit 30718e775cd7cd2a7746608a0c99e86521cf8d7b"
modified:
control_freec.xml
macros.xml
added:
test-data/output/sample.bam_info_woc.txt
test-data/output/sample.bam_ratio_woc.txt
test-data/output/sample.bam_sample_woc.cpn
test-data/output/sample_woc.bam_CNVs
b
diff -r 2c6349fb175c -r f923c54a17ee control_freec.xml
--- a/control_freec.xml Tue Aug 18 08:52:45 2020 -0400
+++ b/control_freec.xml Sat Nov 13 18:28:40 2021 +0000
[
@@ -1,12 +1,10 @@
-<tool id="control_freec" name="Control-FREEC" version="@WRAPPER_VERSION@">
+<tool id="control_freec" name="Control-FREEC" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@">
     <description>detects copy-number changes and allelic imbalances</description>
     <macros>
         <import>macros.xml</import>
     </macros>
-    <requirements>
-        <requirement type="package" version="5.0.1">gawk</requirement>
-        <requirement type="package" version="@TOOL_VERSION@">control-freec</requirement>
-    </requirements>
+    <expand macro="xrefs"/>
+    <expand macro="requirements"/>
     <command detect_errors="exit_code"><![CDATA[
         #if str($reference_source.ref_selector) == 'history':
             ln -s '$reference_source.ref' ./genome.fa &&
@@ -30,7 +28,7 @@
         awk -F " " '/^>/ {close(F); F="chromosomes/"substr($1,2)".fa"} {print >> F}' ./genome.fa &&
 
         ln -s '$WGS_WES.input_sample_file' ./sample.bam &&
-        #if $input_control_file
+        #if $WGS_WES.input_control_file
             ln -s '$WGS_WES.input_control_file' ./control.bam &&
         #end if
 
@@ -97,9 +95,11 @@
 mateOrientation = $WGS_WES.mateOrientation_selector
 
 [control]
+#if $WGS_WES.input_control_file
 mateFile = control.bam
 inputFormat = BAM
 mateOrientation = $WGS_WES.mateOrientation_selector
+#end if
 
 #if int($WGS_WES.advanced_settings.window_section.window) == 0
 [target]
@@ -115,7 +115,7 @@
                 <option value="other">other method</option>
             </param>
             <when value="WGS">
-                <expand macro="def_input_files" />
+                <expand macro="wgs_input_files" />
                 <section name="advanced_settings" title="Advanced WGS settings" expanded="false">
                     <expand macro="WGS" />
                     <expand macro="shared" />
@@ -141,7 +141,9 @@
     </inputs>
     <outputs>
         <data name="out_sample_raw" format="tabular" label="${tool.name} on ${on_string}: Raw copy number profiles (sample)" from_work_dir="output/sample.bam_sample.cpn" />
-        <data name="out_control_raw" format="tabular" label="${tool.name} on ${on_string}: Raw copy number profiles (control)" from_work_dir="output/control.bam_control.cpn" />
+        <data name="out_control_raw" format="tabular" label="${tool.name} on ${on_string}: Raw copy number profiles (control)" from_work_dir="output/control.bam_control.cpn">
+            <filter>WGS_WES['input_control_file']</filter>
+        </data>
         <data name="out_sample_coord" format="tabular" label="${tool.name} on ${on_string}: Coordinates of predicted CN alterations" from_work_dir="output/sample.bam_CNVs" />
         <data name="out_sample_info" format="txt" label="${tool.name} on ${on_string}: Information about FREEC run" from_work_dir="output/sample.bam_info.txt" />
         <data name="out_sample_ratio" format="tabular" label="${tool.name} on ${on_string}: Ratios and predicted CN alterations for each window" from_work_dir="output/sample.bam_ratio.txt" />
@@ -257,6 +259,23 @@
             </output>
             <output name="out_gc_profile" file="output/GC_profile.targetedRegions.cnp" />
         </test>
+        <!-- Test without control file-->
+        <test expect_num_outputs="4">
+            <conditional name="reference_source">
+                <param name="ref_selector" value="history"/>
+                <param name="ref" ftype="fasta" value="genome.fasta" />
+            </conditional>
+            <param name="WGS_WES_selector" value="WGS" />
+            <param name="input_sample_file" ftype="bam" value="tumor.bam" />
+            <param name="mateOrientation_selector" value="0" />
+            <param name="minimalSubclonePresence" value="100" />
+            <param name="window" value="50" />
+            <param name="step" value="10" />
+            <output name="out_sample_raw" file="output/sample.bam_sample_woc.cpn" />
+            <output name="out_sample_coord" file="output/sample_woc.bam_CNVs" />
+            <output name="out_sample_info" file="output/sample.bam_info_woc.txt" />
+            <output name="out_sample_ratio" file="output/sample.bam_ratio_woc.txt" />
+        </test>
     </tests>
     <help><![CDATA[
 Control-FREEC
b
diff -r 2c6349fb175c -r f923c54a17ee macros.xml
--- a/macros.xml Tue Aug 18 08:52:45 2020 -0400
+++ b/macros.xml Sat Nov 13 18:28:40 2021 +0000
b
@@ -1,6 +1,12 @@
 <macros>
-    <token name="@WRAPPER_VERSION@">@TOOL_VERSION@+galaxy1</token>
+    <token name="@VERSION_SUFFIX@">2</token>
     <token name="@TOOL_VERSION@">11.6</token>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="5.0.1">gawk</requirement>
+            <requirement type="package" version="@TOOL_VERSION@">control-freec</requirement>
+        </requirements>
+    </xml>
     <xml name="reference_interface">
         <conditional name="reference_source">
             <param name="ref_selector" type="select" label="Choose the source for the reference genome">
@@ -20,6 +26,11 @@
             </when>
         </conditional>
     </xml>
+    <xml name="xrefs">
+        <xrefs>
+            <xref type="bio.tools">freec</xref>
+        </xrefs>
+    </xml>
     <xml name="WGS">
         <param name="degree" type="select" label="Degree of polynomial" help="">
             <option value="34" selected="True">GC-content based normalization, WGS (3&amp;4)</option>
@@ -115,7 +126,7 @@
             <param name="sex" type="select" label="Sample sex" help="&quot;XX&quot; will exclude chr Y from the analysis. &quot;XY&quot; will not annotate one copy of chr X and Y as a loss.">
                 <option value="XY" selected="True">XY</option>
                 <option value="XX">XX</option>
-                <option value="no value"></option>
+                <option value="no value">Do not expecify sex</option>
             </param>
             <param name="telocentromeric" type="integer" value="50000" label="Length of pre-telomeric and pre-centromeric regions" help="Control-FREEC will not output small CNAs and LOH found within these regions (they are likely to be false because of mappability/genome assembly issues). 50000 is OK for human/mouse genomes. Use smaller values for yeasts and flies. Do not set anything for exome-seq data." />
             <!-- uniqueMatch if true -> gemMappabilityFile is needed-->
@@ -139,13 +150,17 @@
         <param name="input_sample_file" type="data" format="bam" multiple="false" label="Sample file" help="Sample file in .BAM format." />
         <param name="input_control_file" type="data" format="bam" multiple="false" label="Control file" help="Control file in .BAM format." />
     </xml>
+    <xml name="input_control_sample_wgs">
+        <param name="input_sample_file" type="data" format="bam" multiple="false" label="Sample file" help="Sample file in .BAM format." />
+        <param name="input_control_file" type="data" format="bam" multiple="false" optional="true" label="Control file" help="Control file in .BAM format." />
+    </xml>
     <xml name="wes_input_files">
         <expand macro="input_control_sample" />
         <param name="input_capture_file" optional="true" type="data" format="bed" multiple="false" label="BED file with capture regions" help="Capture regions in .BED format." />
         <expand macro="mateOrientation_selector" />
     </xml>
-    <xml name="def_input_files">
-        <expand macro="input_control_sample" />
+    <xml name="wgs_input_files">
+        <expand macro="input_control_sample_wgs" />
         <expand macro="mateOrientation_selector" />
     </xml>
     <xml name="output_section">
b
diff -r 2c6349fb175c -r f923c54a17ee test-data/output/sample.bam_info_woc.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output/sample.bam_info_woc.txt Sat Nov 13 18:28:40 2021 +0000
b
@@ -0,0 +1,13 @@
+Program_Version v11.6
+Sample_Name sample.bam
+Control_Used False
+CGcontent_Used True
+Mappability_Used False
+Looking_For_Subclones False
+Breakpoint_Threshold 0.8
+Window 50
+Number_Of_Reads|Pairs_In_Sample 657
+Number_Of_Reads|Pairs_In_Control 0
+Output_Ploidy 2
+Sample_Purity 1
+Good_Polynomial_Fit True
b
diff -r 2c6349fb175c -r f923c54a17ee test-data/output/sample.bam_ratio_woc.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output/sample.bam_ratio_woc.txt Sat Nov 13 18:28:40 2021 +0000
b
b'@@ -0,0 +1,1659 @@\n+Chromosome\tStart\tRatio\tMedianRatio\tCopyNumber\n+M\t1\t2.45225\t0.489195\t2\n+M\t11\t0.999185\t0.489195\t2\n+M\t21\t0.485727\t0.489195\t2\n+M\t31\t0.681691\t0.489195\t2\n+M\t41\t0\t0.489195\t2\n+M\t51\t0.578857\t0.489195\t2\n+M\t61\t0\t0.489195\t2\n+M\t71\t3.1333\t0.489195\t2\n+M\t81\t1.36338\t0.489195\t2\n+M\t91\t1.15771\t0.489195\t2\n+M\t101\t0.49757\t0.489195\t2\n+M\t111\t0.519271\t0.489195\t2\n+M\t121\t0\t0.489195\t2\n+M\t131\t0\t0.489195\t2\n+M\t141\t0\t0.489195\t2\n+M\t151\t0\t0.489195\t2\n+M\t161\t0\t0.489195\t2\n+M\t171\t0\t0.489195\t2\n+M\t181\t0\t0.489195\t2\n+M\t191\t0\t0.489195\t2\n+M\t201\t0\t0.489195\t2\n+M\t211\t0\t0.489195\t2\n+M\t221\t0\t0.489195\t2\n+M\t231\t0\t0.489195\t2\n+M\t241\t0\t0.489195\t2\n+M\t251\t0\t0.489195\t2\n+M\t261\t0\t0.489195\t2\n+M\t271\t0\t0.489195\t2\n+M\t281\t0\t0.489195\t2\n+M\t291\t0\t0.489195\t2\n+M\t301\t0\t0.489195\t2\n+M\t311\t0\t0.489195\t2\n+M\t321\t0\t0.489195\t2\n+M\t331\t0.499592\t0.489195\t2\n+M\t341\t0.979906\t0.489195\t2\n+M\t351\t1.0361\t0.489195\t2\n+M\t361\t1.36311\t0.489195\t2\n+M\t371\t2.04064\t0.489195\t2\n+M\t381\t1.46531\t0.489195\t2\n+M\t391\t1.55414\t0.489195\t2\n+M\t401\t1.53048\t0.489195\t2\n+M\t411\t2.44976\t0.485727\t1\n+M\t421\t2.93972\t0.485727\t1\n+M\t431\t3.39045\t0.485727\t1\n+M\t441\t3.6349\t0.485727\t1\n+M\t451\t3.99674\t0.485727\t1\n+M\t461\t4.37154\t0.485727\t1\n+M\t471\t5.343\t0.485727\t1\n+M\t481\t7.52514\t0.485727\t1\n+M\t491\t29.7664\t0.485727\t1\n+M\t501\t39.1663\t0.485727\t1\n+M\t511\t25.2611\t0.485727\t1\n+M\t521\t19.1023\t0.485727\t1\n+M\t531\t34.7662\t0.485727\t1\n+M\t541\t26.3712\t0.485727\t1\n+M\t551\t27.6864\t0.485727\t1\n+M\t561\t31.357\t0.485727\t1\n+M\t571\t35.201\t0.485727\t1\n+M\t581\t37.2995\t0.485727\t1\n+M\t591\t36.9698\t0.485727\t1\n+M\t601\t37.2995\t0.485727\t1\n+M\t611\t34.6908\t0.485727\t1\n+M\t621\t32.8269\t0.485727\t1\n+M\t631\t30.0468\t0.485727\t1\n+M\t641\t31.119\t0.485727\t1\n+M\t651\t30.637\t0.485727\t1\n+M\t661\t29.5892\t0.485727\t1\n+M\t671\t29.0791\t0.485727\t1\n+M\t681\t26.6393\t0.485727\t1\n+M\t691\t25.6706\t0.485727\t1\n+M\t701\t24.2175\t0.485727\t1\n+M\t711\t22.0479\t0.485727\t1\n+M\t721\t16.9523\t0.485727\t1\n+M\t731\t13.5618\t0.485727\t1\n+M\t741\t8.74308\t0.485727\t1\n+M\t751\t8.40371\t0.485727\t1\n+M\t761\t9.54367\t0.485727\t1\n+M\t771\t9.54367\t0.485727\t1\n+M\t781\t3.98056\t0.485727\t1\n+M\t791\t4.37154\t0.485727\t1\n+M\t801\t4.49633\t0.485727\t1\n+M\t811\t4.08117\t0.485727\t1\n+M\t821\t2.46466\t0.485727\t1\n+M\t831\t3.11563\t0.485727\t1\n+M\t841\t1.55781\t0.485727\t1\n+M\t851\t0.968701\t0.485727\t1\n+M\t861\t0.578857\t0.485727\t1\n+M\t871\t1.36338\t0.485727\t1\n+M\t881\t0.971454\t0.485727\t1\n+M\t891\t2.04507\t0.485727\t1\n+M\t901\t0.968701\t0.485727\t1\n+M\t911\t0.999185\t0.485727\t1\n+M\t921\t0.968701\t0.485727\t1\n+M\t931\t0.979906\t0.485727\t1\n+M\t941\t0.499592\t0.485727\t1\n+M\t951\t1.0361\t0.485727\t1\n+M\t961\t1.0361\t0.485727\t1\n+M\t971\t0.976874\t0.485727\t1\n+M\t981\t1.03854\t0.485727\t1\n+M\t991\t0.976874\t0.485727\t1\n+M\t1001\t0.519271\t0.485727\t1\n+M\t1011\t0.510147\t0.485727\t1\n+M\t1021\t0\t0.485727\t1\n+M\t1031\t0\t0.485727\t1\n+M\t1041\t0\t0.485727\t1\n+M\t1051\t0\t0.485727\t1\n+M\t1061\t0\t0.485727\t1\n+M\t1071\t0\t0.485727\t1\n+M\t1081\t0\t0.485727\t1\n+M\t1091\t0\t0.485727\t1\n+M\t1101\t0\t0.485727\t1\n+M\t1111\t0\t0.485727\t1\n+M\t1121\t0\t0.485727\t1\n+M\t1131\t0\t0.485727\t1\n+M\t1141\t0\t0.485727\t1\n+M\t1151\t0.485727\t0.485727\t1\n+M\t1161\t0.485727\t0.485727\t1\n+M\t1171\t1.05046\t0.485727\t1\n+M\t1181\t0.968701\t0.485727\t1\n+M\t1191\t0.971454\t0.485727\t1\n+M\t1201\t0.49757\t0.485727\t1\n+M\t1211\t0.48435\t0.485727\t1\n+M\t1221\t0\t0.485727\t1\n+M\t1231\t0\t0.485727\t1\n+M\t1241\t0.971454\t0.485727\t1\n+M\t1251\t0.979906\t0.485727\t1\n+M\t1261\t0.968701\t0.485727\t1\n+M\t1271\t0.979906\t0.485727\t1\n+M\t1281\t0.968701\t0.485727\t1\n+M\t1291\t0\t0.485727\t1\n+M\t1301\t0\t0.485727\t1\n+M\t1311\t0\t0.485727\t1\n+M\t1321\t0\t0.485727\t1\n+M\t1331\t0\t0.485727\t1\n+M\t1341\t0\t0.485727\t1\n+M\t1351\t0\t0.485727\t1\n+M\t1361\t0.518048\t0.485727\t1\n+M\t1371\t0.510159\t0.485727\t1\n+M\t1381\t1.0361\t0.485727\t1\n+M\t1391\t1.53048\t0.485727\t1\n+M\t1401\t1.55414\t0.485727\t1\n+M\t1411\t0.999185\t0.485727\t1\n+M\t1421\t0.979906\t0.485727\t1\n+M\t1431\t0.485727\t0.485727\t1\n+M\t1441\t0.681691\t0.485727\t1\n+M\t1451\t23.3161\t0.485727\t1\n+M\t1461\t1.56665\t0.485727\t1\n+M\t1471\t1.15771\t0.485727\t1\n+M\t1481\t2.44976\t0.485727\t1\n+M\t1491\t2.5508\t0.485727\t1\n+M\t1501\t2.05388\t0.485727\t1\n+M\t1511\t2.44219\t0.485727\t1\n+M\t1521\t1.95375\t0.485727\t1\n+M\t1531\t1.0361\t0.485727\t1\n+M\t1541\t1.03854\t0.485727\t1\n+M\t1551\t0.968701\t0.485727\t1\n+M\t1561\t0.968701\t0.485727\t1\n+M\t1571\t0.979906\t0.485727\t1\n+M\t1581\t0.485727\t0.485727'..b'4991\t0.979906\t0.485727\t1\n+M\t15001\t0.99514\t0.485727\t1\n+M\t15011\t0.999185\t0.485727\t1\n+M\t15021\t0.968701\t0.485727\t1\n+M\t15031\t0.489953\t0.485727\t1\n+M\t15041\t0.489953\t0.485727\t1\n+M\t15051\t0.510147\t0.485727\t1\n+M\t15061\t0\t0.485727\t1\n+M\t15071\t0\t0.485727\t1\n+M\t15081\t0\t0.485727\t1\n+M\t15091\t0\t0.485727\t1\n+M\t15101\t0.49757\t0.485727\t1\n+M\t15111\t0.485727\t0.485727\t1\n+M\t15121\t0.48435\t0.485727\t1\n+M\t15131\t0.525232\t0.485727\t1\n+M\t15141\t0.968701\t0.485727\t1\n+M\t15151\t0.510159\t0.485727\t1\n+M\t15161\t0.499592\t0.485727\t1\n+M\t15171\t0.48435\t0.485727\t1\n+M\t15181\t0.518048\t0.485727\t1\n+M\t15191\t0\t0.485727\t1\n+M\t15201\t0\t0.485727\t1\n+M\t15211\t0\t0.485727\t1\n+M\t15221\t0\t0.485727\t1\n+M\t15231\t0.485727\t0.485727\t1\n+M\t15241\t0.971454\t0.485727\t1\n+M\t15251\t0.979906\t0.485727\t1\n+M\t15261\t0.999185\t0.485727\t1\n+M\t15271\t1.0361\t0.485727\t1\n+M\t15281\t0.499592\t0.485727\t1\n+M\t15291\t0\t0.485727\t1\n+M\t15301\t0\t0.485727\t1\n+M\t15311\t0\t0.485727\t1\n+M\t15321\t0\t0.485727\t1\n+M\t15331\t0\t0.485727\t1\n+M\t15341\t0\t0.485727\t1\n+M\t15351\t0\t0.485727\t1\n+M\t15361\t0\t0.485727\t1\n+M\t15371\t0\t0.485727\t1\n+M\t15381\t0\t0.485727\t1\n+M\t15391\t0\t0.485727\t1\n+M\t15401\t0\t0.485727\t1\n+M\t15411\t0\t0.485727\t1\n+M\t15421\t0.489953\t0.485727\t1\n+M\t15431\t0.499592\t0.485727\t1\n+M\t15441\t1.02029\t0.485727\t1\n+M\t15451\t0.999185\t0.485727\t1\n+M\t15461\t0.979906\t0.485727\t1\n+M\t15471\t0.489953\t0.485727\t1\n+M\t15481\t0.578857\t0.485727\t1\n+M\t15491\t0.525232\t0.485727\t1\n+M\t15501\t0.99514\t0.485727\t1\n+M\t15511\t1.99028\t0.485727\t1\n+M\t15521\t2.48785\t0.485727\t1\n+M\t15531\t3.42967\t0.485727\t1\n+M\t15541\t3.42967\t0.485727\t1\n+M\t15551\t2.9061\t0.485727\t1\n+M\t15561\t2.04064\t0.485727\t1\n+M\t15571\t1.45718\t0.485727\t1\n+M\t15581\t0.49757\t0.485727\t2\n+M\t15591\t0\t0.485727\t2\n+M\t15601\t0\t0.485727\t2\n+M\t15611\t0\t0.485727\t2\n+M\t15621\t0\t0.485727\t2\n+M\t15631\t0\t0.485727\t2\n+M\t15641\t0\t0.485727\t2\n+M\t15651\t0\t0.485727\t2\n+M\t15661\t0\t0.485727\t2\n+M\t15671\t0\t0.485727\t2\n+M\t15681\t0\t0.485727\t2\n+M\t15691\t0.499592\t0.485727\t2\n+M\t15701\t0.485727\t0.485727\t2\n+M\t15711\t0.999185\t0.485727\t2\n+M\t15721\t0.971454\t0.485727\t2\n+M\t15731\t0.99514\t0.485727\t2\n+M\t15741\t0.489953\t0.485727\t2\n+M\t15751\t0.499592\t0.485727\t2\n+M\t15761\t0\t0.485727\t2\n+M\t15771\t0\t0.485727\t2\n+M\t15781\t0\t0.485727\t2\n+M\t15791\t0\t0.485727\t2\n+M\t15801\t0\t0.485727\t2\n+M\t15811\t0\t0.485727\t2\n+M\t15821\t0\t0.485727\t2\n+M\t15831\t0\t0.485727\t2\n+M\t15841\t0\t0.485727\t2\n+M\t15851\t0\t0.485727\t2\n+M\t15861\t0\t0.485727\t2\n+M\t15871\t0\t0.485727\t2\n+M\t15881\t0\t0.485727\t2\n+M\t15891\t0\t0.485727\t2\n+M\t15901\t0\t0.485727\t2\n+M\t15911\t0\t0.485727\t2\n+M\t15921\t0\t0.485727\t2\n+M\t15931\t0\t0.485727\t2\n+M\t15941\t0\t0.485727\t2\n+M\t15951\t0\t0.485727\t2\n+M\t15961\t0\t0.485727\t2\n+M\t15971\t0\t0.485727\t2\n+M\t15981\t0\t0.485727\t2\n+M\t15991\t0\t0.485727\t2\n+M\t16001\t0\t0.485727\t2\n+M\t16011\t0\t0.485727\t2\n+M\t16021\t0\t0.485727\t2\n+M\t16031\t0\t0.485727\t2\n+M\t16041\t0\t0.485727\t2\n+M\t16051\t0\t0.485727\t2\n+M\t16061\t0\t0.485727\t2\n+M\t16071\t0.489953\t0.485727\t2\n+M\t16081\t0.510159\t0.485727\t2\n+M\t16091\t0.519271\t0.485727\t2\n+M\t16101\t0.518048\t0.485727\t2\n+M\t16111\t0.510159\t0.485727\t2\n+M\t16121\t0.510147\t0.485727\t2\n+M\t16131\t0.518048\t0.485727\t2\n+M\t16141\t0.48435\t0.485727\t2\n+M\t16151\t1.05046\t0.485727\t2\n+M\t16161\t0.99514\t0.485727\t2\n+M\t16171\t0.525232\t0.485727\t2\n+M\t16181\t0.99514\t0.485727\t2\n+M\t16191\t1.49878\t0.485727\t2\n+M\t16201\t1.0361\t0.485727\t2\n+M\t16211\t1.02032\t0.485727\t2\n+M\t16221\t0.971454\t0.485727\t2\n+M\t16231\t0.485727\t0.485727\t2\n+M\t16241\t0\t0.485727\t2\n+M\t16251\t0.525232\t0.485727\t2\n+M\t16261\t0.485727\t0.485727\t2\n+M\t16271\t0.519271\t0.485727\t2\n+M\t16281\t0.519271\t0.485727\t2\n+M\t16291\t0.519271\t0.485727\t2\n+M\t16301\t0\t0.485727\t2\n+M\t16311\t0\t0.485727\t2\n+M\t16321\t0\t0.485727\t2\n+M\t16331\t0\t0.485727\t2\n+M\t16341\t0\t0.485727\t2\n+M\t16351\t0\t0.485727\t2\n+M\t16361\t0\t0.485727\t2\n+M\t16371\t0\t0.485727\t2\n+M\t16381\t0\t0.485727\t2\n+M\t16391\t0\t0.485727\t2\n+M\t16401\t0\t0.485727\t2\n+M\t16411\t0\t0.485727\t2\n+M\t16421\t0\t0.485727\t2\n+M\t16431\t0.902181\t0.485727\t2\n+M\t16441\t0.578857\t0.485727\t2\n+M\t16451\t0.578857\t0.485727\t2\n+M\t16461\t0.968701\t0.485727\t2\n+M\t16471\t0.971454\t0.485727\t2\n+M\t16481\t0.489953\t0.485727\t2\n+M\t16491\t0.48435\t0.485727\t2\n+M\t16501\t0.485727\t0.485727\t2\n+M\t16511\t0\t0.485727\t2\n+M\t16521\t0\t0.485727\t2\n+M\t16531\t0\t0.485727\t2\n+M\t16541\t0\t0.485727\t2\n+M\t16551\t0\t0.485727\t2\n+M\t16561\t0\t0.485727\t2\n+M\t16571\t0\t0.485727\t2\n'
b
diff -r 2c6349fb175c -r f923c54a17ee test-data/output/sample.bam_sample_woc.cpn
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output/sample.bam_sample_woc.cpn Sat Nov 13 18:28:40 2021 +0000
b
b'@@ -0,0 +1,1658 @@\n+M\t0\t49\t5\n+M\t10\t59\t2\n+M\t20\t69\t1\n+M\t30\t79\t1\n+M\t40\t89\t0\n+M\t50\t99\t1\n+M\t60\t109\t1\n+M\t70\t119\t2\n+M\t80\t129\t2\n+M\t90\t139\t2\n+M\t100\t149\t1\n+M\t110\t159\t1\n+M\t120\t169\t0\n+M\t130\t179\t0\n+M\t140\t189\t0\n+M\t150\t199\t0\n+M\t160\t209\t0\n+M\t170\t219\t0\n+M\t180\t229\t0\n+M\t190\t239\t0\n+M\t200\t249\t0\n+M\t210\t259\t0\n+M\t220\t269\t0\n+M\t230\t279\t0\n+M\t240\t289\t0\n+M\t250\t299\t0\n+M\t260\t309\t0\n+M\t270\t319\t0\n+M\t280\t329\t0\n+M\t290\t339\t0\n+M\t300\t349\t0\n+M\t310\t359\t0\n+M\t320\t369\t0\n+M\t330\t379\t1\n+M\t340\t389\t2\n+M\t350\t399\t2\n+M\t360\t409\t3\n+M\t370\t419\t4\n+M\t380\t429\t3\n+M\t390\t439\t3\n+M\t400\t449\t3\n+M\t410\t459\t5\n+M\t420\t469\t6\n+M\t430\t479\t7\n+M\t440\t489\t7\n+M\t450\t499\t8\n+M\t460\t509\t9\n+M\t470\t519\t11\n+M\t480\t529\t13\n+M\t490\t539\t19\n+M\t500\t549\t25\n+M\t510\t559\t28\n+M\t520\t569\t33\n+M\t530\t579\t51\n+M\t540\t589\t53\n+M\t550\t599\t57\n+M\t560\t609\t64\n+M\t570\t619\t69\n+M\t580\t629\t72\n+M\t590\t639\t74\n+M\t600\t649\t72\n+M\t610\t659\t68\n+M\t620\t669\t67\n+M\t630\t679\t58\n+M\t640\t689\t61\n+M\t650\t699\t59\n+M\t660\t709\t58\n+M\t670\t719\t57\n+M\t680\t729\t55\n+M\t690\t739\t53\n+M\t700\t749\t50\n+M\t710\t759\t45\n+M\t720\t769\t35\n+M\t730\t779\t28\n+M\t740\t789\t18\n+M\t750\t799\t16\n+M\t760\t809\t14\n+M\t770\t819\t14\n+M\t780\t829\t8\n+M\t790\t839\t9\n+M\t800\t849\t9\n+M\t810\t859\t8\n+M\t820\t869\t6\n+M\t830\t879\t6\n+M\t840\t889\t3\n+M\t850\t899\t2\n+M\t860\t909\t1\n+M\t870\t919\t2\n+M\t880\t929\t2\n+M\t890\t939\t3\n+M\t900\t949\t2\n+M\t910\t959\t2\n+M\t920\t969\t2\n+M\t930\t979\t2\n+M\t940\t989\t1\n+M\t950\t999\t2\n+M\t960\t1009\t2\n+M\t970\t1019\t2\n+M\t980\t1029\t2\n+M\t990\t1039\t2\n+M\t1000\t1049\t1\n+M\t1010\t1059\t1\n+M\t1020\t1069\t0\n+M\t1030\t1079\t0\n+M\t1040\t1089\t0\n+M\t1050\t1099\t0\n+M\t1060\t1109\t0\n+M\t1070\t1119\t0\n+M\t1080\t1129\t0\n+M\t1090\t1139\t0\n+M\t1100\t1149\t0\n+M\t1110\t1159\t0\n+M\t1120\t1169\t0\n+M\t1130\t1179\t0\n+M\t1140\t1189\t0\n+M\t1150\t1199\t1\n+M\t1160\t1209\t1\n+M\t1170\t1219\t2\n+M\t1180\t1229\t2\n+M\t1190\t1239\t2\n+M\t1200\t1249\t1\n+M\t1210\t1259\t1\n+M\t1220\t1269\t0\n+M\t1230\t1279\t0\n+M\t1240\t1289\t2\n+M\t1250\t1299\t2\n+M\t1260\t1309\t2\n+M\t1270\t1319\t2\n+M\t1280\t1329\t2\n+M\t1290\t1339\t0\n+M\t1300\t1349\t0\n+M\t1310\t1359\t0\n+M\t1320\t1369\t0\n+M\t1330\t1379\t0\n+M\t1340\t1389\t0\n+M\t1350\t1399\t0\n+M\t1360\t1409\t1\n+M\t1370\t1419\t1\n+M\t1380\t1429\t2\n+M\t1390\t1439\t3\n+M\t1400\t1449\t3\n+M\t1410\t1459\t2\n+M\t1420\t1469\t2\n+M\t1430\t1479\t1\n+M\t1440\t1489\t1\n+M\t1450\t1499\t1\n+M\t1460\t1509\t1\n+M\t1470\t1519\t2\n+M\t1480\t1529\t5\n+M\t1490\t1539\t5\n+M\t1500\t1549\t5\n+M\t1510\t1559\t5\n+M\t1520\t1569\t4\n+M\t1530\t1579\t2\n+M\t1540\t1589\t2\n+M\t1550\t1599\t2\n+M\t1560\t1609\t2\n+M\t1570\t1619\t2\n+M\t1580\t1629\t1\n+M\t1590\t1639\t0\n+M\t1600\t1649\t0\n+M\t1610\t1659\t0\n+M\t1620\t1669\t0\n+M\t1630\t1679\t0\n+M\t1640\t1689\t2\n+M\t1650\t1699\t3\n+M\t1660\t1709\t3\n+M\t1670\t1719\t3\n+M\t1680\t1729\t3\n+M\t1690\t1739\t1\n+M\t1700\t1749\t0\n+M\t1710\t1759\t0\n+M\t1720\t1769\t0\n+M\t1730\t1779\t0\n+M\t1740\t1789\t0\n+M\t1750\t1799\t0\n+M\t1760\t1809\t0\n+M\t1770\t1819\t0\n+M\t1780\t1829\t0\n+M\t1790\t1839\t0\n+M\t1800\t1849\t1\n+M\t1810\t1859\t2\n+M\t1820\t1869\t2\n+M\t1830\t1879\t2\n+M\t1840\t1889\t2\n+M\t1850\t1899\t1\n+M\t1860\t1909\t0\n+M\t1870\t1919\t0\n+M\t1880\t1929\t0\n+M\t1890\t1939\t0\n+M\t1900\t1949\t0\n+M\t1910\t1959\t0\n+M\t1920\t1969\t0\n+M\t1930\t1979\t0\n+M\t1940\t1989\t0\n+M\t1950\t1999\t0\n+M\t1960\t2009\t0\n+M\t1970\t2019\t0\n+M\t1980\t2029\t0\n+M\t1990\t2039\t0\n+M\t2000\t2049\t0\n+M\t2010\t2059\t0\n+M\t2020\t2069\t0\n+M\t2030\t2079\t0\n+M\t2040\t2089\t0\n+M\t2050\t2099\t0\n+M\t2060\t2109\t0\n+M\t2070\t2119\t0\n+M\t2080\t2129\t0\n+M\t2090\t2139\t0\n+M\t2100\t2149\t0\n+M\t2110\t2159\t0\n+M\t2120\t2169\t0\n+M\t2130\t2179\t0\n+M\t2140\t2189\t0\n+M\t2150\t2199\t0\n+M\t2160\t2209\t0\n+M\t2170\t2219\t0\n+M\t2180\t2229\t0\n+M\t2190\t2239\t0\n+M\t2200\t2249\t0\n+M\t2210\t2259\t1\n+M\t2220\t2269\t1\n+M\t2230\t2279\t1\n+M\t2240\t2289\t1\n+M\t2250\t2299\t1\n+M\t2260\t2309\t0\n+M\t2270\t2319\t1\n+M\t2280\t2329\t1\n+M\t2290\t2339\t1\n+M\t2300\t2349\t2\n+M\t2310\t2359\t2\n+M\t2320\t2369\t1\n+M\t2330\t2379\t1\n+M\t2340\t2389\t1\n+M\t2350\t2399\t1\n+M\t2360\t2409\t1\n+M\t2370\t2419\t1\n+M\t2380\t2429\t2\n+M\t2390\t2439\t2\n+M\t2400\t2449\t1\n+M\t2410\t2459\t1\n+M\t2420\t2469\t1\n+M\t2430\t2479\t0\n+M\t2440\t2489\t1\n+M\t2450\t2499\t1\n+M\t2460\t2509\t1\n+M\t2470\t2519\t1\n+M\t2480\t2529\t3\n+M\t2490\t2539\t3\n+M\t2500\t2549\t4\n+M\t2510\t2559\t5\n+M\t2520\t2569\t6\n+M\t2530\t2579\t4\n+M\t2540\t2589\t3\n+M\t2550\t2599\t3\n+M\t2560\t2609\t2\n+M\t2570\t2619\t1\n+M\t2580\t2629\t3\n+M\t2590\t2639\t4\n+M\t2600\t2649\t3\n+M\t2610\t2659\t3\n+M\t2620\t2669\t3\n+M\t2630\t2679\t1\n+M\t2640\t2689\t0\n+M\t2650\t2699\t0\n+M\t2660\t2709\t0\n+M\t2670\t2719\t0\n+M\t2680\t2729\t0\n+M\t2690\t2739\t0\n+M\t2700\t2749\t0\n+M\t2710\t2759\t0\n+M\t2720\t2769\t0\n+M\t2730\t2779\t0\n+M\t2740\t2789\t0\n+M\t2750\t2799\t1\n+M\t2760\t2809\t1\n+M\t2770'..b'9\t0\n+M\t14230\t14279\t0\n+M\t14240\t14289\t0\n+M\t14250\t14299\t0\n+M\t14260\t14309\t0\n+M\t14270\t14319\t0\n+M\t14280\t14329\t0\n+M\t14290\t14339\t1\n+M\t14300\t14349\t1\n+M\t14310\t14359\t2\n+M\t14320\t14369\t2\n+M\t14330\t14379\t6\n+M\t14340\t14389\t5\n+M\t14350\t14399\t5\n+M\t14360\t14409\t4\n+M\t14370\t14419\t4\n+M\t14380\t14429\t0\n+M\t14390\t14439\t0\n+M\t14400\t14449\t0\n+M\t14410\t14459\t0\n+M\t14420\t14469\t0\n+M\t14430\t14479\t0\n+M\t14440\t14489\t0\n+M\t14450\t14499\t0\n+M\t14460\t14509\t1\n+M\t14470\t14519\t1\n+M\t14480\t14529\t1\n+M\t14490\t14539\t1\n+M\t14500\t14549\t2\n+M\t14510\t14559\t1\n+M\t14520\t14569\t4\n+M\t14530\t14579\t4\n+M\t14540\t14589\t4\n+M\t14550\t14599\t3\n+M\t14560\t14609\t3\n+M\t14570\t14619\t0\n+M\t14580\t14629\t0\n+M\t14590\t14639\t2\n+M\t14600\t14649\t2\n+M\t14610\t14659\t2\n+M\t14620\t14669\t4\n+M\t14630\t14679\t7\n+M\t14640\t14689\t5\n+M\t14650\t14699\t7\n+M\t14660\t14709\t8\n+M\t14670\t14719\t7\n+M\t14680\t14729\t4\n+M\t14690\t14739\t4\n+M\t14700\t14749\t2\n+M\t14710\t14759\t2\n+M\t14720\t14769\t1\n+M\t14730\t14779\t1\n+M\t14740\t14789\t2\n+M\t14750\t14799\t5\n+M\t14760\t14809\t6\n+M\t14770\t14819\t6\n+M\t14780\t14829\t6\n+M\t14790\t14839\t5\n+M\t14800\t14849\t3\n+M\t14810\t14859\t3\n+M\t14820\t14869\t4\n+M\t14830\t14879\t6\n+M\t14840\t14889\t6\n+M\t14850\t14899\t6\n+M\t14860\t14909\t6\n+M\t14870\t14919\t5\n+M\t14880\t14929\t3\n+M\t14890\t14939\t4\n+M\t14900\t14949\t4\n+M\t14910\t14959\t2\n+M\t14920\t14969\t3\n+M\t14930\t14979\t3\n+M\t14940\t14989\t2\n+M\t14950\t14999\t1\n+M\t14960\t15009\t2\n+M\t14970\t15019\t1\n+M\t14980\t15029\t2\n+M\t14990\t15039\t2\n+M\t15000\t15049\t2\n+M\t15010\t15059\t2\n+M\t15020\t15069\t2\n+M\t15030\t15079\t1\n+M\t15040\t15089\t1\n+M\t15050\t15099\t1\n+M\t15060\t15109\t0\n+M\t15070\t15119\t0\n+M\t15080\t15129\t0\n+M\t15090\t15139\t0\n+M\t15100\t15149\t1\n+M\t15110\t15159\t1\n+M\t15120\t15169\t1\n+M\t15130\t15179\t1\n+M\t15140\t15189\t2\n+M\t15150\t15199\t1\n+M\t15160\t15209\t1\n+M\t15170\t15219\t1\n+M\t15180\t15229\t1\n+M\t15190\t15239\t0\n+M\t15200\t15249\t0\n+M\t15210\t15259\t0\n+M\t15220\t15269\t0\n+M\t15230\t15279\t1\n+M\t15240\t15289\t2\n+M\t15250\t15299\t2\n+M\t15260\t15309\t2\n+M\t15270\t15319\t2\n+M\t15280\t15329\t1\n+M\t15290\t15339\t0\n+M\t15300\t15349\t0\n+M\t15310\t15359\t0\n+M\t15320\t15369\t0\n+M\t15330\t15379\t0\n+M\t15340\t15389\t0\n+M\t15350\t15399\t0\n+M\t15360\t15409\t0\n+M\t15370\t15419\t0\n+M\t15380\t15429\t0\n+M\t15390\t15439\t0\n+M\t15400\t15449\t0\n+M\t15410\t15459\t0\n+M\t15420\t15469\t1\n+M\t15430\t15479\t1\n+M\t15440\t15489\t2\n+M\t15450\t15499\t2\n+M\t15460\t15509\t2\n+M\t15470\t15519\t1\n+M\t15480\t15529\t1\n+M\t15490\t15539\t1\n+M\t15500\t15549\t2\n+M\t15510\t15559\t4\n+M\t15520\t15569\t5\n+M\t15530\t15579\t7\n+M\t15540\t15589\t7\n+M\t15550\t15599\t6\n+M\t15560\t15609\t4\n+M\t15570\t15619\t3\n+M\t15580\t15629\t1\n+M\t15590\t15639\t0\n+M\t15600\t15649\t0\n+M\t15610\t15659\t0\n+M\t15620\t15669\t0\n+M\t15630\t15679\t0\n+M\t15640\t15689\t0\n+M\t15650\t15699\t0\n+M\t15660\t15709\t0\n+M\t15670\t15719\t0\n+M\t15680\t15729\t0\n+M\t15690\t15739\t1\n+M\t15700\t15749\t1\n+M\t15710\t15759\t2\n+M\t15720\t15769\t2\n+M\t15730\t15779\t2\n+M\t15740\t15789\t1\n+M\t15750\t15799\t1\n+M\t15760\t15809\t0\n+M\t15770\t15819\t0\n+M\t15780\t15829\t0\n+M\t15790\t15839\t0\n+M\t15800\t15849\t0\n+M\t15810\t15859\t0\n+M\t15820\t15869\t0\n+M\t15830\t15879\t0\n+M\t15840\t15889\t0\n+M\t15850\t15899\t0\n+M\t15860\t15909\t0\n+M\t15870\t15919\t0\n+M\t15880\t15929\t0\n+M\t15890\t15939\t0\n+M\t15900\t15949\t0\n+M\t15910\t15959\t0\n+M\t15920\t15969\t0\n+M\t15930\t15979\t0\n+M\t15940\t15989\t0\n+M\t15950\t15999\t0\n+M\t15960\t16009\t0\n+M\t15970\t16019\t0\n+M\t15980\t16029\t0\n+M\t15990\t16039\t0\n+M\t16000\t16049\t0\n+M\t16010\t16059\t0\n+M\t16020\t16069\t0\n+M\t16030\t16079\t0\n+M\t16040\t16089\t0\n+M\t16050\t16099\t0\n+M\t16060\t16109\t0\n+M\t16070\t16119\t1\n+M\t16080\t16129\t1\n+M\t16090\t16139\t1\n+M\t16100\t16149\t1\n+M\t16110\t16159\t1\n+M\t16120\t16169\t1\n+M\t16130\t16179\t1\n+M\t16140\t16189\t1\n+M\t16150\t16199\t2\n+M\t16160\t16209\t2\n+M\t16170\t16219\t1\n+M\t16180\t16229\t2\n+M\t16190\t16239\t3\n+M\t16200\t16249\t2\n+M\t16210\t16259\t2\n+M\t16220\t16269\t2\n+M\t16230\t16279\t1\n+M\t16240\t16289\t0\n+M\t16250\t16299\t1\n+M\t16260\t16309\t1\n+M\t16270\t16319\t1\n+M\t16280\t16329\t1\n+M\t16290\t16339\t1\n+M\t16300\t16349\t0\n+M\t16310\t16359\t0\n+M\t16320\t16369\t0\n+M\t16330\t16379\t0\n+M\t16340\t16389\t0\n+M\t16350\t16399\t0\n+M\t16360\t16409\t0\n+M\t16370\t16419\t0\n+M\t16380\t16429\t0\n+M\t16390\t16439\t0\n+M\t16400\t16449\t0\n+M\t16410\t16459\t0\n+M\t16420\t16469\t0\n+M\t16430\t16479\t1\n+M\t16440\t16489\t1\n+M\t16450\t16499\t1\n+M\t16460\t16509\t2\n+M\t16470\t16519\t2\n+M\t16480\t16529\t1\n+M\t16490\t16539\t1\n+M\t16500\t16549\t1\n+M\t16510\t16559\t0\n+M\t16520\t16569\t0\n+M\t16530\t16579\t0\n+M\t16540\t16589\t0\n+M\t16550\t16599\t0\n+M\t16560\t16609\t0\n+M\t16570\t16619\t0\n'
b
diff -r 2c6349fb175c -r f923c54a17ee test-data/output/sample_woc.bam_CNVs
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/output/sample_woc.bam_CNVs Sat Nov 13 18:28:40 2021 +0000
b
@@ -0,0 +1,1 @@
+M 410 15619 1 loss