Additional Files:
CMD Executed:
macs14 -t /pasteur/projets/policy01/galaxy-prod/galaxy-dist/database/files/001/122/dataset_1122728.dat -c /pasteur/projets/policy01/galaxy-prod/galaxy-dist/database/files/001/122/dataset_1122715.dat --format='BAM' --name='MACS14_in_Galaxy' --gsize='4639675' --tsize='25' --bw='400' --pvalue='1e-05' --mfold='10,30' --slocal='1000' --llocal='10000' --diag --fe-min='0' --fe-max='32' --fe-step='20'
Messages from MACS:
INFO @ Thu, 06 Feb 2020 13:55:23:
# ARGUMENTS LIST:
# name = MACS14_in_Galaxy
# format = BAM
# ChIP-seq file = /pasteur/projets/policy01/galaxy-prod/galaxy-dist/database/files/001/122/dataset_1122728.dat
# control file = /pasteur/projets/policy01/galaxy-prod/galaxy-dist/database/files/001/122/dataset_1122715.dat
# effective genome size = 4.64e+06
# band width = 400
# model fold = 10,30
# pvalue cutoff = 1.00e-05
# Large dataset will be scaled towards smaller dataset.
# Range for calculating regional lambda is: 1000 bps and 10000 bps
INFO @ Thu, 06 Feb 2020 13:55:23: #1 read tag files...
INFO @ Thu, 06 Feb 2020 13:55:23: #1 read treatment tags...
INFO @ Thu, 06 Feb 2020 13:55:23: tag size: 36
INFO @ Thu, 06 Feb 2020 13:55:29: 1000000
INFO @ Thu, 06 Feb 2020 13:55:35: 2000000
INFO @ Thu, 06 Feb 2020 13:55:40: 3000000
INFO @ Thu, 06 Feb 2020 13:55:43: #1.2 read input tags...
INFO @ Thu, 06 Feb 2020 13:55:49: 1000000
INFO @ Thu, 06 Feb 2020 13:55:55: 2000000
INFO @ Thu, 06 Feb 2020 13:56:01: 3000000
INFO @ Thu, 06 Feb 2020 13:56:07: 4000000
INFO @ Thu, 06 Feb 2020 13:56:13: 5000000
INFO @ Thu, 06 Feb 2020 13:56:18: 6000000
INFO @ Thu, 06 Feb 2020 13:56:23: #1 tag size is determined as 25 bps
INFO @ Thu, 06 Feb 2020 13:56:23: #1 tag size = 25
INFO @ Thu, 06 Feb 2020 13:56:23: #1 total tags in treatment: 2383266
INFO @ Thu, 06 Feb 2020 13:56:23: #1 user defined the maximum tags...
INFO @ Thu, 06 Feb 2020 13:56:23: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s)
INFO @ Thu, 06 Feb 2020 13:56:23: #1 tags after filtering in treatment: 1166947
INFO @ Thu, 06 Feb 2020 13:56:23: #1 Redundant rate of treatment: 0.51
INFO @ Thu, 06 Feb 2020 13:56:23: #1 total tags in control: 6649320
INFO @ Thu, 06 Feb 2020 13:56:23: #1 user defined the maximum tags...
INFO @ Thu, 06 Feb 2020 13:56:23: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s)
INFO @ Thu, 06 Feb 2020 13:56:24: #1 tags after filtering in control: 4400130
INFO @ Thu, 06 Feb 2020 13:56:24: #1 Redundant rate of control: 0.34
INFO @ Thu, 06 Feb 2020 13:56:24: #1 finished!
INFO @ Thu, 06 Feb 2020 13:56:24: #2 Build Peak Model...
INFO @ Thu, 06 Feb 2020 13:56:25: #2 number of paired peaks: 0
WARNING @ Thu, 06 Feb 2020 13:56:25: Too few paired peaks (0) so I can not build the model! Broader your MFOLD range parameter may erase this error. If it still can't build the model, please use --nomodel and --shiftsize 100 instead.
WARNING @ Thu, 06 Feb 2020 13:56:25: Process for pairing-model is terminated!
WARNING @ Thu, 06 Feb 2020 13:56:25: #2 Skipped...
WARNING @ Thu, 06 Feb 2020 13:56:25: #2 Use 100 as shiftsize, 200 as fragment length
INFO @ Thu, 06 Feb 2020 13:56:25: #3 Call peaks...
WARNING @ Thu, 06 Feb 2020 13:56:25: Treatment tags and Control tags are uneven! FDR may be wrong!
INFO @ Thu, 06 Feb 2020 13:56:25: #3 shift treatment data
INFO @ Thu, 06 Feb 2020 13:56:25: #3 merge +/- strand of treatment data
INFO @ Thu, 06 Feb 2020 13:56:25: #3 call peak candidates
INFO @ Thu, 06 Feb 2020 13:56:29: #3 shift control data
INFO @ Thu, 06 Feb 2020 13:56:29: #3 merge +/- strand of control data
INFO @ Thu, 06 Feb 2020 13:56:30: #3 call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:56:44: #3 use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:56:47: #3 Finally, 112 peaks are called!
INFO @ Thu, 06 Feb 2020 13:56:47: #3 find negative peaks by swapping treat and control
INFO @ Thu, 06 Feb 2020 13:56:51: #3 Finally, 0 peaks are called!
INFO @ Thu, 06 Feb 2020 13:56:51: #3 diag: sample 90%
INFO @ Thu, 06 Feb 2020 13:56:59: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:57:03: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:57:15: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:57:17: #3 Finally, 166 peaks are called!
INFO @ Thu, 06 Feb 2020 13:57:17: #3 diag: sample 80%
INFO @ Thu, 06 Feb 2020 13:57:25: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:57:28: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:57:37: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:57:40: #3 Finally, 157 peaks are called!
INFO @ Thu, 06 Feb 2020 13:57:40: #3 diag: sample 70%
INFO @ Thu, 06 Feb 2020 13:57:47: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:57:49: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:57:56: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:57:58: #3 Finally, 154 peaks are called!
INFO @ Thu, 06 Feb 2020 13:57:58: #3 diag: sample 60%
INFO @ Thu, 06 Feb 2020 13:58:04: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:58:06: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:58:12: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:58:13: #3 Finally, 147 peaks are called!
INFO @ Thu, 06 Feb 2020 13:58:13: #3 diag: sample 50%
INFO @ Thu, 06 Feb 2020 13:58:18: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:58:20: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:58:24: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:58:25: #3 Finally, 141 peaks are called!
INFO @ Thu, 06 Feb 2020 13:58:25: #3 diag: sample 40%
INFO @ Thu, 06 Feb 2020 13:58:29: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:58:30: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:58:33: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:58:34: #3 Finally, 135 peaks are called!
INFO @ Thu, 06 Feb 2020 13:58:34: #3 diag: sample 30%
INFO @ Thu, 06 Feb 2020 13:58:37: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:58:38: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:58:39: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:58:40: #3 Finally, 123 peaks are called!
INFO @ Thu, 06 Feb 2020 13:58:40: #3 diag: sample 20%
INFO @ Thu, 06 Feb 2020 13:58:42: #3 diag: call peak candidates
INFO @ Thu, 06 Feb 2020 13:58:42: #3 diag: call negative peak candidates
INFO @ Thu, 06 Feb 2020 13:58:43: #3 diag: use control data to filter peak candidates...
INFO @ Thu, 06 Feb 2020 13:58:43: #3 Finally, 128 peaks are called!
INFO @ Thu, 06 Feb 2020 13:58:43: #4 Write output xls file... MACS14_in_Galaxy_peaks.xls
INFO @ Thu, 06 Feb 2020 13:58:43: #4 Write peak bed file... MACS14_in_Galaxy_peaks.bed
INFO @ Thu, 06 Feb 2020 13:58:43: #4 Write summits bed file... MACS14_in_Galaxy_summits.bed
INFO @ Thu, 06 Feb 2020 13:58:43: #4 Write output xls file for negative peaks... MACS14_in_Galaxy_negative_peaks.xls
INFO @ Thu, 06 Feb 2020 13:58:43: #4 Write diagnosis result ... MACS14_in_Galaxy_diag.xls
INFO @ Thu, 06 Feb 2020 13:58:43: #5 Done! Check the output files!