Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • roschu
    Junior Member
    • Mar 2013
    • 1

    #1

    rsem-calculate-expression takes very long with paired-end reads (using bowtie2)

    I am want to use RSEM for transcript quantification. I have paired-end Illumina 1.5 encoded fastq files, a genome file, and a gtf file from Ensembl.
    First I used rsem-prepare-reference:


    rsem-prepare-reference --gtf myGTF.gtf \
    --bowtie2 \
    --bowtie2-path /my/path/bowtie2-2.2.3 \
    myGenome.fa \
    ./prepREF_bowtie2


    This works well (as far as I can tell) --- computation time: ~5min.
    Then I want to use rsem-calculate-expression:

    rsem-calculate-expression --phred64-quals \
    -p 10 \
    --time \
    --bowtie2 \
    --bowtie2-path /my/path/bowtie2-2.2.3 \
    --paired-end \
    ../RNA_reads/sample1_1.fq \
    ../RNA_reads/sample1_2.fq \
    ./prepREF_bowtie2 \
    ./Sample1_rsem


    I use a linux cluster with 10 CPUs and 10 GB memory. However all the output I get is the bowtie2 command line. I aborted the calculation after 30 hours. Why is this taking so long?
    (Using bowtie instead of bowtie2 seems to work: at least I am getting rsem-parse-alignments outputs. However bowtie was only able to align ~50% of the reads and I got some "Warning: Exhausted best-first chunk memory for read" warnings).
    What am I doing wrong?

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...
    Yesterday, 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
29 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-28-2026, 09:51 AM
0 responses
20 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 10-22-2026, 11:13 AM
0 responses
71 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...