Hi all I am running tophat2 on my proton fastq data reads and have this error which I haven't faced before.
Could you tell me what happened?!
(Sorry for trouble every one , the problem has been sorted! because I didn't have enough space......)
My command:
tophat2 -o folder1 -p 12 --library-type fr-firststrand --keep-fasta-order --GTF genes.gtf Bowtie2Index sample1.fastq
The error is:
[2014-04-30 14:28:28] Beginning TopHat run (v2.0.10)
-----------------------------------------------
[2014-04-30 14:28:28] Checking for Bowtie
Bowtie version: 2.1.0.0
[2014-04-30 14:28:28] Checking for Samtools
Samtools version: 0.1.19.0
[2014-04-30 14:28:28] Checking for Bowtie index files (genome)..
[2014-04-30 14:28:28] Checking for reference FASTA file
[2014-04-30 14:28:28] Generating SAM header for ./Bos_taurus/Ensembl/Btau_4.0/Sequence/Bowtie2Index/genome
[2014-04-30 14:28:30] Reading known junctions from GTF file
[2014-04-30 14:28:33] Preparing reads
left reads: min. length=16, max. length=320, 65134842 kept reads (3312 discarded)
Warning: short reads (<20bp) will make TopHat quite slow and take large amount of memory because they are likely to be mapped in too many places
[2014-04-30 14:43:49] Building transcriptome data files folder1/tmp/genes
[2014-04-30 14:43:57] Building Bowtie index from genes.fa
[2014-04-30 14:45:17] Mapping left_kept_reads to transcriptome genes with Bowtie2
[2014-04-30 15:00:59] Resuming TopHat pipeline with unmapped reads
[2014-04-30 15:00:59] Mapping left_kept_reads.m2g_um to genome genome with Bowtie2
Traceback (most recent call last):
File "/software/bin/tophat", line 4084, in <module>
sys.exit(main())
File "/software/bin/tophat", line 4050, in main
user_supplied_deletions)
File "/software/bin/tophat", line 3541, in spliced_alignment
segment_len)
File "/software/bin/tophat", line 2977, in split_reads
split_record(read_name, read_seq, read_quals, out_segfiles, offsets, params.read_params.color)
File "/software/bin/tophat", line 2925, in split_record
print >> f, "%s|%d:%d:%d" % (read_name,last_seq_offset,seg_num, len(offsets) - 1)
IOError: [Errno 32] Broken pipe
Could you tell me what happened?!
(Sorry for trouble every one , the problem has been sorted! because I didn't have enough space......)
My command:
tophat2 -o folder1 -p 12 --library-type fr-firststrand --keep-fasta-order --GTF genes.gtf Bowtie2Index sample1.fastq
The error is:
[2014-04-30 14:28:28] Beginning TopHat run (v2.0.10)
-----------------------------------------------
[2014-04-30 14:28:28] Checking for Bowtie
Bowtie version: 2.1.0.0
[2014-04-30 14:28:28] Checking for Samtools
Samtools version: 0.1.19.0
[2014-04-30 14:28:28] Checking for Bowtie index files (genome)..
[2014-04-30 14:28:28] Checking for reference FASTA file
[2014-04-30 14:28:28] Generating SAM header for ./Bos_taurus/Ensembl/Btau_4.0/Sequence/Bowtie2Index/genome
[2014-04-30 14:28:30] Reading known junctions from GTF file
[2014-04-30 14:28:33] Preparing reads
left reads: min. length=16, max. length=320, 65134842 kept reads (3312 discarded)
Warning: short reads (<20bp) will make TopHat quite slow and take large amount of memory because they are likely to be mapped in too many places
[2014-04-30 14:43:49] Building transcriptome data files folder1/tmp/genes
[2014-04-30 14:43:57] Building Bowtie index from genes.fa
[2014-04-30 14:45:17] Mapping left_kept_reads to transcriptome genes with Bowtie2
[2014-04-30 15:00:59] Resuming TopHat pipeline with unmapped reads
[2014-04-30 15:00:59] Mapping left_kept_reads.m2g_um to genome genome with Bowtie2
Traceback (most recent call last):
File "/software/bin/tophat", line 4084, in <module>
sys.exit(main())
File "/software/bin/tophat", line 4050, in main
user_supplied_deletions)
File "/software/bin/tophat", line 3541, in spliced_alignment
segment_len)
File "/software/bin/tophat", line 2977, in split_reads
split_record(read_name, read_seq, read_quals, out_segfiles, offsets, params.read_params.color)
File "/software/bin/tophat", line 2925, in split_record
print >> f, "%s|%d:%d:%d" % (read_name,last_seq_offset,seg_num, len(offsets) - 1)
IOError: [Errno 32] Broken pipe
Comment