Job ID = 12265233 SRX = SRX4518003 Genome = dm3 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... Read 53831592 spots for SRR7655607/SRR7655607.sra Written 53831592 spots for SRR7655607/SRR7655607.sra fastq に変換しました。 bowtie でマッピング中... Your job 12265585 ("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:42:07 53831592 reads; of these: 53831592 (100.00%) were paired; of these: 7603123 (14.12%) aligned concordantly 0 times 43969864 (81.68%) aligned concordantly exactly 1 time 2258605 (4.20%) aligned concordantly >1 times ---- 7603123 pairs aligned concordantly 0 times; of these: 3748620 (49.30%) aligned discordantly 1 time ---- 3854503 pairs aligned 0 times concordantly or discordantly; of these: 7709006 mates make up the pairs; of these: 6099643 (79.12%) aligned 0 times 1137482 (14.76%) aligned exactly 1 time 471881 (6.12%) aligned >1 times 94.33% overall alignment rate Time searching: 00:42:07 Overall time: 00:42:07 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 44 files... [bam_rmdup_core] processing reference chr2L... [bam_rmdup_core] processing reference chr2LHet... [bam_rmdup_core] processing reference chr2R... [bam_rmdup_core] processing reference chr2RHet... [bam_rmdup_core] processing reference chr3L... [bam_rmdup_core] processing reference chr3LHet... [bam_rmdup_core] processing reference chr3R... [bam_rmdup_core] processing reference chr3RHet... [bam_rmdup_core] processing reference chr4... [bam_rmdup_core] processing reference chrM... [bam_rmdup_core] processing reference chrU... [bam_rmdup_core] processing reference chrUextra... [bam_rmdup_core] processing reference chrX... [bam_rmdup_core] processing reference chrXHet... [bam_rmdup_core] processing reference chrYHet... [bam_rmdup_core] 5552209 / 49258683 = 0.1127 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:39:11: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.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 @ Sat, 03 Apr 2021 07:39:11: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:39:11: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:39:19: 1000000 INFO @ Sat, 03 Apr 2021 07:39:26: 2000000 INFO @ Sat, 03 Apr 2021 07:39:34: 3000000 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:39:41: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.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 @ Sat, 03 Apr 2021 07:39:41: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:39:41: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:39:42: 4000000 INFO @ Sat, 03 Apr 2021 07:39:45: 1000000 INFO @ Sat, 03 Apr 2021 07:39:49: 5000000 INFO @ Sat, 03 Apr 2021 07:39:51: 2000000 INFO @ Sat, 03 Apr 2021 07:39:55: 3000000 INFO @ Sat, 03 Apr 2021 07:39:57: 6000000 INFO @ Sat, 03 Apr 2021 07:40:00: 4000000 INFO @ Sat, 03 Apr 2021 07:40:04: 7000000 INFO @ Sat, 03 Apr 2021 07:40:05: 5000000 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:40:10: 6000000 INFO @ Sat, 03 Apr 2021 07:40:11: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.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 @ Sat, 03 Apr 2021 07:40:11: #1 read tag files... INFO @ Sat, 03 Apr 2021 07:40:11: #1 read treatment tags... INFO @ Sat, 03 Apr 2021 07:40:11: 8000000 INFO @ Sat, 03 Apr 2021 07:40:15: 7000000 INFO @ Sat, 03 Apr 2021 07:40:18: 1000000 INFO @ Sat, 03 Apr 2021 07:40:18: 9000000 INFO @ Sat, 03 Apr 2021 07:40:20: 8000000 INFO @ Sat, 03 Apr 2021 07:40:25: 9000000 INFO @ Sat, 03 Apr 2021 07:40:26: 2000000 INFO @ Sat, 03 Apr 2021 07:40:26: 10000000 INFO @ Sat, 03 Apr 2021 07:40:30: 10000000 INFO @ Sat, 03 Apr 2021 07:40:33: 3000000 INFO @ Sat, 03 Apr 2021 07:40:33: 11000000 INFO @ Sat, 03 Apr 2021 07:40:35: 11000000 INFO @ Sat, 03 Apr 2021 07:40:41: 12000000 INFO @ Sat, 03 Apr 2021 07:40:41: 4000000 INFO @ Sat, 03 Apr 2021 07:40:41: 12000000 INFO @ Sat, 03 Apr 2021 07:40:46: 13000000 INFO @ Sat, 03 Apr 2021 07:40:48: 5000000 INFO @ Sat, 03 Apr 2021 07:40:48: 13000000 INFO @ Sat, 03 Apr 2021 07:40:51: 14000000 INFO @ Sat, 03 Apr 2021 07:40:56: 14000000 INFO @ Sat, 03 Apr 2021 07:40:56: 6000000 INFO @ Sat, 03 Apr 2021 07:40:57: 15000000 INFO @ Sat, 03 Apr 2021 07:41:02: 16000000 INFO @ Sat, 03 Apr 2021 07:41:03: 15000000 INFO @ Sat, 03 Apr 2021 07:41:03: 7000000 INFO @ Sat, 03 Apr 2021 07:41:07: 17000000 INFO @ Sat, 03 Apr 2021 07:41:10: 16000000 INFO @ Sat, 03 Apr 2021 07:41:11: 8000000 INFO @ Sat, 03 Apr 2021 07:41:12: 18000000 INFO @ Sat, 03 Apr 2021 07:41:17: 17000000 INFO @ Sat, 03 Apr 2021 07:41:17: 19000000 INFO @ Sat, 03 Apr 2021 07:41:18: 9000000 INFO @ Sat, 03 Apr 2021 07:41:23: 20000000 INFO @ Sat, 03 Apr 2021 07:41:25: 18000000 INFO @ Sat, 03 Apr 2021 07:41:25: 10000000 INFO @ Sat, 03 Apr 2021 07:41:28: 21000000 INFO @ Sat, 03 Apr 2021 07:41:32: 19000000 INFO @ Sat, 03 Apr 2021 07:41:33: 11000000 INFO @ Sat, 03 Apr 2021 07:41:33: 22000000 INFO @ Sat, 03 Apr 2021 07:41:38: 23000000 INFO @ Sat, 03 Apr 2021 07:41:40: 20000000 INFO @ Sat, 03 Apr 2021 07:41:40: 12000000 INFO @ Sat, 03 Apr 2021 07:41:43: 24000000 INFO @ Sat, 03 Apr 2021 07:41:47: 21000000 INFO @ Sat, 03 Apr 2021 07:41:48: 13000000 INFO @ Sat, 03 Apr 2021 07:41:48: 25000000 INFO @ Sat, 03 Apr 2021 07:41:54: 26000000 INFO @ Sat, 03 Apr 2021 07:41:55: 22000000 INFO @ Sat, 03 Apr 2021 07:41:55: 14000000 INFO @ Sat, 03 Apr 2021 07:41:59: 27000000 INFO @ Sat, 03 Apr 2021 07:42:02: 23000000 INFO @ Sat, 03 Apr 2021 07:42:03: 15000000 INFO @ Sat, 03 Apr 2021 07:42:05: 28000000 INFO @ Sat, 03 Apr 2021 07:42:10: 24000000 INFO @ Sat, 03 Apr 2021 07:42:10: 16000000 INFO @ Sat, 03 Apr 2021 07:42:10: 29000000 INFO @ Sat, 03 Apr 2021 07:42:16: 30000000 INFO @ Sat, 03 Apr 2021 07:42:17: 25000000 INFO @ Sat, 03 Apr 2021 07:42:18: 17000000 INFO @ Sat, 03 Apr 2021 07:42:21: 31000000 INFO @ Sat, 03 Apr 2021 07:42:25: 26000000 INFO @ Sat, 03 Apr 2021 07:42:25: 18000000 INFO @ Sat, 03 Apr 2021 07:42:26: 32000000 INFO @ Sat, 03 Apr 2021 07:42:31: 33000000 INFO @ Sat, 03 Apr 2021 07:42:32: 27000000 INFO @ Sat, 03 Apr 2021 07:42:33: 19000000 INFO @ Sat, 03 Apr 2021 07:42:36: 34000000 INFO @ Sat, 03 Apr 2021 07:42:40: 28000000 INFO @ Sat, 03 Apr 2021 07:42:40: 20000000 INFO @ Sat, 03 Apr 2021 07:42:41: 35000000 INFO @ Sat, 03 Apr 2021 07:42:47: 36000000 INFO @ Sat, 03 Apr 2021 07:42:48: 29000000 INFO @ Sat, 03 Apr 2021 07:42:48: 21000000 INFO @ Sat, 03 Apr 2021 07:42:52: 37000000 INFO @ Sat, 03 Apr 2021 07:42:55: 30000000 INFO @ Sat, 03 Apr 2021 07:42:55: 22000000 INFO @ Sat, 03 Apr 2021 07:42:57: 38000000 INFO @ Sat, 03 Apr 2021 07:43:02: 39000000 INFO @ Sat, 03 Apr 2021 07:43:03: 31000000 INFO @ Sat, 03 Apr 2021 07:43:03: 23000000 INFO @ Sat, 03 Apr 2021 07:43:07: 40000000 INFO @ Sat, 03 Apr 2021 07:43:10: 32000000 INFO @ Sat, 03 Apr 2021 07:43:11: 24000000 INFO @ Sat, 03 Apr 2021 07:43:12: 41000000 INFO @ Sat, 03 Apr 2021 07:43:17: 42000000 INFO @ Sat, 03 Apr 2021 07:43:18: 33000000 INFO @ Sat, 03 Apr 2021 07:43:18: 25000000 INFO @ Sat, 03 Apr 2021 07:43:22: 43000000 INFO @ Sat, 03 Apr 2021 07:43:25: 34000000 INFO @ Sat, 03 Apr 2021 07:43:26: 26000000 INFO @ Sat, 03 Apr 2021 07:43:27: 44000000 INFO @ Sat, 03 Apr 2021 07:43:32: 45000000 INFO @ Sat, 03 Apr 2021 07:43:33: 35000000 INFO @ Sat, 03 Apr 2021 07:43:33: 27000000 INFO @ Sat, 03 Apr 2021 07:43:37: 46000000 INFO @ Sat, 03 Apr 2021 07:43:41: 36000000 INFO @ Sat, 03 Apr 2021 07:43:41: 28000000 INFO @ Sat, 03 Apr 2021 07:43:42: 47000000 INFO @ Sat, 03 Apr 2021 07:43:47: 48000000 INFO @ Sat, 03 Apr 2021 07:43:48: 37000000 INFO @ Sat, 03 Apr 2021 07:43:49: 29000000 INFO @ Sat, 03 Apr 2021 07:43:52: 49000000 INFO @ Sat, 03 Apr 2021 07:43:56: 38000000 INFO @ Sat, 03 Apr 2021 07:43:56: 30000000 INFO @ Sat, 03 Apr 2021 07:43:57: 50000000 INFO @ Sat, 03 Apr 2021 07:44:03: 51000000 INFO @ Sat, 03 Apr 2021 07:44:04: 39000000 INFO @ Sat, 03 Apr 2021 07:44:04: 31000000 INFO @ Sat, 03 Apr 2021 07:44:08: 52000000 INFO @ Sat, 03 Apr 2021 07:44:12: 40000000 INFO @ Sat, 03 Apr 2021 07:44:12: 32000000 INFO @ Sat, 03 Apr 2021 07:44:13: 53000000 INFO @ Sat, 03 Apr 2021 07:44:18: 54000000 INFO @ Sat, 03 Apr 2021 07:44:19: 41000000 INFO @ Sat, 03 Apr 2021 07:44:20: 33000000 INFO @ Sat, 03 Apr 2021 07:44:23: 55000000 INFO @ Sat, 03 Apr 2021 07:44:27: 42000000 INFO @ Sat, 03 Apr 2021 07:44:27: 34000000 INFO @ Sat, 03 Apr 2021 07:44:28: 56000000 INFO @ Sat, 03 Apr 2021 07:44:33: 57000000 INFO @ Sat, 03 Apr 2021 07:44:34: 43000000 INFO @ Sat, 03 Apr 2021 07:44:35: 35000000 INFO @ Sat, 03 Apr 2021 07:44:38: 58000000 INFO @ Sat, 03 Apr 2021 07:44:41: 44000000 INFO @ Sat, 03 Apr 2021 07:44:42: 36000000 INFO @ Sat, 03 Apr 2021 07:44:43: 59000000 INFO @ Sat, 03 Apr 2021 07:44:48: 60000000 INFO @ Sat, 03 Apr 2021 07:44:49: 45000000 INFO @ Sat, 03 Apr 2021 07:44:49: 37000000 INFO @ Sat, 03 Apr 2021 07:44:53: 61000000 INFO @ Sat, 03 Apr 2021 07:44:56: 46000000 INFO @ Sat, 03 Apr 2021 07:44:56: 38000000 INFO @ Sat, 03 Apr 2021 07:44:57: 62000000 INFO @ Sat, 03 Apr 2021 07:45:02: 63000000 INFO @ Sat, 03 Apr 2021 07:45:03: 47000000 INFO @ Sat, 03 Apr 2021 07:45:04: 39000000 INFO @ Sat, 03 Apr 2021 07:45:07: 64000000 INFO @ Sat, 03 Apr 2021 07:45:11: 48000000 INFO @ Sat, 03 Apr 2021 07:45:12: 40000000 INFO @ Sat, 03 Apr 2021 07:45:12: 65000000 INFO @ Sat, 03 Apr 2021 07:45:17: 66000000 INFO @ Sat, 03 Apr 2021 07:45:18: 49000000 INFO @ Sat, 03 Apr 2021 07:45:19: 41000000 INFO @ Sat, 03 Apr 2021 07:45:22: 67000000 INFO @ Sat, 03 Apr 2021 07:45:26: 50000000 INFO @ Sat, 03 Apr 2021 07:45:27: 42000000 INFO @ Sat, 03 Apr 2021 07:45:27: 68000000 INFO @ Sat, 03 Apr 2021 07:45:32: 69000000 INFO @ Sat, 03 Apr 2021 07:45:34: 51000000 INFO @ Sat, 03 Apr 2021 07:45:35: 43000000 INFO @ Sat, 03 Apr 2021 07:45:38: 70000000 INFO @ Sat, 03 Apr 2021 07:45:41: 52000000 INFO @ Sat, 03 Apr 2021 07:45:42: 44000000 INFO @ Sat, 03 Apr 2021 07:45:43: 71000000 INFO @ Sat, 03 Apr 2021 07:45:48: 72000000 INFO @ Sat, 03 Apr 2021 07:45:49: 53000000 INFO @ Sat, 03 Apr 2021 07:45:50: 45000000 INFO @ Sat, 03 Apr 2021 07:45:53: 73000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Sat, 03 Apr 2021 07:45:56: 54000000 INFO @ Sat, 03 Apr 2021 07:45:57: 46000000 INFO @ Sat, 03 Apr 2021 07:45:59: 74000000 INFO @ Sat, 03 Apr 2021 07:46:04: 55000000 INFO @ Sat, 03 Apr 2021 07:46:04: 75000000 INFO @ Sat, 03 Apr 2021 07:46:05: 47000000 INFO @ Sat, 03 Apr 2021 07:46:09: 76000000 INFO @ Sat, 03 Apr 2021 07:46:11: 56000000 INFO @ Sat, 03 Apr 2021 07:46:13: 48000000 INFO @ Sat, 03 Apr 2021 07:46:15: 77000000 INFO @ Sat, 03 Apr 2021 07:46:19: 57000000 INFO @ Sat, 03 Apr 2021 07:46:20: 78000000 INFO @ Sat, 03 Apr 2021 07:46:20: 49000000 INFO @ Sat, 03 Apr 2021 07:46:25: 79000000 INFO @ Sat, 03 Apr 2021 07:46:27: 58000000 INFO @ Sat, 03 Apr 2021 07:46:28: 50000000 INFO @ Sat, 03 Apr 2021 07:46:30: 80000000 INFO @ Sat, 03 Apr 2021 07:46:34: 59000000 INFO @ Sat, 03 Apr 2021 07:46:35: 81000000 INFO @ Sat, 03 Apr 2021 07:46:36: 51000000 INFO @ Sat, 03 Apr 2021 07:46:40: 82000000 INFO @ Sat, 03 Apr 2021 07:46:42: 60000000 INFO @ Sat, 03 Apr 2021 07:46:43: 52000000 INFO @ Sat, 03 Apr 2021 07:46:45: 83000000 INFO @ Sat, 03 Apr 2021 07:46:49: 61000000 INFO @ Sat, 03 Apr 2021 07:46:50: 84000000 INFO @ Sat, 03 Apr 2021 07:46:50: 53000000 INFO @ Sat, 03 Apr 2021 07:46:55: 85000000 INFO @ Sat, 03 Apr 2021 07:46:57: 62000000 INFO @ Sat, 03 Apr 2021 07:46:58: 54000000 INFO @ Sat, 03 Apr 2021 07:47:01: 86000000 INFO @ Sat, 03 Apr 2021 07:47:04: 63000000 INFO @ Sat, 03 Apr 2021 07:47:06: 55000000 INFO @ Sat, 03 Apr 2021 07:47:06: 87000000 INFO @ Sat, 03 Apr 2021 07:47:11: 88000000 INFO @ Sat, 03 Apr 2021 07:47:12: 64000000 INFO @ Sat, 03 Apr 2021 07:47:13: 56000000 INFO @ Sat, 03 Apr 2021 07:47:18: 89000000 INFO @ Sat, 03 Apr 2021 07:47:20: 65000000 INFO @ Sat, 03 Apr 2021 07:47:21: 57000000 INFO @ Sat, 03 Apr 2021 07:47:23: 90000000 INFO @ Sat, 03 Apr 2021 07:47:26: #1 tag size is determined as 50 bps INFO @ Sat, 03 Apr 2021 07:47:26: #1 tag size = 50 INFO @ Sat, 03 Apr 2021 07:47:26: #1 total tags in treatment: 41372799 INFO @ Sat, 03 Apr 2021 07:47:26: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:47:26: #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:47:27: 66000000 INFO @ Sat, 03 Apr 2021 07:47:27: #1 tags after filtering in treatment: 24695561 INFO @ Sat, 03 Apr 2021 07:47:27: #1 Redundant rate of treatment: 0.40 INFO @ Sat, 03 Apr 2021 07:47:27: #1 finished! INFO @ Sat, 03 Apr 2021 07:47:27: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:47:27: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:47:28: 58000000 INFO @ Sat, 03 Apr 2021 07:47:29: #2 number of paired peaks: 189 WARNING @ Sat, 03 Apr 2021 07:47:29: Fewer paired peaks (189) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 189 pairs to build model! INFO @ Sat, 03 Apr 2021 07:47:29: start model_add_line... INFO @ Sat, 03 Apr 2021 07:47:29: start X-correlation... INFO @ Sat, 03 Apr 2021 07:47:29: end of X-cor INFO @ Sat, 03 Apr 2021 07:47:29: #2 finished! INFO @ Sat, 03 Apr 2021 07:47:29: #2 predicted fragment length is 52 bps INFO @ Sat, 03 Apr 2021 07:47:29: #2 alternative fragment length(s) may be 4,12,32,52 bps INFO @ Sat, 03 Apr 2021 07:47:29: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10_model.r WARNING @ Sat, 03 Apr 2021 07:47:29: #2 Since the d (52) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Sat, 03 Apr 2021 07:47:29: #2 You may need to consider one of the other alternative d(s): 4,12,32,52 WARNING @ Sat, 03 Apr 2021 07:47:29: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Sat, 03 Apr 2021 07:47:29: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:47:29: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:47:34: 67000000 INFO @ Sat, 03 Apr 2021 07:47:36: 59000000 INFO @ Sat, 03 Apr 2021 07:47:42: 68000000 INFO @ Sat, 03 Apr 2021 07:47:43: 60000000 BigWig に変換しました。 INFO @ Sat, 03 Apr 2021 07:47:50: 69000000 INFO @ Sat, 03 Apr 2021 07:47:51: 61000000 INFO @ Sat, 03 Apr 2021 07:47:57: 70000000 INFO @ Sat, 03 Apr 2021 07:47:59: 62000000 INFO @ Sat, 03 Apr 2021 07:48:05: 71000000 INFO @ Sat, 03 Apr 2021 07:48:06: 63000000 INFO @ Sat, 03 Apr 2021 07:48:10: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:48:12: 72000000 INFO @ Sat, 03 Apr 2021 07:48:14: 64000000 INFO @ Sat, 03 Apr 2021 07:48:20: 73000000 INFO @ Sat, 03 Apr 2021 07:48:22: 65000000 INFO @ Sat, 03 Apr 2021 07:48:27: 74000000 INFO @ Sat, 03 Apr 2021 07:48:30: 66000000 INFO @ Sat, 03 Apr 2021 07:48:32: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10_peaks.xls INFO @ Sat, 03 Apr 2021 07:48:32: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:48:32: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.10_summits.bed INFO @ Sat, 03 Apr 2021 07:48:32: Done! pass1 - making usageList (13 chroms): 2 millis pass2 - checking and writing primary data (4166 records, 4 fields): 8 millis CompletedMACS2peakCalling INFO @ Sat, 03 Apr 2021 07:48:35: 75000000 INFO @ Sat, 03 Apr 2021 07:48:37: 67000000 INFO @ Sat, 03 Apr 2021 07:48:43: 76000000 INFO @ Sat, 03 Apr 2021 07:48:45: 68000000 INFO @ Sat, 03 Apr 2021 07:48:50: 77000000 INFO @ Sat, 03 Apr 2021 07:48:53: 69000000 INFO @ Sat, 03 Apr 2021 07:48:58: 78000000 INFO @ Sat, 03 Apr 2021 07:49:00: 70000000 INFO @ Sat, 03 Apr 2021 07:49:05: 79000000 INFO @ Sat, 03 Apr 2021 07:49:08: 71000000 INFO @ Sat, 03 Apr 2021 07:49:12: 80000000 INFO @ Sat, 03 Apr 2021 07:49:14: 72000000 INFO @ Sat, 03 Apr 2021 07:49:19: 81000000 INFO @ Sat, 03 Apr 2021 07:49:22: 73000000 INFO @ Sat, 03 Apr 2021 07:49:27: 82000000 INFO @ Sat, 03 Apr 2021 07:49:29: 74000000 INFO @ Sat, 03 Apr 2021 07:49:34: 83000000 INFO @ Sat, 03 Apr 2021 07:49:37: 75000000 INFO @ Sat, 03 Apr 2021 07:49:42: 84000000 INFO @ Sat, 03 Apr 2021 07:49:45: 76000000 INFO @ Sat, 03 Apr 2021 07:49:50: 85000000 INFO @ Sat, 03 Apr 2021 07:49:52: 77000000 INFO @ Sat, 03 Apr 2021 07:49:57: 86000000 INFO @ Sat, 03 Apr 2021 07:50:00: 78000000 INFO @ Sat, 03 Apr 2021 07:50:04: 87000000 INFO @ Sat, 03 Apr 2021 07:50:07: 79000000 INFO @ Sat, 03 Apr 2021 07:50:12: 88000000 INFO @ Sat, 03 Apr 2021 07:50:15: 80000000 INFO @ Sat, 03 Apr 2021 07:50:20: 89000000 INFO @ Sat, 03 Apr 2021 07:50:22: 81000000 INFO @ Sat, 03 Apr 2021 07:50:27: 90000000 INFO @ Sat, 03 Apr 2021 07:50:30: 82000000 INFO @ Sat, 03 Apr 2021 07:50:32: #1 tag size is determined as 50 bps INFO @ Sat, 03 Apr 2021 07:50:32: #1 tag size = 50 INFO @ Sat, 03 Apr 2021 07:50:32: #1 total tags in treatment: 41372799 INFO @ Sat, 03 Apr 2021 07:50:32: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:50:32: #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:50:33: #1 tags after filtering in treatment: 24695561 INFO @ Sat, 03 Apr 2021 07:50:33: #1 Redundant rate of treatment: 0.40 INFO @ Sat, 03 Apr 2021 07:50:33: #1 finished! INFO @ Sat, 03 Apr 2021 07:50:33: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:50:33: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:50:34: #2 number of paired peaks: 189 WARNING @ Sat, 03 Apr 2021 07:50:34: Fewer paired peaks (189) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 189 pairs to build model! INFO @ Sat, 03 Apr 2021 07:50:34: start model_add_line... INFO @ Sat, 03 Apr 2021 07:50:34: start X-correlation... INFO @ Sat, 03 Apr 2021 07:50:34: end of X-cor INFO @ Sat, 03 Apr 2021 07:50:34: #2 finished! INFO @ Sat, 03 Apr 2021 07:50:34: #2 predicted fragment length is 52 bps INFO @ Sat, 03 Apr 2021 07:50:34: #2 alternative fragment length(s) may be 4,12,32,52 bps INFO @ Sat, 03 Apr 2021 07:50:34: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05_model.r WARNING @ Sat, 03 Apr 2021 07:50:34: #2 Since the d (52) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Sat, 03 Apr 2021 07:50:34: #2 You may need to consider one of the other alternative d(s): 4,12,32,52 WARNING @ Sat, 03 Apr 2021 07:50:34: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Sat, 03 Apr 2021 07:50:34: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:50:34: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:50:38: 83000000 INFO @ Sat, 03 Apr 2021 07:50:46: 84000000 INFO @ Sat, 03 Apr 2021 07:50:53: 85000000 INFO @ Sat, 03 Apr 2021 07:51:01: 86000000 INFO @ Sat, 03 Apr 2021 07:51:08: 87000000 INFO @ Sat, 03 Apr 2021 07:51:16: 88000000 INFO @ Sat, 03 Apr 2021 07:51:20: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:51:23: 89000000 INFO @ Sat, 03 Apr 2021 07:51:31: 90000000 INFO @ Sat, 03 Apr 2021 07:51:35: #1 tag size is determined as 50 bps INFO @ Sat, 03 Apr 2021 07:51:35: #1 tag size = 50 INFO @ Sat, 03 Apr 2021 07:51:35: #1 total tags in treatment: 41372799 INFO @ Sat, 03 Apr 2021 07:51:35: #1 user defined the maximum tags... INFO @ Sat, 03 Apr 2021 07:51:35: #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:51:36: #1 tags after filtering in treatment: 24695561 INFO @ Sat, 03 Apr 2021 07:51:36: #1 Redundant rate of treatment: 0.40 INFO @ Sat, 03 Apr 2021 07:51:36: #1 finished! INFO @ Sat, 03 Apr 2021 07:51:36: #2 Build Peak Model... INFO @ Sat, 03 Apr 2021 07:51:36: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 03 Apr 2021 07:51:38: #2 number of paired peaks: 189 WARNING @ Sat, 03 Apr 2021 07:51:38: Fewer paired peaks (189) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 189 pairs to build model! INFO @ Sat, 03 Apr 2021 07:51:38: start model_add_line... INFO @ Sat, 03 Apr 2021 07:51:38: start X-correlation... INFO @ Sat, 03 Apr 2021 07:51:38: end of X-cor INFO @ Sat, 03 Apr 2021 07:51:38: #2 finished! INFO @ Sat, 03 Apr 2021 07:51:38: #2 predicted fragment length is 52 bps INFO @ Sat, 03 Apr 2021 07:51:38: #2 alternative fragment length(s) may be 4,12,32,52 bps INFO @ Sat, 03 Apr 2021 07:51:38: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20_model.r WARNING @ Sat, 03 Apr 2021 07:51:38: #2 Since the d (52) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Sat, 03 Apr 2021 07:51:38: #2 You may need to consider one of the other alternative d(s): 4,12,32,52 WARNING @ Sat, 03 Apr 2021 07:51:38: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Sat, 03 Apr 2021 07:51:38: #3 Call peaks... INFO @ Sat, 03 Apr 2021 07:51:38: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 03 Apr 2021 07:51:48: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05_peaks.xls INFO @ Sat, 03 Apr 2021 07:51:48: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:51:48: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.05_summits.bed INFO @ Sat, 03 Apr 2021 07:51:48: Done! pass1 - making usageList (14 chroms): 3 millis pass2 - checking and writing primary data (12247 records, 4 fields): 21 millis CompletedMACS2peakCalling INFO @ Sat, 03 Apr 2021 07:52:24: #3 Call peaks for each chromosome... INFO @ Sat, 03 Apr 2021 07:52:51: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20_peaks.xls INFO @ Sat, 03 Apr 2021 07:52:51: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20_peaks.narrowPeak INFO @ Sat, 03 Apr 2021 07:52:51: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4518003/SRX4518003.20_summits.bed INFO @ Sat, 03 Apr 2021 07:52:51: Done! pass1 - making usageList (10 chroms): 0 millis pass2 - checking and writing primary data (307 records, 4 fields): 2 millis CompletedMACS2peakCalling