Question: error in tophat2
0
gravatar for F
22 months ago by
F3.1k
Iran
F3.1k wrote:

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]$ $TOP/tophat -p 8 -G Aspergillus_fumigatus.CADRE.32.gtf -o 1_thout genome.fa clean.fq

[2016-09-03 01:01:10] Beginning TopHat run (v2.1.1)
-----------------------------------------------
[2016-09-03 01:01:10] Checking for Bowtie

          Bowtie version:    2.2.9.0

[2016-09-03 01:01:10] Checking for Bowtie index files (genome)..

Error: Could not find Bowtie 2 index files (genome.fa.*.bt2l)

[izadi@lbox200 bowtie2-2.2.9]$

helpless...

rna-seq software error • 927 views
ADD COMMENTlink modified 22 months ago by geek_y8.6k • written 22 months ago by F3.1k
2
gravatar for Satyajeet Khare
22 months ago by
Satyajeet Khare1.3k
Pune, India
Satyajeet Khare1.3k wrote:

Build genome

bowtie2-build -f genome_bowtie2_build/genome.fa genome_bowtie2_build/genome

That is, keep the name of build files same as that of fasta file and preferably keep them in the same location.

Run tophat. genome_bowtie2_build/genome is the location of genome build files (try removing .fa from your command).

tophat2 -p 10 -o tophat2/sample1/ genome_bowtie2_build/genome sample1.fastq

If you want to provide gtf file, build transcriptome. Genome build location is also provided in the command.

tophat -G basic.annotation.gtf --transcriptome-index=./transcriptome_data/basic.annotation genome_bowtie2_build/genome

and then run tophat with gtf file

tophat2 -p 10 -G basic.annotation.gtf --transcriptome-index transcriptome_data/basic.annotation -o tophat2/Sample1 genome_bowtie2_build/genome sample1.fastq

But if you starting, may be you should start with Hisat.

ADD COMMENTlink modified 22 months ago • written 22 months ago by Satyajeet Khare1.3k

thank uuuuuuuuuuuuuuuuuuuuuuuuuuuu extend to the world u saved me from confusion

ADD REPLYlink written 22 months ago by F3.1k

Hi Satyajeet, From where I can get the basic.annotation.gtf file. Thanks in advance for your help.

ADD REPLYlink written 9 months ago by k.kathirvel9350
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 561 users visited in the last hour