This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How many read files I can use in tophat in the same run?

Good morning friends,

Supposing I have and experiment with 10 samples and each with just sra file. In tophat -r 20 test_ref reads_1.fq reads_2.fq syntax I should mention all of my fastq files? If yes, may you please tell me the order of files. I mean for example I should type reads_1.fq -> reads_10.fq?

Thanks a lot

rna-seq tophat
• 2,726 views • link
updated by Ram Ram • written by zizigolu zizigolu

Run tophat on each sample separately (or write a simple bash script with a list of fastq file names as input and specify different output directory for each of them). In case of paired-end reads_1.fq and reads_2.fq defines for Right and Left.

• 0 views • link
updated by Ram Ram • written by Anushka • Anushka

1 answer

Each sample must be aligned independently. You can not give tophat all of the files for all of the samples at one time and get multiple output files.

• 0 views • link
written by Devon Ryan Devon Ryan

Thanks Devon,

Sorry

I placed the orf_coding.fasta (my reference genome) in the bowtie2 file and indexed the genome in the same bowtie2, and run this syntax but error,

May you please tell me a solution?

[izadi@lbox161 bowtie2-2.2.5]$ $TOP/tophat -r 20 orf_coding.fasta SRR1944914_trimmed.fastq SRR1944926_trimmed.fastq

[2015-09-01 10:27:50] Beginning TopHat run (v2.1.0)
-----------------------------------------------
[2015-09-01 10:27:50] Checking for Bowtie
          Bowtie version:     2.2.5.0
[2015-09-01 10:27:50] Checking for Bowtie index files (genome)..
Error: Could not find Bowtie 2 index files (orf_coding.fasta.*.bt2)
[izadi@lbox161 bowtie2-2.2.5]$ $TOP/tophat -r 20 orf_coding.fasta SRR1944914_trimmed_unmapped.fastq SRR1944926_trimmed_unmapped.fastq

[2015-09-01 10:29:01] Beginning TopHat run (v2.1.0)
-----------------------------------------------
[2015-09-01 10:29:01] Checking for Bowtie
          Bowtie version:     2.2.5.0
[2015-09-01 10:29:01] Checking for Bowtie index files (genome)..
Error: Could not find Bowtie 2 index files (orf_coding.fasta.*.bt2)
[izadi@lbox161 bowtie2-2.2.5]$
• 0 views • link
updated by Ram Ram • written by zizigolu zizigolu

This part of the error message should tell you what to fix:

Error: Could not find Bowtie 2 index files (orf_coding.fasta.*.bt2)
• 0 views • link
updated by Ram Ram • written by Devon Ryan Devon Ryan

Sorry, but after indexing there are coding.1.bt2, coding.2.bt2, coding.3.bt2, coding.4.bt2, coding.rev.1.bt2 and coding.rev.2.bt2.

and orf_coding.fasta in bowtie2 file. Then what I should provide?

• 0 views • link
updated by Ram Ram • written by zizigolu zizigolu

Then you need to write coding rather than orf_coding.fasta.

• 0 views • link
written by Devon Ryan Devon Ryan

Thanks a lot Devon, you saved me I was getting mad since yesterday.

• 0 views • link
updated by Ram Ram • written by zizigolu zizigolu

Log in to answer this question.

More posts like this
  • Tophat and Bowtie installed with Conda not working on test data.
    written by ariel.balter ariel.balter
    I created a python 2.7 environment using conda, and installed tophat and bowtie. I'm trying the test data as specified here: https://ccb.jhu.edu/software/tophat/tutorial.shtml#inst However, I get …
  • run softwares tophat
    written by zh.khodadadi • zh.khodadadi
    i am a begginer user ,i install tophat2,bowtie2,fastqc and cufflink on desktop my 14.04 ubuntu by same these commands: tar -zxvf tophat-2.1.1.Linux_x86_64.tar.gz cd tophat-2.1.1.Linux_x86_64 ./tophat2 …
  • Error: segment-based junction search failed with err =-5
    written by ihabyounis • ihabyounis
    Hello, I am new to using tophat and bowtie for RNA-seq data. I am using MAC OSX10.11.4. I downloaded and unpacked the tophat2 binaries (ophat-2.1.1.OSX_x86_64). …
  • paired-end alignment with tophat2
    written by zizigolu zizigolu
    Hi, I have 12 paired-end fq files. In tophat2 protocol this is the syntax tophat -p 8 -G genes.gtf genome file-1.fq file-2.fq but I read …
  • how I can put bowtie2 with tophat2 in the path?
    written by zizigolu zizigolu
    Hey guys, ``` [izadi@lbox161 ~]$ bash [izadi@lbox161 ~]$ cd /usr/data/nfs6/izadi/gij/test_data [izadi@lbox161 test_data]$ pwd /usr/data/nfs6/izadi/gij/test_data [izadi@lbox161 test_data]$ ls reads_1.fq test_ref.2.bt2 test_ref.fa test_ref.rev.2.ebwt reads_2.fq test_ref.2.ebwt test_ref.rev.1.bt2 test_ref.1.bt2 …
  • Testing Tophat installation
    written by erikfas • erikfas
    I'm just starting out with RNA-seq, and I'm trying to get the installation of Tophat to work. I'm following the "Getting started" instructions on the …
  • tophat Preparing reads [FAILED]
    written by justinfanyi • justinfanyi
    Why is tophat giving this error message? biostat1@biostat1-ubuntu:~/src/tools/test_data$ tophat2 -r 20 test_ref reads_1.fq reads_2.fq [2014-05-24 17:02:09] Beginning TopHat run (v2.0.11) ----------------------------------------------- [2014-05-24 17:02:09] Checking for …
  • Rna-Seq Paired End 76 Bp Read Data From Sra To Fastq To Bam
    written by dfernan dfernan
    <p>Hi,</p> <p>I was wondering if anyone has written a wrapper to process SRA files from SRA to fastq, then do some QC on the fastq …
  • Tophat Prep_Reads Error
    written by Mohammad Reza Bakhtiarizadeh Mohammad Reza Bakhtiarizadeh
    <p>Hi I am running tophat on a test reads (<a href='http://tophat.cbcb.umd.edu/downloads/test_data.tar.gz'>Test data in tophat site</a> and the code was tophat -r 20 test_ref reads_1.fq reads_2.fq) …
  • [Errno 8] Exec Format Error Using Tophat 1.3.1
    written by Marina Manrique Marina Manrique
    <p>Hi all,</p> <p>I've installed Tophat 1.3.1 following <a href='http://tophat.cbcb.umd.edu/tutorial.html#inst'>these instructions</a> </p> <p>and when I run the test command to know if the installation has gone …