Question: rsem-calculate-expression error message
1
gravatar for walker
2.1 years ago by
walker30
walker30 wrote:

I'm trying to calculate exppression from raw fastq, and I am getting error when rsem finish STAR alignment. here is my sh script:

fastp -i 1.fastq.gz -I 2.fastq.gz -o clean_1.fq.gz -O clean_2.fq.gz -j RNA.fastp.json -h RNA.fastp.html
   RSEM-1.3.0/bin/rsem-calculate-expression -p 2 --paired-end --phred33-quals --star --star-path PATH/TO/STAR --star-gzipped-read-file --star-output-genome-bam --estimate-rspd --append-names clean_1.fq.gz  clean_2.fq.gz PATH/TO/reference PATH/TO/2.RSEMstar

and here is my error message:

 ..... started STAR run
 ..... loading genome
 ..... started mapping
..... finished successfully

     "rsem-parse-alignments RSEM/grch37r87_STAR PATH/TO/2.RSEMstar.temp/2.RSEMstar PATH/TO/2.RSEMstar.stat/2.RSEMstar PATH/TO/2.RSEMstar.temp/2.RSEMstar.bam 3 -tag XM" failed! Plase check if you provide correct parameters/options for the pipeline!

anyone knows why and can help solve this questions? Thanks

rna-seq rsem software error • 1.5k views
ADD COMMENTlink modified 2.1 years ago • written 2.1 years ago by walker30
1

Is it normal that you have the same file specified in two arguments? -o clean_2.fq.gz -O clean_2.fq.gz

ADD REPLYlink written 2.1 years ago by Macspider3.2k

Thanks for your answer and sorry for my mistake I write on this question! Acturally, -o clean_1.fq.gz -O clean_2.fq.gz is what i write in my script And I correct it in my question,Thanks!

ADD REPLYlink written 2.1 years ago by walker30
1

You also declare twice the --star argument in RSEM. That shouldn't be a problem, but since it seems that you're not super-solid on your command writing, perhaps check if every file you declare is in the folder you want it to be, check if every tag you added is exactly as it should be, and stuff like that (it's the cause of 99% of the errors in the bioinformatics world, regardless of the experience!).

ADD REPLYlink written 2.1 years ago by Macspider3.2k

Did you find a solution? I am facing the same issue

ADD REPLYlink written 19 months ago by mrpetzel960
1
gravatar for Devon Ryan
2.1 years ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

clean_1.fq.gz isn't being created by fastp, since you specified read 2 twice.

ADD COMMENTlink written 2.1 years ago by Devon Ryan97k

Thanks for your answer and sorry for my mistake I write on this question! Acturally, -o clean_1.fq.gz -O clean_2.fq.gz is what i write in my script

ADD REPLYlink modified 2.1 years ago • written 2.1 years ago by walker30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1083 users visited in the last hour