Job ID = 2008026 sra ファイルのダウンロード中... Read layout: PAIRED fastq に変換中... *** Error in `/home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump': free(): invalid pointer: 0x00007fffe403baef *** ======= Backtrace: ========= /lib64/libc.so.6(+0x81499)[0x7ffff736d499] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x4b5ebc] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x4204b6] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x4204a4] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x42499b] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x425a7b] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x437ae3] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x40fa72] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x411191] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x413b93] /home/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump[0x51a6ea] /lib64/libpthread.so.0(+0x7e25)[0x7ffff79c2e25] /lib64/libc.so.6(clone+0x6d)[0x7ffff73eabad] ======= Memory map: ======== 00400000-006c5000 r-xp 00000000 bd0:3f89c 144116276061639710 /lustre7/home/lustre4/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump.2.9.6 008c5000-008d3000 rw-p 002c5000 bd0:3f89c 144116276061639710 /lustre7/home/lustre4/okishinya/chipatlas/bin/sratoolkit.2.9.6-1-ubuntu64/bin/fasterq-dump.2.9.6 008d3000-011a9000 rw-p 00000000 00:00 0 [heap] 7fff98000000-7fff98023000 rw-p 00000000 00:00 0 7fff98023000-7fff9c000000 ---p 00000000 00:00 0 7fff9d7ff000-7fff9d800000 ---p 00000000 00:00 0 7fff9d800000-7fffa0000000 rw-p 00000000 00:00 0 7fffa0000000-7fffa0026000 rw-p 00000000 00:00 0 7fffa0026000-7fffa4000000 ---p 00000000 00:00 0 7fffa57ff000-7fffa5800000 ---p 00000000 00:00 0 7fffa5800000-7fffa8000000 rw-p 00000000 00:00 0 7fffa8000000-7fffa8034000 rw-p 00000000 00:00 0 7fffa8034000-7fffac000000 ---p 00000000 00:00 0 7fffad7ff000-7fffad800000 ---p 00000000 00:00 0 7fffad800000-7fffb0000000 rw-p 00000000 00:00 0 7fffb0000000-7fffb0025000 rw-p 00000000 00:00 0 7fffb0025000-7fffb4000000 ---p 00000000 00:00 0 7fffb8000000-7fffb8024000 rw-p 00000000 00:00 0 7fffb8024000-7fffbc000000 ---p 00000000 00:00 0 7fffbc000000-7fffbc024000 rw-p 00000000 00:00 0 7fffbc024000-7fffc0000000 ---p 00000000 00:00 0 7fffc0000000-7fffc0b9a000 rw-p 00000000 00:00 0 7fffc0b9a000-7fffc4000000 ---p 00000000 00:00 0 7fffc8000000-7fffcbffd000 rw-p 00000000 00:00 0 7fffcbffd000-7fffcc000000 ---p 00000000 00:00 0 7fffcc000000-7fffcfd95000 rw-p 00000000 00:00 0 7fffcfd95000-7fffd0000000 ---p 00000000 00:00 0 7fffd0000000-7fffd0954000 rw-p 00000000 00:00 0 7fffd0954000-7fffd4000000 ---p 00000000 00:00 0 7fffd47fd000-7fffd47fe000 ---p 00000000 00:00 0 7fffd47fe000-7fffd6ffe000 rw-p 00000000 00:00 0 7fffd6ffe000-7fffd6fff000 ---p 00000000 00:00 0 7fffd6fff000-7fffd97ff000 rw-p 00000000 00:00 0 7fffd97ff000-7fffd9800000 ---p 00000000 00:00 0 7fffd9800000-7fffdc000000 rw-p 00000000 00:00 0 7fffdc000000-7fffdcb40000 rw-p 00000000 00:00 0 7fffdcb40000-7fffe0000000 ---p 00000000 00:00 0 7fffe17ff000-7fffe1800000 ---p 00000000 00:00 0 7fffe1800000-7fffe4000000 rw-p 00000000 00:00 0 7fffe4000000-7fffe7fd8000 rw-p 00000000 00:00 0 7fffe7fd8000-7fffe8000000 ---p 00000000 00:00 0 7fffe97ff000-7fffe9800000 ---p 00000000 00:00 0 7fffe9800000-7fffec000000 rw-p 00000000 00:00 0 7fffec000000-7fffec95e000 rw-p 00000000 00:00 0 7fffec95e000-7ffff0000000 ---p 00000000 00:00 0 7ffff1aa1000-7ffff1ab6000 r-xp 00000000 103:06 6448798287 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 7ffff1ab6000-7ffff1cb5000 ---p 00015000 103:06 6448798287 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 7ffff1cb5000-7ffff1cb6000 r--p 00014000 103:06 6448798287 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 7ffff1cb6000-7ffff1cb7000 rw-p 00015000 103:06 6448798287 /usr/lib64/libgcc_s-4.8.5-20150702.so.1 7ffff1cb7000-7ffff1cb8000 ---p 00000000 00:00 0 7ffff1cb8000-7ffff44b8000 rw-p 00000000 00:00 0 7ffff44b8000-7ffff44b9000 ---p 00000000 00:00 0 7ffff44b9000-7ffff6cb9000 rw-p 00000000 00:00 0 7ffff6cb9000-7ffff6ccf000 r-xp 00000000 103:06 6447972061 /usr/lib64/libresolv-2.17.so 7ffff6ccf000-7ffff6ece000 ---p 00016000 103:06 6447972061 /usr/lib64/libresolv-2.17.so 7ffff6ece000-7ffff6ecf000 r--p 00015000 103:06 6447972061 /usr/lib64/libresolv-2.17.so 7ffff6ecf000-7ffff6ed0000 rw-p 00016000 103:06 6447972061 /usr/lib64/libresolv-2.17.so 7ffff6ed0000-7ffff6ed2000 rw-p 00000000 00:00 0 7ffff6ed2000-7ffff6ed7000 r-xp 00000000 103:06 6447929737 /usr/lib64/libnss_dns-2.17.so 7ffff6ed7000-7ffff70d7000 ---p 00005000 103:06 6447929737 /usr/lib64/libnss_dns-2.17.so 7ffff70d7000-7ffff70d8000 r--p 00005000 103:06 6447929737 /usr/lib64/libnss_dns-2.17.so 7ffff70d8000-7ffff70d9000 rw-p 00006000 103:06 6447929737 /usr/lib64/libnss_dns-2.17.so 7ffff70d9000-7ffff70e5000 r-xp 00000000 103:06 6447929738 /usr/lib64/libnss_files-2.17.so 7ffff70e5000-7ffff72e4000 ---p 0000c000 103:06 6447929738 /usr/lib64/libnss_files-2.17.so 7ffff72e4000-7ffff72e5000 r--p 0000b000 103:06 6447929738 /usr/lib64/libnss_files-2.17.so 7ffff72e5000-7ffff72e6000 rw-p 0000c000 103:06 6447929738 /usr/lib64/libnss_files-2.17.so 7ffff72e6000-7ffff72ec000 rw-p 00000000 00:00 0 7ffff72ec000-7ffff74af000 r-xp 00000000 103:06 6447821271 /usr/lib64/libc-2.17.so 7ffff74af000-7ffff76ae000 ---p 001c3000 103:06 6447821271 /usr/lib64/libc-2.17.so 7ffff76ae000-7ffff76b2000 r--p 001c2000 103:06 6447821271 /usr/lib64/libc-2.17.so 7ffff76b2000-7ffff76b4000 rw-p 001c6000 103:06 6447821271 /usr/lib64/libc-2.17.so 7ffff76b4000-7ffff76b9000 rw-p 00000000 00:00 0 7ffff76b9000-7ffff77ba000 r-xp 00000000 103:06 6447901624 /usr/lib64/libm-2.17.so 7ffff77ba000-7ffff79b9000 ---p 00101000 103:06 6447901624 /usr/lib64/libm-2.17.so 7ffff79b9000-7ffff79ba000 r--p 00100000 103:06 6447901624 /usr/lib64/libm-2.17.so 7ffff79ba000-7ffff79bb000 rw-p 00101000 103:06 6447901624 /usr/lib64/libm-2.17.so 7ffff79bb000-7ffff79d2000 r-xp 00000000 103:06 6447972046 /usr/lib64/libpthread-2.17.so 7ffff79d2000-7ffff7bd1000 ---p 00017000 103:06 6447972046 /usr/lib64/libpthread-2.17.so 7ffff7bd1000-7ffff7bd2000 r--p 00016000 103:06 6447972046 /usr/lib64/libpthread-2.17.so 7ffff7bd2000-7ffff7bd3000 rw-p 00017000 103:06 6447972046 /usr/lib64/libpthread-2.17.so 7ffff7bd3000-7ffff7bd7000 rw-p 00000000 00:00 0 7ffff7bd7000-7ffff7bd9000 r-xp 00000000 103:06 6447847945 /usr/lib64/libdl-2.17.so 7ffff7bd9000-7ffff7dd9000 ---p 00002000 103:06 6447847945 /usr/lib64/libdl-2.17.so 7ffff7dd9000-7ffff7dda000 r--p 00002000 103:06 6447847945 /usr/lib64/libdl-2.17.so 7ffff7dda000-7ffff7ddb000 rw-p 00003000 103:06 6447847945 /usr/lib64/libdl-2.17.so 7ffff7ddb000-7ffff7dfd000 r-xp 00000000 103:06 6447810064 /usr/lib64/ld-2.17.so 7ffff7fdb000-7ffff7fe0000 rw-p 00000000 00:00 0 7ffff7ff8000-7ffff7ffa000 rw-p 00000000 00:00 0 7ffff7ffa000-7ffff7ffc000 r-xp 00000000 00:00 0 [vdso] 7ffff7ffc000-7ffff7ffd000 r--p 00021000 103:06 6447810064 /usr/lib64/ld-2.17.so 7ffff7ffd000-7ffff7ffe000 rw-p 00022000 103:06 6447810064 /usr/lib64/ld-2.17.so 7ffff7ffe000-7ffff7fff000 rw-p 00000000 00:00 0 7ffffffdc000-7ffffffff000 rw-p 00000000 00:00 0 [stack] ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0 [vsyscall] 2019-07-05T08:13:35 fasterq-dump.2.9.6 fatal: SIGNAL - 6 spots read : 1,651,024 reads read : 3,302,048 reads written : 3,302,048 rm: cannot remove ‘[DSE]RR*’: 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: 00:01:34 1651024 reads; of these: 1651024 (100.00%) were paired; of these: 427276 (25.88%) aligned concordantly 0 times 1132017 (68.56%) aligned concordantly exactly 1 time 91731 (5.56%) aligned concordantly >1 times ---- 427276 pairs aligned concordantly 0 times; of these: 27650 (6.47%) aligned discordantly 1 time ---- 399626 pairs aligned 0 times concordantly or discordantly; of these: 799252 mates make up the pairs; of these: 776857 (97.20%) aligned 0 times 15052 (1.88%) aligned exactly 1 time 7343 (0.92%) aligned >1 times 76.47% overall alignment rate Time searching: 00:01:34 Overall time: 00:01:34 マッピングが完了しました。 samtools でBAM に変換中... [samopen] SAM header is present: 17 sequences. [bam_rmdup_core] processing reference chrI... [bam_rmdup_core] processing reference chrII... [bam_rmdup_core] processing reference chrIII... [bam_rmdup_core] processing reference chrIV... [bam_rmdup_core] processing reference chrIX... [bam_rmdup_core] processing reference chrM... [bam_rmdup_core] processing reference chrV... [bam_rmdup_core] processing reference chrVI... [bam_rmdup_core] processing reference chrVII... [bam_rmdup_core] processing reference chrVIII... [bam_rmdup_core] processing reference chrX... [bam_rmdup_core] processing reference chrXI... [bam_rmdup_core] processing reference chrXII... [bam_rmdup_core] processing reference chrXIII... [bam_rmdup_core] processing reference chrXIV... [bam_rmdup_core] processing reference chrXV... [bam_rmdup_core] processing reference chrXVI... [bam_rmdup_core] 346414 / 1248694 = 0.2774 in library ' ' BAM に変換しました。 Bed ファイルを作成中... BedGraph に変換中... INFO @ Fri, 05 Jul 2019 17:22:06: # Command line: callpeak -t /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam -f BAM -g 12100000 -n /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05 -q 1e-05 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam'] # control file = None # effective genome size = 1.21e+07 # 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, 05 Jul 2019 17:22:06: #1 read tag files... INFO @ Fri, 05 Jul 2019 17:22:06: #1 read treatment tags... INFO @ Fri, 05 Jul 2019 17:22:07: # Command line: callpeak -t /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam -f BAM -g 12100000 -n /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10 -q 1e-10 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam'] # control file = None # effective genome size = 1.21e+07 # 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, 05 Jul 2019 17:22:07: #1 read tag files... INFO @ Fri, 05 Jul 2019 17:22:07: #1 read treatment tags... INFO @ Fri, 05 Jul 2019 17:22:08: # Command line: callpeak -t /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam -f BAM -g 12100000 -n /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20 -q 1e-20 # ARGUMENTS LIST: # name = /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20 # format = BAM # ChIP-seq file = ['/home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.bam'] # control file = None # effective genome size = 1.21e+07 # 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, 05 Jul 2019 17:22:08: #1 read tag files... INFO @ Fri, 05 Jul 2019 17:22:08: #1 read treatment tags... INFO @ Fri, 05 Jul 2019 17:22:14: 1000000 INFO @ Fri, 05 Jul 2019 17:22:17: 1000000 INFO @ Fri, 05 Jul 2019 17:22:21: 1000000 INFO @ Fri, 05 Jul 2019 17:22:22: #1 tag size is determined as 100 bps INFO @ Fri, 05 Jul 2019 17:22:22: #1 tag size = 100 INFO @ Fri, 05 Jul 2019 17:22:22: #1 total tags in treatment: 881910 INFO @ Fri, 05 Jul 2019 17:22:22: #1 user defined the maximum tags... INFO @ Fri, 05 Jul 2019 17:22:22: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 05 Jul 2019 17:22:22: #1 tags after filtering in treatment: 850804 INFO @ Fri, 05 Jul 2019 17:22:22: #1 Redundant rate of treatment: 0.04 INFO @ Fri, 05 Jul 2019 17:22:22: #1 finished! INFO @ Fri, 05 Jul 2019 17:22:22: #2 Build Peak Model... INFO @ Fri, 05 Jul 2019 17:22:22: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 05 Jul 2019 17:22:22: #2 number of paired peaks: 180 WARNING @ Fri, 05 Jul 2019 17:22:22: Fewer paired peaks (180) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 180 pairs to build model! INFO @ Fri, 05 Jul 2019 17:22:22: start model_add_line... INFO @ Fri, 05 Jul 2019 17:22:22: start X-correlation... INFO @ Fri, 05 Jul 2019 17:22:22: end of X-cor INFO @ Fri, 05 Jul 2019 17:22:22: #2 finished! INFO @ Fri, 05 Jul 2019 17:22:22: #2 predicted fragment length is 226 bps INFO @ Fri, 05 Jul 2019 17:22:22: #2 alternative fragment length(s) may be 0,169,192,226,480,525,565,583 bps INFO @ Fri, 05 Jul 2019 17:22:22: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05_model.r INFO @ Fri, 05 Jul 2019 17:22:22: #3 Call peaks... INFO @ Fri, 05 Jul 2019 17:22:22: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 05 Jul 2019 17:22:25: #3 Call peaks for each chromosome... INFO @ Fri, 05 Jul 2019 17:22:25: #1 tag size is determined as 100 bps INFO @ Fri, 05 Jul 2019 17:22:25: #1 tag size = 100 INFO @ Fri, 05 Jul 2019 17:22:25: #1 total tags in treatment: 881910 INFO @ Fri, 05 Jul 2019 17:22:25: #1 user defined the maximum tags... INFO @ Fri, 05 Jul 2019 17:22:25: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 05 Jul 2019 17:22:25: #1 tags after filtering in treatment: 850804 INFO @ Fri, 05 Jul 2019 17:22:25: #1 Redundant rate of treatment: 0.04 INFO @ Fri, 05 Jul 2019 17:22:25: #1 finished! INFO @ Fri, 05 Jul 2019 17:22:25: #2 Build Peak Model... INFO @ Fri, 05 Jul 2019 17:22:25: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 05 Jul 2019 17:22:25: #2 number of paired peaks: 180 WARNING @ Fri, 05 Jul 2019 17:22:25: Fewer paired peaks (180) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 180 pairs to build model! INFO @ Fri, 05 Jul 2019 17:22:25: start model_add_line... INFO @ Fri, 05 Jul 2019 17:22:25: start X-correlation... INFO @ Fri, 05 Jul 2019 17:22:25: end of X-cor INFO @ Fri, 05 Jul 2019 17:22:25: #2 finished! INFO @ Fri, 05 Jul 2019 17:22:25: #2 predicted fragment length is 226 bps INFO @ Fri, 05 Jul 2019 17:22:25: #2 alternative fragment length(s) may be 0,169,192,226,480,525,565,583 bps INFO @ Fri, 05 Jul 2019 17:22:25: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10_model.r INFO @ Fri, 05 Jul 2019 17:22:25: #3 Call peaks... INFO @ Fri, 05 Jul 2019 17:22:25: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 05 Jul 2019 17:22:26: #4 Write output xls file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05_peaks.xls INFO @ Fri, 05 Jul 2019 17:22:26: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05_peaks.narrowPeak INFO @ Fri, 05 Jul 2019 17:22:26: #4 Write summits bed file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.05_summits.bed INFO @ Fri, 05 Jul 2019 17:22:26: Done! pass1 - making usageList (15 chroms): 1 millis pass2 - checking and writing primary data (69 records, 4 fields): 3 millis CompletedMACS2peakCalling INFO @ Fri, 05 Jul 2019 17:22:28: #3 Call peaks for each chromosome... INFO @ Fri, 05 Jul 2019 17:22:29: #4 Write output xls file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10_peaks.xls INFO @ Fri, 05 Jul 2019 17:22:29: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10_peaks.narrowPeak INFO @ Fri, 05 Jul 2019 17:22:29: #4 Write summits bed file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.10_summits.bed INFO @ Fri, 05 Jul 2019 17:22:29: Done! pass1 - making usageList (14 chroms): 1 millis pass2 - checking and writing primary data (38 records, 4 fields): 2 millis CompletedMACS2peakCalling INFO @ Fri, 05 Jul 2019 17:22:31: #1 tag size is determined as 100 bps INFO @ Fri, 05 Jul 2019 17:22:31: #1 tag size = 100 INFO @ Fri, 05 Jul 2019 17:22:31: #1 total tags in treatment: 881910 INFO @ Fri, 05 Jul 2019 17:22:31: #1 user defined the maximum tags... INFO @ Fri, 05 Jul 2019 17:22:31: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 05 Jul 2019 17:22:31: #1 tags after filtering in treatment: 850804 INFO @ Fri, 05 Jul 2019 17:22:31: #1 Redundant rate of treatment: 0.04 INFO @ Fri, 05 Jul 2019 17:22:31: #1 finished! INFO @ Fri, 05 Jul 2019 17:22:31: #2 Build Peak Model... INFO @ Fri, 05 Jul 2019 17:22:31: #2 looking for paired plus/minus strand peaks... INFO @ Fri, 05 Jul 2019 17:22:31: #2 number of paired peaks: 180 WARNING @ Fri, 05 Jul 2019 17:22:31: Fewer paired peaks (180) than 1000! Model may not be build well! Lower your MFOLD parameter may erase this warning. Now I will use 180 pairs to build model! INFO @ Fri, 05 Jul 2019 17:22:31: start model_add_line... INFO @ Fri, 05 Jul 2019 17:22:31: start X-correlation... INFO @ Fri, 05 Jul 2019 17:22:31: end of X-cor INFO @ Fri, 05 Jul 2019 17:22:31: #2 finished! INFO @ Fri, 05 Jul 2019 17:22:31: #2 predicted fragment length is 226 bps INFO @ Fri, 05 Jul 2019 17:22:31: #2 alternative fragment length(s) may be 0,169,192,226,480,525,565,583 bps INFO @ Fri, 05 Jul 2019 17:22:31: #2.2 Generate R script for model : /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20_model.r INFO @ Fri, 05 Jul 2019 17:22:31: #3 Call peaks... INFO @ Fri, 05 Jul 2019 17:22:31: #3 Pre-compute pvalue-qvalue table... INFO @ Fri, 05 Jul 2019 17:22:34: #3 Call peaks for each chromosome... INFO @ Fri, 05 Jul 2019 17:22:35: #4 Write output xls file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20_peaks.xls INFO @ Fri, 05 Jul 2019 17:22:35: #4 Write peak in narrowPeak format file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20_peaks.narrowPeak INFO @ Fri, 05 Jul 2019 17:22:35: #4 Write summits bed file... /home/okishinya/chipatlas/results/sacCer3/ERX471812/ERX471812.20_summits.bed INFO @ Fri, 05 Jul 2019 17:22:35: Done! pass1 - making usageList (9 chroms): 1 millis pass2 - checking and writing primary data (13 records, 4 fields): 2 millis CompletedMACS2peakCalling BedGraph に変換しました。 BigWig に変換中... BigWig に変換しました。