changeset 4:11cf21ee4242 draft

Uploaded tool version 2.1.2.0.
author pjbriggs
date Wed, 12 Dec 2018 08:26:16 -0500
parents 4124781932db
children 3c435705aca5
files README.rst macs21_wrapper.xml test-data/test_MACS2.1.1_bw_html_report.zip test-data/test_MACS2.1.1_control_lambda.bdg test-data/test_MACS2.1.1_html_report.zip test-data/test_MACS2.1.1_peaks.xls test-data/test_MACS2.1.1_peaks.xls.re_match test-data/test_MACS2.1.1_peaks_narrowPeak.interval test-data/test_MACS2.1.1_summits.bed test-data/test_MACS2.1.1_treat_pileup.bdg test-data/test_MACS2.1.1_treat_pileup.bw test-data/test_MACS2.1.2_bw_html_report.zip test-data/test_MACS2.1.2_control_lambda.bdg test-data/test_MACS2.1.2_html_report.zip test-data/test_MACS2.1.2_peaks.xls test-data/test_MACS2.1.2_peaks.xls.re_match test-data/test_MACS2.1.2_peaks_narrowPeak.interval test-data/test_MACS2.1.2_summits.bed test-data/test_MACS2.1.2_treat_pileup.bdg test-data/test_MACS2.1.2_treat_pileup.bw
diffstat 20 files changed, 449 insertions(+), 444 deletions(-) [+]
line wrap: on
line diff
--- a/README.rst	Tue Mar 20 11:25:04 2018 -0400
+++ b/README.rst	Wed Dec 12 08:26:16 2018 -0500
@@ -15,7 +15,7 @@
 ======================
 
 Installation via the Galaxy Tool Shed will take of installing the tool wrapper and
-the MACS 2.1 program.
+the MACS 2.1.2 program.
 
 Manual Installation
 ===================
@@ -31,7 +31,7 @@
 
     <tool file="macs21/macs21_wrapper.xml" />
 
-You will also need to install MACS 2.1 and its dependencies:
+You will also need to install MACS 2.1.2 and its dependencies:
 
 - https://pypi.python.org/pypi/MACS2
 
@@ -59,6 +59,7 @@
 ========== ======================================================================
 Version    Changes
 ---------- ----------------------------------------------------------------------
+2.1.2.0    - Update to use MACS 2.1.2
 2.1.1.0    - Update to use MACS 2.1.1 and use conda to resolve dependencies
 2.1.0-6    - Add bedSort step into bigWig file generation; terminate with error
              when MACS finishes with non-zero exit code.
--- a/macs21_wrapper.xml	Tue Mar 20 11:25:04 2018 -0400
+++ b/macs21_wrapper.xml	Wed Dec 12 08:26:16 2018 -0500
@@ -1,9 +1,9 @@
-<tool id="macs2_1_peakcalling" name="MACS2.1.1" version="2.1.1.0">
+<tool id="macs2_1_peakcalling" name="MACS2.1.2" version="2.1.2.0">
   <description>Model-based Analysis of ChIP-Seq: peak calling</description>
   <requirements>
     <requirement type="package" version="2.7">python</requirement>
-    <requirement type="package" version="2.1.1">macs2</requirement>
-    <requirement type="package" version="3.4">R</requirement>
+    <requirement type="package" version="2.1.2">macs2</requirement>
+    <requirement type="package" version="3.5">R</requirement>
     <requirement type="package" version="357">ucsc-fetchchromsizes</requirement>
     <requirement type="package" version="357">ucsc-bedclip</requirement>
     <requirement type="package" version="357">ucsc-bedsort</requirement>
@@ -279,7 +279,7 @@
     <!-- Peak calling without bigwig output -->
     <test>
       <!-- Inputs -->
-      <param name="experiment_name" value="test_MACS2.1.1" />
+      <param name="experiment_name" value="test_MACS2.1.2" />
       <param name="broad_regions" value="" />
       <param name="format" value="BED" />
       <param name="input_chipseq_file1" value="test_region_IP.bed" dbkey="galGal3"
@@ -305,20 +305,20 @@
       <param name="nomodel_type_selector" value="nomodel" />
       <param name="nomodel_type|extsize" value="243" />
       <!-- Outputs -->
-      <output name="output_extra_files" file="test_MACS2.1.1_html_report.zip"
+      <output name="output_extra_files" file="test_MACS2.1.2_html_report.zip"
 	      compare="sim_size" delta="1500" />
-      <output name="output_summits_bed_file" file="test_MACS2.1.1_summits.bed" />
-      <output name="output_narrowpeaks_file" file="test_MACS2.1.1_peaks_narrowPeak.interval" />
+      <output name="output_summits_bed_file" file="test_MACS2.1.2_summits.bed" />
+      <output name="output_narrowpeaks_file" file="test_MACS2.1.2_peaks_narrowPeak.interval" />
       <output name="output_xls_to_interval_peaks_file"
-	      file="test_MACS2.1.1_peaks.xls.re_match"
+	      file="test_MACS2.1.2_peaks.xls.re_match"
 	      compare="re_match" lines_diff="1" />
-      <output name="output_treat_pileup_file" file="test_MACS2.1.1_treat_pileup.bdg" />
-      <output name="output_lambda_bedgraph_file" file="test_MACS2.1.1_control_lambda.bdg" />
+      <output name="output_treat_pileup_file" file="test_MACS2.1.2_treat_pileup.bdg" />
+      <output name="output_lambda_bedgraph_file" file="test_MACS2.1.2_control_lambda.bdg" />
     </test>
     <!-- Peak calling with bigwig output -->
     <test>
       <!-- Inputs -->
-      <param name="experiment_name" value="test_MACS2.1.1" />
+      <param name="experiment_name" value="test_MACS2.1.2" />
       <param name="broad_regions" value="" />
       <param name="format" value="BED" />
       <param name="input_chipseq_file1" value="test_region_IP.bed" dbkey="galGal3"
@@ -344,23 +344,23 @@
       <param name="nomodel_type_selector" value="nomodel" />
       <param name="nomodel_type|extsize" value="243" />
       <!-- Outputs -->
-      <output name="output_extra_files" file="test_MACS2.1.1_bw_html_report.zip"
+      <output name="output_extra_files" file="test_MACS2.1.2_bw_html_report.zip"
 	      compare="sim_size" delta="2500" />
-      <output name="output_summits_bed_file" file="test_MACS2.1.1_summits.bed" />
-      <output name="output_narrowpeaks_file" file="test_MACS2.1.1_peaks_narrowPeak.interval" />
+      <output name="output_summits_bed_file" file="test_MACS2.1.2_summits.bed" />
+      <output name="output_narrowpeaks_file" file="test_MACS2.1.2_peaks_narrowPeak.interval" />
       <output name="output_xls_to_interval_peaks_file"
-	      file="test_MACS2.1.1_peaks.xls.re_match"
+	      file="test_MACS2.1.2_peaks.xls.re_match"
 	      compare="re_match" lines_diff="1" />
-      <output name="output_treat_pileup_file" file="test_MACS2.1.1_treat_pileup.bdg" />
-      <output name="output_lambda_bedgraph_file" file="test_MACS2.1.1_control_lambda.bdg" />
-      <output name="output_bigwig_file" file="test_MACS2.1.1_treat_pileup.bw"
+      <output name="output_treat_pileup_file" file="test_MACS2.1.2_treat_pileup.bdg" />
+      <output name="output_lambda_bedgraph_file" file="test_MACS2.1.2_control_lambda.bdg" />
+      <output name="output_bigwig_file" file="test_MACS2.1.2_treat_pileup.bw"
 	      compare="sim_size" />
     </test>
   </tests>
   <help>
 **What it does**
 
-MACS (Model-based Analysis of ChIP-seq) 2.1.1 provides algorithms for identifying
+MACS (Model-based Analysis of ChIP-seq) 2.1.2 provides algorithms for identifying
 transcript factor binding sites. The program can be used either for ChIP-Seq data alone,
 or with control sample data to improve specificity.
 
--- a/test-data/test_MACS2.1.1_bw_html_report.zip	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,55 +0,0 @@
-<html><head><title>Additional output created by MACS (test_MACS2.1.1)</title></head><body><h3>Additional Files:</h3><p><ul>
-<li><a href="test_MACS2.1.1_peaks.xls">test_MACS2.1.1_peaks.xls</a></li>
-<li><a href="test_MACS2.1.1_treat_pileup.bdg.clipped">test_MACS2.1.1_treat_pileup.bdg.clipped</a></li>
-<li><a href="test_MACS2.1.1_treat_pileup.bdg.clipped.sorted">test_MACS2.1.1_treat_pileup.bdg.clipped.sorted</a></li>
-</ul></p>
-<h3>Messages from MACS:</h3>
-<p><pre>INFO  @ Tue, 20 Mar 2018 14:21:47: 
-# Command line: callpeak -t /tmp/tmpHxmla3/files/000/dataset_9.dat -c /tmp/tmpHxmla3/files/000/dataset_10.dat --format=BED --name=test_MACS2.1.1 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
-# ARGUMENTS LIST:
-# name = test_MACS2.1.1
-# format = BED
-# ChIP-seq file = ['/tmp/tmpHxmla3/files/000/dataset_9.dat']
-# control file = ['/tmp/tmpHxmla3/files/000/dataset_10.dat']
-# effective genome size = 7.75e+08
-# band width = 300
-# model fold = [5, 50]
-# qvalue cutoff = 5.00e-02
-# Larger dataset will be scaled towards smaller dataset.
-# Range for calculating regional lambda is: 1000 bps and 10000 bps
-# Broad region calling is off
-# Paired-End mode is off
-# MACS will save fragment pileup signal per million reads
- 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 read tag files... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 read treatment tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1.2 read input tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 tag size is determined as 50 bps 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 tag size = 50 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  total tags in treatment: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 user defined the maximum tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  tags after filtering in treatment: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  Redundant rate of treatment: 0.00 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  total tags in control: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 user defined the maximum tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  tags after filtering in control: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1  Redundant rate of control: 0.00 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #1 finished! 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Build Peak Model... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Skipped... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Use 243 as fragment length 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Call peaks... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Pre-compute pvalue-qvalue table... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3 In the peak calling step, the following will be performed simultaneously: 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3   Write bedGraph files for treatment pileup (after scaling if necessary)... test_MACS2.1.1_treat_pileup.bdg 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3   Write bedGraph files for control lambda (after scaling if necessary)... test_MACS2.1.1_control_lambda.bdg 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3   --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Call peaks for each chromosome... 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write output xls file... test_MACS2.1.1_peaks.xls 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write peak in narrowPeak format file... test_MACS2.1.1_peaks.narrowPeak 
-INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write summits bed file... test_MACS2.1.1_summits.bed 
-INFO  @ Tue, 20 Mar 2018 14:21:47: Done! 
-</pre></p>
-</body></html>
--- a/test-data/test_MACS2.1.1_control_lambda.bdg	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,153 +0,0 @@
-chr26	0	4102915	0.31355
-chr26	4102915	4103978	486.00000
-chr26	4103978	4103998	972.00000
-chr26	4103998	4104449	1458.00000
-chr26	4104449	4105233	1944.00000
-chr26	4105233	4105326	2430.00000
-chr26	4105326	4105398	2916.00000
-chr26	4105398	4105644	3402.00000
-chr26	4105644	4105855	3888.00000
-chr26	4105855	4105980	4374.00000
-chr26	4105980	4106825	4860.00000
-chr26	4106825	4107471	5346.00000
-chr26	4107471	4107794	5832.00000
-chr26	4107794	4108037	20000.00000
-chr26	4108037	4108498	6804.00000
-chr26	4108498	4108857	9720.00000
-chr26	4108857	4108877	20000.00000
-chr26	4108877	4109100	40000.00000
-chr26	4109100	4109120	20000.00000
-chr26	4109120	4109327	14579.99902
-chr26	4109327	4109570	20000.00000
-chr26	4109570	4109733	8262.00000
-chr26	4109733	4109826	9720.00000
-chr26	4109826	4109898	14579.99902
-chr26	4109898	4109949	19440.00000
-chr26	4109949	4110111	14579.99902
-chr26	4110111	4110204	20000.00000
-chr26	4110204	4110276	40000.00000
-chr26	4110276	4110354	60000.00000
-chr26	4110354	4110447	40000.00000
-chr26	4110447	4110480	24300.00000
-chr26	4110480	4110733	29159.99805
-chr26	4110733	4110765	40000.00000
-chr26	4110765	4110826	24300.00000
-chr26	4110826	4110859	20000.00000
-chr26	4110859	4110976	40000.00000
-chr26	4110976	4111102	20000.00000
-chr26	4111102	4111144	14579.99902
-chr26	4111144	4111325	12150.00000
-chr26	4111325	4111355	14579.99902
-chr26	4111355	4111704	12636.00000
-chr26	4111704	4111947	20000.00000
-chr26	4111947	4112349	13608.00000
-chr26	4112349	4112592	20000.00000
-chr26	4112592	4112706	14579.99902
-chr26	4112706	4112862	20000.00000
-chr26	4112862	4112949	40000.00000
-chr26	4112949	4113105	20000.00000
-chr26	4113105	4113145	13608.00000
-chr26	4113145	4113327	14579.99902
-chr26	4113327	4113524	13608.00000
-chr26	4113524	4113767	20000.00000
-chr26	4113767	4113775	13608.00000
-chr26	4113775	4114145	14579.99902
-chr26	4114145	4114146	12636.00000
-chr26	4114146	4114154	20000.00000
-chr26	4114154	4114389	40000.00000
-chr26	4114389	4114686	24300.00000
-chr26	4114686	4114732	40000.00000
-chr26	4114732	4114855	60000.00000
-chr26	4114855	4114929	40000.00000
-chr26	4114929	4114975	20000.00000
-chr26	4114975	4115116	14579.99902
-chr26	4115116	4115234	19440.00000
-chr26	4115234	4115308	14579.99902
-chr26	4115308	4115326	12150.00000
-chr26	4115326	4115398	11664.00000
-chr26	4115398	4115495	11178.00000
-chr26	4115495	4115635	20000.00000
-chr26	4115635	4115837	24300.00000
-chr26	4115837	4115891	29159.99805
-chr26	4115891	4115990	40000.00000
-chr26	4115990	4116014	60000.00000
-chr26	4116014	4116072	80000.00000
-chr26	4116072	4116134	60000.00000
-chr26	4116134	4116215	40000.00000
-chr26	4116215	4116233	60000.00000
-chr26	4116233	4116257	40000.00000
-chr26	4116257	4116366	29159.99805
-chr26	4116366	4116450	34020.00000
-chr26	4116450	4116512	29159.99805
-chr26	4116512	4116611	24300.00000
-chr26	4116611	4116745	20000.00000
-chr26	4116745	4116831	40000.00000
-chr26	4116831	4116988	20000.00000
-chr26	4116988	4117136	11664.00000
-chr26	4117136	4117471	12150.00000
-chr26	4117471	4117730	11664.00000
-chr26	4117730	4117973	20000.00000
-chr26	4117973	4117983	11178.00000
-chr26	4117983	4118141	10692.00000
-chr26	4118141	4118645	11178.00000
-chr26	4118645	4118886	10692.00000
-chr26	4118886	4118904	11178.00000
-chr26	4118904	4119267	11664.00000
-chr26	4119267	4119275	11178.00000
-chr26	4119275	4119282	10692.00000
-chr26	4119282	4119525	20000.00000
-chr26	4119525	4119734	11178.00000
-chr26	4119734	4119808	10692.00000
-chr26	4119808	4119854	10206.00000
-chr26	4119854	4120286	9720.00000
-chr26	4120286	4120399	10206.00000
-chr26	4120399	4120642	20000.00000
-chr26	4120642	4120714	9720.00000
-chr26	4120714	4120957	20000.00000
-chr26	4120957	4120959	14579.99902
-chr26	4120959	4121021	19440.00000
-chr26	4121021	4121066	14579.99902
-chr26	4121066	4121216	19440.00000
-chr26	4121216	4121338	20000.00000
-chr26	4121338	4121445	40000.00000
-chr26	4121445	4121459	60000.00000
-chr26	4121459	4121581	40000.00000
-chr26	4121581	4121688	20000.00000
-chr26	4121688	4121837	19440.00000
-chr26	4121837	4121959	14579.99902
-chr26	4121959	4122015	9720.00000
-chr26	4122015	4122258	20000.00000
-chr26	4122258	4122851	8262.00000
-chr26	4122851	4123020	7776.00000
-chr26	4123020	4123263	20000.00000
-chr26	4123263	4123386	7776.00000
-chr26	4123386	4123404	9720.00000
-chr26	4123404	4123641	14579.99902
-chr26	4123641	4123764	9720.00000
-chr26	4123764	4123783	20000.00000
-chr26	4123783	4124007	40000.00000
-chr26	4124007	4124026	20000.00000
-chr26	4124026	4124386	14579.99902
-chr26	4124386	4124396	9720.00000
-chr26	4124396	4124639	20000.00000
-chr26	4124639	4124786	9234.00000
-chr26	4124786	4125017	9720.00000
-chr26	4125017	4125164	9234.00000
-chr26	4125164	4125407	20000.00000
-chr26	4125407	4125836	9234.00000
-chr26	4125836	4126337	8748.00000
-chr26	4126337	4126363	8262.00000
-chr26	4126363	4126459	9720.00000
-chr26	4126459	4126509	14579.99902
-chr26	4126509	4126699	19440.00000
-chr26	4126699	4126742	20000.00000
-chr26	4126742	4126837	40000.00000
-chr26	4126837	4126887	60000.00000
-chr26	4126887	4126942	80000.00000
-chr26	4126942	4126985	60000.00000
-chr26	4126985	4127080	40000.00000
-chr26	4127080	4127130	20000.00000
-chr26	4127130	4127320	19440.00000
-chr26	4127320	4127363	14579.99902
-chr26	4127363	4127459	9720.00000
-chr26	4127459	4127761	6804.00000
--- a/test-data/test_MACS2.1.1_html_report.zip	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,53 +0,0 @@
-<html><head><title>Additional output created by MACS (test_MACS2.1.1)</title></head><body><h3>Additional Files:</h3><p><ul>
-<li><a href="test_MACS2.1.1_peaks.xls">test_MACS2.1.1_peaks.xls</a></li>
-</ul></p>
-<h3>Messages from MACS:</h3>
-<p><pre>INFO  @ Tue, 20 Mar 2018 14:21:02: 
-# Command line: callpeak -t /tmp/tmpHxmla3/files/000/dataset_1.dat -c /tmp/tmpHxmla3/files/000/dataset_2.dat --format=BED --name=test_MACS2.1.1 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
-# ARGUMENTS LIST:
-# name = test_MACS2.1.1
-# format = BED
-# ChIP-seq file = ['/tmp/tmpHxmla3/files/000/dataset_1.dat']
-# control file = ['/tmp/tmpHxmla3/files/000/dataset_2.dat']
-# effective genome size = 7.75e+08
-# band width = 300
-# model fold = [5, 50]
-# qvalue cutoff = 5.00e-02
-# Larger dataset will be scaled towards smaller dataset.
-# Range for calculating regional lambda is: 1000 bps and 10000 bps
-# Broad region calling is off
-# Paired-End mode is off
-# MACS will save fragment pileup signal per million reads
- 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 read tag files... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 read treatment tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1.2 read input tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 tag size is determined as 50 bps 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 tag size = 50 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  total tags in treatment: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 user defined the maximum tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  tags after filtering in treatment: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  Redundant rate of treatment: 0.00 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  total tags in control: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 user defined the maximum tags... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  tags after filtering in control: 50 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1  Redundant rate of control: 0.00 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #1 finished! 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Build Peak Model... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Skipped... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Use 243 as fragment length 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Call peaks... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Pre-compute pvalue-qvalue table... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3 In the peak calling step, the following will be performed simultaneously: 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3   Write bedGraph files for treatment pileup (after scaling if necessary)... test_MACS2.1.1_treat_pileup.bdg 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3   Write bedGraph files for control lambda (after scaling if necessary)... test_MACS2.1.1_control_lambda.bdg 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3   --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Call peaks for each chromosome... 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write output xls file... test_MACS2.1.1_peaks.xls 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write peak in narrowPeak format file... test_MACS2.1.1_peaks.narrowPeak 
-INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write summits bed file... test_MACS2.1.1_summits.bed 
-INFO  @ Tue, 20 Mar 2018 14:21:02: Done! 
-</pre></p>
-</body></html>
--- a/test-data/test_MACS2.1.1_peaks.xls	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,30 +0,0 @@
-#peaks file
-# This file is generated by MACS version 2.1.1.20160309
-# Command line: callpeak -t /tmp/tmpHxmla3/files/000/dataset_9.dat -c /tmp/tmpHxmla3/files/000/dataset_10.dat --format=BED --name=test_MACS2.1.1 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
-# ARGUMENTS LIST:
-# name = test_MACS2.1.1
-# format = BED
-# ChIP-seq file = ['/tmp/tmpHxmla3/files/000/dataset_9.dat']
-# control file = ['/tmp/tmpHxmla3/files/000/dataset_10.dat']
-# effective genome size = 7.75e+08
-# band width = 300
-# model fold = [5, 50]
-# qvalue cutoff = 5.00e-02
-# Larger dataset will be scaled towards smaller dataset.
-# Range for calculating regional lambda is: 1000 bps and 10000 bps
-# Broad region calling is off
-# Paired-End mode is off
-# MACS will save fragment pileup signal per million reads
-
-# tag size is determined as 50 bps
-# total tags in treatment: 50
-# tags after filtering in treatment: 50
-# maximum duplicate tags at the same position in treatment = 1
-# Redundant rate in treatment: 0.00
-# total tags in control: 50
-# tags after filtering in control: 50
-# maximum duplicate tags at the same position in control = 1
-# Redundant rate in control: 0.00
-# d = 243
-#chr	start	end	length	abs_summit	pileup	-log10(pvalue)	fold_enrichment	-log10(qvalue)	name
-chr26	4118914	4119282	369	4119130	9.00	9.13132	6.31632	2.51561	test_MACS2.1.1_peak_1
--- a/test-data/test_MACS2.1.1_peaks.xls.re_match	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,30 +0,0 @@
-\#peaks\ file
-\#\ This\ file\ is\ generated\ by\ MACS\ version\ 2\.1\.1\.20160309
-\#\ Command\ line\:\ callpeak\ \-t\ \/tmp\/tmpHxmla3\/files\/000\/dataset\_9\.dat\ \-c\ \/tmp\/tmpHxmla3\/files\/000\/dataset\_10\.dat\ \-\-format\=BED\ \-\-name\=test\_MACS2\.1\.1\ \-\-bw\=300\ \-\-gsize\=775000000\.0\ \-\-nomodel\ \-\-extsize\=243\ \-\-qvalue\=0\.05\ \-B\ \-\-SPMR\ \-\-mfold\ 5\ 50\ \-\-keep\-dup\ 1
-\#\ ARGUMENTS\ LIST\:
-\#\ name\ \=\ test\_MACS2\.1\.1
-\#\ format\ \=\ BED
-\#\ ChIP\-seq\ file\ \=\ \[\'.*\'\]
-\#\ control\ file\ \=\ \[\'.*\'\]
-\#\ effective\ genome\ size\ \=\ 7\.75e\+08
-\#\ band\ width\ \=\ 300
-\#\ model\ fold\ \=\ \[5\,\ 50\]
-\#\ qvalue\ cutoff\ \=\ 5\.00e\-02
-\#\ Larger\ dataset\ will\ be\ scaled\ towards\ smaller\ dataset\.
-\#\ Range\ for\ calculating\ regional\ lambda\ is\:\ 1000\ bps\ and\ 10000\ bps
-\#\ Broad\ region\ calling\ is\ off
-\#\ Paired\-End\ mode\ is\ off
-\#\ MACS\ will\ save\ fragment\ pileup\ signal\ per\ million\ reads
-
-\#\ tag\ size\ is\ determined\ as\ 50\ bps
-\#\ total\ tags\ in\ treatment\:\ 50
-\#\ tags\ after\ filtering\ in\ treatment\:\ 50
-\#\ maximum\ duplicate\ tags\ at\ the\ same\ position\ in\ treatment\ \=\ 1
-\#\ Redundant\ rate\ in\ treatment\:\ 0\.00
-\#\ total\ tags\ in\ control\:\ 50
-\#\ tags\ after\ filtering\ in\ control\:\ 50
-\#\ maximum\ duplicate\ tags\ at\ the\ same\ position\ in\ control\ \=\ 1
-\#\ Redundant\ rate\ in\ control\:\ 0\.00
-\#\ d\ \=\ 243
-\#chr\	start\	end\	length\	abs\_summit\	pileup\	\-log10\(pvalue\)\	fold\_enrichment\	\-log10\(qvalue\)\	name
-chr26\	4118914\	4119282\	369\	4119130\	9\.00\	9\.13132\	6\.31632\	2\.51561\	test\_MACS2\.1\.1\_peak\_1
--- a/test-data/test_MACS2.1.1_peaks_narrowPeak.interval	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,1 +0,0 @@
-chr26	4118913	4119282	test_MACS2.1.1_peak_1	25	.	6.31632	9.13132	2.51561	216
--- a/test-data/test_MACS2.1.1_summits.bed	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,1 +0,0 @@
-chr26	4119129	4119130	test_MACS2.1.1_peak_1	2.51561
--- a/test-data/test_MACS2.1.1_treat_pileup.bdg	Tue Mar 20 11:25:04 2018 -0400
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,100 +0,0 @@
-chr26	0	4107561	0.00000
-chr26	4107561	4107804	20000.00000
-chr26	4107804	4108165	0.00000
-chr26	4108165	4108408	20000.00000
-chr26	4108408	4108531	0.00000
-chr26	4108531	4108774	20000.00000
-chr26	4108774	4109009	0.00000
-chr26	4109009	4109252	20000.00000
-chr26	4109252	4109506	0.00000
-chr26	4109506	4109749	20000.00000
-chr26	4109749	4109782	0.00000
-chr26	4109782	4110025	20000.00000
-chr26	4110025	4110080	0.00000
-chr26	4110080	4110323	20000.00000
-chr26	4110323	4111276	0.00000
-chr26	4111276	4111519	20000.00000
-chr26	4111519	4112110	0.00000
-chr26	4112110	4112353	20000.00000
-chr26	4112353	4112762	0.00000
-chr26	4112762	4113005	20000.00000
-chr26	4113005	4113579	0.00000
-chr26	4113579	4113822	20000.00000
-chr26	4113822	4113899	0.00000
-chr26	4113899	4114142	20000.00000
-chr26	4114142	4115021	0.00000
-chr26	4115021	4115264	20000.00000
-chr26	4115264	4115555	0.00000
-chr26	4115555	4115792	20000.00000
-chr26	4115792	4115798	40000.00000
-chr26	4115798	4116035	20000.00000
-chr26	4116035	4116615	0.00000
-chr26	4116615	4116858	20000.00000
-chr26	4116858	4116931	0.00000
-chr26	4116931	4116988	20000.00000
-chr26	4116988	4117174	40000.00000
-chr26	4117174	4117231	20000.00000
-chr26	4117231	4117308	0.00000
-chr26	4117308	4117551	20000.00000
-chr26	4117551	4117615	0.00000
-chr26	4117615	4117642	20000.00000
-chr26	4117642	4117737	40000.00000
-chr26	4117737	4117858	60000.00000
-chr26	4117858	4117885	40000.00000
-chr26	4117885	4117905	20000.00000
-chr26	4117905	4117980	40000.00000
-chr26	4117980	4118037	20000.00000
-chr26	4118037	4118122	40000.00000
-chr26	4118122	4118148	60000.00000
-chr26	4118148	4118280	40000.00000
-chr26	4118280	4118365	20000.00000
-chr26	4118365	4118516	0.00000
-chr26	4118516	4118531	20000.00000
-chr26	4118531	4118561	40000.00000
-chr26	4118561	4118753	60000.00000
-chr26	4118753	4118759	80000.00000
-chr26	4118759	4118774	60000.00000
-chr26	4118774	4118804	40000.00000
-chr26	4118804	4118812	20000.00000
-chr26	4118812	4118827	40000.00000
-chr26	4118827	4118852	60000.00000
-chr26	4118852	4118898	80000.00000
-chr26	4118898	4118913	100000.00000
-chr26	4118913	4118963	120000.00000
-chr26	4118963	4118967	140000.00000
-chr26	4118967	4118996	160000.00000
-chr26	4118996	4119022	140000.00000
-chr26	4119022	4119047	160000.00000
-chr26	4119047	4119055	180000.00000
-chr26	4119055	4119070	160000.00000
-chr26	4119070	4119077	140000.00000
-chr26	4119077	4119095	160000.00000
-chr26	4119095	4119103	140000.00000
-chr26	4119103	4119118	160000.00000
-chr26	4119118	4119141	180000.00000
-chr26	4119141	4119156	160000.00000
-chr26	4119156	4119163	140000.00000
-chr26	4119163	4119168	160000.00000
-chr26	4119168	4119206	180000.00000
-chr26	4119206	4119210	160000.00000
-chr26	4119210	4119265	140000.00000
-chr26	4119265	4119290	120000.00000
-chr26	4119290	4119320	100000.00000
-chr26	4119320	4119346	80000.00000
-chr26	4119346	4119361	60000.00000
-chr26	4119361	4119406	40000.00000
-chr26	4119406	4119411	20000.00000
-chr26	4119411	4122292	0.00000
-chr26	4122292	4122535	20000.00000
-chr26	4122535	4124351	0.00000
-chr26	4124351	4124452	20000.00000
-chr26	4124452	4124594	40000.00000
-chr26	4124594	4124695	20000.00000
-chr26	4124695	4125766	0.00000
-chr26	4125766	4125809	20000.00000
-chr26	4125809	4126009	40000.00000
-chr26	4126009	4126052	20000.00000
-chr26	4126052	4126452	0.00000
-chr26	4126452	4126695	20000.00000
-chr26	4126695	4127518	0.00000
-chr26	4127518	4127761	20000.00000
Binary file test-data/test_MACS2.1.1_treat_pileup.bw has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_bw_html_report.zip	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,55 @@
+<html><head><title>Additional output created by MACS (test_MACS2.1.2)</title></head><body><h3>Additional Files:</h3><p><ul>
+<li><a href="test_MACS2.1.2_peaks.xls">test_MACS2.1.2_peaks.xls</a></li>
+<li><a href="test_MACS2.1.2_treat_pileup.bdg.clipped">test_MACS2.1.2_treat_pileup.bdg.clipped</a></li>
+<li><a href="test_MACS2.1.2_treat_pileup.bdg.clipped.sorted">test_MACS2.1.2_treat_pileup.bdg.clipped.sorted</a></li>
+</ul></p>
+<h3>Messages from MACS:</h3>
+<p><pre>INFO  @ Tue, 20 Mar 2018 14:21:47: 
+# Command line: callpeak -t /tmp/tmpHxmla3/files/000/dataset_9.dat -c /tmp/tmpHxmla3/files/000/dataset_10.dat --format=BED --name=test_MACS2.1.2 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
+# ARGUMENTS LIST:
+# name = test_MACS2.1.2
+# format = BED
+# ChIP-seq file = ['/tmp/tmpHxmla3/files/000/dataset_9.dat']
+# control file = ['/tmp/tmpHxmla3/files/000/dataset_10.dat']
+# effective genome size = 7.75e+08
+# band width = 300
+# model fold = [5, 50]
+# qvalue cutoff = 5.00e-02
+# Larger dataset will be scaled towards smaller dataset.
+# Range for calculating regional lambda is: 1000 bps and 10000 bps
+# Broad region calling is off
+# Paired-End mode is off
+# MACS will save fragment pileup signal per million reads
+ 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 read tag files... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 read treatment tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1.2 read input tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 tag size is determined as 50 bps 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 tag size = 50 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  total tags in treatment: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 user defined the maximum tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  tags after filtering in treatment: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  Redundant rate of treatment: 0.00 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  total tags in control: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 user defined the maximum tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  tags after filtering in control: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1  Redundant rate of control: 0.00 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #1 finished! 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Build Peak Model... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Skipped... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #2 Use 243 as fragment length 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Call peaks... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Pre-compute pvalue-qvalue table... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3 In the peak calling step, the following will be performed simultaneously: 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3   Write bedGraph files for treatment pileup (after scaling if necessary)... test_MACS2.1.2_treat_pileup.bdg 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3   Write bedGraph files for control lambda (after scaling if necessary)... test_MACS2.1.2_control_lambda.bdg 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3   --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #3 Call peaks for each chromosome... 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write output xls file... test_MACS2.1.2_peaks.xls 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write peak in narrowPeak format file... test_MACS2.1.2_peaks.narrowPeak 
+INFO  @ Tue, 20 Mar 2018 14:21:47: #4 Write summits bed file... test_MACS2.1.2_summits.bed 
+INFO  @ Tue, 20 Mar 2018 14:21:47: Done! 
+</pre></p>
+</body></html>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_control_lambda.bdg	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,153 @@
+chr26	0	4102915	0.31355
+chr26	4102915	4103978	486.00000
+chr26	4103978	4103998	972.00000
+chr26	4103998	4104449	1458.00000
+chr26	4104449	4105233	1944.00000
+chr26	4105233	4105326	2430.00000
+chr26	4105326	4105398	2916.00000
+chr26	4105398	4105644	3402.00000
+chr26	4105644	4105855	3888.00000
+chr26	4105855	4105980	4374.00000
+chr26	4105980	4106825	4860.00000
+chr26	4106825	4107471	5346.00000
+chr26	4107471	4107794	5832.00000
+chr26	4107794	4108037	20000.00000
+chr26	4108037	4108498	6804.00000
+chr26	4108498	4108857	9720.00000
+chr26	4108857	4108877	20000.00000
+chr26	4108877	4109100	40000.00000
+chr26	4109100	4109120	20000.00000
+chr26	4109120	4109327	14579.99902
+chr26	4109327	4109570	20000.00000
+chr26	4109570	4109733	8262.00000
+chr26	4109733	4109826	9720.00000
+chr26	4109826	4109898	14579.99902
+chr26	4109898	4109949	19440.00000
+chr26	4109949	4110111	14579.99902
+chr26	4110111	4110204	20000.00000
+chr26	4110204	4110276	40000.00000
+chr26	4110276	4110354	60000.00000
+chr26	4110354	4110447	40000.00000
+chr26	4110447	4110480	24300.00000
+chr26	4110480	4110733	29159.99805
+chr26	4110733	4110765	40000.00000
+chr26	4110765	4110826	24300.00000
+chr26	4110826	4110859	20000.00000
+chr26	4110859	4110976	40000.00000
+chr26	4110976	4111102	20000.00000
+chr26	4111102	4111144	14579.99902
+chr26	4111144	4111325	12150.00000
+chr26	4111325	4111355	14579.99902
+chr26	4111355	4111704	12636.00000
+chr26	4111704	4111947	20000.00000
+chr26	4111947	4112349	13608.00000
+chr26	4112349	4112592	20000.00000
+chr26	4112592	4112706	14579.99902
+chr26	4112706	4112862	20000.00000
+chr26	4112862	4112949	40000.00000
+chr26	4112949	4113105	20000.00000
+chr26	4113105	4113145	13608.00000
+chr26	4113145	4113327	14579.99902
+chr26	4113327	4113524	13608.00000
+chr26	4113524	4113767	20000.00000
+chr26	4113767	4113775	13608.00000
+chr26	4113775	4114145	14579.99902
+chr26	4114145	4114146	12636.00000
+chr26	4114146	4114154	20000.00000
+chr26	4114154	4114389	40000.00000
+chr26	4114389	4114686	24300.00000
+chr26	4114686	4114732	40000.00000
+chr26	4114732	4114855	60000.00000
+chr26	4114855	4114929	40000.00000
+chr26	4114929	4114975	20000.00000
+chr26	4114975	4115116	14579.99902
+chr26	4115116	4115234	19440.00000
+chr26	4115234	4115308	14579.99902
+chr26	4115308	4115326	12150.00000
+chr26	4115326	4115398	11664.00000
+chr26	4115398	4115495	11178.00000
+chr26	4115495	4115635	20000.00000
+chr26	4115635	4115837	24300.00000
+chr26	4115837	4115891	29159.99805
+chr26	4115891	4115990	40000.00000
+chr26	4115990	4116014	60000.00000
+chr26	4116014	4116072	80000.00000
+chr26	4116072	4116134	60000.00000
+chr26	4116134	4116215	40000.00000
+chr26	4116215	4116233	60000.00000
+chr26	4116233	4116257	40000.00000
+chr26	4116257	4116366	29159.99805
+chr26	4116366	4116450	34020.00000
+chr26	4116450	4116512	29159.99805
+chr26	4116512	4116611	24300.00000
+chr26	4116611	4116745	20000.00000
+chr26	4116745	4116831	40000.00000
+chr26	4116831	4116988	20000.00000
+chr26	4116988	4117136	11664.00000
+chr26	4117136	4117471	12150.00000
+chr26	4117471	4117730	11664.00000
+chr26	4117730	4117973	20000.00000
+chr26	4117973	4117983	11178.00000
+chr26	4117983	4118141	10692.00000
+chr26	4118141	4118645	11178.00000
+chr26	4118645	4118886	10692.00000
+chr26	4118886	4118904	11178.00000
+chr26	4118904	4119267	11664.00000
+chr26	4119267	4119275	11178.00000
+chr26	4119275	4119282	10692.00000
+chr26	4119282	4119525	20000.00000
+chr26	4119525	4119734	11178.00000
+chr26	4119734	4119808	10692.00000
+chr26	4119808	4119854	10206.00000
+chr26	4119854	4120286	9720.00000
+chr26	4120286	4120399	10206.00000
+chr26	4120399	4120642	20000.00000
+chr26	4120642	4120714	9720.00000
+chr26	4120714	4120957	20000.00000
+chr26	4120957	4120959	14579.99902
+chr26	4120959	4121021	19440.00000
+chr26	4121021	4121066	14579.99902
+chr26	4121066	4121216	19440.00000
+chr26	4121216	4121338	20000.00000
+chr26	4121338	4121445	40000.00000
+chr26	4121445	4121459	60000.00000
+chr26	4121459	4121581	40000.00000
+chr26	4121581	4121688	20000.00000
+chr26	4121688	4121837	19440.00000
+chr26	4121837	4121959	14579.99902
+chr26	4121959	4122015	9720.00000
+chr26	4122015	4122258	20000.00000
+chr26	4122258	4122851	8262.00000
+chr26	4122851	4123020	7776.00000
+chr26	4123020	4123263	20000.00000
+chr26	4123263	4123386	7776.00000
+chr26	4123386	4123404	9720.00000
+chr26	4123404	4123641	14579.99902
+chr26	4123641	4123764	9720.00000
+chr26	4123764	4123783	20000.00000
+chr26	4123783	4124007	40000.00000
+chr26	4124007	4124026	20000.00000
+chr26	4124026	4124386	14579.99902
+chr26	4124386	4124396	9720.00000
+chr26	4124396	4124639	20000.00000
+chr26	4124639	4124786	9234.00000
+chr26	4124786	4125017	9720.00000
+chr26	4125017	4125164	9234.00000
+chr26	4125164	4125407	20000.00000
+chr26	4125407	4125836	9234.00000
+chr26	4125836	4126337	8748.00000
+chr26	4126337	4126363	8262.00000
+chr26	4126363	4126459	9720.00000
+chr26	4126459	4126509	14579.99902
+chr26	4126509	4126699	19440.00000
+chr26	4126699	4126742	20000.00000
+chr26	4126742	4126837	40000.00000
+chr26	4126837	4126887	60000.00000
+chr26	4126887	4126942	80000.00000
+chr26	4126942	4126985	60000.00000
+chr26	4126985	4127080	40000.00000
+chr26	4127080	4127130	20000.00000
+chr26	4127130	4127320	19440.00000
+chr26	4127320	4127363	14579.99902
+chr26	4127363	4127459	9720.00000
+chr26	4127459	4127761	6804.00000
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_html_report.zip	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,53 @@
+<html><head><title>Additional output created by MACS (test_MACS2.1.2)</title></head><body><h3>Additional Files:</h3><p><ul>
+<li><a href="test_MACS2.1.2_peaks.xls">test_MACS2.1.2_peaks.xls</a></li>
+</ul></p>
+<h3>Messages from MACS:</h3>
+<p><pre>INFO  @ Tue, 20 Mar 2018 14:21:02: 
+# Command line: callpeak -t /tmp/tmpHxmla3/files/000/dataset_1.dat -c /tmp/tmpHxmla3/files/000/dataset_2.dat --format=BED --name=test_MACS2.1.2 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
+# ARGUMENTS LIST:
+# name = test_MACS2.1.2
+# format = BED
+# ChIP-seq file = ['/tmp/tmpHxmla3/files/000/dataset_1.dat']
+# control file = ['/tmp/tmpHxmla3/files/000/dataset_2.dat']
+# effective genome size = 7.75e+08
+# band width = 300
+# model fold = [5, 50]
+# qvalue cutoff = 5.00e-02
+# Larger dataset will be scaled towards smaller dataset.
+# Range for calculating regional lambda is: 1000 bps and 10000 bps
+# Broad region calling is off
+# Paired-End mode is off
+# MACS will save fragment pileup signal per million reads
+ 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 read tag files... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 read treatment tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1.2 read input tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 tag size is determined as 50 bps 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 tag size = 50 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  total tags in treatment: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 user defined the maximum tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  tags after filtering in treatment: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  Redundant rate of treatment: 0.00 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  total tags in control: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 user defined the maximum tags... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  tags after filtering in control: 50 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1  Redundant rate of control: 0.00 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #1 finished! 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Build Peak Model... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Skipped... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #2 Use 243 as fragment length 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Call peaks... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Pre-compute pvalue-qvalue table... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3 In the peak calling step, the following will be performed simultaneously: 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3   Write bedGraph files for treatment pileup (after scaling if necessary)... test_MACS2.1.2_treat_pileup.bdg 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3   Write bedGraph files for control lambda (after scaling if necessary)... test_MACS2.1.2_control_lambda.bdg 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3   --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #3 Call peaks for each chromosome... 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write output xls file... test_MACS2.1.2_peaks.xls 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write peak in narrowPeak format file... test_MACS2.1.2_peaks.narrowPeak 
+INFO  @ Tue, 20 Mar 2018 14:21:02: #4 Write summits bed file... test_MACS2.1.2_summits.bed 
+INFO  @ Tue, 20 Mar 2018 14:21:02: Done! 
+</pre></p>
+</body></html>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_peaks.xls	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,32 @@
+#peaks file
+# This file is generated by MACS version 2.1.2
+# Command line: callpeak -t /tmp/tmpKNHGvp/files/000/dataset_2.dat -c /tmp/tmpKNHGvp/files/000/dataset_1.dat --format=BED --name=test_MACS2.1.2 --bw=300 --gsize=775000000.0 --nomodel --extsize=243 --qvalue=0.05 -B --SPMR --mfold 5 50 --keep-dup 1
+# ARGUMENTS LIST:
+# name = test_MACS2.1.2
+# format = BED
+# ChIP-seq file = ['/tmp/tmpKNHGvp/files/000/dataset_2.dat']
+# control file = ['/tmp/tmpKNHGvp/files/000/dataset_1.dat']
+# effective genome size = 7.75e+08
+# band width = 300
+# model fold = [5, 50]
+# qvalue cutoff = 5.00e-02
+# The maximum gap between significant sites is assigned as the read length/tag size.
+# The minimum length of peaks is assigned as the predicted fragment length "d".
+# Larger dataset will be scaled towards smaller dataset.
+# Range for calculating regional lambda is: 1000 bps and 10000 bps
+# Broad region calling is off
+# Paired-End mode is off
+# MACS will save fragment pileup signal per million reads
+
+# tag size is determined as 50 bps
+# total tags in treatment: 50
+# tags after filtering in treatment: 50
+# maximum duplicate tags at the same position in treatment = 1
+# Redundant rate in treatment: 0.00
+# total tags in control: 50
+# tags after filtering in control: 50
+# maximum duplicate tags at the same position in control = 1
+# Redundant rate in control: 0.00
+# d = 243
+#chr	start	end	length	abs_summit	pileup	-log10(pvalue)	fold_enrichment	-log10(qvalue)	name
+chr26	4118914	4119282	369	4119130	9.00	9.13132	6.31632	2.51561	test_MACS2.1.2_peak_1
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_peaks.xls.re_match	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,32 @@
+\#peaks\ file
+\#\ This\ file\ is\ generated\ by\ MACS\ version\ 2\.1\.2
+\#\ Command\ line\:\ callpeak\ \-t\ .*\.dat \-c\ .*\.dat\ \-\-format\=BED\ \-\-name\=test\_MACS2\.1\.2\ \-\-bw\=300\ \-\-gsize\=775000000\.0\ \-\-nomodel\ \-\-extsize\=243\ \-\-qvalue\=0\.05\ \-B\ \-\-SPMR\ \-\-mfold\ 5\ 50\ \-\-keep\-dup\ 1
+\#\ ARGUMENTS\ LIST\:
+\#\ name\ \=\ test\_MACS2\.1\.2
+\#\ format\ \=\ BED
+\#\ ChIP\-seq\ file\ \=\ \[\'.*\'\]
+\#\ control\ file\ \=\ \[\'.*\'\]
+\#\ effective\ genome\ size\ \=\ 7\.75e\+08
+\#\ band\ width\ \=\ 300
+\#\ model\ fold\ \=\ \[5\,\ 50\]
+\#\ qvalue\ cutoff\ \=\ 5\.00e\-02
+\#\ The\ maximum\ gap\ between\ significant\ sites\ is\ assigned\ as\ the\ read\ length\/tag\ size\.
+\#\ The\ minimum\ length\ of\ peaks\ is\ assigned\ as\ the\ predicted\ fragment\ length\ \"d\"\.
+\#\ Larger\ dataset\ will\ be\ scaled\ towards\ smaller\ dataset\.
+\#\ Range\ for\ calculating\ regional\ lambda\ is\:\ 1000\ bps\ and\ 10000\ bps
+\#\ Broad\ region\ calling\ is\ off
+\#\ Paired\-End\ mode\ is\ off
+\#\ MACS\ will\ save\ fragment\ pileup\ signal\ per\ million\ reads
+
+\#\ tag\ size\ is\ determined\ as\ 50\ bps
+\#\ total\ tags\ in\ treatment\:\ 50
+\#\ tags\ after\ filtering\ in\ treatment\:\ 50
+\#\ maximum\ duplicate\ tags\ at\ the\ same\ position\ in\ treatment\ \=\ 1
+\#\ Redundant\ rate\ in\ treatment\:\ 0\.00
+\#\ total\ tags\ in\ control\:\ 50
+\#\ tags\ after\ filtering\ in\ control\:\ 50
+\#\ maximum\ duplicate\ tags\ at\ the\ same\ position\ in\ control\ \=\ 1
+\#\ Redundant\ rate\ in\ control\:\ 0\.00
+\#\ d\ \=\ 243
+\#chr\	start\	end\	length\	abs\_summit\	pileup\	\-log10\(pvalue\)\	fold\_enrichment\	\-log10\(qvalue\)\	name
+chr26\	4118914\	4119282\	369\	4119130\	9\.00\	9\.13132\	6\.31632\	2\.51561\	test\_MACS2\.1\.2\_peak\_2
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_peaks_narrowPeak.interval	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,1 @@
+chr26	4118913	4119282	test_MACS2.1.2_peak_1	25	.	6.31632	9.13132	2.51561	216
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_summits.bed	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,1 @@
+chr26	4119129	4119130	test_MACS2.1.2_peak_1	2.51561
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test_MACS2.1.2_treat_pileup.bdg	Wed Dec 12 08:26:16 2018 -0500
@@ -0,0 +1,100 @@
+chr26	0	4107561	0.00000
+chr26	4107561	4107804	20000.00000
+chr26	4107804	4108165	0.00000
+chr26	4108165	4108408	20000.00000
+chr26	4108408	4108531	0.00000
+chr26	4108531	4108774	20000.00000
+chr26	4108774	4109009	0.00000
+chr26	4109009	4109252	20000.00000
+chr26	4109252	4109506	0.00000
+chr26	4109506	4109749	20000.00000
+chr26	4109749	4109782	0.00000
+chr26	4109782	4110025	20000.00000
+chr26	4110025	4110080	0.00000
+chr26	4110080	4110323	20000.00000
+chr26	4110323	4111276	0.00000
+chr26	4111276	4111519	20000.00000
+chr26	4111519	4112110	0.00000
+chr26	4112110	4112353	20000.00000
+chr26	4112353	4112762	0.00000
+chr26	4112762	4113005	20000.00000
+chr26	4113005	4113579	0.00000
+chr26	4113579	4113822	20000.00000
+chr26	4113822	4113899	0.00000
+chr26	4113899	4114142	20000.00000
+chr26	4114142	4115021	0.00000
+chr26	4115021	4115264	20000.00000
+chr26	4115264	4115555	0.00000
+chr26	4115555	4115792	20000.00000
+chr26	4115792	4115798	40000.00000
+chr26	4115798	4116035	20000.00000
+chr26	4116035	4116615	0.00000
+chr26	4116615	4116858	20000.00000
+chr26	4116858	4116931	0.00000
+chr26	4116931	4116988	20000.00000
+chr26	4116988	4117174	40000.00000
+chr26	4117174	4117231	20000.00000
+chr26	4117231	4117308	0.00000
+chr26	4117308	4117551	20000.00000
+chr26	4117551	4117615	0.00000
+chr26	4117615	4117642	20000.00000
+chr26	4117642	4117737	40000.00000
+chr26	4117737	4117858	60000.00000
+chr26	4117858	4117885	40000.00000
+chr26	4117885	4117905	20000.00000
+chr26	4117905	4117980	40000.00000
+chr26	4117980	4118037	20000.00000
+chr26	4118037	4118122	40000.00000
+chr26	4118122	4118148	60000.00000
+chr26	4118148	4118280	40000.00000
+chr26	4118280	4118365	20000.00000
+chr26	4118365	4118516	0.00000
+chr26	4118516	4118531	20000.00000
+chr26	4118531	4118561	40000.00000
+chr26	4118561	4118753	60000.00000
+chr26	4118753	4118759	80000.00000
+chr26	4118759	4118774	60000.00000
+chr26	4118774	4118804	40000.00000
+chr26	4118804	4118812	20000.00000
+chr26	4118812	4118827	40000.00000
+chr26	4118827	4118852	60000.00000
+chr26	4118852	4118898	80000.00000
+chr26	4118898	4118913	100000.00000
+chr26	4118913	4118963	120000.00000
+chr26	4118963	4118967	140000.00000
+chr26	4118967	4118996	160000.00000
+chr26	4118996	4119022	140000.00000
+chr26	4119022	4119047	160000.00000
+chr26	4119047	4119055	180000.00000
+chr26	4119055	4119070	160000.00000
+chr26	4119070	4119077	140000.00000
+chr26	4119077	4119095	160000.00000
+chr26	4119095	4119103	140000.00000
+chr26	4119103	4119118	160000.00000
+chr26	4119118	4119141	180000.00000
+chr26	4119141	4119156	160000.00000
+chr26	4119156	4119163	140000.00000
+chr26	4119163	4119168	160000.00000
+chr26	4119168	4119206	180000.00000
+chr26	4119206	4119210	160000.00000
+chr26	4119210	4119265	140000.00000
+chr26	4119265	4119290	120000.00000
+chr26	4119290	4119320	100000.00000
+chr26	4119320	4119346	80000.00000
+chr26	4119346	4119361	60000.00000
+chr26	4119361	4119406	40000.00000
+chr26	4119406	4119411	20000.00000
+chr26	4119411	4122292	0.00000
+chr26	4122292	4122535	20000.00000
+chr26	4122535	4124351	0.00000
+chr26	4124351	4124452	20000.00000
+chr26	4124452	4124594	40000.00000
+chr26	4124594	4124695	20000.00000
+chr26	4124695	4125766	0.00000
+chr26	4125766	4125809	20000.00000
+chr26	4125809	4126009	40000.00000
+chr26	4126009	4126052	20000.00000
+chr26	4126052	4126452	0.00000
+chr26	4126452	4126695	20000.00000
+chr26	4126695	4127518	0.00000
+chr26	4127518	4127761	20000.00000
Binary file test-data/test_MACS2.1.2_treat_pileup.bw has changed