Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • repinementer
    Member
    • Dec 2009
    • 80

    #1

    Tophat error -Junction database is empty

    I built hg18 indexes in bowtie and I also placed those files in testdata folder
    Could any one please help how to get rid off this???

    Thanx

    Code:
    [Mon Aug  2 21:20:14 2010] Searching for junctions via segment mapping
    Warning: junction database is empty!
    [Mon Aug  2 21:21:59 2010] Joining segment hits
    [Mon Aug  2 21:21:59 2010] Joining segment hits
    [Mon Aug  2 21:21:59 2010] Reporting output tracks
            [FAILED]
    Error: Report generation failed with err = 1
  • repinementer
    Member
    • Dec 2009
    • 80

    #2
    May be I should post the whole error

    I tried with and with out segment length advanced option bylooking at your answers in previous thread. but no luck!

    Code:
    [B]
    $ tophat -r 156 hg18_c SRR015270_1.fastq SRR015271_1.fastq [/B][B]SRR015272_1.fastq SRR015273_1.fastq [/B]
    
    [Mon Aug  2 21:40:27 2010] Beginning TopHat run (v1.0.13)
    -----------------------------------------------
    [Mon Aug  2 21:40:27 2010] Preparing output location ./tophat_out/
    [Mon Aug  2 21:40:27 2010] Checking for Bowtie index files
    [Mon Aug  2 21:40:27 2010] Checking for reference FASTA file
    [Mon Aug  2 21:40:27 2010] Checking for Bowtie
            Bowtie version:          0.12.5.0
    [Mon Aug  2 21:40:27 2010] Checking reads
            seed length:     32bp
            format:          fastq
            quality scale:   phred33 (default)
    [Mon Aug  2 21:41:29 2010] Mapping reads against hg18_c with Bowtie
    [Mon Aug  2 21:41:29 2010] Joining segment hits
    [Mon Aug  2 21:41:29 2010] Mapping reads against hg18_c with Bowtie
    [Mon Aug  2 21:41:29 2010] Joining segment hits
    [Mon Aug  2 21:41:29 2010] Searching for junctions via segment mapping
    Warning: junction database is empty!
    [Mon Aug  2 21:43:13 2010] Joining segment hits
    [Mon Aug  2 21:43:14 2010] Joining segment hits
    [Mon Aug  2 21:43:14 2010] Reporting output tracks
            [FAILED]
    Error: Report generation failed with err = 1
    Traceback (most recent call last):
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1635, in ?
        sys.exit(main())
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1607, in main
        params.gff_annotation)
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1044, in compile_reports
        exit(1)
    TypeError: 'str' object is not callable
    Code:
    $tophat -r 156 --segment-length 32 hg18_c SRR015270_1.fastq SRR015271_1.fastq SRR015272_1.fastq SRR015273_1.fastq 
    
    [Mon Aug  2 22:04:45 2010] Beginning TopHat run (v1.0.13)
    -----------------------------------------------
    [Mon Aug  2 22:04:45 2010] Preparing output location ./tophat_out/
    [Mon Aug  2 22:04:45 2010] Checking for Bowtie index files
    [Mon Aug  2 22:04:45 2010] Checking for reference FASTA file
    [Mon Aug  2 22:04:45 2010] Checking for Bowtie
            Bowtie version:          0.12.5.0
    [Mon Aug  2 22:04:45 2010] Checking reads
            seed length:     32bp
            format:          fastq
            quality scale:   phred33 (default)
    [Mon Aug  2 22:05:47 2010] Mapping reads against hg18_c with Bowtie
    [Mon Aug  2 22:05:47 2010] Joining segment hits
    [Mon Aug  2 22:05:47 2010] Mapping reads against hg18_c with Bowtie
    [Mon Aug  2 22:05:47 2010] Joining segment hits
    [Mon Aug  2 22:05:47 2010] Searching for junctions via segment mapping
    Warning: junction database is empty!
    [Mon Aug  2 22:07:32 2010] Joining segment hits
    [Mon Aug  2 22:07:32 2010] Joining segment hits
    [Mon Aug  2 22:07:32 2010] Reporting output tracks
            [FAILED]
    Error: Report generation failed with err = 1
    Traceback (most recent call last):
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1635, in ?
        sys.exit(main())
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1607, in main
        params.gff_annotation)
      File "/home/bogugk/rnaseq/software/tophat-1.0.13/bin/tophat", line 1044, in compile_reports
        exit(1)
    TypeError: 'str' object is not callable
    segment_junctionLogs
    Code:
    [bogugk@Trivia logs]$ cat segment_juncs.log 
    segment_juncs v1.0.13
    ---------------------------
    Loading reference sequences...
            Loading chr1...done
            Loading chr2...done
            Loading chr3...done
            Loading chr4...done
            Loading chr5...done
            Loading chr6...done
            Loading chr7...done
            Loading chr8...done
            Loading chr9...done
            Loading chr10...done
            Loading chr11...done
            Loading chr12...done
            Loading chr13...done
            Loading chr14...done
            Loading chr15...done
            Loading chr16...done
            Loading chr17...done
            Loading chr18...done
            Loading chr19...done
            Loading chr20...done
            Loading chr21...done
            Loading chr22...done
            Loading chrX...done
            Loading chrY...done
            Loading chrM...done
    Found 0 potential split-segment junctions
    Indexing extensions in ./tophat_out/tmp//left_kept_reads_missing.fq
    Can't open file ./tophat_out/tmp//left_kept_reads_missing.fq for reading, skipping...
    Indexing extensions in ./tophat_out/tmp//right_kept_reads_missing.fq
    Can't open file ./tophat_out/tmp//right_kept_reads_missing.fq for reading, skipping...
    Looking for junctions by island end pairings
    Adding hits from segment file 0 to coverage map
    Adding hits from segment file 1 to coverage map
    Map covers 0 bases
    Map covers 0 bases in sufficiently long segments
    Map contains 1 good islands
    0 are left looking bases
    0 are right looking bases
    Collecting potential splice sites in islands
    reporting synthetic splice junctions...
    Found 0 potential island-end pairing junctions
    done
    Reporting potential splice junctions...done
    Reported 0 total possible splices
    Last edited by repinementer; 08-02-2010, 06:10 AM. Reason: Added another command output

    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, Yesterday, 11:13 AM
    0 responses
    21 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, Yesterday, 10:48 AM
    0 responses
    18 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    36 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    42 views
    0 reactions
    Last Post SEQadmin2  
    Working...