Comments (3)
@sopenaml, what command did you use when you got this error?
from rsem.
It's a bit convoluted, but I'm running the nfcore/rnaseq pipeline using STAR-rsem to align/quantify counts per gene and my fastq files are quite big. Below is the error obtained that also includes the step where the error occurred.
Error executing process > 'NFCORE_RNASEQ:RNASEQ:QUANTIFY_RSEM:RSEM_CALCULATEEXPRESSION (HepG2_Rep3)'
Caused by:
Process `NFCORE_RNASEQ:RNASEQ:QUANTIFY_RSEM:RSEM_CALCULATEEXPRESSION (HepG2_Rep3)` terminated with an error exit status (255)
Command executed:
INDEX=`find -L ./ -name "*.grp" | sed 's/\.grp$//'`
rsem-calculate-expression \
--num-threads 12 \
--temporary-folder ./tmp/ \
--strandedness forward \
--paired-end \
--star --star-output-genome-bam --star-gzipped-read-file --estimate-rspd --seed 1 \
HepG2_Rep3_1_val_1.fq.gz HepG2_Rep3_2_val_2.fq.gz \
$INDEX \
HepG2_Rep3
cat <<-END_VERSIONS > versions.yml
"NFCORE_RNASEQ:RNASEQ:QUANTIFY_RSEM:RSEM_CALCULATEEXPRESSION":
rsem: $(rsem-calculate-expression --version | sed -e "s/Current version: RSEM v//g")
star: $(STAR --version | sed -e "s/STAR_//g")
END_VERSIONS
Command exit status:
255
Command output:
STAR --genomeDir ./rsem --outSAMunmapped Within --outFilterType BySJout --outSAMattributes NH HI AS NM MD --outFilterMultimapNmax 20 --outFilterMismatchNmax 999 --outFilterMismatchNoverLmax 0.04 --alignIntronMin 20 --alignIntronMax 1000000 --alignMatesGapMax 1000000 --alignSJoverhangMin 8 --alignSJDBoverhangMin 1 --sjdbScore 1 --runThreadN 12 --genomeLoad NoSharedMemory --outSAMtype BAM Unsorted --quantMode TranscriptomeSAM --outSAMheaderHD @HD VN:1.4 SO:unsorted --outFileNamePrefix ./tmp//HepG2_Rep3 --readFilesCommand zcat --readFilesIn HepG2_Rep3_1_val_1.fq.gz HepG2_Rep3_2_val_2.fq.gz
STAR --genomeDir ./rsem --outSAMunmapped Within --outFilterType BySJout --outSAMattributes NH HI AS NM MD --outFilterMultimapNmax 20 --outFilterMismatchNmax 999 --outFilterMismatchNoverLmax 0.04 --alignIntronMin 20 --alignIntronMax 1000000 --alignMatesGapMax 1000000 --alignSJoverhangMin 8 --alignSJDBoverhangMin 1 --sjdbScore 1 --runThreadN 12 --genomeLoad NoSharedMemory --outSAMtype BAM Unsorted --quantMode TranscriptomeSAM --outSAMheaderHD @HD VN:1.4 SO:unsorted --outFileNamePrefix ./tmp//HepG2_Rep3 --readFilesCommand zcat --readFilesIn HepG2_Rep3_1_val_1.fq.gz HepG2_Rep3_2_val_2.fq.gz
STAR version: 2.7.10a compiled: 2022-01-14T18:50:00-05:00 :/home/dobin/data/STAR/STARcode/STAR.master/source
Jul 12 17:07:07 ..... started STAR run
Jul 12 17:07:07 ..... loading genome
Jul 12 17:07:43 ..... started mapping
"STAR --genomeDir ./rsem --outSAMunmapped Within --outFilterType BySJout --outSAMattributes NH HI AS NM MD --outFilterMultimapNmax 20 --outFilterMismatchNmax 999 --outFilterMismatchNoverLmax 0.04 --alignIntronMin 20 --alignIntronMax 1000000 --alignMatesGapMax 1000000 --alignSJoverhangMin 8 --alignSJDBoverhangMin 1 --sjdbScore 1 --runThreadN 12 --genomeLoad NoSharedMemory --outSAMtype BAM Unsorted --quantMode TranscriptomeSAM --outSAMheaderHD @HD VN:1.4 SO:unsorted --outFileNamePrefix ./tmp//HepG2_Rep3 --readFilesCommand zcat --readFilesIn HepG2_Rep3_1_val_1.fq.gz HepG2_Rep3_2_val_2.fq.gz" failed! Plase check if you provide correct parameters/options for the pipeline!
Command error:
EXITING because of fatal error: buffer size for SJ output is too small
Solution: increase input parameter --limitOutSJcollapsed
Jul 12 17:35:08 ...... FATAL ERROR, exiting
from rsem.
STAR's default buffer size for SJ output is usually fine for human data. Just for testing and perhaps as a temporary fix, would you like to replace rsem-calculate-expression's line 489 by the following?
" --outSAMheaderHD \@HD VN:1.4 SO:unsorted --limitOutSJcollapsed 2000000 ".
assuming you want to increase it to 2000000.
from rsem.
Related Issues (20)
- The SAM/BAM file declares less than one reference sequence! HOT 2
- genomeSAsparseD option for RSEM Star?
- Problem implementing a reverse-strand specific analysis HOT 2
- about gffcompare
- RSEM deal with multimapping? HOT 1
- Pearl error
- error when using convert-sam-for-rsem to a RSEM compatible BAM file.
- rsem-prepare-reference arguments
- Using STAR BAMs from genomic alignment HOT 2
- rsem-parse-alignments and rsem-build-read-index
- Gtf file
- cleaned.fq.gz: command not found HOT 3
- Two mates are aligned to two different transcripts HOT 1
- rsem-generate-ngvector command not found HOT 1
- bam only containing aligned reads produced different expression profile
- Not sure if rsem is using --num-threads option...
- Single-end fragment length mean and SD
- rsem-extract-reference-transcripts fails with "Error Message: Strand is neither '+' nor '-'!"
- Different read counts from STAR and RSEM HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from rsem.