Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • smleighton
    Junior Member
    • Feb 2012
    • 3

    #1

    tophat error with large file

    I am getting errors running large rna files. If the file is larger than about 140 GB then I get the following error. If I split the files up and run each individually then it works fine. Does anyone know if there is a limit to the size of input to the tophat application?

    [Wed Mar 21 00:07:31 2012] Beginning TopHat run for step PreProcess (v1.3.1)
    -----------------------------------------------
    [Wed Mar 21 00:07:31 2012] Preparing output location ./
    [Wed Mar 21 00:07:32 2012] Checking for Bowtie index files
    [Wed Mar 21 00:07:32 2012] Checking for reference FASTA file
    [Wed Mar 21 00:07:32 2012] Checking for Bowtie
    Bowtie version: 0.12.7.0
    [Wed Mar 21 00:07:32 2012] Checking for Samtools
    Samtools Version: 0.1.18
    [Wed Mar 21 00:07:32 2012] Generating SAM header for /hg19/data/hg19
    format: fastq
    quality scale: solexa33 (reads generated with GA pipeline version < 1.3)
    [Wed Mar 21 00:07:54 2012] Preparing reads
    [FAILED]
    Error running 'prep_reads'
    Error: beginning of quality values record not found! (@CCFFFDDHHHHHJJJJJJHIJJIJIGIIIGJIJJIIIIIGJDHIJJJJIJ>F=FHHIIIIHHHGHHFFDDDEEECEEDCDEDDDDCDDDDBDDDDDCCC)
  • aoifemcm
    Junior Member
    • Oct 2012
    • 3

    #2
    I'm having exactly the same problem. Anyone have a suggestion?

    I'm using tophat v2.0.10.

    Comment

    • WhatsOEver
      Senior Member
      • Apr 2012
      • 215

      #3
      Did you concatenate the files somehow? The fastq format is normally easy to parse as a new entry has to have 4 lines (header, seq, def, qual). If this order gets somehow corrupted, such errors can appear (It looks to me like tophat identifies the "@CCFFF..." as the header, but its actually the quality line...).

      what do you get from the following command:

      Code:
      grep "@CCFFFDDHHHHHJJJJJJHIJJIJIGIIIGJIJJIIIIIGJDHIJJJJIJ>F=FHHIIIIHHHGHHFFDDDEEECEEDCDEDDDDCDDDDBDDDDDCCC" -A 4 -B 4 ./yourFastqFile.fastq
      [EDIT]
      @aoifemcm: you need to exchange "@CCF..." with what you have in the error, of course
      Last edited by WhatsOEver; 10-02-2014, 12:19 AM.

      Comment

      • aoifemcm
        Junior Member
        • Oct 2012
        • 3

        #4
        Yes, these files were concatenated. Tophat ran perfectly well on each constituent file separately, and when I redid the concatenation it solved the problem.
        I guess corruption was at fault.

        Thanks!

        Comment

        • WhatsOEver
          Senior Member
          • Apr 2012
          • 215

          #5
          Though its good to hear you solved your problem, I would be very interested to hear how you concatenated the files and if the error is reproducible.

          Comment

          Latest Articles

          Collapse

          • SEQadmin2
            Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
            by SEQadmin2



            Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

            Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
            ...
            10-14-2026, 09:55 AM
          • SEQadmin2
            New Genomics Technologies Take Aim at Long-Standing Limits
            by SEQadmin2


            Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

            We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
            ...
            09-28-2026, 10:25 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 10-22-2026, 11:13 AM
          0 responses
          35 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 10:48 AM
          0 responses
          25 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-15-2026, 10:52 AM
          0 responses
          41 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-14-2026, 05:39 AM
          0 responses
          45 views
          0 reactions
          Last Post SEQadmin2  
          Working...