Job ID = 6627198 SRX = SRX8521316 Genome = dm3 sra ファイルのダウンロード中... Read layout: SINGLE fastq に変換中... 2020-07-14T01:32:46 prefetch.2.10.7: 1) Downloading 'SRR11978749'... 2020-07-14T01:32:46 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:33:01 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:33:01 prefetch.2.10.7: 'SRR11978749' is valid 2020-07-14T01:33:01 prefetch.2.10.7: 1) 'SRR11978749' was downloaded successfully 2020-07-14T01:33:01 prefetch.2.10.7: 'SRR11978749' has 0 unresolved dependencies Read 177263 spots for SRR11978749/SRR11978749.sra Written 177263 spots for SRR11978749/SRR11978749.sra 2020-07-14T01:33:12 prefetch.2.10.7: 1) Downloading 'SRR11978750'... 2020-07-14T01:33:12 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:34:11 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:34:11 prefetch.2.10.7: 'SRR11978750' is valid 2020-07-14T01:34:11 prefetch.2.10.7: 1) 'SRR11978750' was downloaded successfully 2020-07-14T01:34:11 prefetch.2.10.7: 'SRR11978750' has 0 unresolved dependencies Read 3114721 spots for SRR11978750/SRR11978750.sra Written 3114721 spots for SRR11978750/SRR11978750.sra 2020-07-14T01:34:33 prefetch.2.10.7: 1) Downloading 'SRR11978751'... 2020-07-14T01:34:33 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:35:15 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:35:15 prefetch.2.10.7: 'SRR11978751' is valid 2020-07-14T01:35:15 prefetch.2.10.7: 1) 'SRR11978751' was downloaded successfully 2020-07-14T01:35:15 prefetch.2.10.7: 'SRR11978751' has 0 unresolved dependencies Read 3316781 spots for SRR11978751/SRR11978751.sra Written 3316781 spots for SRR11978751/SRR11978751.sra 2020-07-14T01:35:36 prefetch.2.10.7: 1) Downloading 'SRR11978752'... 2020-07-14T01:35:36 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:36:22 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:36:22 prefetch.2.10.7: 'SRR11978752' is valid 2020-07-14T01:36:22 prefetch.2.10.7: 1) 'SRR11978752' was downloaded successfully 2020-07-14T01:36:22 prefetch.2.10.7: 'SRR11978752' has 0 unresolved dependencies Read 3237458 spots for SRR11978752/SRR11978752.sra Written 3237458 spots for SRR11978752/SRR11978752.sra 2020-07-14T01:36:48 prefetch.2.10.7: 1) Downloading 'SRR11978753'... 2020-07-14T01:36:48 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:37:56 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:37:57 prefetch.2.10.7: 'SRR11978753' is valid 2020-07-14T01:37:57 prefetch.2.10.7: 1) 'SRR11978753' was downloaded successfully 2020-07-14T01:37:57 prefetch.2.10.7: 'SRR11978753' has 0 unresolved dependencies Read 7895143 spots for SRR11978753/SRR11978753.sra Written 7895143 spots for SRR11978753/SRR11978753.sra 2020-07-14T01:38:34 prefetch.2.10.7: 1) Downloading 'SRR11978754'... 2020-07-14T01:38:34 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:41:02 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:41:03 prefetch.2.10.7: 'SRR11978754' is valid 2020-07-14T01:41:03 prefetch.2.10.7: 1) 'SRR11978754' was downloaded successfully 2020-07-14T01:41:03 prefetch.2.10.7: 'SRR11978754' has 0 unresolved dependencies Read 7946031 spots for SRR11978754/SRR11978754.sra Written 7946031 spots for SRR11978754/SRR11978754.sra 2020-07-14T01:41:40 prefetch.2.10.7: 1) Downloading 'SRR11978755'... 2020-07-14T01:41:40 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:43:16 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:43:17 prefetch.2.10.7: 'SRR11978755' is valid 2020-07-14T01:43:17 prefetch.2.10.7: 1) 'SRR11978755' was downloaded successfully 2020-07-14T01:43:17 prefetch.2.10.7: 'SRR11978755' has 0 unresolved dependencies Read 7975508 spots for SRR11978755/SRR11978755.sra Written 7975508 spots for SRR11978755/SRR11978755.sra 2020-07-14T01:43:53 prefetch.2.10.7: 1) Downloading 'SRR11978756'... 2020-07-14T01:43:53 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:44:50 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:44:50 prefetch.2.10.7: 'SRR11978756' is valid 2020-07-14T01:44:50 prefetch.2.10.7: 1) 'SRR11978756' was downloaded successfully 2020-07-14T01:44:50 prefetch.2.10.7: 'SRR11978756' has 0 unresolved dependencies Read 7845684 spots for SRR11978756/SRR11978756.sra Written 7845684 spots for SRR11978756/SRR11978756.sra 2020-07-14T01:45:27 prefetch.2.10.7: 1) Downloading 'SRR11978757'... 2020-07-14T01:45:27 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:45:39 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:45:39 prefetch.2.10.7: 'SRR11978757' is valid 2020-07-14T01:45:39 prefetch.2.10.7: 1) 'SRR11978757' was downloaded successfully 2020-07-14T01:45:39 prefetch.2.10.7: 'SRR11978757' has 0 unresolved dependencies Read 167993 spots for SRR11978757/SRR11978757.sra Written 167993 spots for SRR11978757/SRR11978757.sra 2020-07-14T01:45:50 prefetch.2.10.7: 1) Downloading 'SRR11978758'... 2020-07-14T01:45:50 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:46:07 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:46:07 prefetch.2.10.7: 'SRR11978758' is valid 2020-07-14T01:46:07 prefetch.2.10.7: 1) 'SRR11978758' was downloaded successfully 2020-07-14T01:46:07 prefetch.2.10.7: 'SRR11978758' has 0 unresolved dependencies Read 179130 spots for SRR11978758/SRR11978758.sra Written 179130 spots for SRR11978758/SRR11978758.sra 2020-07-14T01:46:20 prefetch.2.10.7: 1) Downloading 'SRR11978759'... 2020-07-14T01:46:20 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:46:34 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:46:34 prefetch.2.10.7: 'SRR11978759' is valid 2020-07-14T01:46:34 prefetch.2.10.7: 1) 'SRR11978759' was downloaded successfully 2020-07-14T01:46:34 prefetch.2.10.7: 'SRR11978759' has 0 unresolved dependencies Read 175885 spots for SRR11978759/SRR11978759.sra Written 175885 spots for SRR11978759/SRR11978759.sra 2020-07-14T01:46:49 prefetch.2.10.7: 1) Downloading 'SRR11978760'... 2020-07-14T01:46:49 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:47:36 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:47:36 prefetch.2.10.7: 'SRR11978760' is valid 2020-07-14T01:47:36 prefetch.2.10.7: 1) 'SRR11978760' was downloaded successfully 2020-07-14T01:47:36 prefetch.2.10.7: 'SRR11978760' has 0 unresolved dependencies Read 6590397 spots for SRR11978760/SRR11978760.sra Written 6590397 spots for SRR11978760/SRR11978760.sra 2020-07-14T01:48:12 prefetch.2.10.7: 1) Downloading 'SRR11978761'... 2020-07-14T01:48:12 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:49:03 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:49:03 prefetch.2.10.7: 'SRR11978761' is valid 2020-07-14T01:49:03 prefetch.2.10.7: 1) 'SRR11978761' was downloaded successfully 2020-07-14T01:49:03 prefetch.2.10.7: 'SRR11978761' has 0 unresolved dependencies Read 6427417 spots for SRR11978761/SRR11978761.sra Written 6427417 spots for SRR11978761/SRR11978761.sra 2020-07-14T01:49:35 prefetch.2.10.7: 1) Downloading 'SRR11978762'... 2020-07-14T01:49:35 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:50:15 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:50:15 prefetch.2.10.7: 'SRR11978762' is valid 2020-07-14T01:50:15 prefetch.2.10.7: 1) 'SRR11978762' was downloaded successfully 2020-07-14T01:50:15 prefetch.2.10.7: 'SRR11978762' has 0 unresolved dependencies Read 6619091 spots for SRR11978762/SRR11978762.sra Written 6619091 spots for SRR11978762/SRR11978762.sra 2020-07-14T01:50:47 prefetch.2.10.7: 1) Downloading 'SRR11978763'... 2020-07-14T01:50:47 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:51:59 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:51:59 prefetch.2.10.7: 'SRR11978763' is valid 2020-07-14T01:51:59 prefetch.2.10.7: 1) 'SRR11978763' was downloaded successfully 2020-07-14T01:51:59 prefetch.2.10.7: 'SRR11978763' has 0 unresolved dependencies Read 6440667 spots for SRR11978763/SRR11978763.sra Written 6440667 spots for SRR11978763/SRR11978763.sra 2020-07-14T01:52:31 prefetch.2.10.7: 1) Downloading 'SRR11978764'... 2020-07-14T01:52:31 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T01:53:04 prefetch.2.10.7: HTTPS download succeed 2020-07-14T01:53:04 prefetch.2.10.7: 'SRR11978764' is valid 2020-07-14T01:53:04 prefetch.2.10.7: 1) 'SRR11978764' was downloaded successfully 2020-07-14T01:53:04 prefetch.2.10.7: 'SRR11978764' has 0 unresolved dependencies Read 3264541 spots for SRR11978764/SRR11978764.sra Written 3264541 spots for SRR11978764/SRR11978764.sra fastq に変換しました。 bowtie でマッピング中... Your job 6627509 ("srTdm6") 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:10:34 71373710 reads; of these: 71373710 (100.00%) were unpaired; of these: 59773410 (83.75%) aligned 0 times 9057528 (12.69%) aligned exactly 1 time 2542772 (3.56%) aligned >1 times 16.25% overall alignment rate Time searching: 00:10:34 Overall time: 00:10:34 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 8 files... [bam_rmdupse_core] 3017295 / 11600300 = 0.2601 in library ' ' BAM に変換しました。 Bed ファイルを作成中... WARNING: Skipping mount /opt/pkg/singularity/3.5.2/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Tue, 14 Jul 2020 11:07:43: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam'] # control file = None # effective genome size = 1.20e+08 # 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 @ Tue, 14 Jul 2020 11:07:43: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:07:43: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:07:49: 1000000 INFO @ Tue, 14 Jul 2020 11:07:56: 2000000 INFO @ Tue, 14 Jul 2020 11:08:03: 3000000 INFO @ Tue, 14 Jul 2020 11:08:10: 4000000 WARNING: Skipping mount /opt/pkg/singularity/3.5.2/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Tue, 14 Jul 2020 11:08:13: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam'] # control file = None # effective genome size = 1.20e+08 # 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 @ Tue, 14 Jul 2020 11:08:13: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:08:13: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:08:18: 5000000 INFO @ Tue, 14 Jul 2020 11:08:21: 1000000 INFO @ Tue, 14 Jul 2020 11:08:26: 6000000 INFO @ Tue, 14 Jul 2020 11:08:30: 2000000 INFO @ Tue, 14 Jul 2020 11:08:34: 7000000 INFO @ Tue, 14 Jul 2020 11:08:38: 3000000 BedGraph に変換中... WARNING: Skipping mount /opt/pkg/singularity/3.5.2/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Tue, 14 Jul 2020 11:08:42: 8000000 INFO @ Tue, 14 Jul 2020 11:08:43: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.bam'] # control file = None # effective genome size = 1.20e+08 # 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 @ Tue, 14 Jul 2020 11:08:43: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:08:43: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:08:45: 4000000 INFO @ Tue, 14 Jul 2020 11:08:47: #1 tag size is determined as 62 bps INFO @ Tue, 14 Jul 2020 11:08:47: #1 tag size = 62 INFO @ Tue, 14 Jul 2020 11:08:47: #1 total tags in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:08:47: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:08:47: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 14 Jul 2020 11:08:47: #1 tags after filtering in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:08:47: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:08:47: #1 finished! INFO @ Tue, 14 Jul 2020 11:08:47: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:08:47: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:08:48: #2 number of paired peaks: 288 WARNING @ Tue, 14 Jul 2020 11:08:48: Fewer paired peaks (288) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 288 pairs to build model! INFO @ Tue, 14 Jul 2020 11:08:48: start model_add_line... INFO @ Tue, 14 Jul 2020 11:08:48: start X-correlation... INFO @ Tue, 14 Jul 2020 11:08:48: end of X-cor INFO @ Tue, 14 Jul 2020 11:08:48: #2 finished! INFO @ Tue, 14 Jul 2020 11:08:48: #2 predicted fragment length is 48 bps INFO @ Tue, 14 Jul 2020 11:08:48: #2 alternative fragment length(s) may be 48 bps INFO @ Tue, 14 Jul 2020 11:08:48: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05_model.r WARNING @ Tue, 14 Jul 2020 11:08:48: #2 Since the d (48) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:08:48: #2 You may need to consider one of the other alternative d(s): 48 WARNING @ Tue, 14 Jul 2020 11:08:48: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Tue, 14 Jul 2020 11:08:48: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:08:48: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 14 Jul 2020 11:08:51: 1000000 INFO @ Tue, 14 Jul 2020 11:08:53: 5000000 INFO @ Tue, 14 Jul 2020 11:08:59: 2000000 INFO @ Tue, 14 Jul 2020 11:09:00: 6000000 INFO @ Tue, 14 Jul 2020 11:09:05: #3 Call peaks for each chromosome... INFO @ Tue, 14 Jul 2020 11:09:07: 3000000 INFO @ Tue, 14 Jul 2020 11:09:07: 7000000 INFO @ Tue, 14 Jul 2020 11:09:13: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05_peaks.xls INFO @ Tue, 14 Jul 2020 11:09:13: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:09:13: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.05_summits.bed INFO @ Tue, 14 Jul 2020 11:09:13: Done! pass1 - making usageList (15 chroms): 1 millis pass2 - checking and writing primary data (2145 records, 4 fields): 4 millis CompletedMACS2peakCalling INFO @ Tue, 14 Jul 2020 11:09:14: 4000000 INFO @ Tue, 14 Jul 2020 11:09:15: 8000000 INFO @ Tue, 14 Jul 2020 11:09:19: #1 tag size is determined as 62 bps INFO @ Tue, 14 Jul 2020 11:09:19: #1 tag size = 62 INFO @ Tue, 14 Jul 2020 11:09:19: #1 total tags in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:09:19: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:09:19: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 14 Jul 2020 11:09:19: #1 tags after filtering in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:09:19: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:09:19: #1 finished! INFO @ Tue, 14 Jul 2020 11:09:19: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:09:19: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:09:20: #2 number of paired peaks: 288 WARNING @ Tue, 14 Jul 2020 11:09:20: Fewer paired peaks (288) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 288 pairs to build model! INFO @ Tue, 14 Jul 2020 11:09:20: start model_add_line... INFO @ Tue, 14 Jul 2020 11:09:20: start X-correlation... INFO @ Tue, 14 Jul 2020 11:09:20: end of X-cor INFO @ Tue, 14 Jul 2020 11:09:20: #2 finished! INFO @ Tue, 14 Jul 2020 11:09:20: #2 predicted fragment length is 48 bps INFO @ Tue, 14 Jul 2020 11:09:20: #2 alternative fragment length(s) may be 48 bps INFO @ Tue, 14 Jul 2020 11:09:20: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10_model.r WARNING @ Tue, 14 Jul 2020 11:09:20: #2 Since the d (48) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:09:20: #2 You may need to consider one of the other alternative d(s): 48 WARNING @ Tue, 14 Jul 2020 11:09:20: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Tue, 14 Jul 2020 11:09:20: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:09:20: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 14 Jul 2020 11:09:21: 5000000 INFO @ Tue, 14 Jul 2020 11:09:27: 6000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Tue, 14 Jul 2020 11:09:33: 7000000 INFO @ Tue, 14 Jul 2020 11:09:36: #3 Call peaks for each chromosome... INFO @ Tue, 14 Jul 2020 11:09:39: 8000000 INFO @ Tue, 14 Jul 2020 11:09:43: #1 tag size is determined as 62 bps INFO @ Tue, 14 Jul 2020 11:09:43: #1 tag size = 62 INFO @ Tue, 14 Jul 2020 11:09:43: #1 total tags in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:09:43: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:09:43: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 14 Jul 2020 11:09:43: #1 tags after filtering in treatment: 8583005 INFO @ Tue, 14 Jul 2020 11:09:43: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:09:43: #1 finished! INFO @ Tue, 14 Jul 2020 11:09:43: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:09:43: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:09:43: #2 number of paired peaks: 288 WARNING @ Tue, 14 Jul 2020 11:09:43: Fewer paired peaks (288) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 288 pairs to build model! INFO @ Tue, 14 Jul 2020 11:09:43: start model_add_line... INFO @ Tue, 14 Jul 2020 11:09:43: start X-correlation... INFO @ Tue, 14 Jul 2020 11:09:43: end of X-cor INFO @ Tue, 14 Jul 2020 11:09:43: #2 finished! INFO @ Tue, 14 Jul 2020 11:09:43: #2 predicted fragment length is 48 bps INFO @ Tue, 14 Jul 2020 11:09:43: #2 alternative fragment length(s) may be 48 bps INFO @ Tue, 14 Jul 2020 11:09:43: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20_model.r WARNING @ Tue, 14 Jul 2020 11:09:44: #2 Since the d (48) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:09:44: #2 You may need to consider one of the other alternative d(s): 48 WARNING @ Tue, 14 Jul 2020 11:09:44: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Tue, 14 Jul 2020 11:09:44: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:09:44: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 14 Jul 2020 11:09:45: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10_peaks.xls INFO @ Tue, 14 Jul 2020 11:09:45: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:09:45: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.10_summits.bed INFO @ Tue, 14 Jul 2020 11:09:45: Done! pass1 - making usageList (14 chroms): 1 millis pass2 - checking and writing primary data (1002 records, 4 fields): 2 millis CompletedMACS2peakCalling BigWig に変換しました。 INFO @ Tue, 14 Jul 2020 11:10:00: #3 Call peaks for each chromosome... INFO @ Tue, 14 Jul 2020 11:10:09: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20_peaks.xls INFO @ Tue, 14 Jul 2020 11:10:09: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:10:09: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521316/SRX8521316.20_summits.bed INFO @ Tue, 14 Jul 2020 11:10:09: Done! pass1 - making usageList (7 chroms): 0 millis pass2 - checking and writing primary data (157 records, 4 fields): 1 millis CompletedMACS2peakCalling