The following have been reloaded with a version change: 1) GCCcore/9.3.0 => GCCcore/8.3.0 2) XZ/5.2.5-GCCcore-9.3.0 => XZ/5.2.4-GCCcore-8.3.0 3) binutils/.2.34-GCCcore-9.3.0 => binutils/.2.32-GCCcore-8.3.0 4) bzip2/1.0.8-GCCcore-9.3.0 => bzip2/1.0.8-GCCcore-8.3.0 5) ncurses/.6.2-GCCcore-9.3.0 => ncurses/.6.1-GCCcore-8.3.0 6) zlib/1.2.11-GCCcore-9.3.0 => zlib/1.2.11-GCCcore-8.3.0 INFO @ Tue, 15 Jun 2021 21:58:19: # Command line: callpeak -t H526_ASCL2_S14_R1_001.fastq.gz_bt2_snstv.sorted.rd.bam -c H526_input_S13_R1_001.fastq.gz_bt2_snstv.sorted.rd.bam -g hs -q 0.05 -n H526_ASCL2_S14_R1_001.fastq.gz # ARGUMENTS LIST: # name = H526_ASCL2_S14_R1_001.fastq.gz # format = AUTO # ChIP-seq file = ['H526_ASCL2_S14_R1_001.fastq.gz_bt2_snstv.sorted.rd.bam'] # control file = ['H526_input_S13_R1_001.fastq.gz_bt2_snstv.sorted.rd.bam'] # effective genome size = 2.70e+09 # 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 INFO @ Tue, 15 Jun 2021 21:58:19: #1 read tag files... INFO @ Tue, 15 Jun 2021 21:58:19: #1 read treatment tags... INFO @ Tue, 15 Jun 2021 21:58:19: Detected format is: BAM INFO @ Tue, 15 Jun 2021 21:58:19: * Input file is gzipped. INFO @ Tue, 15 Jun 2021 21:58:22: 1000000 INFO @ Tue, 15 Jun 2021 21:58:25: 2000000 INFO @ Tue, 15 Jun 2021 21:58:28: 3000000 INFO @ Tue, 15 Jun 2021 21:58:31: 4000000 INFO @ Tue, 15 Jun 2021 21:58:34: 5000000 INFO @ Tue, 15 Jun 2021 21:58:37: 6000000 INFO @ Tue, 15 Jun 2021 21:58:40: 7000000 INFO @ Tue, 15 Jun 2021 21:58:42: 8000000 INFO @ Tue, 15 Jun 2021 21:58:45: 9000000 INFO @ Tue, 15 Jun 2021 21:58:53: 9630680 reads have been read. INFO @ Tue, 15 Jun 2021 21:58:53: #1.2 read input tags... INFO @ Tue, 15 Jun 2021 21:58:53: Detected format is: BAM INFO @ Tue, 15 Jun 2021 21:58:53: * Input file is gzipped. INFO @ Tue, 15 Jun 2021 21:58:56: 1000000 INFO @ Tue, 15 Jun 2021 21:58:59: 2000000 INFO @ Tue, 15 Jun 2021 21:59:01: 3000000 INFO @ Tue, 15 Jun 2021 21:59:04: 4000000 INFO @ Tue, 15 Jun 2021 21:59:07: 5000000 INFO @ Tue, 15 Jun 2021 21:59:10: 6000000 INFO @ Tue, 15 Jun 2021 21:59:12: 7000000 INFO @ Tue, 15 Jun 2021 21:59:15: 8000000 INFO @ Tue, 15 Jun 2021 21:59:18: 9000000 INFO @ Tue, 15 Jun 2021 21:59:21: 10000000 INFO @ Tue, 15 Jun 2021 21:59:23: 11000000 INFO @ Tue, 15 Jun 2021 21:59:26: 12000000 INFO @ Tue, 15 Jun 2021 21:59:29: 13000000 INFO @ Tue, 15 Jun 2021 21:59:32: 14000000 INFO @ Tue, 15 Jun 2021 21:59:35: 15000000 INFO @ Tue, 15 Jun 2021 21:59:37: 16000000 INFO @ Tue, 15 Jun 2021 21:59:40: 17000000 INFO @ Tue, 15 Jun 2021 21:59:43: 18000000 INFO @ Tue, 15 Jun 2021 21:59:46: 19000000 INFO @ Tue, 15 Jun 2021 21:59:49: 20000000 INFO @ Tue, 15 Jun 2021 21:59:51: 21000000 INFO @ Tue, 15 Jun 2021 21:59:54: 22000000 INFO @ Tue, 15 Jun 2021 21:59:57: 23000000 INFO @ Tue, 15 Jun 2021 22:00:00: 24000000 INFO @ Tue, 15 Jun 2021 22:00:07: 24660004 reads have been read. INFO @ Tue, 15 Jun 2021 22:00:08: #1 tag size is determined as 76 bps INFO @ Tue, 15 Jun 2021 22:00:08: #1 tag size = 76.0 INFO @ Tue, 15 Jun 2021 22:00:08: #1 total tags in treatment: 9630680 INFO @ Tue, 15 Jun 2021 22:00:08: #1 user defined the maximum tags... INFO @ Tue, 15 Jun 2021 22:00:08: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 15 Jun 2021 22:00:08: #1 tags after filtering in treatment: 9630658 INFO @ Tue, 15 Jun 2021 22:00:08: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 15 Jun 2021 22:00:08: #1 total tags in control: 24660004 INFO @ Tue, 15 Jun 2021 22:00:08: #1 user defined the maximum tags... INFO @ Tue, 15 Jun 2021 22:00:08: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 15 Jun 2021 22:00:08: #1 tags after filtering in control: 24660003 INFO @ Tue, 15 Jun 2021 22:00:08: #1 Redundant rate of control: 0.00 INFO @ Tue, 15 Jun 2021 22:00:08: #1 finished! INFO @ Tue, 15 Jun 2021 22:00:08: #2 Build Peak Model... INFO @ Tue, 15 Jun 2021 22:00:08: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 15 Jun 2021 22:00:11: #2 number of paired peaks: 27670 INFO @ Tue, 15 Jun 2021 22:00:11: start model_add_line... INFO @ Tue, 15 Jun 2021 22:00:11: start X-correlation... INFO @ Tue, 15 Jun 2021 22:00:11: end of X-cor INFO @ Tue, 15 Jun 2021 22:00:11: #2 finished! INFO @ Tue, 15 Jun 2021 22:00:11: #2 predicted fragment length is 88 bps INFO @ Tue, 15 Jun 2021 22:00:11: #2 alternative fragment length(s) may be 88 bps INFO @ Tue, 15 Jun 2021 22:00:11: #2.2 Generate R script for model : H526_ASCL2_S14_R1_001.fastq.gz_model.r WARNING @ Tue, 15 Jun 2021 22:00:11: #2 Since the d (88) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 15 Jun 2021 22:00:11: #2 You may need to consider one of the other alternative d(s): 88 WARNING @ Tue, 15 Jun 2021 22:00:11: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Tue, 15 Jun 2021 22:00:11: #3 Call peaks... INFO @ Tue, 15 Jun 2021 22:00:11: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 15 Jun 2021 22:01:58: #3 Call peaks for each chromosome... INFO @ Tue, 15 Jun 2021 22:02:54: #4 Write output xls file... H526_ASCL2_S14_R1_001.fastq.gz_peaks.xls INFO @ Tue, 15 Jun 2021 22:02:54: #4 Write peak in narrowPeak format file... H526_ASCL2_S14_R1_001.fastq.gz_peaks.narrowPeak INFO @ Tue, 15 Jun 2021 22:02:54: #4 Write summits bed file... H526_ASCL2_S14_R1_001.fastq.gz_summits.bed INFO @ Tue, 15 Jun 2021 22:02:54: Done!