Job ID = 12264781 SRX = SRX4082387 Genome = ce10 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... Read 23213305 spots for SRR7164205/SRR7164205.sra Written 23213305 spots for SRR7164205/SRR7164205.sra Read 47115866 spots for SRR7164206/SRR7164206.sra Written 47115866 spots for SRR7164206/SRR7164206.sra fastq に変換しました。 bowtie でマッピング中... Your job 12265505 ("srTce11") has been submitted Time loading reference: 00:00:00 Time loading forward index: 00:00:00 Time loading mirror index: 00:00:00 Multiseed full-index search: 00:51:26 70329171 reads; of these: 70329171 (100.00%) were paired; of these: 3188725 (4.53%) aligned concordantly 0 times 61179946 (86.99%) aligned concordantly exactly 1 time 5960500 (8.48%) aligned concordantly >1 times ---- 3188725 pairs aligned concordantly 0 times; of these: 1407010 (44.12%) aligned discordantly 1 time ---- 1781715 pairs aligned 0 times concordantly or discordantly; of these: 3563430 mates make up the pairs; of these: 2271086 (63.73%) aligned 0 times 863392 (24.23%) aligned exactly 1 time 428952 (12.04%) aligned >1 times 98.39% overall alignment rate Time searching: 00:51:26 Overall time: 00:51:26 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 7 sequences. [bam_sort_core] merging from 36 files... [bam_rmdup_core] processing reference chrI... [bam_rmdup_core] processing reference chrII... [bam_rmdup_core] processing reference chrIII... [bam_rmdup_core] processing reference chrIV... [bam_rmdup_core] processing reference chrV... [bam_rmdup_core] processing reference chrX... [bam_rmdup_core] processing reference chrM... [bam_rmdup_core] 30516902 / 68464877 = 0.4457 in library ' ' BAM に変換しました。 Bed ファイルを作成中... WARNING: Skipping mount /opt/pkg/singularity/3.7.1/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Sat, 03 Apr 2021 07:21:11: # Command line: callpeak -t /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam -f BAM -g ce -n /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam'] # control file = None # effective genome size = 9.00e+07 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-05 # 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 INFO @ Sat, 03 Apr 2021 07:21:11: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:21:11: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:21:16: 1000000 INFO @ Sat, 03 Apr 2021 07:21:21: 2000000 INFO @ Sat, 03 Apr 2021 07:21:28: 3000000 INFO @ Sat, 03 Apr 2021 07:21:36: 4000000 WARNING: Skipping mount /opt/pkg/singularity/3.7.1/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Sat, 03 Apr 2021 07:21:40: # Command line: callpeak -t /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam -f BAM -g ce -n /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam'] # control file = None # effective genome size = 9.00e+07 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-10 # 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 INFO @ Sat, 03 Apr 2021 07:21:40: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:21:40: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:21:41: 5000000 INFO @ Sat, 03 Apr 2021 07:21:45: 1000000 INFO @ Sat, 03 Apr 2021 07:21:46: 6000000 INFO @ Sat, 03 Apr 2021 07:21:52: 2000000 INFO @ Sat, 03 Apr 2021 07:21:52: 7000000 INFO @ Sat, 03 Apr 2021 07:21:58: 3000000 INFO @ Sat, 03 Apr 2021 07:21:59: 8000000 INFO @ Sat, 03 Apr 2021 07:22:05: 9000000 INFO @ Sat, 03 Apr 2021 07:22:05: 4000000 BedGraph に変換中... WARNING: Skipping mount /opt/pkg/singularity/3.7.1/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Sat, 03 Apr 2021 07:22:11: # Command line: callpeak -t /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam -f BAM -g ce -n /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.bam'] # control file = None # effective genome size = 9.00e+07 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-20 # 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 INFO @ Sat, 03 Apr 2021 07:22:11: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:22:11: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:22:12: 10000000 INFO @ Sat, 03 Apr 2021 07:22:12: 5000000 INFO @ Sat, 03 Apr 2021 07:22:17: 1000000 INFO @ Sat, 03 Apr 2021 07:22:20: 11000000 INFO @ Sat, 03 Apr 2021 07:22:20: 6000000 INFO @ Sat, 03 Apr 2021 07:22:24: 2000000 INFO @ Sat, 03 Apr 2021 07:22:26: 12000000 INFO @ Sat, 03 Apr 2021 07:22:27: 7000000 INFO @ Sat, 03 Apr 2021 07:22:32: 3000000 INFO @ Sat, 03 Apr 2021 07:22:34: 8000000 INFO @ Sat, 03 Apr 2021 07:22:34: 13000000 INFO @ Sat, 03 Apr 2021 07:22:39: 4000000 INFO @ Sat, 03 Apr 2021 07:22:41: 9000000 INFO @ Sat, 03 Apr 2021 07:22:42: 14000000 INFO @ Sat, 03 Apr 2021 07:22:46: 5000000 INFO @ Sat, 03 Apr 2021 07:22:47: 10000000 INFO @ Sat, 03 Apr 2021 07:22:49: 15000000 INFO @ Sat, 03 Apr 2021 07:22:53: 6000000 INFO @ Sat, 03 Apr 2021 07:22:55: 11000000 INFO @ Sat, 03 Apr 2021 07:22:56: 16000000 INFO @ Sat, 03 Apr 2021 07:23:01: 7000000 INFO @ Sat, 03 Apr 2021 07:23:01: 12000000 INFO @ Sat, 03 Apr 2021 07:23:04: 17000000 INFO @ Sat, 03 Apr 2021 07:23:08: 13000000 INFO @ Sat, 03 Apr 2021 07:23:08: 8000000 INFO @ Sat, 03 Apr 2021 07:23:11: 18000000 INFO @ Sat, 03 Apr 2021 07:23:15: 9000000 INFO @ Sat, 03 Apr 2021 07:23:15: 14000000 INFO @ Sat, 03 Apr 2021 07:23:20: 19000000 INFO @ Sat, 03 Apr 2021 07:23:22: 15000000 INFO @ Sat, 03 Apr 2021 07:23:23: 10000000 INFO @ Sat, 03 Apr 2021 07:23:27: 20000000 INFO @ Sat, 03 Apr 2021 07:23:29: 16000000 INFO @ Sat, 03 Apr 2021 07:23:32: 11000000 INFO @ Sat, 03 Apr 2021 07:23:34: 21000000 INFO @ Sat, 03 Apr 2021 07:23:36: 17000000 INFO @ Sat, 03 Apr 2021 07:23:41: 22000000 INFO @ Sat, 03 Apr 2021 07:23:43: 12000000 INFO @ Sat, 03 Apr 2021 07:23:43: 18000000 INFO @ Sat, 03 Apr 2021 07:23:47: 23000000 INFO @ Sat, 03 Apr 2021 07:23:50: 19000000 INFO @ Sat, 03 Apr 2021 07:23:53: 24000000 INFO @ Sat, 03 Apr 2021 07:23:54: 13000000 INFO @ Sat, 03 Apr 2021 07:23:56: 20000000 INFO @ Sat, 03 Apr 2021 07:24:00: 25000000 INFO @ Sat, 03 Apr 2021 07:24:03: 21000000 INFO @ Sat, 03 Apr 2021 07:24:03: 14000000 INFO @ Sat, 03 Apr 2021 07:24:06: 26000000 INFO @ Sat, 03 Apr 2021 07:24:12: 22000000 INFO @ Sat, 03 Apr 2021 07:24:13: 15000000 INFO @ Sat, 03 Apr 2021 07:24:14: 27000000 INFO @ Sat, 03 Apr 2021 07:24:19: 23000000 INFO @ Sat, 03 Apr 2021 07:24:20: 28000000 INFO @ Sat, 03 Apr 2021 07:24:21: 16000000 INFO @ Sat, 03 Apr 2021 07:24:26: 24000000 INFO @ Sat, 03 Apr 2021 07:24:27: 29000000 INFO @ Sat, 03 Apr 2021 07:24:29: 17000000 INFO @ Sat, 03 Apr 2021 07:24:33: 25000000 INFO @ Sat, 03 Apr 2021 07:24:34: 30000000 INFO @ Sat, 03 Apr 2021 07:24:36: 18000000 INFO @ Sat, 03 Apr 2021 07:24:39: 26000000 INFO @ Sat, 03 Apr 2021 07:24:40: 31000000 INFO @ Sat, 03 Apr 2021 07:24:44: 19000000 INFO @ Sat, 03 Apr 2021 07:24:46: 27000000 INFO @ Sat, 03 Apr 2021 07:24:47: 32000000 INFO @ Sat, 03 Apr 2021 07:24:51: 20000000 INFO @ Sat, 03 Apr 2021 07:24:53: 28000000 INFO @ Sat, 03 Apr 2021 07:24:54: 33000000 INFO @ Sat, 03 Apr 2021 07:24:59: 21000000 INFO @ Sat, 03 Apr 2021 07:25:00: 29000000 INFO @ Sat, 03 Apr 2021 07:25:01: 34000000 INFO @ Sat, 03 Apr 2021 07:25:06: 30000000 INFO @ Sat, 03 Apr 2021 07:25:07: 22000000 INFO @ Sat, 03 Apr 2021 07:25:07: 35000000 INFO @ Sat, 03 Apr 2021 07:25:13: 31000000 INFO @ Sat, 03 Apr 2021 07:25:13: 36000000 INFO @ Sat, 03 Apr 2021 07:25:14: 23000000 INFO @ Sat, 03 Apr 2021 07:25:20: 37000000 INFO @ Sat, 03 Apr 2021 07:25:20: 32000000 INFO @ Sat, 03 Apr 2021 07:25:21: 24000000 INFO @ Sat, 03 Apr 2021 07:25:27: 38000000 INFO @ Sat, 03 Apr 2021 07:25:27: 33000000 INFO @ Sat, 03 Apr 2021 07:25:29: 25000000 INFO @ Sat, 03 Apr 2021 07:25:33: 39000000 INFO @ Sat, 03 Apr 2021 07:25:35: 34000000 INFO @ Sat, 03 Apr 2021 07:25:35: 26000000 INFO @ Sat, 03 Apr 2021 07:25:40: 40000000 INFO @ Sat, 03 Apr 2021 07:25:43: 27000000 INFO @ Sat, 03 Apr 2021 07:25:44: 35000000 INFO @ Sat, 03 Apr 2021 07:25:47: 41000000 INFO @ Sat, 03 Apr 2021 07:25:50: 28000000 INFO @ Sat, 03 Apr 2021 07:25:53: 42000000 INFO @ Sat, 03 Apr 2021 07:25:54: 36000000 INFO @ Sat, 03 Apr 2021 07:25:57: 29000000 INFO @ Sat, 03 Apr 2021 07:26:00: 43000000 INFO @ Sat, 03 Apr 2021 07:26:03: 37000000 INFO @ Sat, 03 Apr 2021 07:26:05: 30000000 INFO @ Sat, 03 Apr 2021 07:26:07: 44000000 INFO @ Sat, 03 Apr 2021 07:26:12: 31000000 INFO @ Sat, 03 Apr 2021 07:26:12: 38000000 INFO @ Sat, 03 Apr 2021 07:26:15: 45000000 INFO @ Sat, 03 Apr 2021 07:26:19: 32000000 INFO @ Sat, 03 Apr 2021 07:26:21: 39000000 INFO @ Sat, 03 Apr 2021 07:26:23: 46000000 INFO @ Sat, 03 Apr 2021 07:26:26: 33000000 INFO @ Sat, 03 Apr 2021 07:26:30: 40000000 INFO @ Sat, 03 Apr 2021 07:26:32: 47000000 INFO @ Sat, 03 Apr 2021 07:26:33: 34000000 INFO @ Sat, 03 Apr 2021 07:26:39: 41000000 INFO @ Sat, 03 Apr 2021 07:26:40: 35000000 INFO @ Sat, 03 Apr 2021 07:26:41: 48000000 INFO @ Sat, 03 Apr 2021 07:26:48: 36000000 INFO @ Sat, 03 Apr 2021 07:26:48: 42000000 INFO @ Sat, 03 Apr 2021 07:26:49: 49000000 INFO @ Sat, 03 Apr 2021 07:26:55: 37000000 INFO @ Sat, 03 Apr 2021 07:26:56: 50000000 INFO @ Sat, 03 Apr 2021 07:26:57: 43000000 INFO @ Sat, 03 Apr 2021 07:27:02: 51000000 INFO @ Sat, 03 Apr 2021 07:27:03: 38000000 INFO @ Sat, 03 Apr 2021 07:27:06: 44000000 INFO @ Sat, 03 Apr 2021 07:27:08: 52000000 INFO @ Sat, 03 Apr 2021 07:27:10: 39000000 INFO @ Sat, 03 Apr 2021 07:27:14: 53000000 INFO @ Sat, 03 Apr 2021 07:27:15: 45000000 INFO @ Sat, 03 Apr 2021 07:27:18: 40000000 INFO @ Sat, 03 Apr 2021 07:27:20: 54000000 INFO @ Sat, 03 Apr 2021 07:27:22: 46000000 INFO @ Sat, 03 Apr 2021 07:27:26: 41000000 INFO @ Sat, 03 Apr 2021 07:27:27: 55000000 INFO @ Sat, 03 Apr 2021 07:27:31: 47000000 INFO @ Sat, 03 Apr 2021 07:27:33: 56000000 INFO @ Sat, 03 Apr 2021 07:27:34: 42000000 INFO @ Sat, 03 Apr 2021 07:27:38: 48000000 INFO @ Sat, 03 Apr 2021 07:27:39: 57000000 INFO @ Sat, 03 Apr 2021 07:27:41: 43000000 INFO @ Sat, 03 Apr 2021 07:27:46: 49000000 INFO @ Sat, 03 Apr 2021 07:27:46: 58000000 INFO @ Sat, 03 Apr 2021 07:27:48: 44000000 INFO @ Sat, 03 Apr 2021 07:27:53: 59000000 INFO @ Sat, 03 Apr 2021 07:27:53: 50000000 INFO @ Sat, 03 Apr 2021 07:27:55: 45000000 INFO @ Sat, 03 Apr 2021 07:27:59: 60000000 INFO @ Sat, 03 Apr 2021 07:28:01: 51000000 INFO @ Sat, 03 Apr 2021 07:28:03: 46000000 INFO @ Sat, 03 Apr 2021 07:28:05: 61000000 INFO @ Sat, 03 Apr 2021 07:28:08: 52000000 INFO @ Sat, 03 Apr 2021 07:28:10: 47000000 INFO @ Sat, 03 Apr 2021 07:28:11: 62000000 INFO @ Sat, 03 Apr 2021 07:28:15: 53000000 INFO @ Sat, 03 Apr 2021 07:28:17: 63000000 INFO @ Sat, 03 Apr 2021 07:28:17: 48000000 INFO @ Sat, 03 Apr 2021 07:28:22: 54000000 INFO @ Sat, 03 Apr 2021 07:28:23: 64000000 INFO @ Sat, 03 Apr 2021 07:28:25: 49000000 INFO @ Sat, 03 Apr 2021 07:28:30: 65000000 INFO @ Sat, 03 Apr 2021 07:28:31: 55000000 INFO @ Sat, 03 Apr 2021 07:28:33: 50000000 INFO @ Sat, 03 Apr 2021 07:28:36: 66000000 INFO @ Sat, 03 Apr 2021 07:28:39: 56000000 INFO @ Sat, 03 Apr 2021 07:28:41: 51000000 INFO @ Sat, 03 Apr 2021 07:28:42: 67000000 INFO @ Sat, 03 Apr 2021 07:28:46: 57000000 INFO @ Sat, 03 Apr 2021 07:28:48: 52000000 INFO @ Sat, 03 Apr 2021 07:28:49: 68000000 INFO @ Sat, 03 Apr 2021 07:28:53: 58000000 INFO @ Sat, 03 Apr 2021 07:28:55: 69000000 INFO @ Sat, 03 Apr 2021 07:28:56: 53000000 INFO @ Sat, 03 Apr 2021 07:29:01: 59000000 INFO @ Sat, 03 Apr 2021 07:29:02: 70000000 INFO @ Sat, 03 Apr 2021 07:29:04: 54000000 INFO @ Sat, 03 Apr 2021 07:29:08: 60000000 INFO @ Sat, 03 Apr 2021 07:29:08: 71000000 INFO @ Sat, 03 Apr 2021 07:29:10: 55000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Sat, 03 Apr 2021 07:29:14: 72000000 INFO @ Sat, 03 Apr 2021 07:29:15: 61000000 INFO @ Sat, 03 Apr 2021 07:29:17: 56000000 INFO @ Sat, 03 Apr 2021 07:29:21: 73000000 INFO @ Sat, 03 Apr 2021 07:29:22: 62000000 INFO @ Sat, 03 Apr 2021 07:29:24: 57000000 INFO @ Sat, 03 Apr 2021 07:29:27: 74000000 INFO @ Sat, 03 Apr 2021 07:29:29: 63000000 INFO @ Sat, 03 Apr 2021 07:29:31: 58000000 INFO @ Sat, 03 Apr 2021 07:29:32: 75000000 INFO @ Sat, 03 Apr 2021 07:29:36: 64000000 INFO @ Sat, 03 Apr 2021 07:29:38: 59000000 INFO @ Sat, 03 Apr 2021 07:29:38: 76000000 INFO @ Sat, 03 Apr 2021 07:29:43: 65000000 INFO @ Sat, 03 Apr 2021 07:29:43: 77000000 INFO @ Sat, 03 Apr 2021 07:29:44: 60000000 INFO @ Sat, 03 Apr 2021 07:29:46: #1 tag size is determined as 32 bps INFO @ Sat, 03 Apr 2021 07:29:46: #1 tag size = 32 INFO @ Sat, 03 Apr 2021 07:29:46: #1 total tags in treatment: 36762111 INFO @ Sat, 03 Apr 2021 07:29:46: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:29:46: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 03 Apr 2021 07:29:47: #1 tags after filtering in treatment: 19919695 INFO @ Sat, 03 Apr 2021 07:29:47: #1 Redundant rate of treatment: 0.46 INFO @ Sat, 03 Apr 2021 07:29:47: #1 finished! INFO @ Sat, 03 Apr 2021 07:29:47: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:29:47: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:29:48: #2 number of paired peaks: 234 WARNING @ Sat, 03 Apr 2021 07:29:48: Fewer paired peaks (234) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 234 pairs to build model! INFO @ Sat, 03 Apr 2021 07:29:48: start model_add_line... INFO @ Sat, 03 Apr 2021 07:29:48: start X-correlation... INFO @ Sat, 03 Apr 2021 07:29:48: end of X-cor INFO @ Sat, 03 Apr 2021 07:29:48: #2 finished! INFO @ Sat, 03 Apr 2021 07:29:48: #2 predicted fragment length is 100 bps INFO @ Sat, 03 Apr 2021 07:29:48: #2 alternative fragment length(s) may be 100 bps INFO @ Sat, 03 Apr 2021 07:29:48: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05_model.r INFO @ Sat, 03 Apr 2021 07:29:48: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:29:48: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:29:50: 66000000 INFO @ Sat, 03 Apr 2021 07:29:51: 61000000 INFO @ Sat, 03 Apr 2021 07:29:57: 67000000 INFO @ Sat, 03 Apr 2021 07:29:58: 62000000 INFO @ Sat, 03 Apr 2021 07:30:04: 68000000 INFO @ Sat, 03 Apr 2021 07:30:05: 63000000 INFO @ Sat, 03 Apr 2021 07:30:10: 69000000 INFO @ Sat, 03 Apr 2021 07:30:12: 64000000 INFO @ Sat, 03 Apr 2021 07:30:17: 70000000 INFO @ Sat, 03 Apr 2021 07:30:18: 65000000 INFO @ Sat, 03 Apr 2021 07:30:24: 71000000 INFO @ Sat, 03 Apr 2021 07:30:25: 66000000 INFO @ Sat, 03 Apr 2021 07:30:25: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:30:31: 72000000 INFO @ Sat, 03 Apr 2021 07:30:32: 67000000 INFO @ Sat, 03 Apr 2021 07:30:38: 73000000 INFO @ Sat, 03 Apr 2021 07:30:38: 68000000 BigWig に変換しました。 INFO @ Sat, 03 Apr 2021 07:30:44: #4 Write output xls file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05_peaks.xls INFO @ Sat, 03 Apr 2021 07:30:44: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:30:44: #4 Write summits bed file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.05_summits.bed INFO @ Sat, 03 Apr 2021 07:30:44: Done! pass1 - making usageList (7 chroms): 4 millis pass2 - checking and writing primary data (10852 records, 4 fields): 15 millis CompletedMACS2peakCalling INFO @ Sat, 03 Apr 2021 07:30:44: 74000000 INFO @ Sat, 03 Apr 2021 07:30:45: 69000000 INFO @ Sat, 03 Apr 2021 07:30:51: 75000000 INFO @ Sat, 03 Apr 2021 07:30:52: 70000000 INFO @ Sat, 03 Apr 2021 07:30:57: 76000000 INFO @ Sat, 03 Apr 2021 07:30:58: 71000000 INFO @ Sat, 03 Apr 2021 07:31:03: 77000000 INFO @ Sat, 03 Apr 2021 07:31:05: 72000000 INFO @ Sat, 03 Apr 2021 07:31:06: #1 tag size is determined as 32 bps INFO @ Sat, 03 Apr 2021 07:31:06: #1 tag size = 32 INFO @ Sat, 03 Apr 2021 07:31:06: #1 total tags in treatment: 36762111 INFO @ Sat, 03 Apr 2021 07:31:06: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:31:06: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 03 Apr 2021 07:31:07: #1 tags after filtering in treatment: 19919695 INFO @ Sat, 03 Apr 2021 07:31:07: #1 Redundant rate of treatment: 0.46 INFO @ Sat, 03 Apr 2021 07:31:07: #1 finished! INFO @ Sat, 03 Apr 2021 07:31:07: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:31:07: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:31:08: #2 number of paired peaks: 234 WARNING @ Sat, 03 Apr 2021 07:31:08: Fewer paired peaks (234) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 234 pairs to build model! INFO @ Sat, 03 Apr 2021 07:31:08: start model_add_line... INFO @ Sat, 03 Apr 2021 07:31:08: start X-correlation... INFO @ Sat, 03 Apr 2021 07:31:08: end of X-cor INFO @ Sat, 03 Apr 2021 07:31:08: #2 finished! INFO @ Sat, 03 Apr 2021 07:31:08: #2 predicted fragment length is 100 bps INFO @ Sat, 03 Apr 2021 07:31:08: #2 alternative fragment length(s) may be 100 bps INFO @ Sat, 03 Apr 2021 07:31:08: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10_model.r INFO @ Sat, 03 Apr 2021 07:31:08: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:31:08: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:31:11: 73000000 INFO @ Sat, 03 Apr 2021 07:31:18: 74000000 INFO @ Sat, 03 Apr 2021 07:31:24: 75000000 INFO @ Sat, 03 Apr 2021 07:31:30: 76000000 INFO @ Sat, 03 Apr 2021 07:31:36: 77000000 INFO @ Sat, 03 Apr 2021 07:31:38: #1 tag size is determined as 32 bps INFO @ Sat, 03 Apr 2021 07:31:38: #1 tag size = 32 INFO @ Sat, 03 Apr 2021 07:31:38: #1 total tags in treatment: 36762111 INFO @ Sat, 03 Apr 2021 07:31:38: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:31:38: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 03 Apr 2021 07:31:39: #1 tags after filtering in treatment: 19919695 INFO @ Sat, 03 Apr 2021 07:31:39: #1 Redundant rate of treatment: 0.46 INFO @ Sat, 03 Apr 2021 07:31:39: #1 finished! INFO @ Sat, 03 Apr 2021 07:31:39: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:31:39: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:31:40: #2 number of paired peaks: 234 WARNING @ Sat, 03 Apr 2021 07:31:40: Fewer paired peaks (234) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 234 pairs to build model! INFO @ Sat, 03 Apr 2021 07:31:40: start model_add_line... INFO @ Sat, 03 Apr 2021 07:31:41: start X-correlation... INFO @ Sat, 03 Apr 2021 07:31:41: end of X-cor INFO @ Sat, 03 Apr 2021 07:31:41: #2 finished! INFO @ Sat, 03 Apr 2021 07:31:41: #2 predicted fragment length is 100 bps INFO @ Sat, 03 Apr 2021 07:31:41: #2 alternative fragment length(s) may be 100 bps INFO @ Sat, 03 Apr 2021 07:31:41: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20_model.r INFO @ Sat, 03 Apr 2021 07:31:41: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:31:41: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:31:45: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:32:04: #4 Write output xls file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10_peaks.xls INFO @ Sat, 03 Apr 2021 07:32:04: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:32:04: #4 Write summits bed file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.10_summits.bed INFO @ Sat, 03 Apr 2021 07:32:04: Done! pass1 - making usageList (7 chroms): 2 millis pass2 - checking and writing primary data (6494 records, 4 fields): 9 millis CompletedMACS2peakCalling INFO @ Sat, 03 Apr 2021 07:32:17: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:32:35: #4 Write output xls file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20_peaks.xls INFO @ Sat, 03 Apr 2021 07:32:35: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:32:35: #4 Write summits bed file... /home/okishinya/chipatlas/results/ce10/SRX4082387/SRX4082387.20_summits.bed INFO @ Sat, 03 Apr 2021 07:32:35: Done! pass1 - making usageList (6 chroms): 1 millis pass2 - checking and writing primary data (3107 records, 4 fields): 6 millis CompletedMACS2peakCalling