Job ID = 14167113 SRX = SRX4712935 Genome = dm3 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... Read 57043990 spots for SRR7874065/SRR7874065.sra Written 57043990 spots for SRR7874065/SRR7874065.sra fastq に変換しました。 bowtie でマッピング中... Your job 14168202 ("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: 02:53:13 57043990 reads; of these: 57043990 (100.00%) were paired; of these: 17671815 (30.98%) aligned concordantly 0 times 32372667 (56.75%) aligned concordantly exactly 1 time 6999508 (12.27%) aligned concordantly >1 times ---- 17671815 pairs aligned concordantly 0 times; of these: 7344931 (41.56%) aligned discordantly 1 time ---- 10326884 pairs aligned 0 times concordantly or discordantly; of these: 20653768 mates make up the pairs; of these: 16110383 (78.00%) aligned 0 times 1259477 (6.10%) aligned exactly 1 time 3283908 (15.90%) aligned >1 times 85.88% overall alignment rate Time searching: 02:53:14 Overall time: 02:53:14 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 40 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] 5096680 / 45983172 = 0.1108 in library ' ' BAM に変換しました。 Bed ファイルを作成中... WARNING: Skipping mount /opt/pkg/singularity/3.8.3/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Fri, 10 Dec 2021 13:38:00: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.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 @ Fri, 10 Dec 2021 13:38:00: #1 read tag files... INFO @ Fri, 10 Dec 2021 13:38:00: #1 read treatment tags... INFO @ Fri, 10 Dec 2021 13:38:07: 1000000 INFO @ Fri, 10 Dec 2021 13:38:13: 2000000 INFO @ Fri, 10 Dec 2021 13:38:19: 3000000 INFO @ Fri, 10 Dec 2021 13:38:26: 4000000 WARNING: Skipping mount /opt/pkg/singularity/3.8.3/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Fri, 10 Dec 2021 13:38:30: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.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 @ Fri, 10 Dec 2021 13:38:30: #1 read tag files... INFO @ Fri, 10 Dec 2021 13:38:30: #1 read treatment tags... INFO @ Fri, 10 Dec 2021 13:38:32: 5000000 INFO @ Fri, 10 Dec 2021 13:38:37: 1000000 INFO @ Fri, 10 Dec 2021 13:38:40: 6000000 INFO @ Fri, 10 Dec 2021 13:38:44: 2000000 INFO @ Fri, 10 Dec 2021 13:38:47: 7000000 INFO @ Fri, 10 Dec 2021 13:38:51: 3000000 INFO @ Fri, 10 Dec 2021 13:38:53: 8000000 BedGraph に変換中... WARNING: Skipping mount /opt/pkg/singularity/3.8.3/var/singularity/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container INFO @ Fri, 10 Dec 2021 13:38:58: 4000000 INFO @ Fri, 10 Dec 2021 13:39:00: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.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 @ Fri, 10 Dec 2021 13:39:00: #1 read tag files... INFO @ Fri, 10 Dec 2021 13:39:00: #1 read treatment tags... INFO @ Fri, 10 Dec 2021 13:39:00: 9000000 INFO @ Fri, 10 Dec 2021 13:39:05: 5000000 INFO @ Fri, 10 Dec 2021 13:39:08: 10000000 INFO @ Fri, 10 Dec 2021 13:39:08: 1000000 INFO @ Fri, 10 Dec 2021 13:39:13: 6000000 INFO @ Fri, 10 Dec 2021 13:39:14: 11000000 INFO @ Fri, 10 Dec 2021 13:39:18: 2000000 INFO @ Fri, 10 Dec 2021 13:39:21: 7000000 INFO @ Fri, 10 Dec 2021 13:39:22: 12000000 INFO @ Fri, 10 Dec 2021 13:39:25: 3000000 INFO @ Fri, 10 Dec 2021 13:39:28: 8000000 INFO @ Fri, 10 Dec 2021 13:39:29: 13000000 INFO @ Fri, 10 Dec 2021 13:39:32: 4000000 INFO @ Fri, 10 Dec 2021 13:39:34: 9000000 INFO @ Fri, 10 Dec 2021 13:39:36: 14000000 INFO @ Fri, 10 Dec 2021 13:39:39: 5000000 INFO @ Fri, 10 Dec 2021 13:39:41: 10000000 INFO @ Fri, 10 Dec 2021 13:39:44: 15000000 INFO @ Fri, 10 Dec 2021 13:39:47: 6000000 INFO @ Fri, 10 Dec 2021 13:39:47: 11000000 INFO @ Fri, 10 Dec 2021 13:39:51: 16000000 INFO @ Fri, 10 Dec 2021 13:39:54: 12000000 INFO @ Fri, 10 Dec 2021 13:39:55: 7000000 INFO @ Fri, 10 Dec 2021 13:39:59: 17000000 INFO @ Fri, 10 Dec 2021 13:40:01: 13000000 INFO @ Fri, 10 Dec 2021 13:40:03: 8000000 INFO @ Fri, 10 Dec 2021 13:40:07: 18000000 INFO @ Fri, 10 Dec 2021 13:40:08: 14000000 INFO @ Fri, 10 Dec 2021 13:40:11: 9000000 INFO @ Fri, 10 Dec 2021 13:40:15: 15000000 INFO @ Fri, 10 Dec 2021 13:40:16: 19000000 INFO @ Fri, 10 Dec 2021 13:40:19: 10000000 INFO @ Fri, 10 Dec 2021 13:40:22: 16000000 INFO @ Fri, 10 Dec 2021 13:40:23: 20000000 INFO @ Fri, 10 Dec 2021 13:40:26: 11000000 INFO @ Fri, 10 Dec 2021 13:40:29: 17000000 INFO @ Fri, 10 Dec 2021 13:40:31: 21000000 INFO @ Fri, 10 Dec 2021 13:40:34: 12000000 INFO @ Fri, 10 Dec 2021 13:40:36: 18000000 INFO @ Fri, 10 Dec 2021 13:40:38: 22000000 INFO @ Fri, 10 Dec 2021 13:40:41: 13000000 INFO @ Fri, 10 Dec 2021 13:40:43: 19000000 INFO @ Fri, 10 Dec 2021 13:40:46: 23000000 INFO @ Fri, 10 Dec 2021 13:40:49: 14000000 INFO @ Fri, 10 Dec 2021 13:40:53: 20000000 INFO @ Fri, 10 Dec 2021 13:40:54: 24000000 INFO @ Fri, 10 Dec 2021 13:40:57: 15000000 INFO @ Fri, 10 Dec 2021 13:41:02: 25000000 INFO @ Fri, 10 Dec 2021 13:41:03: 21000000 INFO @ Fri, 10 Dec 2021 13:41:04: 16000000 INFO @ Fri, 10 Dec 2021 13:41:11: 26000000 INFO @ Fri, 10 Dec 2021 13:41:11: 22000000 INFO @ Fri, 10 Dec 2021 13:41:12: 17000000 INFO @ Fri, 10 Dec 2021 13:41:19: 23000000 INFO @ Fri, 10 Dec 2021 13:41:20: 27000000 INFO @ Fri, 10 Dec 2021 13:41:20: 18000000 INFO @ Fri, 10 Dec 2021 13:41:26: 24000000 INFO @ Fri, 10 Dec 2021 13:41:30: 28000000 INFO @ Fri, 10 Dec 2021 13:41:30: 19000000 INFO @ Fri, 10 Dec 2021 13:41:36: 25000000 INFO @ Fri, 10 Dec 2021 13:41:40: 20000000 INFO @ Fri, 10 Dec 2021 13:41:40: 29000000 INFO @ Fri, 10 Dec 2021 13:41:45: 26000000 INFO @ Fri, 10 Dec 2021 13:41:48: 21000000 INFO @ Fri, 10 Dec 2021 13:41:49: 30000000 INFO @ Fri, 10 Dec 2021 13:41:53: 27000000 INFO @ Fri, 10 Dec 2021 13:41:56: 22000000 INFO @ Fri, 10 Dec 2021 13:41:56: 31000000 INFO @ Fri, 10 Dec 2021 13:42:00: 28000000 INFO @ Fri, 10 Dec 2021 13:42:04: 32000000 INFO @ Fri, 10 Dec 2021 13:42:04: 23000000 INFO @ Fri, 10 Dec 2021 13:42:07: 29000000 INFO @ Fri, 10 Dec 2021 13:42:12: 24000000 INFO @ Fri, 10 Dec 2021 13:42:13: 33000000 INFO @ Fri, 10 Dec 2021 13:42:15: 30000000 INFO @ Fri, 10 Dec 2021 13:42:19: 25000000 INFO @ Fri, 10 Dec 2021 13:42:20: 34000000 INFO @ Fri, 10 Dec 2021 13:42:22: 31000000 INFO @ Fri, 10 Dec 2021 13:42:27: 26000000 INFO @ Fri, 10 Dec 2021 13:42:27: 35000000 INFO @ Fri, 10 Dec 2021 13:42:29: 32000000 INFO @ Fri, 10 Dec 2021 13:42:34: 27000000 INFO @ Fri, 10 Dec 2021 13:42:35: 36000000 INFO @ Fri, 10 Dec 2021 13:42:36: 33000000 INFO @ Fri, 10 Dec 2021 13:42:41: 28000000 INFO @ Fri, 10 Dec 2021 13:42:42: 34000000 INFO @ Fri, 10 Dec 2021 13:42:42: 37000000 INFO @ Fri, 10 Dec 2021 13:42:48: 29000000 INFO @ Fri, 10 Dec 2021 13:42:49: 35000000 INFO @ Fri, 10 Dec 2021 13:42:50: 38000000 INFO @ Fri, 10 Dec 2021 13:42:56: 30000000 INFO @ Fri, 10 Dec 2021 13:42:56: 36000000 INFO @ Fri, 10 Dec 2021 13:42:58: 39000000 INFO @ Fri, 10 Dec 2021 13:43:03: 37000000 INFO @ Fri, 10 Dec 2021 13:43:03: 31000000 INFO @ Fri, 10 Dec 2021 13:43:05: 40000000 INFO @ Fri, 10 Dec 2021 13:43:09: 38000000 INFO @ Fri, 10 Dec 2021 13:43:10: 32000000 INFO @ Fri, 10 Dec 2021 13:43:12: 41000000 INFO @ Fri, 10 Dec 2021 13:43:16: 39000000 INFO @ Fri, 10 Dec 2021 13:43:17: 33000000 INFO @ Fri, 10 Dec 2021 13:43:20: 42000000 INFO @ Fri, 10 Dec 2021 13:43:23: 40000000 INFO @ Fri, 10 Dec 2021 13:43:24: 34000000 INFO @ Fri, 10 Dec 2021 13:43:28: 43000000 INFO @ Fri, 10 Dec 2021 13:43:30: 41000000 INFO @ Fri, 10 Dec 2021 13:43:30: 35000000 INFO @ Fri, 10 Dec 2021 13:43:36: 44000000 INFO @ Fri, 10 Dec 2021 13:43:36: 42000000 INFO @ Fri, 10 Dec 2021 13:43:38: 36000000 INFO @ Fri, 10 Dec 2021 13:43:42: 43000000 INFO @ Fri, 10 Dec 2021 13:43:43: 45000000 INFO @ Fri, 10 Dec 2021 13:43:45: 37000000 INFO @ Fri, 10 Dec 2021 13:43:50: 44000000 INFO @ Fri, 10 Dec 2021 13:43:51: 46000000 INFO @ Fri, 10 Dec 2021 13:43:53: 38000000 INFO @ Fri, 10 Dec 2021 13:43:56: 45000000 INFO @ Fri, 10 Dec 2021 13:43:58: 47000000 INFO @ Fri, 10 Dec 2021 13:44:00: 39000000 INFO @ Fri, 10 Dec 2021 13:44:02: 46000000 INFO @ Fri, 10 Dec 2021 13:44:06: 48000000 INFO @ Fri, 10 Dec 2021 13:44:08: 40000000 INFO @ Fri, 10 Dec 2021 13:44:09: 47000000 INFO @ Fri, 10 Dec 2021 13:44:14: 49000000 INFO @ Fri, 10 Dec 2021 13:44:15: 48000000 INFO @ Fri, 10 Dec 2021 13:44:15: 41000000 INFO @ Fri, 10 Dec 2021 13:44:21: 50000000 INFO @ Fri, 10 Dec 2021 13:44:22: 49000000 INFO @ Fri, 10 Dec 2021 13:44:23: 42000000 INFO @ Fri, 10 Dec 2021 13:44:28: 50000000 INFO @ Fri, 10 Dec 2021 13:44:29: 51000000 INFO @ Fri, 10 Dec 2021 13:44:31: 43000000 INFO @ Fri, 10 Dec 2021 13:44:35: 51000000 INFO @ Fri, 10 Dec 2021 13:44:37: 52000000 INFO @ Fri, 10 Dec 2021 13:44:38: 44000000 INFO @ Fri, 10 Dec 2021 13:44:42: 52000000 INFO @ Fri, 10 Dec 2021 13:44:45: 53000000 INFO @ Fri, 10 Dec 2021 13:44:46: 45000000 INFO @ Fri, 10 Dec 2021 13:44:50: 53000000 INFO @ Fri, 10 Dec 2021 13:44:53: 54000000 INFO @ Fri, 10 Dec 2021 13:44:53: 46000000 INFO @ Fri, 10 Dec 2021 13:44:57: 54000000 INFO @ Fri, 10 Dec 2021 13:45:00: 55000000 INFO @ Fri, 10 Dec 2021 13:45:01: 47000000 INFO @ Fri, 10 Dec 2021 13:45:03: 55000000 INFO @ Fri, 10 Dec 2021 13:45:08: 56000000 INFO @ Fri, 10 Dec 2021 13:45:08: 48000000 INFO @ Fri, 10 Dec 2021 13:45:10: 56000000 INFO @ Fri, 10 Dec 2021 13:45:15: 57000000 INFO @ Fri, 10 Dec 2021 13:45:16: 49000000 INFO @ Fri, 10 Dec 2021 13:45:16: 57000000 INFO @ Fri, 10 Dec 2021 13:45:23: 58000000 INFO @ Fri, 10 Dec 2021 13:45:23: 58000000 INFO @ Fri, 10 Dec 2021 13:45:23: 50000000 INFO @ Fri, 10 Dec 2021 13:45:30: 59000000 INFO @ Fri, 10 Dec 2021 13:45:31: 59000000 INFO @ Fri, 10 Dec 2021 13:45:31: 51000000 INFO @ Fri, 10 Dec 2021 13:45:37: 60000000 INFO @ Fri, 10 Dec 2021 13:45:37: 60000000 INFO @ Fri, 10 Dec 2021 13:45:40: 52000000 INFO @ Fri, 10 Dec 2021 13:45:44: 61000000 INFO @ Fri, 10 Dec 2021 13:45:46: 61000000 INFO @ Fri, 10 Dec 2021 13:45:47: 53000000 INFO @ Fri, 10 Dec 2021 13:45:51: 62000000 INFO @ Fri, 10 Dec 2021 13:45:54: 54000000 INFO @ Fri, 10 Dec 2021 13:45:54: 62000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Fri, 10 Dec 2021 13:45:58: 63000000 INFO @ Fri, 10 Dec 2021 13:46:01: 55000000 INFO @ Fri, 10 Dec 2021 13:46:02: 63000000 INFO @ Fri, 10 Dec 2021 13:46:05: 64000000 INFO @ Fri, 10 Dec 2021 13:46:09: 56000000 INFO @ Fri, 10 Dec 2021 13:46:09: 64000000 INFO @ Fri, 10 Dec 2021 13:46:12: 65000000 INFO @ Fri, 10 Dec 2021 13:46:17: 57000000 INFO @ Fri, 10 Dec 2021 13:46:17: 65000000 INFO @ Fri, 10 Dec 2021 13:46:18: 66000000 INFO @ Fri, 10 Dec 2021 13:46:24: 58000000 INFO @ Fri, 10 Dec 2021 13:46:25: 67000000 INFO @ Fri, 10 Dec 2021 13:46:25: 66000000 INFO @ Fri, 10 Dec 2021 13:46:32: 68000000 INFO @ Fri, 10 Dec 2021 13:46:32: 59000000 INFO @ Fri, 10 Dec 2021 13:46:32: 67000000 INFO @ Fri, 10 Dec 2021 13:46:39: 60000000 INFO @ Fri, 10 Dec 2021 13:46:40: 69000000 INFO @ Fri, 10 Dec 2021 13:46:41: 68000000 INFO @ Fri, 10 Dec 2021 13:46:47: 61000000 INFO @ Fri, 10 Dec 2021 13:46:47: 70000000 INFO @ Fri, 10 Dec 2021 13:46:49: 69000000 INFO @ Fri, 10 Dec 2021 13:46:55: 62000000 INFO @ Fri, 10 Dec 2021 13:46:56: 71000000 INFO @ Fri, 10 Dec 2021 13:46:57: 70000000 INFO @ Fri, 10 Dec 2021 13:47:02: 63000000 INFO @ Fri, 10 Dec 2021 13:47:03: 72000000 INFO @ Fri, 10 Dec 2021 13:47:06: 71000000 INFO @ Fri, 10 Dec 2021 13:47:10: 73000000 INFO @ Fri, 10 Dec 2021 13:47:13: 72000000 INFO @ Fri, 10 Dec 2021 13:47:13: 64000000 INFO @ Fri, 10 Dec 2021 13:47:17: 74000000 INFO @ Fri, 10 Dec 2021 13:47:20: 73000000 INFO @ Fri, 10 Dec 2021 13:47:21: 65000000 INFO @ Fri, 10 Dec 2021 13:47:23: 75000000 INFO @ Fri, 10 Dec 2021 13:47:27: 74000000 BigWig に変換しました。 INFO @ Fri, 10 Dec 2021 13:47:28: 66000000 INFO @ Fri, 10 Dec 2021 13:47:30: 76000000 INFO @ Fri, 10 Dec 2021 13:47:34: 75000000 INFO @ Fri, 10 Dec 2021 13:47:37: 67000000 INFO @ Fri, 10 Dec 2021 13:47:37: 77000000 INFO @ Fri, 10 Dec 2021 13:47:43: 76000000 INFO @ Fri, 10 Dec 2021 13:47:45: 78000000 INFO @ Fri, 10 Dec 2021 13:47:46: 68000000 INFO @ Fri, 10 Dec 2021 13:47:50: 77000000 INFO @ Fri, 10 Dec 2021 13:47:54: 79000000 INFO @ Fri, 10 Dec 2021 13:47:55: 69000000 INFO @ Fri, 10 Dec 2021 13:47:58: 78000000 INFO @ Fri, 10 Dec 2021 13:48:01: 80000000 INFO @ Fri, 10 Dec 2021 13:48:02: 70000000 INFO @ Fri, 10 Dec 2021 13:48:06: 79000000 INFO @ Fri, 10 Dec 2021 13:48:08: 81000000 INFO @ Fri, 10 Dec 2021 13:48:11: 71000000 INFO @ Fri, 10 Dec 2021 13:48:13: 80000000 INFO @ Fri, 10 Dec 2021 13:48:15: 82000000 INFO @ Fri, 10 Dec 2021 13:48:18: 72000000 INFO @ Fri, 10 Dec 2021 13:48:20: 81000000 INFO @ Fri, 10 Dec 2021 13:48:22: 83000000 INFO @ Fri, 10 Dec 2021 13:48:25: 73000000 INFO @ Fri, 10 Dec 2021 13:48:27: 82000000 INFO @ Fri, 10 Dec 2021 13:48:29: 84000000 INFO @ Fri, 10 Dec 2021 13:48:33: 74000000 INFO @ Fri, 10 Dec 2021 13:48:35: 83000000 INFO @ Fri, 10 Dec 2021 13:48:35: 85000000 INFO @ Fri, 10 Dec 2021 13:48:40: 75000000 INFO @ Fri, 10 Dec 2021 13:48:42: 84000000 INFO @ Fri, 10 Dec 2021 13:48:43: 86000000 INFO @ Fri, 10 Dec 2021 13:48:47: 76000000 INFO @ Fri, 10 Dec 2021 13:48:50: 87000000 INFO @ Fri, 10 Dec 2021 13:48:51: 85000000 INFO @ Fri, 10 Dec 2021 13:48:55: 77000000 INFO @ Fri, 10 Dec 2021 13:48:56: #1 tag size is determined as 76 bps INFO @ Fri, 10 Dec 2021 13:48:56: #1 tag size = 76 INFO @ Fri, 10 Dec 2021 13:48:56: #1 total tags in treatment: 35469307 INFO @ Fri, 10 Dec 2021 13:48:56: #1 user defined the maximum tags... INFO @ Fri, 10 Dec 2021 13:48:56: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 10 Dec 2021 13:48:57: #1 tags after filtering in treatment: 22696791 INFO @ Fri, 10 Dec 2021 13:48:57: #1 Redundant rate of treatment: 0.36 INFO @ Fri, 10 Dec 2021 13:48:57: #1 finished! INFO @ Fri, 10 Dec 2021 13:48:57: #2 Build Peak Model... INFO @ Fri, 10 Dec 2021 13:48:57: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 10 Dec 2021 13:48:59: 86000000 INFO @ Fri, 10 Dec 2021 13:48:59: #2 number of paired peaks: 5335 INFO @ Fri, 10 Dec 2021 13:48:59: start model_add_line... INFO @ Fri, 10 Dec 2021 13:48:59: start X-correlation... INFO @ Fri, 10 Dec 2021 13:49:00: end of X-cor INFO @ Fri, 10 Dec 2021 13:49:00: #2 finished! INFO @ Fri, 10 Dec 2021 13:49:00: #2 predicted fragment length is 96 bps INFO @ Fri, 10 Dec 2021 13:49:00: #2 alternative fragment length(s) may be 4,96 bps INFO @ Fri, 10 Dec 2021 13:49:00: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10_model.r WARNING @ Fri, 10 Dec 2021 13:49:00: #2 Since the d (96) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Fri, 10 Dec 2021 13:49:00: #2 You may need to consider one of the other alternative d(s): 4,96 WARNING @ Fri, 10 Dec 2021 13:49:00: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Fri, 10 Dec 2021 13:49:00: #3 Call peaks... INFO @ Fri, 10 Dec 2021 13:49:00: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 10 Dec 2021 13:49:02: 78000000 INFO @ Fri, 10 Dec 2021 13:49:07: 87000000 INFO @ Fri, 10 Dec 2021 13:49:10: 79000000 INFO @ Fri, 10 Dec 2021 13:49:14: #1 tag size is determined as 76 bps INFO @ Fri, 10 Dec 2021 13:49:14: #1 tag size = 76 INFO @ Fri, 10 Dec 2021 13:49:14: #1 total tags in treatment: 35469307 INFO @ Fri, 10 Dec 2021 13:49:14: #1 user defined the maximum tags... INFO @ Fri, 10 Dec 2021 13:49:14: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 10 Dec 2021 13:49:15: #1 tags after filtering in treatment: 22696791 INFO @ Fri, 10 Dec 2021 13:49:15: #1 Redundant rate of treatment: 0.36 INFO @ Fri, 10 Dec 2021 13:49:15: #1 finished! INFO @ Fri, 10 Dec 2021 13:49:15: #2 Build Peak Model... INFO @ Fri, 10 Dec 2021 13:49:15: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 10 Dec 2021 13:49:17: 80000000 INFO @ Fri, 10 Dec 2021 13:49:17: #2 number of paired peaks: 5335 INFO @ Fri, 10 Dec 2021 13:49:17: start model_add_line... INFO @ Fri, 10 Dec 2021 13:49:17: start X-correlation... INFO @ Fri, 10 Dec 2021 13:49:18: end of X-cor INFO @ Fri, 10 Dec 2021 13:49:18: #2 finished! INFO @ Fri, 10 Dec 2021 13:49:18: #2 predicted fragment length is 96 bps INFO @ Fri, 10 Dec 2021 13:49:18: #2 alternative fragment length(s) may be 4,96 bps INFO @ Fri, 10 Dec 2021 13:49:18: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05_model.r WARNING @ Fri, 10 Dec 2021 13:49:18: #2 Since the d (96) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Fri, 10 Dec 2021 13:49:18: #2 You may need to consider one of the other alternative d(s): 4,96 WARNING @ Fri, 10 Dec 2021 13:49:18: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Fri, 10 Dec 2021 13:49:18: #3 Call peaks... INFO @ Fri, 10 Dec 2021 13:49:18: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 10 Dec 2021 13:49:23: 81000000 INFO @ Fri, 10 Dec 2021 13:49:30: 82000000 INFO @ Fri, 10 Dec 2021 13:49:37: 83000000 INFO @ Fri, 10 Dec 2021 13:49:43: 84000000 INFO @ Fri, 10 Dec 2021 13:49:50: 85000000 INFO @ Fri, 10 Dec 2021 13:49:53: #3 Call peaks for each chromosome... INFO @ Fri, 10 Dec 2021 13:49:57: 86000000 INFO @ Fri, 10 Dec 2021 13:50:04: 87000000 INFO @ Fri, 10 Dec 2021 13:50:10: #1 tag size is determined as 76 bps INFO @ Fri, 10 Dec 2021 13:50:10: #1 tag size = 76 INFO @ Fri, 10 Dec 2021 13:50:10: #1 total tags in treatment: 35469307 INFO @ Fri, 10 Dec 2021 13:50:10: #1 user defined the maximum tags... INFO @ Fri, 10 Dec 2021 13:50:10: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 10 Dec 2021 13:50:11: #1 tags after filtering in treatment: 22696791 INFO @ Fri, 10 Dec 2021 13:50:11: #1 Redundant rate of treatment: 0.36 INFO @ Fri, 10 Dec 2021 13:50:11: #1 finished! INFO @ Fri, 10 Dec 2021 13:50:11: #2 Build Peak Model... INFO @ Fri, 10 Dec 2021 13:50:11: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 10 Dec 2021 13:50:11: #3 Call peaks for each chromosome... INFO @ Fri, 10 Dec 2021 13:50:13: #2 number of paired peaks: 5335 INFO @ Fri, 10 Dec 2021 13:50:13: start model_add_line... INFO @ Fri, 10 Dec 2021 13:50:14: start X-correlation... INFO @ Fri, 10 Dec 2021 13:50:14: end of X-cor INFO @ Fri, 10 Dec 2021 13:50:14: #2 finished! INFO @ Fri, 10 Dec 2021 13:50:14: #2 predicted fragment length is 96 bps INFO @ Fri, 10 Dec 2021 13:50:14: #2 alternative fragment length(s) may be 4,96 bps INFO @ Fri, 10 Dec 2021 13:50:14: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20_model.r WARNING @ Fri, 10 Dec 2021 13:50:14: #2 Since the d (96) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Fri, 10 Dec 2021 13:50:14: #2 You may need to consider one of the other alternative d(s): 4,96 WARNING @ Fri, 10 Dec 2021 13:50:14: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Fri, 10 Dec 2021 13:50:14: #3 Call peaks... INFO @ Fri, 10 Dec 2021 13:50:14: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 10 Dec 2021 13:50:19: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10_peaks.xls INFO @ Fri, 10 Dec 2021 13:50:19: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10_peaks.narrowPeak INFO @ Fri, 10 Dec 2021 13:50:19: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.10_summits.bed INFO @ Fri, 10 Dec 2021 13:50:19: Done! pass1 - making usageList (15 chroms): 2 millis pass2 - checking and writing primary data (10056 records, 4 fields): 15 millis CompletedMACS2peakCalling INFO @ Fri, 10 Dec 2021 13:50:36: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05_peaks.xls INFO @ Fri, 10 Dec 2021 13:50:36: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05_peaks.narrowPeak INFO @ Fri, 10 Dec 2021 13:50:36: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.05_summits.bed INFO @ Fri, 10 Dec 2021 13:50:36: Done! pass1 - making usageList (15 chroms): 16 millis pass2 - checking and writing primary data (15358 records, 4 fields): 20 millis CompletedMACS2peakCalling INFO @ Fri, 10 Dec 2021 13:51:06: #3 Call peaks for each chromosome... INFO @ Fri, 10 Dec 2021 13:51:30: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20_peaks.xls INFO @ Fri, 10 Dec 2021 13:51:30: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20_peaks.narrowPeak INFO @ Fri, 10 Dec 2021 13:51:30: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX4712935/SRX4712935.20_summits.bed INFO @ Fri, 10 Dec 2021 13:51:30: Done! pass1 - making usageList (15 chroms): 2 millis pass2 - checking and writing primary data (5755 records, 4 fields): 10 millis CompletedMACS2peakCalling