Job ID = 14167453 SRX = ERX3978911 Genome = dm3 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... Read 110150411 spots for ERR3975980/ERR3975980.sra Written 110150411 spots for ERR3975980/ERR3975980.sra Read 55217729 spots for ERR3976020/ERR3976020.sra Written 55217729 spots for ERR3976020/ERR3976020.sra fastq に変換しました。 bowtie でマッピング中... Your job 14171265 ("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: 20:06:30 165368140 reads; of these: 165368140 (100.00%) were paired; of these: 51658564 (31.24%) aligned concordantly 0 times 87072011 (52.65%) aligned concordantly exactly 1 time 26637565 (16.11%) aligned concordantly >1 times ---- 51658564 pairs aligned concordantly 0 times; of these: 4306932 (8.34%) aligned discordantly 1 time ---- 47351632 pairs aligned 0 times concordantly or discordantly; of these: 94703264 mates make up the pairs; of these: 29878397 (31.55%) aligned 0 times 30887127 (32.61%) aligned exactly 1 time 33937740 (35.84%) aligned >1 times 90.97% overall alignment rate Time searching: 20:06:31 Overall time: 20:06:31 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 124 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] 106052580 / 116881887 = 0.9073 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 @ Sat, 11 Dec 2021 09:53:04: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.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, 11 Dec 2021 09:53:04: #1 read tag files... INFO @ Sat, 11 Dec 2021 09:53:04: #1 read treatment tags... INFO @ Sat, 11 Dec 2021 09:53:10: 1000000 INFO @ Sat, 11 Dec 2021 09:53:16: 2000000 INFO @ Sat, 11 Dec 2021 09:53:21: 3000000 INFO @ Sat, 11 Dec 2021 09:53:27: 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 @ Sat, 11 Dec 2021 09:53:33: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.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, 11 Dec 2021 09:53:33: #1 read tag files... INFO @ Sat, 11 Dec 2021 09:53:33: #1 read treatment tags... INFO @ Sat, 11 Dec 2021 09:53:34: 5000000 INFO @ Sat, 11 Dec 2021 09:53:41: 6000000 INFO @ Sat, 11 Dec 2021 09:53:42: 1000000 INFO @ Sat, 11 Dec 2021 09:53:48: 7000000 INFO @ Sat, 11 Dec 2021 09:53:50: 2000000 INFO @ Sat, 11 Dec 2021 09:53:55: 8000000 INFO @ Sat, 11 Dec 2021 09:53:59: 3000000 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 @ Sat, 11 Dec 2021 09:54:03: 9000000 INFO @ Sat, 11 Dec 2021 09:54:03: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.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, 11 Dec 2021 09:54:03: #1 read tag files... INFO @ Sat, 11 Dec 2021 09:54:03: #1 read treatment tags... INFO @ Sat, 11 Dec 2021 09:54:07: 4000000 INFO @ Sat, 11 Dec 2021 09:54:10: 10000000 INFO @ Sat, 11 Dec 2021 09:54:12: 1000000 INFO @ Sat, 11 Dec 2021 09:54:16: 5000000 INFO @ Sat, 11 Dec 2021 09:54:17: 11000000 INFO @ Sat, 11 Dec 2021 09:54:20: 2000000 INFO @ Sat, 11 Dec 2021 09:54:25: 6000000 INFO @ Sat, 11 Dec 2021 09:54:25: 12000000 INFO @ Sat, 11 Dec 2021 09:54:29: 3000000 INFO @ Sat, 11 Dec 2021 09:54:32: 13000000 INFO @ Sat, 11 Dec 2021 09:54:33: 7000000 INFO @ Sat, 11 Dec 2021 09:54:38: 4000000 INFO @ Sat, 11 Dec 2021 09:54:40: 14000000 INFO @ Sat, 11 Dec 2021 09:54:42: 8000000 INFO @ Sat, 11 Dec 2021 09:54:47: 5000000 INFO @ Sat, 11 Dec 2021 09:54:47: 15000000 INFO @ Sat, 11 Dec 2021 09:54:51: 9000000 INFO @ Sat, 11 Dec 2021 09:54:55: 16000000 INFO @ Sat, 11 Dec 2021 09:54:55: 6000000 INFO @ Sat, 11 Dec 2021 09:55:00: 10000000 INFO @ Sat, 11 Dec 2021 09:55:03: 17000000 INFO @ Sat, 11 Dec 2021 09:55:04: 7000000 INFO @ Sat, 11 Dec 2021 09:55:08: 11000000 INFO @ Sat, 11 Dec 2021 09:55:10: 18000000 INFO @ Sat, 11 Dec 2021 09:55:13: 8000000 INFO @ Sat, 11 Dec 2021 09:55:17: 12000000 INFO @ Sat, 11 Dec 2021 09:55:18: 19000000 INFO @ Sat, 11 Dec 2021 09:55:22: 9000000 INFO @ Sat, 11 Dec 2021 09:55:26: 13000000 INFO @ Sat, 11 Dec 2021 09:55:26: 20000000 INFO @ Sat, 11 Dec 2021 09:55:31: 10000000 INFO @ Sat, 11 Dec 2021 09:55:34: 21000000 INFO @ Sat, 11 Dec 2021 09:55:34: 14000000 INFO @ Sat, 11 Dec 2021 09:55:40: 11000000 INFO @ Sat, 11 Dec 2021 09:55:42: 22000000 INFO @ Sat, 11 Dec 2021 09:55:43: 15000000 INFO @ Sat, 11 Dec 2021 09:55:49: 12000000 INFO @ Sat, 11 Dec 2021 09:55:50: 23000000 INFO @ Sat, 11 Dec 2021 09:55:52: 16000000 INFO @ Sat, 11 Dec 2021 09:55:57: 24000000 INFO @ Sat, 11 Dec 2021 09:55:58: 13000000 INFO @ Sat, 11 Dec 2021 09:56:01: 17000000 INFO @ Sat, 11 Dec 2021 09:56:05: 25000000 INFO @ Sat, 11 Dec 2021 09:56:07: 14000000 INFO @ Sat, 11 Dec 2021 09:56:09: 18000000 INFO @ Sat, 11 Dec 2021 09:56:13: 26000000 INFO @ Sat, 11 Dec 2021 09:56:16: 15000000 INFO @ Sat, 11 Dec 2021 09:56:18: 19000000 INFO @ Sat, 11 Dec 2021 09:56:21: 27000000 INFO @ Sat, 11 Dec 2021 09:56:25: 16000000 INFO @ Sat, 11 Dec 2021 09:56:27: 20000000 INFO @ Sat, 11 Dec 2021 09:56:29: 28000000 INFO @ Sat, 11 Dec 2021 09:56:34: 17000000 INFO @ Sat, 11 Dec 2021 09:56:36: 21000000 INFO @ Sat, 11 Dec 2021 09:56:37: 29000000 INFO @ Sat, 11 Dec 2021 09:56:43: 18000000 INFO @ Sat, 11 Dec 2021 09:56:45: 22000000 INFO @ Sat, 11 Dec 2021 09:56:45: 30000000 INFO @ Sat, 11 Dec 2021 09:56:52: 19000000 INFO @ Sat, 11 Dec 2021 09:56:53: 31000000 INFO @ Sat, 11 Dec 2021 09:56:54: 23000000 INFO @ Sat, 11 Dec 2021 09:57:00: 32000000 INFO @ Sat, 11 Dec 2021 09:57:01: 20000000 INFO @ Sat, 11 Dec 2021 09:57:03: 24000000 INFO @ Sat, 11 Dec 2021 09:57:08: 33000000 INFO @ Sat, 11 Dec 2021 09:57:10: 21000000 INFO @ Sat, 11 Dec 2021 09:57:12: 25000000 INFO @ Sat, 11 Dec 2021 09:57:16: 34000000 INFO @ Sat, 11 Dec 2021 09:57:19: 22000000 INFO @ Sat, 11 Dec 2021 09:57:20: 26000000 INFO @ Sat, 11 Dec 2021 09:57:24: 35000000 INFO @ Sat, 11 Dec 2021 09:57:28: 23000000 INFO @ Sat, 11 Dec 2021 09:57:29: 27000000 INFO @ Sat, 11 Dec 2021 09:57:32: 36000000 INFO @ Sat, 11 Dec 2021 09:57:37: 24000000 INFO @ Sat, 11 Dec 2021 09:57:38: 28000000 INFO @ Sat, 11 Dec 2021 09:57:40: 37000000 INFO @ Sat, 11 Dec 2021 09:57:46: 25000000 INFO @ Sat, 11 Dec 2021 09:57:47: 29000000 INFO @ Sat, 11 Dec 2021 09:57:48: 38000000 INFO @ Sat, 11 Dec 2021 09:57:55: 26000000 INFO @ Sat, 11 Dec 2021 09:57:55: 30000000 INFO @ Sat, 11 Dec 2021 09:57:56: 39000000 INFO @ Sat, 11 Dec 2021 09:58:04: 40000000 INFO @ Sat, 11 Dec 2021 09:58:04: 27000000 INFO @ Sat, 11 Dec 2021 09:58:04: 31000000 INFO @ Sat, 11 Dec 2021 09:58:12: 41000000 INFO @ Sat, 11 Dec 2021 09:58:13: 32000000 INFO @ Sat, 11 Dec 2021 09:58:13: 28000000 INFO @ Sat, 11 Dec 2021 09:58:20: 42000000 INFO @ Sat, 11 Dec 2021 09:58:22: 33000000 INFO @ Sat, 11 Dec 2021 09:58:22: 29000000 INFO @ Sat, 11 Dec 2021 09:58:28: 43000000 INFO @ Sat, 11 Dec 2021 09:58:30: 34000000 INFO @ Sat, 11 Dec 2021 09:58:31: 30000000 INFO @ Sat, 11 Dec 2021 09:58:36: 44000000 INFO @ Sat, 11 Dec 2021 09:58:39: 35000000 INFO @ Sat, 11 Dec 2021 09:58:40: 31000000 INFO @ Sat, 11 Dec 2021 09:58:44: 45000000 INFO @ Sat, 11 Dec 2021 09:58:48: 36000000 INFO @ Sat, 11 Dec 2021 09:58:49: 32000000 BedGraph に変換しました。 BigWig に変換中... INFO @ Sat, 11 Dec 2021 09:58:52: 46000000 INFO @ Sat, 11 Dec 2021 09:58:57: 33000000 INFO @ Sat, 11 Dec 2021 09:58:57: 37000000 INFO @ Sat, 11 Dec 2021 09:58:59: 47000000 INFO @ Sat, 11 Dec 2021 09:59:06: 34000000 INFO @ Sat, 11 Dec 2021 09:59:07: 38000000 INFO @ Sat, 11 Dec 2021 09:59:07: 48000000 INFO @ Sat, 11 Dec 2021 09:59:14: 35000000 INFO @ Sat, 11 Dec 2021 09:59:15: 49000000 INFO @ Sat, 11 Dec 2021 09:59:16: 39000000 INFO @ Sat, 11 Dec 2021 09:59:23: 36000000 INFO @ Sat, 11 Dec 2021 09:59:23: 50000000 INFO @ Sat, 11 Dec 2021 09:59:25: 40000000 BigWig に変換しました。 INFO @ Sat, 11 Dec 2021 09:59:31: 37000000 INFO @ Sat, 11 Dec 2021 09:59:32: 51000000 INFO @ Sat, 11 Dec 2021 09:59:34: 41000000 INFO @ Sat, 11 Dec 2021 09:59:39: 38000000 INFO @ Sat, 11 Dec 2021 09:59:40: 52000000 INFO @ Sat, 11 Dec 2021 09:59:44: 42000000 INFO @ Sat, 11 Dec 2021 09:59:48: 39000000 INFO @ Sat, 11 Dec 2021 09:59:48: 53000000 INFO @ Sat, 11 Dec 2021 09:59:53: 43000000 INFO @ Sat, 11 Dec 2021 09:59:56: 40000000 INFO @ Sat, 11 Dec 2021 09:59:57: 54000000 INFO @ Sat, 11 Dec 2021 10:00:02: 44000000 INFO @ Sat, 11 Dec 2021 10:00:05: 41000000 INFO @ Sat, 11 Dec 2021 10:00:05: 55000000 INFO @ Sat, 11 Dec 2021 10:00:11: 45000000 INFO @ Sat, 11 Dec 2021 10:00:13: 56000000 INFO @ Sat, 11 Dec 2021 10:00:13: 42000000 INFO @ Sat, 11 Dec 2021 10:00:20: 46000000 INFO @ Sat, 11 Dec 2021 10:00:21: 57000000 INFO @ Sat, 11 Dec 2021 10:00:22: 43000000 INFO @ Sat, 11 Dec 2021 10:00:29: 47000000 INFO @ Sat, 11 Dec 2021 10:00:30: 58000000 INFO @ Sat, 11 Dec 2021 10:00:31: 44000000 INFO @ Sat, 11 Dec 2021 10:00:38: 48000000 INFO @ Sat, 11 Dec 2021 10:00:38: 59000000 INFO @ Sat, 11 Dec 2021 10:00:39: 45000000 INFO @ Sat, 11 Dec 2021 10:00:47: 49000000 INFO @ Sat, 11 Dec 2021 10:00:47: 60000000 INFO @ Sat, 11 Dec 2021 10:00:48: 46000000 INFO @ Sat, 11 Dec 2021 10:00:55: 61000000 INFO @ Sat, 11 Dec 2021 10:00:56: 50000000 INFO @ Sat, 11 Dec 2021 10:00:57: 47000000 INFO @ Sat, 11 Dec 2021 10:01:03: 62000000 INFO @ Sat, 11 Dec 2021 10:01:06: 51000000 INFO @ Sat, 11 Dec 2021 10:01:06: 48000000 INFO @ Sat, 11 Dec 2021 10:01:12: 63000000 INFO @ Sat, 11 Dec 2021 10:01:15: 52000000 INFO @ Sat, 11 Dec 2021 10:01:15: 49000000 INFO @ Sat, 11 Dec 2021 10:01:20: 64000000 INFO @ Sat, 11 Dec 2021 10:01:24: 50000000 INFO @ Sat, 11 Dec 2021 10:01:25: 53000000 INFO @ Sat, 11 Dec 2021 10:01:30: 65000000 INFO @ Sat, 11 Dec 2021 10:01:34: 51000000 INFO @ Sat, 11 Dec 2021 10:01:35: 54000000 INFO @ Sat, 11 Dec 2021 10:01:39: 66000000 INFO @ Sat, 11 Dec 2021 10:01:43: 52000000 INFO @ Sat, 11 Dec 2021 10:01:44: 55000000 INFO @ Sat, 11 Dec 2021 10:01:48: 67000000 INFO @ Sat, 11 Dec 2021 10:01:53: 53000000 INFO @ Sat, 11 Dec 2021 10:01:54: 56000000 INFO @ Sat, 11 Dec 2021 10:01:57: 68000000 INFO @ Sat, 11 Dec 2021 10:02:02: 54000000 INFO @ Sat, 11 Dec 2021 10:02:03: 57000000 INFO @ Sat, 11 Dec 2021 10:02:05: 69000000 INFO @ Sat, 11 Dec 2021 10:02:11: 55000000 INFO @ Sat, 11 Dec 2021 10:02:13: 58000000 INFO @ Sat, 11 Dec 2021 10:02:14: 70000000 INFO @ Sat, 11 Dec 2021 10:02:21: 56000000 INFO @ Sat, 11 Dec 2021 10:02:22: 71000000 INFO @ Sat, 11 Dec 2021 10:02:22: 59000000 INFO @ Sat, 11 Dec 2021 10:02:30: 72000000 INFO @ Sat, 11 Dec 2021 10:02:30: 57000000 INFO @ Sat, 11 Dec 2021 10:02:32: 60000000 INFO @ Sat, 11 Dec 2021 10:02:38: 73000000 INFO @ Sat, 11 Dec 2021 10:02:39: 58000000 INFO @ Sat, 11 Dec 2021 10:02:41: 61000000 INFO @ Sat, 11 Dec 2021 10:02:47: 74000000 INFO @ Sat, 11 Dec 2021 10:02:48: 59000000 INFO @ Sat, 11 Dec 2021 10:02:51: 62000000 INFO @ Sat, 11 Dec 2021 10:02:56: 75000000 INFO @ Sat, 11 Dec 2021 10:02:58: 60000000 INFO @ Sat, 11 Dec 2021 10:03:01: 63000000 INFO @ Sat, 11 Dec 2021 10:03:06: 76000000 INFO @ Sat, 11 Dec 2021 10:03:07: 61000000 INFO @ Sat, 11 Dec 2021 10:03:10: 64000000 INFO @ Sat, 11 Dec 2021 10:03:14: 77000000 INFO @ Sat, 11 Dec 2021 10:03:16: 62000000 INFO @ Sat, 11 Dec 2021 10:03:20: 65000000 INFO @ Sat, 11 Dec 2021 10:03:22: 78000000 INFO @ Sat, 11 Dec 2021 10:03:26: 63000000 INFO @ Sat, 11 Dec 2021 10:03:30: 66000000 INFO @ Sat, 11 Dec 2021 10:03:30: 79000000 INFO @ Sat, 11 Dec 2021 10:03:35: 64000000 INFO @ Sat, 11 Dec 2021 10:03:39: 80000000 INFO @ Sat, 11 Dec 2021 10:03:40: 67000000 INFO @ Sat, 11 Dec 2021 10:03:45: 65000000 INFO @ Sat, 11 Dec 2021 10:03:47: 81000000 INFO @ Sat, 11 Dec 2021 10:03:49: 68000000 INFO @ Sat, 11 Dec 2021 10:03:55: 66000000 INFO @ Sat, 11 Dec 2021 10:03:56: 82000000 INFO @ Sat, 11 Dec 2021 10:03:59: 69000000 INFO @ Sat, 11 Dec 2021 10:04:04: 83000000 INFO @ Sat, 11 Dec 2021 10:04:04: 67000000 INFO @ Sat, 11 Dec 2021 10:04:08: 70000000 INFO @ Sat, 11 Dec 2021 10:04:12: 84000000 INFO @ Sat, 11 Dec 2021 10:04:14: 68000000 INFO @ Sat, 11 Dec 2021 10:04:17: 71000000 INFO @ Sat, 11 Dec 2021 10:04:20: 85000000 INFO @ Sat, 11 Dec 2021 10:04:23: 69000000 INFO @ Sat, 11 Dec 2021 10:04:26: 72000000 INFO @ Sat, 11 Dec 2021 10:04:28: 86000000 INFO @ Sat, 11 Dec 2021 10:04:32: 70000000 INFO @ Sat, 11 Dec 2021 10:04:35: 73000000 INFO @ Sat, 11 Dec 2021 10:04:37: 87000000 INFO @ Sat, 11 Dec 2021 10:04:40: 71000000 INFO @ Sat, 11 Dec 2021 10:04:44: 74000000 INFO @ Sat, 11 Dec 2021 10:04:45: 88000000 INFO @ Sat, 11 Dec 2021 10:04:49: 72000000 INFO @ Sat, 11 Dec 2021 10:04:51: #1 tag size is determined as 80 bps INFO @ Sat, 11 Dec 2021 10:04:51: #1 tag size = 80 INFO @ Sat, 11 Dec 2021 10:04:51: #1 total tags in treatment: 10811656 INFO @ Sat, 11 Dec 2021 10:04:51: #1 user defined the maximum tags... INFO @ Sat, 11 Dec 2021 10:04:51: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 11 Dec 2021 10:04:51: #1 tags after filtering in treatment: 7600108 INFO @ Sat, 11 Dec 2021 10:04:51: #1 Redundant rate of treatment: 0.30 INFO @ Sat, 11 Dec 2021 10:04:51: #1 finished! INFO @ Sat, 11 Dec 2021 10:04:51: #2 Build Peak Model... INFO @ Sat, 11 Dec 2021 10:04:51: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 11 Dec 2021 10:04:52: #2 number of paired peaks: 6698 INFO @ Sat, 11 Dec 2021 10:04:52: start model_add_line... INFO @ Sat, 11 Dec 2021 10:04:52: start X-correlation... INFO @ Sat, 11 Dec 2021 10:04:52: end of X-cor INFO @ Sat, 11 Dec 2021 10:04:52: #2 finished! INFO @ Sat, 11 Dec 2021 10:04:52: #2 predicted fragment length is 196 bps INFO @ Sat, 11 Dec 2021 10:04:52: #2 alternative fragment length(s) may be 196 bps INFO @ Sat, 11 Dec 2021 10:04:52: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05_model.r INFO @ Sat, 11 Dec 2021 10:04:52: #3 Call peaks... INFO @ Sat, 11 Dec 2021 10:04:52: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 11 Dec 2021 10:04:53: 75000000 INFO @ Sat, 11 Dec 2021 10:04:58: 73000000 INFO @ Sat, 11 Dec 2021 10:05:02: 76000000 INFO @ Sat, 11 Dec 2021 10:05:06: 74000000 INFO @ Sat, 11 Dec 2021 10:05:11: 77000000 INFO @ Sat, 11 Dec 2021 10:05:13: #3 Call peaks for each chromosome... INFO @ Sat, 11 Dec 2021 10:05:15: 75000000 INFO @ Sat, 11 Dec 2021 10:05:20: 78000000 INFO @ Sat, 11 Dec 2021 10:05:22: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05_peaks.xls INFO @ Sat, 11 Dec 2021 10:05:22: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05_peaks.narrowPeak INFO @ Sat, 11 Dec 2021 10:05:22: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.05_summits.bed INFO @ Sat, 11 Dec 2021 10:05:22: Done! pass1 - making usageList (14 chroms): 2 millis pass2 - checking and writing primary data (12285 records, 4 fields): 12 millis CompletedMACS2peakCalling INFO @ Sat, 11 Dec 2021 10:05:23: 76000000 INFO @ Sat, 11 Dec 2021 10:05:29: 79000000 INFO @ Sat, 11 Dec 2021 10:05:32: 77000000 INFO @ Sat, 11 Dec 2021 10:05:38: 80000000 INFO @ Sat, 11 Dec 2021 10:05:40: 78000000 INFO @ Sat, 11 Dec 2021 10:05:48: 81000000 INFO @ Sat, 11 Dec 2021 10:05:49: 79000000 INFO @ Sat, 11 Dec 2021 10:05:57: 82000000 INFO @ Sat, 11 Dec 2021 10:05:58: 80000000 INFO @ Sat, 11 Dec 2021 10:06:06: 83000000 INFO @ Sat, 11 Dec 2021 10:06:07: 81000000 INFO @ Sat, 11 Dec 2021 10:06:15: 84000000 INFO @ Sat, 11 Dec 2021 10:06:17: 82000000 INFO @ Sat, 11 Dec 2021 10:06:24: 85000000 INFO @ Sat, 11 Dec 2021 10:06:25: 83000000 INFO @ Sat, 11 Dec 2021 10:06:32: 86000000 INFO @ Sat, 11 Dec 2021 10:06:34: 84000000 INFO @ Sat, 11 Dec 2021 10:06:41: 87000000 INFO @ Sat, 11 Dec 2021 10:06:42: 85000000 INFO @ Sat, 11 Dec 2021 10:06:50: 88000000 INFO @ Sat, 11 Dec 2021 10:06:51: 86000000 INFO @ Sat, 11 Dec 2021 10:06:56: #1 tag size is determined as 80 bps INFO @ Sat, 11 Dec 2021 10:06:56: #1 tag size = 80 INFO @ Sat, 11 Dec 2021 10:06:56: #1 total tags in treatment: 10811656 INFO @ Sat, 11 Dec 2021 10:06:56: #1 user defined the maximum tags... INFO @ Sat, 11 Dec 2021 10:06:56: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 11 Dec 2021 10:06:56: #1 tags after filtering in treatment: 7600108 INFO @ Sat, 11 Dec 2021 10:06:56: #1 Redundant rate of treatment: 0.30 INFO @ Sat, 11 Dec 2021 10:06:56: #1 finished! INFO @ Sat, 11 Dec 2021 10:06:56: #2 Build Peak Model... INFO @ Sat, 11 Dec 2021 10:06:56: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 11 Dec 2021 10:06:57: #2 number of paired peaks: 6698 INFO @ Sat, 11 Dec 2021 10:06:57: start model_add_line... INFO @ Sat, 11 Dec 2021 10:06:57: start X-correlation... INFO @ Sat, 11 Dec 2021 10:06:57: end of X-cor INFO @ Sat, 11 Dec 2021 10:06:57: #2 finished! INFO @ Sat, 11 Dec 2021 10:06:57: #2 predicted fragment length is 196 bps INFO @ Sat, 11 Dec 2021 10:06:57: #2 alternative fragment length(s) may be 196 bps INFO @ Sat, 11 Dec 2021 10:06:57: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10_model.r INFO @ Sat, 11 Dec 2021 10:06:57: #3 Call peaks... INFO @ Sat, 11 Dec 2021 10:06:57: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 11 Dec 2021 10:06:59: 87000000 INFO @ Sat, 11 Dec 2021 10:07:06: 88000000 INFO @ Sat, 11 Dec 2021 10:07:12: #1 tag size is determined as 80 bps INFO @ Sat, 11 Dec 2021 10:07:12: #1 tag size = 80 INFO @ Sat, 11 Dec 2021 10:07:12: #1 total tags in treatment: 10811656 INFO @ Sat, 11 Dec 2021 10:07:12: #1 user defined the maximum tags... INFO @ Sat, 11 Dec 2021 10:07:12: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Sat, 11 Dec 2021 10:07:12: #1 tags after filtering in treatment: 7600108 INFO @ Sat, 11 Dec 2021 10:07:12: #1 Redundant rate of treatment: 0.30 INFO @ Sat, 11 Dec 2021 10:07:12: #1 finished! INFO @ Sat, 11 Dec 2021 10:07:12: #2 Build Peak Model... INFO @ Sat, 11 Dec 2021 10:07:12: #2 looking for paired plus/minus strand peaks... INFO @ Sat, 11 Dec 2021 10:07:13: #2 number of paired peaks: 6698 INFO @ Sat, 11 Dec 2021 10:07:13: start model_add_line... INFO @ Sat, 11 Dec 2021 10:07:13: start X-correlation... INFO @ Sat, 11 Dec 2021 10:07:13: end of X-cor INFO @ Sat, 11 Dec 2021 10:07:13: #2 finished! INFO @ Sat, 11 Dec 2021 10:07:13: #2 predicted fragment length is 196 bps INFO @ Sat, 11 Dec 2021 10:07:13: #2 alternative fragment length(s) may be 196 bps INFO @ Sat, 11 Dec 2021 10:07:13: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20_model.r INFO @ Sat, 11 Dec 2021 10:07:13: #3 Call peaks... INFO @ Sat, 11 Dec 2021 10:07:13: #3 Pre-compute pvalue-qvalue table... INFO @ Sat, 11 Dec 2021 10:07:17: #3 Call peaks for each chromosome... INFO @ Sat, 11 Dec 2021 10:07:26: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10_peaks.xls INFO @ Sat, 11 Dec 2021 10:07:26: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10_peaks.narrowPeak INFO @ Sat, 11 Dec 2021 10:07:26: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.10_summits.bed INFO @ Sat, 11 Dec 2021 10:07:26: Done! pass1 - making usageList (14 chroms): 2 millis pass2 - checking and writing primary data (7969 records, 4 fields): 10 millis CompletedMACS2peakCalling INFO @ Sat, 11 Dec 2021 10:07:33: #3 Call peaks for each chromosome... INFO @ Sat, 11 Dec 2021 10:07:43: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20_peaks.xls INFO @ Sat, 11 Dec 2021 10:07:43: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20_peaks.narrowPeak INFO @ Sat, 11 Dec 2021 10:07:43: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/ERX3978911/ERX3978911.20_summits.bed INFO @ Sat, 11 Dec 2021 10:07:43: Done! pass1 - making usageList (14 chroms): 1 millis pass2 - checking and writing primary data (4585 records, 4 fields): 6 millis CompletedMACS2peakCalling