Entering edit mode
6.1 years ago
amitunited0532
▴
40
[2018-03-17 18:25:31] Beginning TopHat run (v2.1.1)
-----------------------------------------------
[2018-03-17 18:25:31] Checking for Bowtie
Bowtie version: 2.3.2.0
[2018-03-17 18:25:31] Checking for Bowtie index files (genome)..
[2018-03-17 18:25:31] Checking for reference FASTA file
[2018-03-17 18:25:31] Generating SAM header for /media/amit/New Volume/genome/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome
[2018-03-17 18:25:31] Reading known junctions from GTF file
[2018-03-17 18:25:36] Preparing reads
left reads: min. length=50, max. length=2469, 230262 kept reads (8215 discarded)
[2018-03-17 18:26:08] Building transcriptome data files ./tophat_out/tmp/genes
[2018-03-17 18:27:51] Building Bowtie index from genes.fa
[2018-03-17 18:40:12] Mapping left_kept_reads to transcriptome genes with Bowtie2
[2018-03-17 19:31:47] Resuming TopHat pipeline with unmapped reads
[2018-03-17 19:31:47] Mapping left_kept_reads.m2g_um to genome genome with Bowtie2
[FAILED]
Error running bowtie:
Warning: Output file '-' was specified without -S. This will not work in future Bowtie 2 versions. Please use -S instead.
stat: No such file or directory
Warning: Could not open read file "Volume/genome/Homo_sapiens/UCSC/hg19/Sequence/Bowtie2Index/genome" for reading; skipping...
Error: No input read files were valid
(ERR): bowtie2-align exited with value 1
Here we go again:
You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. (If you can't get access to that publication, let me know and I'll -cough- help you.) There are also other alternatives, including alignment with STAR and bbmap, or pseudo-alignment using salmon.
That was already conveyed half-way up this thread: C: Error running bowtie
Oh right, I see. While I understand you and OP like to solve this issue, there is no much use in spending time on it, is there?
After running HiSAT2 its showing error like: (ERR): hisat2-align exited with value 1
Then wouldn't you rather try to solve the HISAT error? When my car breaks down I'm not getting a horse, I fix the car.
Please provide the full command you are using for this run.
There seem to be multiple issues. You may be missing a leading
/before the Volume in your genome index path specification. You are not providing correct path for the fastq files and/or they may not be in right format.I have all file in external disk except fastq file, Is it rit for run this command?
Before we go forward consider using some other program that is more current compared to TopHat. STAR/bbmap etc. TopHat is no longer recommended for use even by its authors.
You have a space in the file paths that is causing this problem. You could try following command (or eliminate the space in the path name and try again).
It iS not working.
Rename
New VolumetoNew_Volumeand try with the new name minus the space in name. You know for sure that the file is present in that directory?Again it showing same error.
what does this command show?
and try this next
Witthout "_" showing path of directory but not with "_"
With New_Volume
Without New Volume
You did not rename the directory like I suggested you to. Do
they try tophat command I gave above with the new name without the space in New Volume..