INFO @ 20 Dec 2024 06:52:56: [31 MB] # Command line: callpeak -t 2_celltype_fragments/pseudorep2/C7_sorted.tsv -f BED -n C7_pseudoreplicate2 -g hs --outdir 3_peaks -p 0.01 --shift -75 --extsize 150 --nomodel -B --SPMR --keep-dup all --call-summits # ARGUMENTS LIST: # name = C7_pseudoreplicate2 # format = BED # ChIP-seq file = ['2_celltype_fragments/pseudorep2/C7_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 @ 20 Dec 2024 06:52:56: [31 MB] #1 read tag files... INFO @ 20 Dec 2024 06:52:56: [31 MB] #1 read treatment tags... INFO @ 20 Dec 2024 06:52:58: [98 MB] 1000000 reads parsed INFO @ 20 Dec 2024 06:52:59: [101 MB] 2000000 reads parsed INFO @ 20 Dec 2024 06:53:01: [103 MB] 3000000 reads parsed INFO @ 20 Dec 2024 06:53:03: [105 MB] 4000000 reads parsed INFO @ 20 Dec 2024 06:53:05: [114 MB] 5000000 reads parsed INFO @ 20 Dec 2024 06:53:06: [117 MB] 6000000 reads parsed INFO @ 20 Dec 2024 06:53:08: [124 MB] 7000000 reads parsed INFO @ 20 Dec 2024 06:53:10: [127 MB] 8000000 reads parsed INFO @ 20 Dec 2024 06:53:12: [135 MB] 9000000 reads parsed INFO @ 20 Dec 2024 06:53:14: [139 MB] 10000000 reads parsed INFO @ 20 Dec 2024 06:53:16: [142 MB] 11000000 reads parsed INFO @ 20 Dec 2024 06:53:18: [147 MB] 12000000 reads parsed INFO @ 20 Dec 2024 06:53:19: [151 MB] 13000000 reads parsed INFO @ 20 Dec 2024 06:53:21: [156 MB] 14000000 reads parsed INFO @ 20 Dec 2024 06:53:23: [160 MB] 15000000 reads parsed INFO @ 20 Dec 2024 06:53:25: [162 MB] 16000000 reads parsed INFO @ 20 Dec 2024 06:53:27: [168 MB] 17000000 reads parsed INFO @ 20 Dec 2024 06:53:29: [171 MB] 18000000 reads parsed INFO @ 20 Dec 2024 06:53:31: [178 MB] 19000000 reads parsed INFO @ 20 Dec 2024 06:53:32: [180 MB] 20000000 reads parsed INFO @ 20 Dec 2024 06:53:34: [186 MB] 21000000 reads parsed INFO @ 20 Dec 2024 06:53:36: [190 MB] 22000000 reads parsed INFO @ 20 Dec 2024 06:53:38: [192 MB] 23000000 reads parsed INFO @ 20 Dec 2024 06:53:40: [194 MB] 24000000 reads parsed INFO @ 20 Dec 2024 06:53:42: [198 MB] 25000000 reads parsed INFO @ 20 Dec 2024 06:53:43: [206 MB] 26000000 reads parsed INFO @ 20 Dec 2024 06:53:45: [211 MB] 27000000 reads parsed INFO @ 20 Dec 2024 06:53:47: [215 MB] 28000000 reads parsed INFO @ 20 Dec 2024 06:53:48: [217 MB] 29000000 reads parsed INFO @ 20 Dec 2024 06:53:50: [222 MB] 30000000 reads parsed INFO @ 20 Dec 2024 06:53:52: [225 MB] 31000000 reads parsed INFO @ 20 Dec 2024 06:53:54: [231 MB] 32000000 reads parsed INFO @ 20 Dec 2024 06:53:56: [234 MB] 33000000 reads parsed INFO @ 20 Dec 2024 06:53:57: [238 MB] 34000000 reads parsed INFO @ 20 Dec 2024 06:53:59: [242 MB] 35000000 reads parsed INFO @ 20 Dec 2024 06:54:01: [246 MB] 36000000 reads parsed INFO @ 20 Dec 2024 06:54:03: [250 MB] 37000000 reads parsed INFO @ 20 Dec 2024 06:54:05: [254 MB] 38000000 reads parsed INFO @ 20 Dec 2024 06:54:07: [258 MB] 39000000 reads parsed INFO @ 20 Dec 2024 06:54:09: [264 MB] 40000000 reads parsed INFO @ 20 Dec 2024 06:54:11: [267 MB] 41000000 reads parsed INFO @ 20 Dec 2024 06:54:13: [271 MB] 42000000 reads parsed INFO @ 20 Dec 2024 06:54:13: [271 MB] #1 tag size is determined as 1 bps INFO @ 20 Dec 2024 06:54:13: [271 MB] #1 tag size = 1.0 INFO @ 20 Dec 2024 06:54:13: [271 MB] #1 total tags in treatment: 42127722 INFO @ 20 Dec 2024 06:54:13: [271 MB] #1 finished! INFO @ 20 Dec 2024 06:54:13: [271 MB] #2 Build Peak Model... INFO @ 20 Dec 2024 06:54:13: [271 MB] #2 Skipped... INFO @ 20 Dec 2024 06:54:13: [271 MB] #2 Use 150 as fragment length INFO @ 20 Dec 2024 06:54:13: [271 MB] #2 Sequencing ends will be shifted towards 5' by 75 bp(s) INFO @ 20 Dec 2024 06:54:13: [271 MB] #3 Call peaks... INFO @ 20 Dec 2024 06:54:13: [271 MB] #3 Going to call summits inside each peak ... INFO @ 20 Dec 2024 06:54:13: [271 MB] #3 Call peaks with given -log10pvalue cutoff: 2.00000 ... INFO @ 20 Dec 2024 06:54:13: [271 MB] #3 Pre-compute pvalue-qvalue table... INFO @ 20 Dec 2024 06:55:28: [1143 MB] #3 In the peak calling step, the following will be performed simultaneously: INFO @ 20 Dec 2024 06:55:28: [1143 MB] #3 Write bedGraph files for treatment pileup (after scaling if necessary)... C7_pseudoreplicate2_treat_pileup.bdg INFO @ 20 Dec 2024 06:55:28: [1143 MB] #3 Write bedGraph files for control lambda (after scaling if necessary)... C7_pseudoreplicate2_control_lambda.bdg INFO @ 20 Dec 2024 06:55:28: [1143 MB] #3 --SPMR is requested, so pileup will be normalized by sequencing depth in million reads. INFO @ 20 Dec 2024 06:55:28: [1143 MB] #3 Call peaks for each chromosome... INFO @ 20 Dec 2024 06:57:00: [1143 MB] #4 Write output xls file... 3_peaks/C7_pseudoreplicate2_peaks.xls INFO @ 20 Dec 2024 06:57:01: [1143 MB] #4 Write peak in narrowPeak format file... 3_peaks/C7_pseudoreplicate2_peaks.narrowPeak INFO @ 20 Dec 2024 06:57:02: [1143 MB] #4 Write summits bed file... 3_peaks/C7_pseudoreplicate2_summits.bed INFO @ 20 Dec 2024 06:57:02: [1143 MB] Done!