Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • wokai001
    Member
    • Nov 2010
    • 20

    #1

    tophat rna-seq crash

    I'm using tophat for aligning Illumina rna-seq data. It crashed in one of the last steps:

    [Tue Jan 10 08:26:56 2012] Beginning TopHat run (v1.4.0)
    -----------------------------------------------
    [Tue Jan 10 08:26:56 2012] Preparing output location ./tophat_out_bwa/
    [Tue Jan 10 08:26:56 2012] Checking for Bowtie index files
    [Tue Jan 10 08:26:56 2012] Checking for reference FASTA file
    [Tue Jan 10 08:26:56 2012] Checking for Bowtie
    Bowtie version: 0.12.7.0
    [Tue Jan 10 08:26:56 2012] Checking for Samtools
    Samtools Version: 0.1.17
    [Tue Jan 10 08:26:56 2012] Generating SAM header for hshiv
    format: fastq
    quality scale: phred33 (default)
    [Tue Jan 10 08:27:16 2012] Preparing reads
    left reads: min. length=101, count=86332313
    [Tue Jan 10 08:55:48 2012] Mapping left_kept_reads against hshiv with Bowtie
    [Tue Jan 10 09:35:06 2012] Processing bowtie hits
    [Tue Jan 10 10:30:04 2012] Mapping left_kept_reads_seg1 against hshiv with Bowtie (1/4)
    [Tue Jan 10 12:41:18 2012] Mapping left_kept_reads_seg2 against hshiv with Bowtie (2/4)
    [Tue Jan 10 14:54:35 2012] Mapping left_kept_reads_seg3 against hshiv with Bowtie (3/4)
    [Tue Jan 10 17:08:07 2012] Mapping left_kept_reads_seg4 against hshiv with Bowtie (4/4)
    [Tue Jan 10 18:56:53 2012] Searching for junctions via segment mapping
    [FAILED]
    Error: segment-based junction search failed with err =1
    Error: could not get read# 839911731 from stream!

    I used tophat on Ubuntu Linux (linux kernel 3.0.0-14-generic). Actually, I first aligned my reads with bwa (bwa aln, then bwa samse) and extracted the unmapped reads with samtools into bam. From there, I used the (tophat) bam2fastx tool to convert the unmapped reads back into fastq. The unmapped reads (about 30%) then go as Input into tophat.

    Has someone any idea what's going on there?

    Thanks for your help
    Wolfgang

Latest Articles

Collapse

  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2




    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification, detection...
    10-30-2026, 06:46 AM
  • 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

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
33 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-28-2026, 09:51 AM
0 responses
21 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
72 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 10:48 AM
0 responses
48 views
0 reactions
Last Post SEQadmin2  
Working...