Job ID = 10609137 sra ファイルのダウンロード中... Completed: 667727K bytes transferred in 37 seconds (146805K bits/sec), in 1 file. sra ファイルのダウンロードが完了しました。 Read layout: PAIRED fastq に変換中... 2018-05-03T22:04:30 fastq-dump.2.9.0 sys: timeout exhausted while reading file within network system module - mbedtls_ssl_read returned -76 ( NET - Reading information from the socket failed ) 2018-05-03T22:05:10 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - mbedtls_ssl_handshake returned -76 ( NET - Reading information from the socket failed ) 2018-05-03T22:05:10 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - ktls_handshake failed while accessing '130.14.250.27' from '172.20.0.116' 2018-05-03T22:05:10 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - Failed to create TLS stream for 'sra-download.ncbi.nlm.nih.gov' (130.14.250.27) from '172.20.0.116' 2018-05-03T22:05:10 fastq-dump.2.9.0 err: connection failed while opening file within cryptographic module - error with https open 'https://sra-download.ncbi.nlm.nih.gov/traces/refseq/AE014297.2' 2018-05-03T22:06:23 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - mbedtls_ssl_handshake returned -76 ( NET - Reading information from the socket failed ) 2018-05-03T22:06:23 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - ktls_handshake failed while accessing '130.14.250.25' from '172.20.0.116' 2018-05-03T22:06:23 fastq-dump.2.9.0 sys: connection failed while opening file within cryptographic module - Failed to create TLS stream for 'sra-download.ncbi.nlm.nih.gov' (130.14.250.25) from '172.20.0.116' 2018-05-03T22:06:23 fastq-dump.2.9.0 err: connection failed while opening file within cryptographic module - error with https open 'https://sra-download.ncbi.nlm.nih.gov/traces/refseq/AE014297.2' 2018-05-03T22:06:34 fastq-dump.2.9.0 err: data corrupt while selecting function within transform module - failed /home/okishinya/chipatlas/results/dm3/SRX3145638/SRR5989886.sra ============================================================= An error occurred during processing. A report was generated into the file '/home/okishinya/ncbi_error_report.xml'. If the problem persists, you may consider sending the file to 'sra@ncbi.nlm.nih.gov' for assistance. ============================================================= rm: cannot remove `[DSE]RX*': そのようなファイルやディレクトリはありません fastq に変換しました。 bowtie でマッピング中... Warning: Could not open read file "SRX3145638_2.fq" for reading; skipping... Error: No input read files were valid (ERR): bowtie2-align exited with value 1 マッピングが完了しました。 samtools でBAM に変換中... [samopen] no @SQ lines in the header. [sam_read1] missing header? Abort! [bam_header_read] EOF marker is absent. The input is probably truncated. awk: コマンドライン:1: (FILENAME=- FNR=1) 致命的: ゼロによる除算が試みられました awk: コマンドライン:1: (FILENAME=- FNR=1) 致命的: ゼロによる除算が試みられました awk: コマンドライン:1: (FILENAME=- FNR=1) 致命的: ゼロによる除算が試みられました BAM に変換しました。 Bed ファイルを作成中... BedGraph に変換中... *****ERROR: Unrecognized parameter: SRX3145638.bam ***** Tool: bedtools genomecov (aka genomeCoverageBed) Version: v2.17.0 Summary: Compute the coverage of a feature file among a genome. Usage: bedtools genomecov [OPTIONS] -i -g Options: -ibam The input file is in BAM format. Note: BAM _must_ be sorted by position -d Report the depth at each genome position (with one-based coordinates). Default behavior is to report a histogram. -dz Report the depth at each genome position (with zero-based coordinates). Reports only non-zero positions. Default behavior is to report a histogram. -bg Report depth in BedGraph format. For details, see: genome.ucsc.edu/goldenPath/help/bedgraph.html -bga Report depth in BedGraph format, as above (-bg). However with this option, regions with zero coverage are also reported. This allows one to quickly extract all regions of a genome with 0 coverage by applying: "grep -w 0$" to the output. -split Treat "split" BAM or BED12 entries as distinct BED intervals. when computing coverage. For BAM files, this uses the CIGAR "N" and "D" operations to infer the blocks for computing coverage. For BED12 files, this uses the BlockCount, BlockStarts, and BlockEnds fields (i.e., columns 10,11,12). -strand Calculate coverage of intervals from a specific strand. With BED files, requires at least 6 columns (strand is column 6). - (STRING): can be + or - -5 Calculate coverage of 5" positions (instead of entire interval). -3 Calculate coverage of 3" positions (instead of entire interval). -max Combine all positions with a depth >= max into a single bin in the histogram. Irrelevant for -d and -bedGraph - (INTEGER) -scale Scale the coverage by a constant factor. Each coverage value is multiplied by this factor before being reported. Useful for normalizing coverage by, e.g., reads per million (RPM). - Default is 1.0; i.e., unscaled. - (FLOAT) -trackline Adds a UCSC/Genome-Browser track line definition in the first line of the output. - See here for more details about track line definition: http://genome.ucsc.edu/goldenPath/help/bedgraph.html - NOTE: When adding a trackline definition, the output BedGraph can be easily uploaded to the Genome Browser as a custom track, BUT CAN NOT be converted into a BigWig file (w/o removing the first line). -trackopts Writes additional track line definition parameters in the first line. - Example: -trackopts 'name="My Track" visibility=2 color=255,30,30' Note the use of single-quotes if you have spaces in your parameters. - (TEXT) Notes: (1) The genome file should tab delimited and structured as follows: For example, Human (hg19): chr1 249250621 chr2 243199373 ... chr18_gl000207_random 4262 (2) The input BED (-i) file must be grouped by chromosome. A simple "sort -k 1,1 > .sorted" will suffice. (3) The input BAM (-ibam) file must be sorted by position. A "samtools sort " should suffice. Tips: One can use the UCSC Genome Browser's MySQL database to extract chromosome sizes. For example, H. sapiens: mysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -e \ "select chrom, size from hg19.chromInfo" > hg19.genome BedGraph に変換しました。 BigWig に変換中... needLargeMem: trying to allocate 0 bytes (limit: 100000000000) BigWig に変換しました。 INFO @ Fri, 04 May 2018 07:07:26: # Command line: callpeak -t SRX3145638.bam -f BAM -g dm -n SRX3145638.10 -q 1e-10 # ARGUMENTS LIST: # name = SRX3145638.10 # format = BAM # ChIP-seq file = ['SRX3145638.bam'] # control file = None # effective genome size = 1.20e+08 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-10 # Larger dataset will be scaled towards smaller dataset. # Range for calculating regional lambda is: 10000 bps # Broad region calling is off # Paired-End mode is off INFO @ Fri, 04 May 2018 07:07:26: # Command line: callpeak -t SRX3145638.bam -f BAM -g dm -n SRX3145638.05 -q 1e-05 # ARGUMENTS LIST: # name = SRX3145638.05 # format = BAM # ChIP-seq file = ['SRX3145638.bam'] # control file = None # effective genome size = 1.20e+08 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-05 # Larger dataset will be scaled towards smaller dataset. # Range for calculating regional lambda is: 10000 bps # Broad region calling is off # Paired-End mode is off INFO @ Fri, 04 May 2018 07:07:26: # Command line: callpeak -t SRX3145638.bam -f BAM -g dm -n SRX3145638.20 -q 1e-20 # ARGUMENTS LIST: # name = SRX3145638.20 # format = BAM # ChIP-seq file = ['SRX3145638.bam'] # control file = None # effective genome size = 1.20e+08 # band width = 300 # model fold = [5, 50] # qvalue cutoff = 1.00e-20 # Larger dataset will be scaled towards smaller dataset. # Range for calculating regional lambda is: 10000 bps # Broad region calling is off # Paired-End mode is off INFO @ Fri, 04 May 2018 07:07:26: #1 read tag files... INFO @ Fri, 04 May 2018 07:07:26: #1 read tag files... INFO @ Fri, 04 May 2018 07:07:26: #1 read treatment tags... INFO @ Fri, 04 May 2018 07:07:26: #1 read tag files... INFO @ Fri, 04 May 2018 07:07:26: #1 read treatment tags... INFO @ Fri, 04 May 2018 07:07:26: #1 read treatment tags... Exception Exception Exception structstructstruct..errorerror: .: error'': uunnp'paacucnkk p arreecqquuikrie rrse esa qsatur isitrnesg aa rsgturmienngt aorfg ulmeenngtt ho f4 'l in e'nMgAtChS 24.'I in O'.MPAaCrSs2e.rI.OB.APMaPrasresre.rB.AtMsPiazres'e ignored r.tsizrei'n ignored g argument of length 4' in 'MACS2.IO.Parser.BAMParser.tsize' ignored INFO @ Fri, 04 May 2018 07:07:26: #1 tag size is determined as 0 bps INFO @ Fri, 04 May 2018 07:07:26: #1 tag size is determined as 0 bps INFO @ Fri, 04 May 2018 07:07:26: #1 tag size is determined as 0 bps INFO @ Fri, 04 May 2018 07:07:26: #1 tag size = 0 INFO @ Fri, 04 May 2018 07:07:26: #1 tag size = 0 INFO @ Fri, 04 May 2018 07:07:26: #1 tag size = 0 INFO @ Fri, 04 May 2018 07:07:26: #1 total tags in treatment: 0 INFO @ Fri, 04 May 2018 07:07:26: #1 total tags in treatment: 0 INFO @ Fri, 04 May 2018 07:07:26: #1 total tags in treatment: 0 INFO @ Fri, 04 May 2018 07:07:26: #1 user defined the maximum tags... INFO @ Fri, 04 May 2018 07:07:26: #1 user defined the maximum tags... INFO @ Fri, 04 May 2018 07:07:26: #1 user defined the maximum tags... INFO @ Fri, 04 May 2018 07:07:26: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 04 May 2018 07:07:26: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 04 May 2018 07:07:26: #1 filter out redundant tags at the same location and the same strand by allowing at most 1 tag(s) INFO @ Fri, 04 May 2018 07:07:26: #1 tags after filtering in treatment: 0 INFO @ Fri, 04 May 2018 07:07:26: #1 tags after filtering in treatment: 0 INFO @ Fri, 04 May 2018 07:07:26: #1 tags after filtering in treatment: 0 Traceback (most recent call last): Traceback (most recent call last): File "/usr/local/bin/macs2", line 5, in File "/usr/local/bin/macs2", line 5, in Traceback (most recent call last): pkg_resources.run_script('MACS2==2.1.1.20160309', 'macs2') pkg_resources.run_script('MACS2==2.1.1.20160309', 'macs2') File "/usr/local/bin/macs2", line 5, in File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 489, in run_script File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 489, in run_script pkg_resources.run_script('MACS2==2.1.1.20160309', 'macs2') File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 489, in run_script File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 1214, in run_script File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 1214, in run_script File "build/bdist.linux-x86_64/egg/pkg_resources.py", line 1214, in run_script File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 617, in File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 617, in File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 617, in File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 57, in main File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 57, in main File "/usr/local/lib/python2.7/site-packages/MACS2-2.1.1.20160309-py2.7-linux-x86_64.egg/EGG-INFO/scripts/macs2", line 57, in main File "build/bdist.linux-x86_64/egg/MACS2/callpeak_cmd.py", line 112, in run File "build/bdist.linux-x86_64/egg/MACS2/callpeak_cmd.py", line 112, in run File "build/bdist.linux-x86_64/egg/MACS2/callpeak_cmd.py", line 112, in run ZeroDivisionError: ZeroDivisionErrorfloat division by zero: float division by zeroZeroDivisionError : float division by zero cat: SRX3145638.05_peaks.narrowPeakcat: SRX3145638.20_peaks.narrowPeak: そのようなファイルやディレクトリはありません : そのようなファイルやディレクトリはありません cat: SRX3145638.10_peaks.narrowPeak: そのようなファイルやディレクトリはありません pass1 - making usageList (0 chroms)pass1 - making usageList (0 chroms)pass1 - making usageList (0 chroms): 2 millis : 2 millis : 2 millis needLargeMem: trying to allocate 0 bytes (limit: 17179869184) needLargeMem: trying to allocate 0 bytes (limit: 17179869184) needLargeMem: trying to allocate 0 bytes (limit: 17179869184) rm: cannot remove `SRX3145638.10_model.r'rm: cannot remove `SRX3145638.20_model.r'rm: cannot remove `SRX3145638.05_model.r': そのようなファイルやディレクトリはありません : そのようなファイルやディレクトリはありません : そのようなファイルやディレクトリはありません rm: rm: cannot remove `SRX3145638.05_*.xls'cannot remove `SRX3145638.10_*.xls': そのようなファイルやディレクトリはありません: そのようなファイルやディレクトリはありません rm: rm: cannot remove `SRX3145638.10_peaks.narrowPeak'cannot remove `SRX3145638.05_peaks.narrowPeak'rm: : そのようなファイルやディレクトリはありません: そのようなファイルやディレクトリはありません cannot remove `SRX3145638.20_*.xls' : そのようなファイルやディレクトリはありません rm: cannot remove `SRX3145638.20_peaks.narrowPeak': そのようなファイルやディレクトリはありません CompletedMACS2peakCalling CompletedMACS2peakCalling CompletedMACS2peakCalling