Job ID = 1293885 sra ファイルのダウンロード中... Read layout: SINGLE fastq に変換中... spots read : 124,117,236 reads read : 124,117,236 reads written : 124,117,236 rm: cannot remove ‘[DSE]RR*’: No such file or directory rm: cannot remove ‘fastqDump_tmp*’: No such file or directory fastq に変換しました。 bowtie でマッピング中... 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:53:59 124117236 reads; of these: 124117236 (100.00%) were unpaired; of these: 9221849 (7.43%) aligned 0 times 75344934 (60.70%) aligned exactly 1 time 39550453 (31.87%) aligned >1 times 92.57% overall alignment rate Time searching: 00:53:59 Overall time: 00:53:59 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 48 files... [bam_rmdupse_core] 62292122 / 114895387 = 0.5422 in library ' ' BAM に変換しました。 Bed ファイルを作成中... BedGraph に変換中... INFO @ Mon, 03 Jun 2019 04:28:55: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.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 @ Mon, 03 Jun 2019 04:28:55: #1 read tag files... INFO @ Mon, 03 Jun 2019 04:28:55: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 04:28:55: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.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 @ Mon, 03 Jun 2019 04:28:55: #1 read tag files... INFO @ Mon, 03 Jun 2019 04:28:55: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 04:28:55: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.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 @ Mon, 03 Jun 2019 04:28:55: #1 read tag files... INFO @ Mon, 03 Jun 2019 04:28:55: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 04:29:03: 1000000 INFO @ Mon, 03 Jun 2019 04:29:04: 1000000 INFO @ Mon, 03 Jun 2019 04:29:04: 1000000 INFO @ Mon, 03 Jun 2019 04:29:11: 2000000 INFO @ Mon, 03 Jun 2019 04:29:12: 2000000 INFO @ Mon, 03 Jun 2019 04:29:13: 2000000 INFO @ Mon, 03 Jun 2019 04:29:19: 3000000 INFO @ Mon, 03 Jun 2019 04:29:21: 3000000 INFO @ Mon, 03 Jun 2019 04:29:22: 3000000 INFO @ Mon, 03 Jun 2019 04:29:26: 4000000 INFO @ Mon, 03 Jun 2019 04:29:29: 4000000 INFO @ Mon, 03 Jun 2019 04:29:31: 4000000 INFO @ Mon, 03 Jun 2019 04:29:34: 5000000 INFO @ Mon, 03 Jun 2019 04:29:37: 5000000 INFO @ Mon, 03 Jun 2019 04:29:40: 5000000 INFO @ Mon, 03 Jun 2019 04:29:41: 6000000 INFO @ Mon, 03 Jun 2019 04:29:45: 6000000 INFO @ Mon, 03 Jun 2019 04:29:48: 6000000 INFO @ Mon, 03 Jun 2019 04:29:49: 7000000 INFO @ Mon, 03 Jun 2019 04:29:53: 7000000 INFO @ Mon, 03 Jun 2019 04:29:56: 8000000 INFO @ Mon, 03 Jun 2019 04:29:57: 7000000 INFO @ Mon, 03 Jun 2019 04:30:01: 8000000 INFO @ Mon, 03 Jun 2019 04:30:04: 9000000 INFO @ Mon, 03 Jun 2019 04:30:05: 8000000 INFO @ Mon, 03 Jun 2019 04:30:08: 9000000 INFO @ Mon, 03 Jun 2019 04:30:11: 10000000 INFO @ Mon, 03 Jun 2019 04:30:13: 9000000 INFO @ Mon, 03 Jun 2019 04:30:16: 10000000 INFO @ Mon, 03 Jun 2019 04:30:19: 11000000 INFO @ Mon, 03 Jun 2019 04:30:21: 10000000 INFO @ Mon, 03 Jun 2019 04:30:24: 11000000 INFO @ Mon, 03 Jun 2019 04:30:26: 12000000 INFO @ Mon, 03 Jun 2019 04:30:29: 11000000 INFO @ Mon, 03 Jun 2019 04:30:32: 12000000 INFO @ Mon, 03 Jun 2019 04:30:34: 13000000 INFO @ Mon, 03 Jun 2019 04:30:37: 12000000 INFO @ Mon, 03 Jun 2019 04:30:40: 13000000 INFO @ Mon, 03 Jun 2019 04:30:41: 14000000 INFO @ Mon, 03 Jun 2019 04:30:45: 13000000 INFO @ Mon, 03 Jun 2019 04:30:48: 14000000 INFO @ Mon, 03 Jun 2019 04:30:49: 15000000 INFO @ Mon, 03 Jun 2019 04:30:53: 14000000 INFO @ Mon, 03 Jun 2019 04:30:55: 15000000 INFO @ Mon, 03 Jun 2019 04:30:56: 16000000 INFO @ Mon, 03 Jun 2019 04:31:01: 15000000 INFO @ Mon, 03 Jun 2019 04:31:03: 16000000 INFO @ Mon, 03 Jun 2019 04:31:04: 17000000 INFO @ Mon, 03 Jun 2019 04:31:10: 16000000 INFO @ Mon, 03 Jun 2019 04:31:11: 17000000 INFO @ Mon, 03 Jun 2019 04:31:11: 18000000 INFO @ Mon, 03 Jun 2019 04:31:18: 17000000 INFO @ Mon, 03 Jun 2019 04:31:19: 19000000 INFO @ Mon, 03 Jun 2019 04:31:19: 18000000 INFO @ Mon, 03 Jun 2019 04:31:26: 18000000 INFO @ Mon, 03 Jun 2019 04:31:26: 20000000 INFO @ Mon, 03 Jun 2019 04:31:27: 19000000 INFO @ Mon, 03 Jun 2019 04:31:34: 21000000 INFO @ Mon, 03 Jun 2019 04:31:34: 19000000 INFO @ Mon, 03 Jun 2019 04:31:34: 20000000 INFO @ Mon, 03 Jun 2019 04:31:41: 22000000 INFO @ Mon, 03 Jun 2019 04:31:42: 21000000 INFO @ Mon, 03 Jun 2019 04:31:42: 20000000 INFO @ Mon, 03 Jun 2019 04:31:49: 23000000 INFO @ Mon, 03 Jun 2019 04:31:50: 22000000 INFO @ Mon, 03 Jun 2019 04:31:50: 21000000 INFO @ Mon, 03 Jun 2019 04:31:56: 24000000 INFO @ Mon, 03 Jun 2019 04:31:58: 23000000 INFO @ Mon, 03 Jun 2019 04:31:58: 22000000 INFO @ Mon, 03 Jun 2019 04:32:04: 25000000 INFO @ Mon, 03 Jun 2019 04:32:06: 24000000 INFO @ Mon, 03 Jun 2019 04:32:06: 23000000 INFO @ Mon, 03 Jun 2019 04:32:11: 26000000 INFO @ Mon, 03 Jun 2019 04:32:13: 25000000 INFO @ Mon, 03 Jun 2019 04:32:14: 24000000 INFO @ Mon, 03 Jun 2019 04:32:19: 27000000 INFO @ Mon, 03 Jun 2019 04:32:21: 26000000 INFO @ Mon, 03 Jun 2019 04:32:22: 25000000 INFO @ Mon, 03 Jun 2019 04:32:26: 28000000 INFO @ Mon, 03 Jun 2019 04:32:29: 27000000 INFO @ Mon, 03 Jun 2019 04:32:30: 26000000 INFO @ Mon, 03 Jun 2019 04:32:34: 29000000 INFO @ Mon, 03 Jun 2019 04:32:37: 28000000 INFO @ Mon, 03 Jun 2019 04:32:38: 27000000 INFO @ Mon, 03 Jun 2019 04:32:42: 30000000 INFO @ Mon, 03 Jun 2019 04:32:45: 29000000 INFO @ Mon, 03 Jun 2019 04:32:47: 28000000 INFO @ Mon, 03 Jun 2019 04:32:49: 31000000 INFO @ Mon, 03 Jun 2019 04:32:53: 30000000 INFO @ Mon, 03 Jun 2019 04:32:55: 29000000 INFO @ Mon, 03 Jun 2019 04:32:56: 32000000 INFO @ Mon, 03 Jun 2019 04:33:00: 31000000 INFO @ Mon, 03 Jun 2019 04:33:03: 30000000 INFO @ Mon, 03 Jun 2019 04:33:04: 33000000 INFO @ Mon, 03 Jun 2019 04:33:08: 32000000 INFO @ Mon, 03 Jun 2019 04:33:11: 31000000 INFO @ Mon, 03 Jun 2019 04:33:11: 34000000 INFO @ Mon, 03 Jun 2019 04:33:16: 33000000 INFO @ Mon, 03 Jun 2019 04:33:19: 35000000 INFO @ Mon, 03 Jun 2019 04:33:19: 32000000 INFO @ Mon, 03 Jun 2019 04:33:24: 34000000 INFO @ Mon, 03 Jun 2019 04:33:26: 36000000 INFO @ Mon, 03 Jun 2019 04:33:29: 33000000 INFO @ Mon, 03 Jun 2019 04:33:31: 35000000 INFO @ Mon, 03 Jun 2019 04:33:34: 37000000 INFO @ Mon, 03 Jun 2019 04:33:37: 34000000 INFO @ Mon, 03 Jun 2019 04:33:39: 36000000 INFO @ Mon, 03 Jun 2019 04:33:42: 38000000 INFO @ Mon, 03 Jun 2019 04:33:45: 35000000 INFO @ Mon, 03 Jun 2019 04:33:47: 37000000 INFO @ Mon, 03 Jun 2019 04:33:49: 39000000 INFO @ Mon, 03 Jun 2019 04:33:53: 36000000 INFO @ Mon, 03 Jun 2019 04:33:55: 38000000 INFO @ Mon, 03 Jun 2019 04:33:57: 40000000 INFO @ Mon, 03 Jun 2019 04:34:01: 37000000 INFO @ Mon, 03 Jun 2019 04:34:03: 39000000 INFO @ Mon, 03 Jun 2019 04:34:04: 41000000 INFO @ Mon, 03 Jun 2019 04:34:09: 38000000 INFO @ Mon, 03 Jun 2019 04:34:11: 40000000 INFO @ Mon, 03 Jun 2019 04:34:12: 42000000 INFO @ Mon, 03 Jun 2019 04:34:17: 39000000 INFO @ Mon, 03 Jun 2019 04:34:19: 41000000 INFO @ Mon, 03 Jun 2019 04:34:20: 43000000 INFO @ Mon, 03 Jun 2019 04:34:25: 40000000 INFO @ Mon, 03 Jun 2019 04:34:27: 42000000 INFO @ Mon, 03 Jun 2019 04:34:27: 44000000 INFO @ Mon, 03 Jun 2019 04:34:33: 41000000 INFO @ Mon, 03 Jun 2019 04:34:34: 43000000 INFO @ Mon, 03 Jun 2019 04:34:35: 45000000 INFO @ Mon, 03 Jun 2019 04:34:41: 42000000 INFO @ Mon, 03 Jun 2019 04:34:42: 44000000 INFO @ Mon, 03 Jun 2019 04:34:43: 46000000 INFO @ Mon, 03 Jun 2019 04:34:49: 43000000 INFO @ Mon, 03 Jun 2019 04:34:50: 47000000 INFO @ Mon, 03 Jun 2019 04:34:51: 45000000 INFO @ Mon, 03 Jun 2019 04:34:57: 44000000 INFO @ Mon, 03 Jun 2019 04:34:58: 48000000 INFO @ Mon, 03 Jun 2019 04:34:58: 46000000 INFO @ Mon, 03 Jun 2019 04:35:05: 45000000 INFO @ Mon, 03 Jun 2019 04:35:05: 49000000 INFO @ Mon, 03 Jun 2019 04:35:06: 47000000 INFO @ Mon, 03 Jun 2019 04:35:13: 46000000 INFO @ Mon, 03 Jun 2019 04:35:13: 50000000 INFO @ Mon, 03 Jun 2019 04:35:14: 48000000 INFO @ Mon, 03 Jun 2019 04:35:20: 51000000 INFO @ Mon, 03 Jun 2019 04:35:20: 47000000 INFO @ Mon, 03 Jun 2019 04:35:22: 49000000 INFO @ Mon, 03 Jun 2019 04:35:28: 52000000 INFO @ Mon, 03 Jun 2019 04:35:28: 48000000 INFO @ Mon, 03 Jun 2019 04:35:29: 50000000 INFO @ Mon, 03 Jun 2019 04:35:33: #1 tag size is determined as 50 bps INFO @ Mon, 03 Jun 2019 04:35:33: #1 tag size = 50 INFO @ Mon, 03 Jun 2019 04:35:33: #1 total tags in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:35:33: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 04:35:33: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 04:35:34: #1 tags after filtering in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:35:34: #1 Redundant rate of treatment: 0.00 INFO @ Mon, 03 Jun 2019 04:35:34: #1 finished! INFO @ Mon, 03 Jun 2019 04:35:34: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 04:35:34: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 04:35:36: 49000000 INFO @ Mon, 03 Jun 2019 04:35:37: 51000000 INFO @ Mon, 03 Jun 2019 04:35:38: #2 number of paired peaks: 0 WARNING @ Mon, 03 Jun 2019 04:35:38: Too few paired peaks (0) 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 @ Mon, 03 Jun 2019 04:35:38: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.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/SRX135527/SRX135527.20_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.20_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.20_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling INFO @ Mon, 03 Jun 2019 04:35:43: 50000000 INFO @ Mon, 03 Jun 2019 04:35:45: 52000000 INFO @ Mon, 03 Jun 2019 04:35:50: #1 tag size is determined as 50 bps INFO @ Mon, 03 Jun 2019 04:35:50: #1 tag size = 50 INFO @ Mon, 03 Jun 2019 04:35:50: #1 total tags in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:35:50: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 04:35:50: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 04:35:51: 51000000 INFO @ Mon, 03 Jun 2019 04:35:51: #1 tags after filtering in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:35:51: #1 Redundant rate of treatment: 0.00 INFO @ Mon, 03 Jun 2019 04:35:51: #1 finished! INFO @ Mon, 03 Jun 2019 04:35:51: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 04:35:51: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 04:35:56: #2 number of paired peaks: 0 WARNING @ Mon, 03 Jun 2019 04:35:56: Too few paired peaks (0) 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 @ Mon, 03 Jun 2019 04:35:56: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.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/SRX135527/SRX135527.05_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.05_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.05_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling INFO @ Mon, 03 Jun 2019 04:35:59: 52000000 INFO @ Mon, 03 Jun 2019 04:36:04: #1 tag size is determined as 50 bps INFO @ Mon, 03 Jun 2019 04:36:04: #1 tag size = 50 INFO @ Mon, 03 Jun 2019 04:36:04: #1 total tags in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:36:04: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 04:36:04: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 04:36:05: #1 tags after filtering in treatment: 52603265 INFO @ Mon, 03 Jun 2019 04:36:05: #1 Redundant rate of treatment: 0.00 INFO @ Mon, 03 Jun 2019 04:36:05: #1 finished! INFO @ Mon, 03 Jun 2019 04:36:05: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 04:36:05: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 04:36:09: #2 number of paired peaks: 0 WARNING @ Mon, 03 Jun 2019 04:36:09: Too few paired peaks (0) 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 @ Mon, 03 Jun 2019 04:36:09: Process for pairing-model is terminated! cut: /home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10_peaks.narrowPeak: No such file or directory pass1 - making usageList (0 chroms): 2 millis needLargeMem: trying to allocate 0 bytes (limit: 17179869184) rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10_model.r’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10_*.xls’: No such file or directory rm: cannot remove ‘/home/okishinya/chipatlas/results/dm3/SRX135527/SRX135527.10_peaks.narrowPeak’: No such file or directory CompletedMACS2peakCalling BedGraph に変換しました。 BigWig に変換中... BigWig に変換しました。