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