Hello, galaxy users and developers, I wrote a perl script conducting saturation analysis which called macs several times.The tool's output was right but ,errorAn error occurred with this dataset:INFO @ Sat, 18 Jan 2014 21:46:33: # ARGUMENTS LIST: # name = /usr/local/galaxy-dist/database/job_working_directory/000/0.output # format = BED # ChIP-seq file = 0.test.bed # control file = 0.control.bed # effective genome size = 2.70e+09 # tag size = 2 Tool execution generated the following error message: INFO @ Sat, 18 Jan 2014 21:46:33: # ARGUMENTS LIST: # name = /usr/local/galaxy-dist/database/job_working_directory/000/0.output # format = BED # ChIP-seq file = 0.test.bed # control file = 0.control.bed # effective genome size = 2.70e+09 # tag size = 28 # band width = 300 # model fold = 16 # pvalue cutoff = 1.00e-05 # Ranges for calculating regional lambda are : peak_region,1000,5000,10000 INFO @ Sat, 18 Jan 2014 21:46:33: #1 read tag files... INFO @ Sat, 18 Jan 2014 21:46:33: #1 read treatment tags... INFO @ Sat, 18 Jan 2014 21:46:37: #1.2 read input tags... INFO @ Sat, 18 Jan 2014 21:46:39: #1 Background Redundant rate: 0.04 INFO @ Sat, 18 Jan 2014 21:46:39: #1 finished! INFO @ Sat, 18 Jan 2014 21:46:39: #2 Build Peak Model... INFO @ Sat, 18 Jan 2014 21:46:47: #2 number of paired peaks: 11875 INFO @ Sat, 18 Jan 2014 21:46:47: #2 finished! INFO @ Sat, 18 Jan 2014 21:46:47: #2.2 Generate R script for model : /usr/local/galaxy-dist/database/job_working_directory/000/0.output_model.r INFO @ Sat, 18 Jan 2014 21:46:47: #3 Call peaks... INFO @ Sat, 18 Jan 2014 21:46:47: #3 shift treatment data INFO @ Sat, 18 Jan 2014 21:46:47: #3 merge +/- strand of treatment data INFO @ Sat, 18 Jan 2014 21:46:47: #3 call peak candidates INFO @ Sat, 18 Jan 2014 21:46:48: #3 shift control data INFO @ Sat, 18 Jan 2014 21:46:48: #3 merge +/- strand of control data INFO @ Sat, 18 Jan 2014 21:46:48: #3 call negative peak candidates INFO @ Sat, 18 Jan 2014 21:46:48: #3 use control data to filter peak candidates... INFO @ Sat, 18 Jan 2014 21:46:48: #3 Finally, 2003 peaks are called! INFO @ Sat, 18 Jan 2014 21:46:48: #3 find negative peaks by reversing treat and control INFO @ Sat, 18 Jan 2014 21:46:49: #3 Finally, 120 peaks are called! INFO @ Sat, 18 Jan 2014 21:46:49: #4 Write output xls file... /usr/local/galaxy-dist/database/job_working_directory/000/0.output_peaks.xls INFO @ Sat, 18 Jan 2014 21:46:49: #4 Write output bed file... /usr/local/galaxy-dist/database/job_working_directory/000/0.output_peaks.bed INFO @ Sat, 18 Jan 2014 21:46:49: #4 Write output xls file for negative peaks... /usr/local/galaxy-dist/database/job_working_directory/000/0.output_negative_peaks.xls INFO @ Sat, 18 Jan 2014 21:46:49: #5 Done! Check the output files!
Just ordinary macs output, we see it in galaxy-macs report. The only problem is that they're in stderr. With this error I cannot do other works with the output of this tool.How can I deal with it?Thanks! -Han --------------------------------