Error in Bowtie2

Hi The script I have is

tophat -p 4 --b2-very-sensitive -G '/home/rasoul/Desktop/Data/ngs/genome/annotation.gtf' /Home/rasoul/Desktop/Data/ngs/genome/genome '/Home/rasoul/Desktop/Data/ngs/SRR941816/keep16.gz' 



Error
[2020-04-05 18:15:48] Beginning TopHat run (v2.1.1)
-----------------------------------------------
[2020-04-05 18:15:48] Checking for Bowtie
          Bowtie version:    2.2.5.0
[2020-04-05 18:15:48] Checking for Bowtie index files (genome)..
Error: Could not find Bowtie 2 index files (/Home/rasoul/Desktop/Data/ngs/genome/genome.*.bt2l)

rna-seq

• 844 views • link

updated 6.3 years ago by Kevin Blighe

• written 6.3 years ago by rasoulfarzaneh •

If you are trying to run a transcriptome only run with -G option then you are not supposed to provide reads for that run.

For example, in order to just prepare the transcriptome index files for a specific annotation, an initial, single TopHat run could be invoked like this:
tophat -G known_genes.gtf \
    --transcriptome-index=transcriptome_data/known \
    hg19

It also looks like you need to use /Home/rasoul/Desktop/Data/ngs/genome/ as the index base.

You should also consider switching to a newer aligner than TopHat.

• 0 views • link

written 6.3 years ago by GenoMax

0 answers

No answers yet.

More posts like this

tophat2 runs with error:: segment-based junction search failed with err =-6
written 4.0 years ago by pxd •

The command I use is : tophat -p 8 -r 101 -G gencode.v41.annotation.gtf -o GM12878_1 pxd ENCFF001REJ.fastq ENCFF001REK.fastq [2022-08-18 17:30:19] Beginning TopHat run (v2.1.1) ----------------------------------------------- …
Error: Could not find Bowtie 2 index files (genome_index_base.*.bt2l)
written 4.3 years ago by rishav513 •

when trying to map my reads to the index file, it shows an error. please help > $tophat2 genome_index_base reads1.fq reads2.fq [2022-04-14 12:41:18] Beginning TopHat …
Bowtie 2 index files exist but return Error: Could not find Bowtie 2 index files (humanGRCh38.*.bt2…
written 7.1 years ago by skurkure •

As I have stated above, I have received an error in running tophat, with the output below. The command I supplied was tophat humanGRCh38 gerald_C2DBEACXX_3.fq,gerald_C1TD1ACXX_8_ACAGTG.fq …
Couldn't Build Bowtie index with err = 1
written 8.4 years ago by nathan.yeung •

I'm trying to run tophat2 on RNA-seq data but it can't be accomplished as it continues to give the error message of `Couldn't Build Bowtie …
Error: Could not find Bowtie 2 index files (GRCh38.*.bt2l)
written 9.8 years ago by ashMC •

I am running tophat/2.1.1 and using bowtie/2.2.9. I run the following commands, and get the following errors. tophat -G reference_files/Homo_sapiens.GRCh38.86.gtf -p4 -o tophat_results/ reference_files/GRCh38 data/S_01_L003_R1.fastq …
error in tophat2
written 9.9 years ago by zizigolu

hi, I have already built index for genome fasta by **bowtie2-build genome.fa aspergillus.fumigatus** but I don't know why tophat can't find them :( [izadi@lbox200 bowtie2-2.2.9]$ …
Unable to find Bowtie2 index
written 10.2 years ago by EpiExplorer

Hi, I am unable to run Tophat2 as I get an error. Here is the command I run: tophat2 -p 5 -r 62 –library-type fr-firststrand …
tophat2 Error: segment-based junction search failed with err =-9
written 10.2 years ago by marongiu.luigi

Dear all, I am aligning HISeq RNA-Seq data to the human reference genome. I downloaded the Homo_sapiens.GRCh38.dna.toplevel.fa.gz and the Homo_sapiens.GRCh38.84.gtf.gz reference files and I built …
Tophat syntax to align query sequences to reference index
written 10.3 years ago by marongiu.luigi

Dear all,<br> I need to align some fastq files to the human genome. I prepared the reference index files and transcriptomes with bowtie and tophat …
Error: Could not find Bowtie 2 index files (p.*.bt2)
written 10.6 years ago by sangita_b

Hi, I am repeating alignments using TopHat, I am using the same Bowtie index as I used previously (i.e I have successfully aligned reads using …