INFO @ 10 Jun 2025 23:28:13: [32 MB] # Command line: callpeak -t 2_celltype_fragments/pseudorep1/1_sorted.tsv -f BED -n 1_pseudoreplicate1 -g hs --outdir 3_peaks -p 0.01 --shift -75 --extsize 150 --nomodel -B --SPMR --keep-dup all --call-summits # ARGUMENTS LIST: # name = 1_pseudoreplicate1 # format = BED # ChIP-seq file = ['2_celltype_fragments/pseudorep1/1_sorted.tsv'] # control file = None # effective genome size = 2.91e+09 # band width = 300 # model fold = [5, 50] # pvalue cutoff = 1.00e-02 # qvalue will not be calculated and reported as -1 in the final output. # 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: 10000 bps # Broad region calling is off # Paired-End mode is off # Searching for subpeak summits is on # MACS will save fragment pileup signal per million reads INFO @ 10 Jun 2025 23:28:13: [32 MB] #1 read tag files... INFO @ 10 Jun 2025 23:28:13: [32 MB] #1 read treatment tags... INFO @ 10 Jun 2025 23:28:15: [99 MB] 1000000 reads parsed INFO @ 10 Jun 2025 23:28:17: [103 MB] 2000000 reads parsed INFO @ 10 Jun 2025 23:28:18: [105 MB] 3000000 reads parsed INFO @ 10 Jun 2025 23:28:20: [109 MB] 4000000 reads parsed INFO @ 10 Jun 2025 23:28:21: [111 MB] 5000000 reads parsed INFO @ 10 Jun 2025 23:28:23: [112 MB] 6000000 reads parsed INFO @ 10 Jun 2025 23:28:25: [115 MB] 7000000 reads parsed INFO @ 10 Jun 2025 23:28:26: [121 MB] 8000000 reads parsed INFO @ 10 Jun 2025 23:28:28: [124 MB] 9000000 reads parsed INFO @ 10 Jun 2025 23:28:30: [129 MB] 10000000 reads parsed INFO @ 10 Jun 2025 23:28:31: [138 MB] 11000000 reads parsed INFO @ 10 Jun 2025 23:28:32: [141 MB] 12000000 reads parsed INFO @ 10 Jun 2025 23:28:34: [144 MB] 13000000 reads parsed INFO @ 10 Jun 2025 23:28:35: [147 MB] 14000000 reads parsed INFO @ 10 Jun 2025 23:28:37: [150 MB] 15000000 reads parsed INFO @ 10 Jun 2025 23:28:38: [165 MB] 16000000 reads parsed INFO @ 10 Jun 2025 23:28:40: [165 MB] 17000000 reads parsed INFO @ 10 Jun 2025 23:28:41: [167 MB] 18000000 reads parsed INFO @ 10 Jun 2025 23:28:42: [170 MB] 19000000 reads parsed INFO @ 10 Jun 2025 23:28:44: [172 MB] 20000000 reads parsed INFO @ 10 Jun 2025 23:28:46: [179 MB] 21000000 reads parsed INFO @ 10 Jun 2025 23:28:47: [184 MB] 22000000 reads parsed INFO @ 10 Jun 2025 23:28:48: [187 MB] 23000000 reads parsed INFO @ 10 Jun 2025 23:28:50: [187 MB] 24000000 reads parsed INFO @ 10 Jun 2025 23:28:52: [190 MB] 25000000 reads parsed INFO @ 10 Jun 2025 23:28:54: [200 MB] 26000000 reads parsed INFO @ 10 Jun 2025 23:28:55: [203 MB] 27000000 reads parsed INFO @ 10 Jun 2025 23:28:57: [208 MB] 28000000 reads parsed INFO @ 10 Jun 2025 23:28:59: [215 MB] 29000000 reads parsed INFO @ 10 Jun 2025 23:29:00: [218 MB] 30000000 reads parsed INFO @ 10 Jun 2025 23:29:02: [218 MB] 31000000 reads parsed INFO @ 10 Jun 2025 23:29:04: [222 MB] 32000000 reads parsed INFO @ 10 Jun 2025 23:29:05: [226 MB] 33000000 reads parsed INFO @ 10 Jun 2025 23:29:07: [228 MB] 34000000 reads parsed INFO @ 10 Jun 2025 23:29:09: [236 MB] 35000000 reads parsed INFO @ 10 Jun 2025 23:29:11: [240 MB] 36000000 reads parsed INFO @ 10 Jun 2025 23:29:12: [242 MB] 37000000 reads parsed INFO @ 10 Jun 2025 23:29:14: [244 MB] 38000000 reads parsed INFO @ 10 Jun 2025 23:29:15: [249 MB] 39000000 reads parsed INFO @ 10 Jun 2025 23:29:17: [253 MB] 40000000 reads parsed INFO @ 10 Jun 2025 23:29:19: [257 MB] 41000000 reads parsed INFO @ 10 Jun 2025 23:29:20: [260 MB] 42000000 reads parsed INFO @ 10 Jun 2025 23:29:22: [264 MB] 43000000 reads parsed INFO @ 10 Jun 2025 23:29:24: [271 MB] 44000000 reads parsed INFO @ 10 Jun 2025 23:29:25: [274 MB] 45000000 reads parsed INFO @ 10 Jun 2025 23:29:27: [278 MB] 46000000 reads parsed INFO @ 10 Jun 2025 23:29:29: [283 MB] 47000000 reads parsed INFO @ 10 Jun 2025 23:29:31: [286 MB] 48000000 reads parsed INFO @ 10 Jun 2025 23:29:32: [290 MB] 49000000 reads parsed INFO @ 10 Jun 2025 23:29:34: [293 MB] 50000000 reads parsed INFO @ 10 Jun 2025 23:29:35: [297 MB] 51000000 reads parsed INFO @ 10 Jun 2025 23:29:36: [302 MB] 52000000 reads parsed INFO @ 10 Jun 2025 23:29:38: [305 MB] 53000000 reads parsed INFO @ 10 Jun 2025 23:29:39: [309 MB] 54000000 reads parsed INFO @ 10 Jun 2025 23:29:41: [313 MB] 55000000 reads parsed INFO @ 10 Jun 2025 23:29:42: [317 MB] 56000000 reads parsed INFO @ 10 Jun 2025 23:29:44: [320 MB] 57000000 reads parsed INFO @ 10 Jun 2025 23:29:45: [325 MB] 58000000 reads parsed INFO @ 10 Jun 2025 23:29:47: [330 MB] 59000000 reads parsed INFO @ 10 Jun 2025 23:29:48: [334 MB] 60000000 reads parsed INFO @ 10 Jun 2025 23:29:50: [338 MB] 61000000 reads parsed INFO @ 10 Jun 2025 23:29:52: [343 MB] 62000000 reads parsed INFO @ 10 Jun 2025 23:29:53: [348 MB] 63000000 reads parsed INFO @ 10 Jun 2025 23:29:55: [352 MB] 64000000 reads parsed INFO @ 10 Jun 2025 23:29:57: [356 MB] 65000000 reads parsed INFO @ 10 Jun 2025 23:29:58: [360 MB] 66000000 reads parsed INFO @ 10 Jun 2025 23:29:59: [362 MB] 67000000 reads parsed INFO @ 10 Jun 2025 23:30:01: [367 MB] 68000000 reads parsed INFO @ 10 Jun 2025 23:30:03: [371 MB] 69000000 reads parsed INFO @ 10 Jun 2025 23:30:05: [375 MB] 70000000 reads parsed INFO @ 10 Jun 2025 23:30:06: [378 MB] 71000000 reads parsed INFO @ 10 Jun 2025 23:30:08: [385 MB] 72000000 reads parsed INFO @ 10 Jun 2025 23:30:09: [389 MB] 73000000 reads parsed INFO @ 10 Jun 2025 23:30:11: [390 MB] 74000000 reads parsed INFO @ 10 Jun 2025 23:30:12: [394 MB] 75000000 reads parsed INFO @ 10 Jun 2025 23:30:14: [398 MB] 76000000 reads parsed INFO @ 10 Jun 2025 23:30:16: [402 MB] 77000000 reads parsed INFO @ 10 Jun 2025 23:30:17: [406 MB] 78000000 reads parsed INFO @ 10 Jun 2025 23:30:19: [410 MB] 79000000 reads parsed INFO @ 10 Jun 2025 23:30:21: [413 MB] 80000000 reads parsed INFO @ 10 Jun 2025 23:30:22: [417 MB] 81000000 reads parsed INFO @ 10 Jun 2025 23:30:24: [420 MB] 82000000 reads parsed INFO @ 10 Jun 2025 23:30:25: [425 MB] 83000000 reads parsed INFO @ 10 Jun 2025 23:30:27: [429 MB] 84000000 reads parsed INFO @ 10 Jun 2025 23:30:28: [433 MB] 85000000 reads parsed INFO @ 10 Jun 2025 23:30:30: [436 MB] 86000000 reads parsed INFO @ 10 Jun 2025 23:30:31: [441 MB] 87000000 reads parsed INFO @ 10 Jun 2025 23:30:33: [444 MB] 88000000 reads parsed INFO @ 10 Jun 2025 23:30:35: [448 MB] 89000000 reads parsed INFO @ 10 Jun 2025 23:30:36: [452 MB] 90000000 reads parsed INFO @ 10 Jun 2025 23:30:38: [456 MB] 91000000 reads parsed INFO @ 10 Jun 2025 23:30:40: [460 MB] 92000000 reads parsed INFO @ 10 Jun 2025 23:30:41: [464 MB] 93000000 reads parsed INFO @ 10 Jun 2025 23:30:43: [467 MB] 94000000 reads parsed INFO @ 10 Jun 2025 23:30:44: [470 MB] 95000000 reads parsed INFO @ 10 Jun 2025 23:30:46: [475 MB] 96000000 reads parsed INFO @ 10 Jun 2025 23:30:47: [479 MB] 97000000 reads parsed INFO @ 10 Jun 2025 23:30:49: [483 MB] 98000000 reads parsed INFO @ 10 Jun 2025 23:30:50: [487 MB] 99000000 reads parsed INFO @ 10 Jun 2025 23:30:51: [491 MB] 100000000 reads parsed INFO @ 10 Jun 2025 23:30:53: [492 MB] #1 tag size is determined as 1 bps INFO @ 10 Jun 2025 23:30:53: [492 MB] #1 tag size = 1.0 INFO @ 10 Jun 2025 23:30:53: [492 MB] #1 total tags in treatment: 100706523 INFO @ 10 Jun 2025 23:30:53: [492 MB] #1 finished! INFO @ 10 Jun 2025 23:30:53: [492 MB] #2 Build Peak Model... INFO @ 10 Jun 2025 23:30:53: [492 MB] #2 Skipped... INFO @ 10 Jun 2025 23:30:53: [492 MB] #2 Use 150 as fragment length INFO @ 10 Jun 2025 23:30:53: [492 MB] #2 Sequencing ends will be shifted towards 5' by 75 bp(s) INFO @ 10 Jun 2025 23:30:53: [492 MB] #3 Call peaks... INFO @ 10 Jun 2025 23:30:53: [492 MB] #3 Going to call summits inside each peak ... INFO @ 10 Jun 2025 23:30:53: [492 MB] #3 Call peaks with given -log10pvalue cutoff: 2.00000 ... INFO @ 10 Jun 2025 23:30:53: [492 MB] #3 Pre-compute pvalue-qvalue table... INFO @ 10 Jun 2025 23:34:43: [2425 MB] #3 In the peak calling step, the following will be performed simultaneously: INFO @ 10 Jun 2025 23:34:43: [2425 MB] #3 Write bedGraph files for treatment pileup (after scaling if necessary)... 1_pseudoreplicate1_treat_pileup.bdg INFO @ 10 Jun 2025 23:34:43: [2425 MB] #3 Write bedGraph files for control lambda (after scaling if necessary)... 1_pseudoreplicate1_control_lambda.bdg INFO @ 10 Jun 2025 23:34:43: [2425 MB] #3 --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. INFO @ 10 Jun 2025 23:34:43: [2425 MB] #3 Call peaks for each chromosome... INFO @ 10 Jun 2025 23:51:20: [2425 MB] #4 Write output xls file... 3_peaks/1_pseudoreplicate1_peaks.xls INFO @ 10 Jun 2025 23:51:23: [2425 MB] #4 Write peak in narrowPeak format file... 3_peaks/1_pseudoreplicate1_peaks.narrowPeak INFO @ 10 Jun 2025 23:51:25: [2425 MB] #4 Write summits bed file... 3_peaks/1_pseudoreplicate1_summits.bed INFO @ 10 Jun 2025 23:51:27: [2425 MB] Done!