Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • James
    Member
    • Mar 2010
    • 23

    another tophat "could not execute prep_reads" error

    Hi,

    Im getting this error:
    Code:
    tophat -o RNA-seq-14.11.10/tophat_out_0H -G maygtf.gtf --no-novel-juncs --solexa1.3-quals indexes/_May_2009-masked reads/RNA-seq/R43s_1_sequence_0H.fastq
    
    
    [Mon Nov 15 12:52:59 2010] Beginning TopHat run (v1.1.0)
    -----------------------------------------------
    [Mon Nov 15 12:52:59 2010] Preparing output location RNA-seq-14.11.10/tophat_out_0H//
    [Mon Nov 15 12:52:59 2010] Checking for Bowtie index files
    [Mon Nov 15 12:52:59 2010] Checking for reference FASTA file
    [Mon Nov 15 12:52:59 2010] Checking for Bowtie
    	Bowtie version:			 0.12.7.0
    [Mon Nov 15 12:52:59 2010] Checking for Samtools
    	Samtools version:		 0.1.9.0
    [Mon Nov 15 12:52:59 2010] Checking reads
    	min read length: 36bp, max read length: 36bp
    	format:		 fastq
    	quality scale:	 phred64 (reads generated with GA pipeline version >= 1.3)
    [Mon Nov 15 12:55:49 2010] Reading known junctions from GTF file
    	[FAILED]
    Error: could not execute prep_reads
    I check the log file for prep_reads it has:

    prep_reads v1.1.3 (1680)
    ---------------------------
    prep_reads: unrecognized option `--gtf-annotations'

    Anybody had this error? or can help? Thanks James
  • SongLi
    Member
    • Oct 2010
    • 19

    #2
    I got the same error for my dataset today, without those options:

    tophat -o Mytophat_out -G mygtf.gtf indexes/myindex reads/myreads.fasta

    The error message is :

    [FAILED]
    Error: could not execute prep_reads


    and in the prep_reads.log is:

    prep_reads v1.1.3 (1680)
    ---------------------------
    somedir/prep_reads: unrecognized option `--gtf-annotations'


    Hope someone can help.


    Originally posted by James View Post
    Hi,

    Im getting this error:
    Code:
    tophat -o RNA-seq-14.11.10/tophat_out_0H -G maygtf.gtf --no-novel-juncs --solexa1.3-quals indexes/_May_2009-masked reads/RNA-seq/R43s_1_sequence_0H.fastq
    
    
    [Mon Nov 15 12:52:59 2010] Beginning TopHat run (v1.1.0)
    -----------------------------------------------
    [Mon Nov 15 12:52:59 2010] Preparing output location RNA-seq-14.11.10/tophat_out_0H//
    [Mon Nov 15 12:52:59 2010] Checking for Bowtie index files
    [Mon Nov 15 12:52:59 2010] Checking for reference FASTA file
    [Mon Nov 15 12:52:59 2010] Checking for Bowtie
    	Bowtie version:			 0.12.7.0
    [Mon Nov 15 12:52:59 2010] Checking for Samtools
    	Samtools version:		 0.1.9.0
    [Mon Nov 15 12:52:59 2010] Checking reads
    	min read length: 36bp, max read length: 36bp
    	format:		 fastq
    	quality scale:	 phred64 (reads generated with GA pipeline version >= 1.3)
    [Mon Nov 15 12:55:49 2010] Reading known junctions from GTF file
    	[FAILED]
    Error: could not execute prep_reads
    I check the log file for prep_reads it has:

    prep_reads v1.1.3 (1680)
    ---------------------------
    prep_reads: unrecognized option `--gtf-annotations'

    Anybody had this error? or can help? Thanks James
    Last edited by SongLi; 11-16-2010, 11:21 AM.

    Comment

    • James
      Member
      • Mar 2010
      • 23

      #3
      I fixed this problem. It may have been two things:

      It may have been calling an older prep_reads I had lurking in my path.

      It said I was calling tophat v1.1.0 even though I downloaded v1.1.3. I just re-downloaded the binaries and replaced them. Now it tells me I'm running v1.1.3.

      I'm running Mac OSX 10.5.8. by the way.

      Comment

      • SongLi
        Member
        • Oct 2010
        • 19

        #4
        Thank you James,

        However, it may not be the reason, because in the log file you posted, it says v1.1.3 there.

        I still have the problem on my side.



        Originally posted by James View Post
        I fixed this problem. It may have been two things:

        It may have been calling an older prep_reads I had lurking in my path.

        It said I was calling tophat v1.1.0 even though I downloaded v1.1.3. I just re-downloaded the binaries and replaced them. Now it tells me I'm running v1.1.3.

        I'm running Mac OSX 10.5.8. by the way.

        Comment

        • SongLi
          Member
          • Oct 2010
          • 19

          #5
          Hi James,

          Since you have solved the problem, would you mind post what the messages you get when the read RTF file passed?

          On my side the error is:
          [Tue Nov 16 17:02:17 2010] Reading known junctions from GTF file
          [FAILED]
          Error: could not execute prep_reads


          That may help me to figure out what's wrong in my case.
          Thanks,



          Originally posted by James View Post
          I fixed this problem. It may have been two things:

          It may have been calling an older prep_reads I had lurking in my path.

          It said I was calling tophat v1.1.0 even though I downloaded v1.1.3. I just re-downloaded the binaries and replaced them. Now it tells me I'm running v1.1.3.

          I'm running Mac OSX 10.5.8. by the way.

          Comment

          • oliviera
            Member
            • Apr 2010
            • 31

            #6
            I have the same problem with tophat 1.1.3, and do not manage to solve the problem.

            I used :
            ./tophat --solexa1.3-quals --no-novel-juncs -G Danio_rerio.Zv9.60.gtf -o /5cq15 Danio_rerio.Zv9.60 5.fastq

            the log file is the same as yours:
            prep_reads v1.1.3 (1680)
            ---------------------------
            /home/olivier/src/illumina_software_version/tophat_113/bin/prep_reads: unrecognized option '--gtf-annotations'

            Funny thing however in the gtf_juncs log file: seems to do something their
            gtf_juncs v1.1.3 (1680)
            ---------------------------
            Extracted 218405 junctions from /Danio_rerio.Zv9.60.gtf

            I tried gtf and gff format and both failed.
            When I run without -G option it works fine...
            Any feedback would be much appreciated

            Cheers

            Oliviera

            Comment

            • oliviera
              Member
              • Apr 2010
              • 31

              #7
              I just tested the new release 1.1.4 and now the problem seems to be fixed!
              Thanx a lot


              Oliviera

              Comment

              • SongLi
                Member
                • Oct 2010
                • 19

                #8
                Thank you TopHat team for the timingly fix!

                Don't know where to post, so I just post here.

                Originally posted by oliviera View Post
                I just tested the new release 1.1.4 and now the problem seems to be fixed!
                Thanx a lot


                Oliviera

                Comment

                Latest Articles

                Collapse

                • SEQadmin2
                  Cancer Drug Resistance: The Lingering Barrier to Rising Survival
                  by SEQadmin2



                  Cancer survival rates have significantly increased in the last few decades in the United States, reaching a combined 70% 5-year survival rate by 2021. Behind this number, there are years of research to find new therapies, drug targets, and early detection methods. But there is one core challenge that keeps slowing down these advances, and it’s about drug resistance.

                  There is no single reason why many patients don’t respond to treatment as expected. Cancer is...
                  Today, 05:17 AM
                • GATTACAT
                  Reply to Nine Things a Sample Prep Scientist Thinks About Before Sequencing
                  by GATTACAT
                  Love this - good data definitely starts from good input, and poor input can only give relatively poor data. I particularly like the mention of Nanodrop/absorbance based methods for quantification. It's such a toss up if you'll get an accurate reading or what amounts to a randomly generated number, and a lot of library/sequencing related issues can be traced back to poor quant.
                  07-01-2026, 11:43 AM
                • SEQadmin2
                  Nine Things a Sample Prep Scientist Thinks About Before Sequencing
                  by SEQadmin2


                  I’m not a sequencing expert. I’m a purification scientist who uses NGS to evaluate workflows my group develops. With this perspective, we think about the sample first and the NGS workflow second. The sequencer is an exceptionally honest reporter, but it can only report on what you give it, so whether you get clean, interpretable data from an NGS workflow is largely determined before you begin.

                  Here are nine questions we think about, in roughly the order they matter, before...
                  06-18-2026, 07:11 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by SEQadmin2, Today, 10:08 AM
                0 responses
                6 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, Yesterday, 11:05 AM
                0 responses
                8 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 07-02-2026, 11:08 AM
                0 responses
                31 views
                0 reactions
                Last Post SEQadmin2  
                Started by SEQadmin2, 06-30-2026, 05:37 AM
                0 responses
                29 views
                0 reactions
                Last Post SEQadmin2  
                Working...