Job ID = 1293912 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... 2019-06-02T18:05:29 fasterq-dump.2.9.6 sys: connection not found while validating within network system module - Failed to Make Connection in KClientHttpOpen to 'sra-download-internal.ncbi.nlm.nih.gov:443' 2019-06-02T18:05:29 fasterq-dump.2.9.6 err: connection not found while validating within network system module - error with https open 'https://sra-download-internal.ncbi.nlm.nih.gov/sos/sra-pub-run-2/SRR495380/SRR495380.2' 2019-06-02T18:05:40 fasterq-dump.2.9.6 err: cmn_iter.c cmn_iter_open_tbl().VDBManagerOpenTableRead( 'SRR495380' ) -> RC(rcDB,rcMgr,rcOpening,rcTable,rcIncorrect) 2019-06-02T18:05:40 fasterq-dump.2.9.6 err: make_fastq_iter() -> RC(rcDB,rcMgr,rcOpening,rcTable,rcIncorrect) 2019-06-02T18:16:02 fasterq-dump.2.9.6 sys: connection not found while validating within network system module - Failed to Make Connection in KClientHttpOpen to 'sra-download-internal.ncbi.nlm.nih.gov:443' 2019-06-02T18:16:02 fasterq-dump.2.9.6 err: connection not found while validating within network system module - error with https open 'https://sra-download-internal.ncbi.nlm.nih.gov/sos/sra-pub-run-2/SRR495380/SRR495380.2' 2019-06-02T18:16:13 fasterq-dump.2.9.6 err: cmn_iter.c cmn_iter_open_tbl().VDBManagerOpenTableRead( 'SRR495380' ) -> RC(rcDB,rcMgr,rcOpening,rcTable,rcIncorrect) 2019-06-02T18:16:13 fasterq-dump.2.9.6 err: make_fastq_iter() -> RC(rcDB,rcMgr,rcOpening,rcTable,rcIncorrect) 2019-06-02T18:18:08 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:18:08 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:18:08 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:18:08 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:28:28 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:28:28 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:28:28 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2019-06-02T18:28:28 fasterq-dump.2.9.6 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) spots read : 60,541,438 reads read : 121,082,876 reads written : 121,082,876 rm: cannot remove ‘[DSE]RR*’: No such file or directory rm: cannot remove ‘fastqDump_tmp*’: No such file or directory fastq に変換しました。 bowtie でマッピング中... Time loading reference: 00:00:00 Time loading forward index: 00:00:00 Time loading mirror index: 00:00:00 Multiseed full-index search: 06:06:39 60541438 reads; of these: 60541438 (100.00%) were paired; of these: 3613732 (5.97%) aligned concordantly 0 times 31076526 (51.33%) aligned concordantly exactly 1 time 25851180 (42.70%) aligned concordantly >1 times ---- 3613732 pairs aligned concordantly 0 times; of these: 727029 (20.12%) aligned discordantly 1 time ---- 2886703 pairs aligned 0 times concordantly or discordantly; of these: 5773406 mates make up the pairs; of these: 3825703 (66.26%) aligned 0 times 762286 (13.20%) aligned exactly 1 time 1185417 (20.53%) aligned >1 times 96.84% overall alignment rate Time searching: 06:06:39 Overall time: 06:06:39 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 15 sequences. [bam_sort_core] merging from 52 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] 9401602 / 57261992 = 0.1642 in library ' ' BAM に変換しました。 Bed ファイルを作成中... BedGraph に変換中... INFO @ Mon, 03 Jun 2019 10:46:59: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.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 @ Mon, 03 Jun 2019 10:46:59: #1 read tag files... INFO @ Mon, 03 Jun 2019 10:46:59: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 10:46:59: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.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 @ Mon, 03 Jun 2019 10:46:59: #1 read tag files... INFO @ Mon, 03 Jun 2019 10:46:59: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 10:46:59: # Command line: callpeak -t /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.bam -f BAM -g dm -n /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.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 @ Mon, 03 Jun 2019 10:46:59: #1 read tag files... INFO @ Mon, 03 Jun 2019 10:46:59: #1 read treatment tags... INFO @ Mon, 03 Jun 2019 10:47:10: 1000000 INFO @ Mon, 03 Jun 2019 10:47:11: 1000000 INFO @ Mon, 03 Jun 2019 10:47:12: 1000000 INFO @ Mon, 03 Jun 2019 10:47:22: 2000000 INFO @ Mon, 03 Jun 2019 10:47:23: 2000000 INFO @ Mon, 03 Jun 2019 10:47:26: 2000000 INFO @ Mon, 03 Jun 2019 10:47:33: 3000000 INFO @ Mon, 03 Jun 2019 10:47:34: 3000000 INFO @ Mon, 03 Jun 2019 10:47:38: 3000000 INFO @ Mon, 03 Jun 2019 10:47:43: 4000000 INFO @ Mon, 03 Jun 2019 10:47:45: 4000000 INFO @ Mon, 03 Jun 2019 10:47:50: 4000000 INFO @ Mon, 03 Jun 2019 10:47:53: 5000000 INFO @ Mon, 03 Jun 2019 10:47:56: 5000000 INFO @ Mon, 03 Jun 2019 10:48:03: 5000000 INFO @ Mon, 03 Jun 2019 10:48:04: 6000000 INFO @ Mon, 03 Jun 2019 10:48:07: 6000000 INFO @ Mon, 03 Jun 2019 10:48:14: 7000000 INFO @ Mon, 03 Jun 2019 10:48:15: 6000000 INFO @ Mon, 03 Jun 2019 10:48:18: 7000000 INFO @ Mon, 03 Jun 2019 10:48:25: 8000000 INFO @ Mon, 03 Jun 2019 10:48:27: 7000000 INFO @ Mon, 03 Jun 2019 10:48:30: 8000000 INFO @ Mon, 03 Jun 2019 10:48:35: 9000000 INFO @ Mon, 03 Jun 2019 10:48:39: 8000000 INFO @ Mon, 03 Jun 2019 10:48:40: 9000000 INFO @ Mon, 03 Jun 2019 10:48:45: 10000000 INFO @ Mon, 03 Jun 2019 10:48:52: 10000000 INFO @ Mon, 03 Jun 2019 10:48:52: 9000000 INFO @ Mon, 03 Jun 2019 10:48:57: 11000000 INFO @ Mon, 03 Jun 2019 10:49:03: 11000000 INFO @ Mon, 03 Jun 2019 10:49:05: 10000000 INFO @ Mon, 03 Jun 2019 10:49:07: 12000000 INFO @ Mon, 03 Jun 2019 10:49:14: 12000000 INFO @ Mon, 03 Jun 2019 10:49:17: 11000000 INFO @ Mon, 03 Jun 2019 10:49:18: 13000000 INFO @ Mon, 03 Jun 2019 10:49:25: 13000000 INFO @ Mon, 03 Jun 2019 10:49:28: 14000000 INFO @ Mon, 03 Jun 2019 10:49:29: 12000000 INFO @ Mon, 03 Jun 2019 10:49:36: 14000000 INFO @ Mon, 03 Jun 2019 10:49:38: 15000000 INFO @ Mon, 03 Jun 2019 10:49:42: 13000000 INFO @ Mon, 03 Jun 2019 10:49:47: 15000000 INFO @ Mon, 03 Jun 2019 10:49:48: 16000000 INFO @ Mon, 03 Jun 2019 10:49:54: 14000000 INFO @ Mon, 03 Jun 2019 10:49:58: 16000000 INFO @ Mon, 03 Jun 2019 10:49:59: 17000000 INFO @ Mon, 03 Jun 2019 10:50:06: 15000000 INFO @ Mon, 03 Jun 2019 10:50:09: 18000000 INFO @ Mon, 03 Jun 2019 10:50:09: 17000000 INFO @ Mon, 03 Jun 2019 10:50:19: 16000000 INFO @ Mon, 03 Jun 2019 10:50:19: 19000000 INFO @ Mon, 03 Jun 2019 10:50:20: 18000000 INFO @ Mon, 03 Jun 2019 10:50:30: 20000000 INFO @ Mon, 03 Jun 2019 10:50:31: 17000000 INFO @ Mon, 03 Jun 2019 10:50:31: 19000000 INFO @ Mon, 03 Jun 2019 10:50:41: 21000000 INFO @ Mon, 03 Jun 2019 10:50:43: 20000000 INFO @ Mon, 03 Jun 2019 10:50:44: 18000000 INFO @ Mon, 03 Jun 2019 10:50:51: 22000000 INFO @ Mon, 03 Jun 2019 10:50:53: 21000000 INFO @ Mon, 03 Jun 2019 10:50:56: 19000000 INFO @ Mon, 03 Jun 2019 10:51:02: 23000000 INFO @ Mon, 03 Jun 2019 10:51:04: 22000000 INFO @ Mon, 03 Jun 2019 10:51:08: 20000000 INFO @ Mon, 03 Jun 2019 10:51:13: 24000000 INFO @ Mon, 03 Jun 2019 10:51:15: 23000000 INFO @ Mon, 03 Jun 2019 10:51:20: 21000000 INFO @ Mon, 03 Jun 2019 10:51:23: 25000000 INFO @ Mon, 03 Jun 2019 10:51:26: 24000000 INFO @ Mon, 03 Jun 2019 10:51:32: 22000000 INFO @ Mon, 03 Jun 2019 10:51:34: 26000000 INFO @ Mon, 03 Jun 2019 10:51:37: 25000000 INFO @ Mon, 03 Jun 2019 10:51:45: 23000000 INFO @ Mon, 03 Jun 2019 10:51:45: 27000000 INFO @ Mon, 03 Jun 2019 10:51:48: 26000000 INFO @ Mon, 03 Jun 2019 10:51:57: 28000000 INFO @ Mon, 03 Jun 2019 10:51:57: 24000000 INFO @ Mon, 03 Jun 2019 10:52:00: 27000000 INFO @ Mon, 03 Jun 2019 10:52:09: 29000000 INFO @ Mon, 03 Jun 2019 10:52:10: 25000000 INFO @ Mon, 03 Jun 2019 10:52:11: 28000000 INFO @ Mon, 03 Jun 2019 10:52:19: 30000000 INFO @ Mon, 03 Jun 2019 10:52:22: 29000000 INFO @ Mon, 03 Jun 2019 10:52:22: 26000000 INFO @ Mon, 03 Jun 2019 10:52:29: 31000000 INFO @ Mon, 03 Jun 2019 10:52:33: 30000000 INFO @ Mon, 03 Jun 2019 10:52:34: 27000000 INFO @ Mon, 03 Jun 2019 10:52:40: 32000000 INFO @ Mon, 03 Jun 2019 10:52:44: 31000000 INFO @ Mon, 03 Jun 2019 10:52:46: 28000000 INFO @ Mon, 03 Jun 2019 10:52:50: 33000000 INFO @ Mon, 03 Jun 2019 10:52:55: 32000000 INFO @ Mon, 03 Jun 2019 10:52:59: 29000000 INFO @ Mon, 03 Jun 2019 10:53:00: 34000000 INFO @ Mon, 03 Jun 2019 10:53:05: 33000000 INFO @ Mon, 03 Jun 2019 10:53:11: 35000000 INFO @ Mon, 03 Jun 2019 10:53:11: 30000000 INFO @ Mon, 03 Jun 2019 10:53:15: 34000000 INFO @ Mon, 03 Jun 2019 10:53:21: 36000000 INFO @ Mon, 03 Jun 2019 10:53:23: 31000000 INFO @ Mon, 03 Jun 2019 10:53:26: 35000000 INFO @ Mon, 03 Jun 2019 10:53:32: 37000000 INFO @ Mon, 03 Jun 2019 10:53:36: 32000000 INFO @ Mon, 03 Jun 2019 10:53:37: 36000000 INFO @ Mon, 03 Jun 2019 10:53:43: 38000000 INFO @ Mon, 03 Jun 2019 10:53:48: 37000000 INFO @ Mon, 03 Jun 2019 10:53:49: 33000000 INFO @ Mon, 03 Jun 2019 10:53:54: 39000000 INFO @ Mon, 03 Jun 2019 10:53:58: 38000000 INFO @ Mon, 03 Jun 2019 10:54:01: 34000000 INFO @ Mon, 03 Jun 2019 10:54:04: 40000000 INFO @ Mon, 03 Jun 2019 10:54:08: 39000000 INFO @ Mon, 03 Jun 2019 10:54:13: 35000000 INFO @ Mon, 03 Jun 2019 10:54:14: 41000000 INFO @ Mon, 03 Jun 2019 10:54:18: 40000000 INFO @ Mon, 03 Jun 2019 10:54:25: 42000000 INFO @ Mon, 03 Jun 2019 10:54:25: 36000000 INFO @ Mon, 03 Jun 2019 10:54:29: 41000000 INFO @ Mon, 03 Jun 2019 10:54:35: 43000000 INFO @ Mon, 03 Jun 2019 10:54:37: 37000000 INFO @ Mon, 03 Jun 2019 10:54:39: 42000000 INFO @ Mon, 03 Jun 2019 10:54:46: 44000000 INFO @ Mon, 03 Jun 2019 10:54:50: 43000000 INFO @ Mon, 03 Jun 2019 10:54:50: 38000000 INFO @ Mon, 03 Jun 2019 10:54:56: 45000000 INFO @ Mon, 03 Jun 2019 10:55:00: 44000000 INFO @ Mon, 03 Jun 2019 10:55:02: 39000000 INFO @ Mon, 03 Jun 2019 10:55:07: 46000000 INFO @ Mon, 03 Jun 2019 10:55:11: 45000000 INFO @ Mon, 03 Jun 2019 10:55:15: 40000000 INFO @ Mon, 03 Jun 2019 10:55:17: 47000000 INFO @ Mon, 03 Jun 2019 10:55:21: 46000000 INFO @ Mon, 03 Jun 2019 10:55:27: 41000000 INFO @ Mon, 03 Jun 2019 10:55:28: 48000000 INFO @ Mon, 03 Jun 2019 10:55:32: 47000000 INFO @ Mon, 03 Jun 2019 10:55:38: 49000000 INFO @ Mon, 03 Jun 2019 10:55:39: 42000000 INFO @ Mon, 03 Jun 2019 10:55:42: 48000000 INFO @ Mon, 03 Jun 2019 10:55:49: 50000000 INFO @ Mon, 03 Jun 2019 10:55:51: 43000000 INFO @ Mon, 03 Jun 2019 10:55:52: 49000000 INFO @ Mon, 03 Jun 2019 10:55:59: 51000000 INFO @ Mon, 03 Jun 2019 10:56:02: 50000000 INFO @ Mon, 03 Jun 2019 10:56:03: 44000000 INFO @ Mon, 03 Jun 2019 10:56:09: 52000000 INFO @ Mon, 03 Jun 2019 10:56:13: 51000000 INFO @ Mon, 03 Jun 2019 10:56:15: 45000000 INFO @ Mon, 03 Jun 2019 10:56:20: 53000000 INFO @ Mon, 03 Jun 2019 10:56:23: 52000000 INFO @ Mon, 03 Jun 2019 10:56:28: 46000000 INFO @ Mon, 03 Jun 2019 10:56:30: 54000000 INFO @ Mon, 03 Jun 2019 10:56:34: 53000000 INFO @ Mon, 03 Jun 2019 10:56:41: 47000000 INFO @ Mon, 03 Jun 2019 10:56:42: 55000000 INFO @ Mon, 03 Jun 2019 10:56:46: 54000000 INFO @ Mon, 03 Jun 2019 10:56:53: 56000000 INFO @ Mon, 03 Jun 2019 10:56:53: 48000000 INFO @ Mon, 03 Jun 2019 10:56:56: 55000000 INFO @ Mon, 03 Jun 2019 10:57:02: 57000000 INFO @ Mon, 03 Jun 2019 10:57:05: 49000000 INFO @ Mon, 03 Jun 2019 10:57:06: 56000000 INFO @ Mon, 03 Jun 2019 10:57:13: 58000000 INFO @ Mon, 03 Jun 2019 10:57:16: 57000000 INFO @ Mon, 03 Jun 2019 10:57:17: 50000000 INFO @ Mon, 03 Jun 2019 10:57:23: 59000000 INFO @ Mon, 03 Jun 2019 10:57:26: 58000000 INFO @ Mon, 03 Jun 2019 10:57:29: 51000000 INFO @ Mon, 03 Jun 2019 10:57:32: 60000000 INFO @ Mon, 03 Jun 2019 10:57:36: 59000000 INFO @ Mon, 03 Jun 2019 10:57:42: 52000000 INFO @ Mon, 03 Jun 2019 10:57:42: 61000000 INFO @ Mon, 03 Jun 2019 10:57:46: 60000000 INFO @ Mon, 03 Jun 2019 10:57:52: 62000000 INFO @ Mon, 03 Jun 2019 10:57:55: 53000000 INFO @ Mon, 03 Jun 2019 10:57:55: 61000000 INFO @ Mon, 03 Jun 2019 10:58:01: 63000000 INFO @ Mon, 03 Jun 2019 10:58:05: 62000000 INFO @ Mon, 03 Jun 2019 10:58:08: 54000000 INFO @ Mon, 03 Jun 2019 10:58:11: 64000000 INFO @ Mon, 03 Jun 2019 10:58:15: 63000000 INFO @ Mon, 03 Jun 2019 10:58:21: 65000000 INFO @ Mon, 03 Jun 2019 10:58:22: 55000000 INFO @ Mon, 03 Jun 2019 10:58:25: 64000000 INFO @ Mon, 03 Jun 2019 10:58:31: 66000000 INFO @ Mon, 03 Jun 2019 10:58:35: 65000000 INFO @ Mon, 03 Jun 2019 10:58:35: 56000000 INFO @ Mon, 03 Jun 2019 10:58:41: 67000000 INFO @ Mon, 03 Jun 2019 10:58:45: 66000000 INFO @ Mon, 03 Jun 2019 10:58:48: 57000000 INFO @ Mon, 03 Jun 2019 10:58:51: 68000000 INFO @ Mon, 03 Jun 2019 10:58:54: 67000000 INFO @ Mon, 03 Jun 2019 10:59:00: 69000000 INFO @ Mon, 03 Jun 2019 10:59:00: 58000000 INFO @ Mon, 03 Jun 2019 10:59:04: 68000000 INFO @ Mon, 03 Jun 2019 10:59:10: 70000000 INFO @ Mon, 03 Jun 2019 10:59:12: 59000000 INFO @ Mon, 03 Jun 2019 10:59:13: 69000000 INFO @ Mon, 03 Jun 2019 10:59:20: 71000000 INFO @ Mon, 03 Jun 2019 10:59:23: 70000000 INFO @ Mon, 03 Jun 2019 10:59:24: 60000000 INFO @ Mon, 03 Jun 2019 10:59:30: 72000000 INFO @ Mon, 03 Jun 2019 10:59:33: 71000000 INFO @ Mon, 03 Jun 2019 10:59:36: 61000000 INFO @ Mon, 03 Jun 2019 10:59:40: 73000000 INFO @ Mon, 03 Jun 2019 10:59:43: 72000000 INFO @ Mon, 03 Jun 2019 10:59:47: 62000000 INFO @ Mon, 03 Jun 2019 10:59:50: 74000000 INFO @ Mon, 03 Jun 2019 10:59:53: 73000000 INFO @ Mon, 03 Jun 2019 10:59:59: 63000000 INFO @ Mon, 03 Jun 2019 11:00:00: 75000000 INFO @ Mon, 03 Jun 2019 11:00:03: 74000000 INFO @ Mon, 03 Jun 2019 11:00:10: 76000000 INFO @ Mon, 03 Jun 2019 11:00:10: 64000000 INFO @ Mon, 03 Jun 2019 11:00:13: 75000000 INFO @ Mon, 03 Jun 2019 11:00:20: 77000000 INFO @ Mon, 03 Jun 2019 11:00:22: 65000000 INFO @ Mon, 03 Jun 2019 11:00:23: 76000000 INFO @ Mon, 03 Jun 2019 11:00:31: 78000000 INFO @ Mon, 03 Jun 2019 11:00:32: 77000000 INFO @ Mon, 03 Jun 2019 11:00:34: 66000000 INFO @ Mon, 03 Jun 2019 11:00:41: 79000000 INFO @ Mon, 03 Jun 2019 11:00:42: 78000000 INFO @ Mon, 03 Jun 2019 11:00:45: 67000000 INFO @ Mon, 03 Jun 2019 11:00:51: 80000000 INFO @ Mon, 03 Jun 2019 11:00:52: 79000000 INFO @ Mon, 03 Jun 2019 11:00:56: 68000000 INFO @ Mon, 03 Jun 2019 11:01:01: 80000000 INFO @ Mon, 03 Jun 2019 11:01:01: 81000000 INFO @ Mon, 03 Jun 2019 11:01:08: 69000000 INFO @ Mon, 03 Jun 2019 11:01:11: 81000000 INFO @ Mon, 03 Jun 2019 11:01:11: 82000000 INFO @ Mon, 03 Jun 2019 11:01:19: 70000000 INFO @ Mon, 03 Jun 2019 11:01:21: 82000000 INFO @ Mon, 03 Jun 2019 11:01:21: 83000000 INFO @ Mon, 03 Jun 2019 11:01:31: 71000000 INFO @ Mon, 03 Jun 2019 11:01:31: 83000000 INFO @ Mon, 03 Jun 2019 11:01:34: 84000000 INFO @ Mon, 03 Jun 2019 11:01:41: 84000000 INFO @ Mon, 03 Jun 2019 11:01:43: 72000000 INFO @ Mon, 03 Jun 2019 11:01:44: 85000000 INFO @ Mon, 03 Jun 2019 11:01:51: 85000000 INFO @ Mon, 03 Jun 2019 11:01:53: 86000000 INFO @ Mon, 03 Jun 2019 11:01:55: 73000000 INFO @ Mon, 03 Jun 2019 11:02:00: 86000000 INFO @ Mon, 03 Jun 2019 11:02:03: 87000000 INFO @ Mon, 03 Jun 2019 11:02:06: 74000000 INFO @ Mon, 03 Jun 2019 11:02:10: 87000000 INFO @ Mon, 03 Jun 2019 11:02:14: 88000000 INFO @ Mon, 03 Jun 2019 11:02:18: 75000000 INFO @ Mon, 03 Jun 2019 11:02:21: 88000000 INFO @ Mon, 03 Jun 2019 11:02:24: 89000000 INFO @ Mon, 03 Jun 2019 11:02:29: 76000000 INFO @ Mon, 03 Jun 2019 11:02:31: 89000000 INFO @ Mon, 03 Jun 2019 11:02:34: 90000000 INFO @ Mon, 03 Jun 2019 11:02:41: 77000000 INFO @ Mon, 03 Jun 2019 11:02:41: 90000000 INFO @ Mon, 03 Jun 2019 11:02:44: 91000000 INFO @ Mon, 03 Jun 2019 11:02:50: 91000000 INFO @ Mon, 03 Jun 2019 11:02:52: 78000000 INFO @ Mon, 03 Jun 2019 11:02:53: 92000000 INFO @ Mon, 03 Jun 2019 11:03:00: 92000000 INFO @ Mon, 03 Jun 2019 11:03:03: 93000000 INFO @ Mon, 03 Jun 2019 11:03:03: 79000000 INFO @ Mon, 03 Jun 2019 11:03:10: 93000000 INFO @ Mon, 03 Jun 2019 11:03:14: 94000000 INFO @ Mon, 03 Jun 2019 11:03:16: 80000000 INFO @ Mon, 03 Jun 2019 11:03:20: 94000000 INFO @ Mon, 03 Jun 2019 11:03:24: 95000000 INFO @ Mon, 03 Jun 2019 11:03:27: 81000000 INFO @ Mon, 03 Jun 2019 11:03:30: 95000000 INFO @ Mon, 03 Jun 2019 11:03:34: 96000000 INFO @ Mon, 03 Jun 2019 11:03:38: 82000000 INFO @ Mon, 03 Jun 2019 11:03:41: 96000000 INFO @ Mon, 03 Jun 2019 11:03:43: 97000000 INFO @ Mon, 03 Jun 2019 11:03:50: 83000000 INFO @ Mon, 03 Jun 2019 11:03:51: 97000000 INFO @ Mon, 03 Jun 2019 11:03:53: 98000000 INFO @ Mon, 03 Jun 2019 11:03:57: #1 tag size is determined as 100 bps INFO @ Mon, 03 Jun 2019 11:03:57: #1 tag size = 100 INFO @ Mon, 03 Jun 2019 11:03:57: #1 total tags in treatment: 47585890 INFO @ Mon, 03 Jun 2019 11:03:57: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 11:03:57: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 11:03:59: #1 tags after filtering in treatment: 38151387 INFO @ Mon, 03 Jun 2019 11:03:59: #1 Redundant rate of treatment: 0.20 INFO @ Mon, 03 Jun 2019 11:03:59: #1 finished! INFO @ Mon, 03 Jun 2019 11:03:59: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 11:03:59: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 11:04:00: 98000000 INFO @ Mon, 03 Jun 2019 11:04:01: 84000000 INFO @ Mon, 03 Jun 2019 11:04:03: #2 number of paired peaks: 364 WARNING @ Mon, 03 Jun 2019 11:04:03: Fewer paired peaks (364) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 364 pairs to build model! INFO @ Mon, 03 Jun 2019 11:04:03: start model_add_line... INFO @ Mon, 03 Jun 2019 11:04:04: start X-correlation... INFO @ Mon, 03 Jun 2019 11:04:04: end of X-cor INFO @ Mon, 03 Jun 2019 11:04:04: #2 finished! INFO @ Mon, 03 Jun 2019 11:04:04: #2 predicted fragment length is 104 bps INFO @ Mon, 03 Jun 2019 11:04:04: #2 alternative fragment length(s) may be 2,104 bps INFO @ Mon, 03 Jun 2019 11:04:04: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20_model.r WARNING @ Mon, 03 Jun 2019 11:04:04: #2 Since the d (104) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Mon, 03 Jun 2019 11:04:04: #2 You may need to consider one of the other alternative d(s): 2,104 WARNING @ Mon, 03 Jun 2019 11:04:04: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Mon, 03 Jun 2019 11:04:04: #3 Call peaks... INFO @ Mon, 03 Jun 2019 11:04:04: #3 Pre-compute pvalue-qvalue table... INFO @ Mon, 03 Jun 2019 11:04:05: #1 tag size is determined as 100 bps INFO @ Mon, 03 Jun 2019 11:04:05: #1 tag size = 100 INFO @ Mon, 03 Jun 2019 11:04:05: #1 total tags in treatment: 47585890 INFO @ Mon, 03 Jun 2019 11:04:05: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 11:04:05: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 11:04:07: #1 tags after filtering in treatment: 38151387 INFO @ Mon, 03 Jun 2019 11:04:07: #1 Redundant rate of treatment: 0.20 INFO @ Mon, 03 Jun 2019 11:04:07: #1 finished! INFO @ Mon, 03 Jun 2019 11:04:07: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 11:04:07: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 11:04:11: #2 number of paired peaks: 364 WARNING @ Mon, 03 Jun 2019 11:04:11: Fewer paired peaks (364) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 364 pairs to build model! INFO @ Mon, 03 Jun 2019 11:04:11: start model_add_line... INFO @ Mon, 03 Jun 2019 11:04:12: start X-correlation... INFO @ Mon, 03 Jun 2019 11:04:12: end of X-cor INFO @ Mon, 03 Jun 2019 11:04:12: #2 finished! INFO @ Mon, 03 Jun 2019 11:04:12: #2 predicted fragment length is 104 bps INFO @ Mon, 03 Jun 2019 11:04:12: #2 alternative fragment length(s) may be 2,104 bps INFO @ Mon, 03 Jun 2019 11:04:12: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10_model.r WARNING @ Mon, 03 Jun 2019 11:04:12: #2 Since the d (104) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Mon, 03 Jun 2019 11:04:12: #2 You may need to consider one of the other alternative d(s): 2,104 WARNING @ Mon, 03 Jun 2019 11:04:12: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Mon, 03 Jun 2019 11:04:12: #3 Call peaks... INFO @ Mon, 03 Jun 2019 11:04:12: #3 Pre-compute pvalue-qvalue table... INFO @ Mon, 03 Jun 2019 11:04:12: 85000000 INFO @ Mon, 03 Jun 2019 11:04:23: 86000000 INFO @ Mon, 03 Jun 2019 11:04:35: 87000000 INFO @ Mon, 03 Jun 2019 11:04:46: 88000000 INFO @ Mon, 03 Jun 2019 11:04:58: 89000000 INFO @ Mon, 03 Jun 2019 11:05:09: 90000000 INFO @ Mon, 03 Jun 2019 11:05:20: 91000000 INFO @ Mon, 03 Jun 2019 11:05:31: 92000000 INFO @ Mon, 03 Jun 2019 11:05:35: #3 Call peaks for each chromosome... INFO @ Mon, 03 Jun 2019 11:05:42: 93000000 INFO @ Mon, 03 Jun 2019 11:05:44: #3 Call peaks for each chromosome... INFO @ Mon, 03 Jun 2019 11:05:53: 94000000 INFO @ Mon, 03 Jun 2019 11:06:04: 95000000 INFO @ Mon, 03 Jun 2019 11:06:15: 96000000 INFO @ Mon, 03 Jun 2019 11:06:19: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20_peaks.xls INFO @ Mon, 03 Jun 2019 11:06:19: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20_peaks.narrowPeak INFO @ Mon, 03 Jun 2019 11:06:19: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.20_summits.bed INFO @ Mon, 03 Jun 2019 11:06:19: Done! pass1 - making usageList (14 chroms): 1 millis pass2 - checking and writing primary data (2520 records, 4 fields): 8 millis CompletedMACS2peakCalling INFO @ Mon, 03 Jun 2019 11:06:26: 97000000 INFO @ Mon, 03 Jun 2019 11:06:28: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10_peaks.xls INFO @ Mon, 03 Jun 2019 11:06:28: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10_peaks.narrowPeak INFO @ Mon, 03 Jun 2019 11:06:28: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.10_summits.bed INFO @ Mon, 03 Jun 2019 11:06:28: Done! pass1 - making usageList (15 chroms): 2 millis pass2 - checking and writing primary data (6849 records, 4 fields): 10 millis CompletedMACS2peakCalling INFO @ Mon, 03 Jun 2019 11:06:37: 98000000 INFO @ Mon, 03 Jun 2019 11:06:43: #1 tag size is determined as 100 bps INFO @ Mon, 03 Jun 2019 11:06:43: #1 tag size = 100 INFO @ Mon, 03 Jun 2019 11:06:43: #1 total tags in treatment: 47585890 INFO @ Mon, 03 Jun 2019 11:06:43: #1 user defined the maximum tags... INFO @ Mon, 03 Jun 2019 11:06:43: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Mon, 03 Jun 2019 11:06:45: #1 tags after filtering in treatment: 38151387 INFO @ Mon, 03 Jun 2019 11:06:45: #1 Redundant rate of treatment: 0.20 INFO @ Mon, 03 Jun 2019 11:06:45: #1 finished! INFO @ Mon, 03 Jun 2019 11:06:45: #2 Build Peak Model... INFO @ Mon, 03 Jun 2019 11:06:45: #2 looking for paired plus/minus strand peaks... INFO @ Mon, 03 Jun 2019 11:06:49: #2 number of paired peaks: 364 WARNING @ Mon, 03 Jun 2019 11:06:49: Fewer paired peaks (364) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 364 pairs to build model! INFO @ Mon, 03 Jun 2019 11:06:49: start model_add_line... INFO @ Mon, 03 Jun 2019 11:06:49: start X-correlation... INFO @ Mon, 03 Jun 2019 11:06:49: end of X-cor INFO @ Mon, 03 Jun 2019 11:06:49: #2 finished! INFO @ Mon, 03 Jun 2019 11:06:49: #2 predicted fragment length is 104 bps INFO @ Mon, 03 Jun 2019 11:06:49: #2 alternative fragment length(s) may be 2,104 bps INFO @ Mon, 03 Jun 2019 11:06:49: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05_model.r WARNING @ Mon, 03 Jun 2019 11:06:49: #2 Since the d (104) calculated from paired-peaks are smaller than 2*tag length, it may be influenced by unknown sequencing problem! WARNING @ Mon, 03 Jun 2019 11:06:49: #2 You may need to consider one of the other alternative d(s): 2,104 WARNING @ Mon, 03 Jun 2019 11:06:49: #2 You can restart the process with --nomodel --extsize XXX with your choice or an arbitrary number. Nontheless, MACS will continute computing. INFO @ Mon, 03 Jun 2019 11:06:49: #3 Call peaks... INFO @ Mon, 03 Jun 2019 11:06:49: #3 Pre-compute pvalue-qvalue table... BedGraph に変換しました。 BigWig に変換中... INFO @ Mon, 03 Jun 2019 11:08:20: #3 Call peaks for each chromosome... INFO @ Mon, 03 Jun 2019 11:09:05: #4 Write output xls file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05_peaks.xls INFO @ Mon, 03 Jun 2019 11:09:05: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05_peaks.narrowPeak INFO @ Mon, 03 Jun 2019 11:09:05: #4 Write summits bed file... /home/okishinya/chipatlas/results/dm3/SRX147006/SRX147006.05_summits.bed INFO @ Mon, 03 Jun 2019 11:09:05: Done! pass1 - making usageList (15 chroms): 4 millis pass2 - checking and writing primary data (11948 records, 4 fields): 16 millis CompletedMACS2peakCalling BigWig に変換しました。