Job ID = 6627397 SRX = SRX8521379 Genome = dm3 sra ファイルのダウンロード中... Read layout: SINGLE fastq に変換中... 2020-07-14T02:18:44 prefetch.2.10.7: 1) Downloading 'SRR11978234'... 2020-07-14T02:18:44 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:19:00 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:19:00 prefetch.2.10.7: 'SRR11978234' is valid 2020-07-14T02:19:00 prefetch.2.10.7: 1) 'SRR11978234' was downloaded successfully 2020-07-14T02:19:00 prefetch.2.10.7: 'SRR11978234' has 0 unresolved dependencies Read 356305 spots for SRR11978234/SRR11978234.sra Written 356305 spots for SRR11978234/SRR11978234.sra 2020-07-14T02:19:12 prefetch.2.10.7: 1) Downloading 'SRR11978235'... 2020-07-14T02:19:12 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:21:31 prefetch.2.10.7: HTTPS download failed 2020-07-14T02:21:31 prefetch.2.10.7: 1) failed to download SRR11978235 2020-07-14T02:21:45 prefetch.2.10.7: 1) Downloading 'SRR11978235'... 2020-07-14T02:21:45 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:21:45 prefetch.2.10.7: Continue download of 'SRR11978235' from 51886318 2020-07-14T02:21:53 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:21:53 prefetch.2.10.7: 'SRR11978235' is valid 2020-07-14T02:21:53 prefetch.2.10.7: 1) 'SRR11978235' was downloaded successfully 2020-07-14T02:21:53 prefetch.2.10.7: 'SRR11978235' has 0 unresolved dependencies Read 1698219 spots for SRR11978235/SRR11978235.sra Written 1698219 spots for SRR11978235/SRR11978235.sra 2020-07-14T02:22:13 prefetch.2.10.7: 1) Downloading 'SRR11978236'... 2020-07-14T02:22:13 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:22:42 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:22:42 prefetch.2.10.7: 'SRR11978236' is valid 2020-07-14T02:22:42 prefetch.2.10.7: 1) 'SRR11978236' was downloaded successfully 2020-07-14T02:22:42 prefetch.2.10.7: 'SRR11978236' has 0 unresolved dependencies Read 1802063 spots for SRR11978236/SRR11978236.sra Written 1802063 spots for SRR11978236/SRR11978236.sra 2020-07-14T02:22:59 prefetch.2.10.7: 1) Downloading 'SRR11978237'... 2020-07-14T02:22:59 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:23:22 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:23:22 prefetch.2.10.7: 'SRR11978237' is valid 2020-07-14T02:23:22 prefetch.2.10.7: 1) 'SRR11978237' was downloaded successfully 2020-07-14T02:23:22 prefetch.2.10.7: 'SRR11978237' has 0 unresolved dependencies Read 1764510 spots for SRR11978237/SRR11978237.sra Written 1764510 spots for SRR11978237/SRR11978237.sra 2020-07-14T02:23:48 prefetch.2.10.7: 1) Downloading 'SRR11978238'... 2020-07-14T02:23:48 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:24:20 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:24:20 prefetch.2.10.7: 'SRR11978238' is valid 2020-07-14T02:24:20 prefetch.2.10.7: 1) 'SRR11978238' was downloaded successfully 2020-07-14T02:24:20 prefetch.2.10.7: 'SRR11978238' has 0 unresolved dependencies Read 3974048 spots for SRR11978238/SRR11978238.sra Written 3974048 spots for SRR11978238/SRR11978238.sra 2020-07-14T02:24:44 prefetch.2.10.7: 1) Downloading 'SRR11978239'... 2020-07-14T02:24:44 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:25:24 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:25:25 prefetch.2.10.7: 'SRR11978239' is valid 2020-07-14T02:25:25 prefetch.2.10.7: 1) 'SRR11978239' was downloaded successfully 2020-07-14T02:25:25 prefetch.2.10.7: 'SRR11978239' has 0 unresolved dependencies Read 4043353 spots for SRR11978239/SRR11978239.sra Written 4043353 spots for SRR11978239/SRR11978239.sra 2020-07-14T02:25:48 prefetch.2.10.7: 1) Downloading 'SRR11978240'... 2020-07-14T02:25:48 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:26:57 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:26:57 prefetch.2.10.7: 'SRR11978240' is valid 2020-07-14T02:26:57 prefetch.2.10.7: 1) 'SRR11978240' was downloaded successfully 2020-07-14T02:26:57 prefetch.2.10.7: 'SRR11978240' has 0 unresolved dependencies Read 4037216 spots for SRR11978240/SRR11978240.sra Written 4037216 spots for SRR11978240/SRR11978240.sra 2020-07-14T02:27:25 prefetch.2.10.7: 1) Downloading 'SRR11978241'... 2020-07-14T02:27:25 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:28:14 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:28:14 prefetch.2.10.7: 'SRR11978241' is valid 2020-07-14T02:28:14 prefetch.2.10.7: 1) 'SRR11978241' was downloaded successfully 2020-07-14T02:28:14 prefetch.2.10.7: 'SRR11978241' has 0 unresolved dependencies Read 3998767 spots for SRR11978241/SRR11978241.sra Written 3998767 spots for SRR11978241/SRR11978241.sra 2020-07-14T02:28:41 prefetch.2.10.7: 1) Downloading 'SRR11978242'... 2020-07-14T02:28:41 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:29:00 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:29:00 prefetch.2.10.7: 'SRR11978242' is valid 2020-07-14T02:29:00 prefetch.2.10.7: 1) 'SRR11978242' was downloaded successfully 2020-07-14T02:29:00 prefetch.2.10.7: 'SRR11978242' has 0 unresolved dependencies Read 344382 spots for SRR11978242/SRR11978242.sra Written 344382 spots for SRR11978242/SRR11978242.sra 2020-07-14T02:29:15 prefetch.2.10.7: 1) Downloading 'SRR11978243'... 2020-07-14T02:29:15 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:29:54 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:29:54 prefetch.2.10.7: 'SRR11978243' is valid 2020-07-14T02:29:54 prefetch.2.10.7: 1) 'SRR11978243' was downloaded successfully 2020-07-14T02:29:54 prefetch.2.10.7: 'SRR11978243' has 0 unresolved dependencies Read 361015 spots for SRR11978243/SRR11978243.sra Written 361015 spots for SRR11978243/SRR11978243.sra 2020-07-14T02:30:10 prefetch.2.10.7: 1) Downloading 'SRR11978244'... 2020-07-14T02:30:10 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:30:24 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:30:24 prefetch.2.10.7: 'SRR11978244' is valid 2020-07-14T02:30:24 prefetch.2.10.7: 1) 'SRR11978244' was downloaded successfully 2020-07-14T02:30:24 prefetch.2.10.7: 'SRR11978244' has 0 unresolved dependencies Read 362038 spots for SRR11978244/SRR11978244.sra Written 362038 spots for SRR11978244/SRR11978244.sra 2020-07-14T02:30:40 prefetch.2.10.7: 1) Downloading 'SRR11978245'... 2020-07-14T02:30:40 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:31:57 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:31:57 prefetch.2.10.7: 'SRR11978245' is valid 2020-07-14T02:31:57 prefetch.2.10.7: 1) 'SRR11978245' was downloaded successfully 2020-07-14T02:31:57 prefetch.2.10.7: 'SRR11978245' has 0 unresolved dependencies Read 7519134 spots for SRR11978245/SRR11978245.sra Written 7519134 spots for SRR11978245/SRR11978245.sra 2020-07-14T02:32:36 prefetch.2.10.7: 1) Downloading 'SRR11978246'... 2020-07-14T02:32:36 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:34:09 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:34:10 prefetch.2.10.7: 'SRR11978246' is valid 2020-07-14T02:34:10 prefetch.2.10.7: 1) 'SRR11978246' was downloaded successfully 2020-07-14T02:34:10 prefetch.2.10.7: 'SRR11978246' has 0 unresolved dependencies Read 7373120 spots for SRR11978246/SRR11978246.sra Written 7373120 spots for SRR11978246/SRR11978246.sra 2020-07-14T02:34:48 prefetch.2.10.7: 1) Downloading 'SRR11978247'... 2020-07-14T02:34:48 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:35:29 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:35:30 prefetch.2.10.7: 'SRR11978247' is valid 2020-07-14T02:35:30 prefetch.2.10.7: 1) 'SRR11978247' was downloaded successfully 2020-07-14T02:35:30 prefetch.2.10.7: 'SRR11978247' has 0 unresolved dependencies Read 7527770 spots for SRR11978247/SRR11978247.sra Written 7527770 spots for SRR11978247/SRR11978247.sra 2020-07-14T02:36:09 prefetch.2.10.7: 1) Downloading 'SRR11978248'... 2020-07-14T02:36:09 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:37:10 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:37:10 prefetch.2.10.7: 'SRR11978248' is valid 2020-07-14T02:37:10 prefetch.2.10.7: 1) 'SRR11978248' was downloaded successfully 2020-07-14T02:37:10 prefetch.2.10.7: 'SRR11978248' has 0 unresolved dependencies Read 7457624 spots for SRR11978248/SRR11978248.sra Written 7457624 spots for SRR11978248/SRR11978248.sra 2020-07-14T02:37:45 prefetch.2.10.7: 1) Downloading 'SRR11978249'... 2020-07-14T02:37:45 prefetch.2.10.7: Downloading via HTTPS... 2020-07-14T02:38:09 prefetch.2.10.7: HTTPS download succeed 2020-07-14T02:38:09 prefetch.2.10.7: 'SRR11978249' is valid 2020-07-14T02:38:09 prefetch.2.10.7: 1) 'SRR11978249' was downloaded successfully 2020-07-14T02:38:09 prefetch.2.10.7: 'SRR11978249' has 0 unresolved dependencies Read 1752665 spots for SRR11978249/SRR11978249.sra Written 1752665 spots for SRR11978249/SRR11978249.sra fastq に変換しました。 bowtie でマッピング中... Your job 6627640 ("srTdm6") has been submitted Time loading reference: 00:00:00 Time loading forward index: 00:00:01 Time loading mirror index: 00:00:00 Multiseed full-index search: 00:09:21 54372229 reads; of these: 54372229 (100.00%) were unpaired; of these: 45359621 (83.42%) aligned 0 times 6826868 (12.56%) aligned exactly 1 time 2185740 (4.02%) aligned >1 times 16.58% overall alignment rate Time searching: 00:09:22 Overall time: 00:09:22 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_rmdupse_core] 1994085 / 9012608 = 0.2213 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:51:02: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.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:51:02: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:51:02: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:51:09: 1000000 INFO @ Tue, 14 Jul 2020 11:51:14: 2000000 INFO @ Tue, 14 Jul 2020 11:51:19: 3000000 INFO @ Tue, 14 Jul 2020 11:51:25: 4000000 INFO @ Tue, 14 Jul 2020 11:51:30: 5000000 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:51:32: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.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:51:32: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:51:32: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:51:36: 6000000 INFO @ Tue, 14 Jul 2020 11:51:38: 1000000 INFO @ Tue, 14 Jul 2020 11:51:41: 7000000 INFO @ Tue, 14 Jul 2020 11:51:42: #1 tag size is determined as 65 bps INFO @ Tue, 14 Jul 2020 11:51:42: #1 tag size = 65 INFO @ Tue, 14 Jul 2020 11:51:42: #1 total tags in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:51:42: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:51:42: #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:51:42: #1 tags after filtering in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:51:42: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:51:42: #1 finished! INFO @ Tue, 14 Jul 2020 11:51:42: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:51:42: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:51:42: #2 number of paired peaks: 335 WARNING @ Tue, 14 Jul 2020 11:51:42: Fewer paired peaks (335) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 335 pairs to build model! INFO @ Tue, 14 Jul 2020 11:51:42: start model_add_line... INFO @ Tue, 14 Jul 2020 11:51:42: start X-correlation... INFO @ Tue, 14 Jul 2020 11:51:42: end of X-cor INFO @ Tue, 14 Jul 2020 11:51:42: #2 finished! INFO @ Tue, 14 Jul 2020 11:51:42: #2 predicted fragment length is 59 bps INFO @ Tue, 14 Jul 2020 11:51:42: #2 alternative fragment length(s) may be 59 bps INFO @ Tue, 14 Jul 2020 11:51:42: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05_model.r WARNING @ Tue, 14 Jul 2020 11:51:42: #2 Since the d (59) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:51:42: #2 You may need to consider one of the other alternative d(s): 59 WARNING @ Tue, 14 Jul 2020 11:51:42: #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:51:42: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:51:42: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 14 Jul 2020 11:51:44: 2000000 INFO @ Tue, 14 Jul 2020 11:51:49: 3000000 INFO @ Tue, 14 Jul 2020 11:51:55: 4000000 INFO @ Tue, 14 Jul 2020 11:51:57: #3 Call peaks for each chromosome... INFO @ Tue, 14 Jul 2020 11:52:00: 5000000 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:52:02: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.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:52:02: #1 read tag files... INFO @ Tue, 14 Jul 2020 11:52:02: #1 read treatment tags... INFO @ Tue, 14 Jul 2020 11:52:04: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05_peaks.xls INFO @ Tue, 14 Jul 2020 11:52:04: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:52:04: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.05_summits.bed INFO @ Tue, 14 Jul 2020 11:52:04: Done! pass1 - making usageList (15 chroms): 1 millis pass2 - checking and writing primary data (1865 records, 4 fields): 3 millis CompletedMACS2peakCalling INFO @ Tue, 14 Jul 2020 11:52:06: 6000000 INFO @ Tue, 14 Jul 2020 11:52:09: 1000000 INFO @ Tue, 14 Jul 2020 11:52:12: 7000000 INFO @ Tue, 14 Jul 2020 11:52:12: #1 tag size is determined as 65 bps INFO @ Tue, 14 Jul 2020 11:52:12: #1 tag size = 65 INFO @ Tue, 14 Jul 2020 11:52:12: #1 total tags in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:52:12: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:52:12: #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:52:12: #1 tags after filtering in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:52:12: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:52:12: #1 finished! INFO @ Tue, 14 Jul 2020 11:52:12: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:52:12: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:52:12: #2 number of paired peaks: 335 WARNING @ Tue, 14 Jul 2020 11:52:12: Fewer paired peaks (335) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 335 pairs to build model! INFO @ Tue, 14 Jul 2020 11:52:12: start model_add_line... INFO @ Tue, 14 Jul 2020 11:52:13: start X-correlation... INFO @ Tue, 14 Jul 2020 11:52:13: end of X-cor INFO @ Tue, 14 Jul 2020 11:52:13: #2 finished! INFO @ Tue, 14 Jul 2020 11:52:13: #2 predicted fragment length is 59 bps INFO @ Tue, 14 Jul 2020 11:52:13: #2 alternative fragment length(s) may be 59 bps INFO @ Tue, 14 Jul 2020 11:52:13: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10_model.r WARNING @ Tue, 14 Jul 2020 11:52:13: #2 Since the d (59) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:52:13: #2 You may need to consider one of the other alternative d(s): 59 WARNING @ Tue, 14 Jul 2020 11:52:13: #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:52:13: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:52:13: #3 Pre-compute pvalue-qvalue table... INFO @ Tue, 14 Jul 2020 11:52:14: 2000000 INFO @ Tue, 14 Jul 2020 11:52:20: 3000000 INFO @ Tue, 14 Jul 2020 11:52:26: 4000000 INFO @ Tue, 14 Jul 2020 11:52:27: #3 Call peaks for each chromosome... INFO @ Tue, 14 Jul 2020 11:52:31: 5000000 INFO @ Tue, 14 Jul 2020 11:52:35: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10_peaks.xls INFO @ Tue, 14 Jul 2020 11:52:35: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:52:35: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.10_summits.bed INFO @ Tue, 14 Jul 2020 11:52:35: Done! pass1 - making usageList (14 chroms): 0 millis pass2 - checking and writing primary data (884 records, 4 fields): 3 millis CompletedMACS2peakCalling INFO @ Tue, 14 Jul 2020 11:52:37: 6000000 INFO @ Tue, 14 Jul 2020 11:52:42: 7000000 INFO @ Tue, 14 Jul 2020 11:52:42: #1 tag size is determined as 65 bps INFO @ Tue, 14 Jul 2020 11:52:42: #1 tag size = 65 INFO @ Tue, 14 Jul 2020 11:52:42: #1 total tags in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:52:42: #1 user defined the maximum tags... INFO @ Tue, 14 Jul 2020 11:52:42: #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:52:42: #1 tags after filtering in treatment: 7018523 INFO @ Tue, 14 Jul 2020 11:52:42: #1 Redundant rate of treatment: 0.00 INFO @ Tue, 14 Jul 2020 11:52:42: #1 finished! INFO @ Tue, 14 Jul 2020 11:52:42: #2 Build Peak Model... INFO @ Tue, 14 Jul 2020 11:52:42: #2 looking for paired plus/minus strand peaks... INFO @ Tue, 14 Jul 2020 11:52:43: #2 number of paired peaks: 335 WARNING @ Tue, 14 Jul 2020 11:52:43: Fewer paired peaks (335) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 335 pairs to build model! INFO @ Tue, 14 Jul 2020 11:52:43: start model_add_line... INFO @ Tue, 14 Jul 2020 11:52:43: start X-correlation... INFO @ Tue, 14 Jul 2020 11:52:43: end of X-cor INFO @ Tue, 14 Jul 2020 11:52:43: #2 finished! INFO @ Tue, 14 Jul 2020 11:52:43: #2 predicted fragment length is 59 bps INFO @ Tue, 14 Jul 2020 11:52:43: #2 alternative fragment length(s) may be 59 bps INFO @ Tue, 14 Jul 2020 11:52:43: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20_model.r WARNING @ Tue, 14 Jul 2020 11:52:43: #2 Since the d (59) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Tue, 14 Jul 2020 11:52:43: #2 You may need to consider one of the other alternative d(s): 59 WARNING @ Tue, 14 Jul 2020 11:52:43: #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:52:43: #3 Call peaks... INFO @ Tue, 14 Jul 2020 11:52:43: #3 Pre-compute pvalue-qvalue table... BedGraph に変換しました。 BigWig に変換中... INFO @ Tue, 14 Jul 2020 11:52:57: #3 Call peaks for each chromosome... BigWig に変換しました。 INFO @ Tue, 14 Jul 2020 11:53:05: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20_peaks.xls INFO @ Tue, 14 Jul 2020 11:53:05: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20_peaks.narrowPeak INFO @ Tue, 14 Jul 2020 11:53:05: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX8521379/SRX8521379.20_summits.bed INFO @ Tue, 14 Jul 2020 11:53:05: Done! pass1 - making usageList (10 chroms): 1 millis pass2 - checking and writing primary data (190 records, 4 fields): 1 millis CompletedMACS2peakCalling