Job ID = 6528007 SRX = SRX3723041 Genome = dm3 sra ファイルのダウンロード中... Read layout: SINGLE fastq に変換中... 2020-06-29T14:38:22 prefetch.2.10.7: 1) Downloading 'SRR6750408'... 2020-06-29T14:38:22 prefetch.2.10.7: Downloading via HTTPS... 2020-06-29T14:40:40 prefetch.2.10.7: HTTPS download succeed 2020-06-29T14:40:41 prefetch.2.10.7: 'SRR6750408' is valid 2020-06-29T14:40:41 prefetch.2.10.7: 1) 'SRR6750408' was downloaded successfully 2020-06-29T14:40:41 prefetch.2.10.7: 'SRR6750408' has 0 unresolved dependencies Read 12603035 spots for SRR6750408/SRR6750408.sra Written 12603035 spots for SRR6750408/SRR6750408.sra 2020-06-29T14:41:49 prefetch.2.10.7: 1) Downloading 'SRR6750409'... 2020-06-29T14:41:49 prefetch.2.10.7: Downloading via HTTPS... 2020-06-29T14:43:22 prefetch.2.10.7: HTTPS download succeed 2020-06-29T14:43:23 prefetch.2.10.7: 'SRR6750409' is valid 2020-06-29T14:43:23 prefetch.2.10.7: 1) 'SRR6750409' was downloaded successfully 2020-06-29T14:43:23 prefetch.2.10.7: 'SRR6750409' has 0 unresolved dependencies Read 12471928 spots for SRR6750409/SRR6750409.sra Written 12471928 spots for SRR6750409/SRR6750409.sra 2020-06-29T14:44:31 prefetch.2.10.7: 1) Downloading 'SRR6750410'... 2020-06-29T14:44:31 prefetch.2.10.7: Downloading via HTTPS... 2020-06-29T14:46:30 prefetch.2.10.7: HTTPS download succeed 2020-06-29T14:46:31 prefetch.2.10.7: 'SRR6750410' is valid 2020-06-29T14:46:31 prefetch.2.10.7: 1) 'SRR6750410' was downloaded successfully 2020-06-29T14:46:31 prefetch.2.10.7: 'SRR6750410' has 0 unresolved dependencies Read 13154751 spots for SRR6750410/SRR6750410.sra Written 13154751 spots for SRR6750410/SRR6750410.sra 2020-06-29T14:47:41 prefetch.2.10.7: 1) Downloading 'SRR6750411'... 2020-06-29T14:47:41 prefetch.2.10.7: Downloading via HTTPS... 2020-06-29T14:48:59 prefetch.2.10.7: HTTPS download succeed 2020-06-29T14:49:00 prefetch.2.10.7: 'SRR6750411' is valid 2020-06-29T14:49:00 prefetch.2.10.7: 1) 'SRR6750411' was downloaded successfully 2020-06-29T14:49:00 prefetch.2.10.7: 'SRR6750411' has 0 unresolved dependencies Read 13015888 spots for SRR6750411/SRR6750411.sra Written 13015888 spots for SRR6750411/SRR6750411.sra fastq に変換しました。 bowtie でマッピング中... Time loading reference: 00:00:00 Time loading forward index: 00:00:00 Time loading mirror index: 00:00:01 Multiseed full-index search: 00:22:21 51245602 reads; of these: 51245602 (100.00%) were unpaired; of these: 4657636 (9.09%) aligned 0 times 34706850 (67.73%) aligned exactly 1 time 11881116 (23.18%) aligned >1 times 90.91% overall alignment rate Time searching: 00:22:22 Overall time: 00:22:22 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 20 files... [bam_rmdupse_core] 12537301 / 46587966 = 0.2691 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, 30 Jun 2020 00:55:27: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.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, 30 Jun 2020 00:55:27: #1 read tag files... INFO @ Tue, 30 Jun 2020 00:55:27: #1 read treatment tags... INFO @ Tue, 30 Jun 2020 00:55:33: 1000000 INFO @ Tue, 30 Jun 2020 00:55:40: 2000000 INFO @ Tue, 30 Jun 2020 00:55:47: 3000000 INFO @ Tue, 30 Jun 2020 00:55:54: 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, 30 Jun 2020 00:55:57: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.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, 30 Jun 2020 00:55:57: #1 read tag files... INFO @ Tue, 30 Jun 2020 00:55:57: #1 read treatment tags... INFO @ Tue, 30 Jun 2020 00:56:01: 5000000 INFO @ Tue, 30 Jun 2020 00:56:04: 1000000 INFO @ Tue, 30 Jun 2020 00:56:08: 6000000 INFO @ Tue, 30 Jun 2020 00:56:10: 2000000 INFO @ Tue, 30 Jun 2020 00:56:16: 7000000 INFO @ Tue, 30 Jun 2020 00:56:17: 3000000 INFO @ Tue, 30 Jun 2020 00:56:23: 8000000 INFO @ Tue, 30 Jun 2020 00:56:24: 4000000 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, 30 Jun 2020 00:56:27: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.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, 30 Jun 2020 00:56:27: #1 read tag files... INFO @ Tue, 30 Jun 2020 00:56:27: #1 read treatment tags... INFO @ Tue, 30 Jun 2020 00:56:31: 9000000 INFO @ Tue, 30 Jun 2020 00:56:31: 5000000 INFO @ Tue, 30 Jun 2020 00:56:35: 1000000 INFO @ Tue, 30 Jun 2020 00:56:38: 6000000 INFO @ Tue, 30 Jun 2020 00:56:38: 10000000 INFO @ Tue, 30 Jun 2020 00:56:43: 2000000 INFO @ Tue, 30 Jun 2020 00:56:45: 7000000 INFO @ Tue, 30 Jun 2020 00:56:46: 11000000 INFO @ Tue, 30 Jun 2020 00:56:51: 3000000 INFO @ Tue, 30 Jun 2020 00:56:53: 8000000 INFO @ Tue, 30 Jun 2020 00:56:53: 12000000 INFO @ Tue, 30 Jun 2020 00:56:58: 4000000 INFO @ Tue, 30 Jun 2020 00:57:00: 9000000 INFO @ Tue, 30 Jun 2020 00:57:01: 13000000 INFO @ Tue, 30 Jun 2020 00:57:06: 5000000 INFO @ Tue, 30 Jun 2020 00:57:08: 10000000 INFO @ Tue, 30 Jun 2020 00:57:08: 14000000 INFO @ Tue, 30 Jun 2020 00:57:14: 6000000 INFO @ Tue, 30 Jun 2020 00:57:16: 15000000 INFO @ Tue, 30 Jun 2020 00:57:16: 11000000 INFO @ Tue, 30 Jun 2020 00:57:21: 7000000 INFO @ Tue, 30 Jun 2020 00:57:23: 16000000 INFO @ Tue, 30 Jun 2020 00:57:23: 12000000 INFO @ Tue, 30 Jun 2020 00:57:28: 8000000 INFO @ Tue, 30 Jun 2020 00:57:30: 17000000 INFO @ Tue, 30 Jun 2020 00:57:31: 13000000 INFO @ Tue, 30 Jun 2020 00:57:35: 9000000 INFO @ Tue, 30 Jun 2020 00:57:37: 18000000 INFO @ Tue, 30 Jun 2020 00:57:38: 14000000 INFO @ Tue, 30 Jun 2020 00:57:43: 10000000 INFO @ Tue, 30 Jun 2020 00:57:44: 19000000 INFO @ Tue, 30 Jun 2020 00:57:45: 15000000 INFO @ Tue, 30 Jun 2020 00:57:50: 11000000 INFO @ Tue, 30 Jun 2020 00:57:51: 20000000 INFO @ Tue, 30 Jun 2020 00:57:52: 16000000 INFO @ Tue, 30 Jun 2020 00:57:57: 12000000 INFO @ Tue, 30 Jun 2020 00:57:58: 21000000 INFO @ Tue, 30 Jun 2020 00:57:58: 17000000 INFO @ Tue, 30 Jun 2020 00:58:04: 13000000 INFO @ Tue, 30 Jun 2020 00:58:05: 18000000 INFO @ Tue, 30 Jun 2020 00:58:06: 22000000 INFO @ Tue, 30 Jun 2020 00:58:12: 19000000 INFO @ Tue, 30 Jun 2020 00:58:12: 14000000 INFO @ Tue, 30 Jun 2020 00:58:13: 23000000 INFO @ Tue, 30 Jun 2020 00:58:19: 20000000 INFO @ Tue, 30 Jun 2020 00:58:20: 15000000 INFO @ Tue, 30 Jun 2020 00:58:20: 24000000 INFO @ Tue, 30 Jun 2020 00:58:25: 21000000 INFO @ Tue, 30 Jun 2020 00:58:27: 16000000 INFO @ Tue, 30 Jun 2020 00:58:28: 25000000 INFO @ Tue, 30 Jun 2020 00:58:32: 22000000 INFO @ Tue, 30 Jun 2020 00:58:34: 17000000 INFO @ Tue, 30 Jun 2020 00:58:35: 26000000 INFO @ Tue, 30 Jun 2020 00:58:39: 23000000 INFO @ Tue, 30 Jun 2020 00:58:41: 18000000 INFO @ Tue, 30 Jun 2020 00:58:43: 27000000 INFO @ Tue, 30 Jun 2020 00:58:46: 24000000 INFO @ Tue, 30 Jun 2020 00:58:48: 19000000 INFO @ Tue, 30 Jun 2020 00:58:50: 28000000 INFO @ Tue, 30 Jun 2020 00:58:53: 25000000 INFO @ Tue, 30 Jun 2020 00:58:55: 20000000 INFO @ Tue, 30 Jun 2020 00:58:57: 29000000 INFO @ Tue, 30 Jun 2020 00:58:59: 26000000 INFO @ Tue, 30 Jun 2020 00:59:02: 21000000 INFO @ Tue, 30 Jun 2020 00:59:04: 30000000 INFO @ Tue, 30 Jun 2020 00:59:06: 27000000 INFO @ Tue, 30 Jun 2020 00:59:10: 22000000 INFO @ Tue, 30 Jun 2020 00:59:11: 31000000 INFO @ Tue, 30 Jun 2020 00:59:13: 28000000 INFO @ Tue, 30 Jun 2020 00:59:17: 23000000 INFO @ Tue, 30 Jun 2020 00:59:19: 32000000 INFO @ Tue, 30 Jun 2020 00:59:19: 29000000 INFO @ Tue, 30 Jun 2020 00:59:24: 24000000 INFO @ Tue, 30 Jun 2020 00:59:26: 33000000 INFO @ Tue, 30 Jun 2020 00:59:26: 30000000 INFO @ Tue, 30 Jun 2020 00:59:32: 25000000 INFO @ Tue, 30 Jun 2020 00:59:33: 31000000 INFO @ Tue, 30 Jun 2020 00:59:33: 34000000 INFO @ Tue, 30 Jun 2020 00:59:34: #1 tag size is determined as 75 bps INFO @ Tue, 30 Jun 2020 00:59:34: #1 tag size = 75 INFO @ Tue, 30 Jun 2020 00:59:34: #1 total tags in treatment: 34050665 INFO @ Tue, 30 Jun 2020 00:59:34: #1 user defined the maximum tags... INFO @ Tue, 30 Jun 2020 00:59:34: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 30 Jun 2020 00:59:34: #1 tags after filtering in treatment: 34050665 INFO @ Tue, 30 Jun 2020 00:59:34: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 30 Jun 2020 00:59:34: #1 finished! INFO @ Tue, 30 Jun 2020 00:59:34: #2 Build Peak Model... INFO @ Tue, 30 Jun 2020 00:59:34: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 30 Jun 2020 00:59:37: #2 number of paired peaks: 3 WARNING @ Tue, 30 Jun 2020 00:59:37: Too few paired peaks (3) so I can not build the model! Broader your MFOLD range parameter may erase this error. If it still can't build the model, we suggest to use --nomodel and --extsize 147 or other fixed number instead. WARNING @ Tue, 30 Jun 2020 00:59:37: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05_peaks.narrowPeak: No such file or directory pass1 - making usageList (0 chroms): 1 millis needLargeMem: trying to allocate 0 bytes (limit: 17179869184) rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.05_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling INFO @ Tue, 30 Jun 2020 00:59:39: 26000000 INFO @ Tue, 30 Jun 2020 00:59:42: 32000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Tue, 30 Jun 2020 00:59:46: 27000000 INFO @ Tue, 30 Jun 2020 00:59:48: 33000000 INFO @ Tue, 30 Jun 2020 00:59:53: 28000000 INFO @ Tue, 30 Jun 2020 00:59:55: 34000000 INFO @ Tue, 30 Jun 2020 00:59:56: #1 tag size is determined as 75 bps INFO @ Tue, 30 Jun 2020 00:59:56: #1 tag size = 75 INFO @ Tue, 30 Jun 2020 00:59:56: #1 total tags in treatment: 34050665 INFO @ Tue, 30 Jun 2020 00:59:56: #1 user defined the maximum tags... INFO @ Tue, 30 Jun 2020 00:59:56: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 30 Jun 2020 00:59:56: #1 tags after filtering in treatment: 34050665 INFO @ Tue, 30 Jun 2020 00:59:56: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 30 Jun 2020 00:59:56: #1 finished! INFO @ Tue, 30 Jun 2020 00:59:56: #2 Build Peak Model... INFO @ Tue, 30 Jun 2020 00:59:56: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 30 Jun 2020 00:59:58: #2 number of paired peaks: 3 WARNING @ Tue, 30 Jun 2020 00:59:58: Too few paired peaks (3) so I can not build the model! Broader your MFOLD range parameter may erase this error. If it still can't build the model, we suggest to use --nomodel and --extsize 147 or other fixed number instead. WARNING @ Tue, 30 Jun 2020 00:59:58: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10_peaks.narrowPeak: No such file or directory pass1 - making usageList (0 chroms): 1 millis needLargeMem: trying to allocate 0 bytes (limit: 17179869184) rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.10_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling INFO @ Tue, 30 Jun 2020 00:59:59: 29000000 INFO @ Tue, 30 Jun 2020 01:00:06: 30000000 INFO @ Tue, 30 Jun 2020 01:00:12: 31000000 INFO @ Tue, 30 Jun 2020 01:00:19: 32000000 INFO @ Tue, 30 Jun 2020 01:00:25: 33000000 INFO @ Tue, 30 Jun 2020 01:00:32: 34000000 INFO @ Tue, 30 Jun 2020 01:00:33: #1 tag size is determined as 75 bps INFO @ Tue, 30 Jun 2020 01:00:33: #1 tag size = 75 INFO @ Tue, 30 Jun 2020 01:00:33: #1 total tags in treatment: 34050665 INFO @ Tue, 30 Jun 2020 01:00:33: #1 user defined the maximum tags... INFO @ Tue, 30 Jun 2020 01:00:33: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Tue, 30 Jun 2020 01:00:33: #1 tags after filtering in treatment: 34050665 INFO @ Tue, 30 Jun 2020 01:00:33: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 30 Jun 2020 01:00:33: #1 finished! INFO @ Tue, 30 Jun 2020 01:00:33: #2 Build Peak Model... INFO @ Tue, 30 Jun 2020 01:00:33: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 30 Jun 2020 01:00:36: #2 number of paired peaks: 3 WARNING @ Tue, 30 Jun 2020 01:00:36: Too few paired peaks (3) so I can not build the model! Broader your MFOLD range parameter may erase this error. If it still can't build the model, we suggest to use --nomodel and --extsize 147 or other fixed number instead. WARNING @ Tue, 30 Jun 2020 01:00:36: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20_peaks.narrowPeak: No such file or directory pass1 - making usageList (0 chroms): 1 millis needLargeMem: trying to allocate 0 bytes (limit: 17179869184) rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX3723041/SRX3723041.20_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling BigWig に変換しました。