rsem-calculate-expression error

[Crystal Gowda]

I have prepared reference using the command:

rsem-prepare-reference -p 5 --bowtie2 cd_hit_est.fasta ref_output_12

And the files i obtained are:

ref_12.1.bt2

ref_12.2.bt2

ref_12.3.bt2 

ref_12.4.bt2 

ref_12.grp

ref_12.idx.fa

ref_12.n2g.idx.fa 

ref_12.rev.1.bt2

ref_12.rev.2.bt2

ref_12.seq

ref_12.ti 

ref_12.transcripts.fa

then i tried running the command:

rsem-calculate-expression -p 5 --bowtie2 --ci-memory 5000 --paired-end SRR5183402_1.fastq SRR5183402_2.fastq ref_output_12 SSE_fm

and this is the error I am getting:

bowtie2 -q --phred33 --sensitive --dpad 0 --gbar 99999999 --mp 1,1 --np 1 --score-min L,0,-0.1 -I 1 -X 1000 --no-mixed --no-discordant -p 5 -k 200 -x ref_output_12 -1 SRR5183402_1.fastq -2 SRR5183402_2.fastq 2> SSE_fm.log | samtools view -b -o SSE_fm.temp/SSE_fm.bam -
[main_samview] fail to read the header from "-".
"bowtie2 -q --phred33 --sensitive --dpad 0 --gbar 99999999 --mp 1,1 --np 1 --score-min L,0,-0.1 -I 1 -X 1000 --no-mixed --no-discordant -p 5 -k 200 -x ref_output_12 -1 SRR5183402_1.fastq -2 SRR5183402_2.fastq 2> SSE_fm.log | samtools view -b -o SSE_fm.temp/SSE_fm.bam -" failed! Plase check if you provide correct parameters/options for the pipeline!

Previosly I had used trinity output, but i got same error. This time i am using clustered data, still the same error.

I have used bowtie2 option while creating both reference and rsem-calculate-expression.

Any help on this would be highly appreciated. Thank you very much.

Best regards,

Bhoomika

[Herman Xin Yang Leong]

Hi Bhoomika and other RSEM users,

I ran into the exact same issue, so I am wondering if you managed to solve this particular error. It would be of tremendous help if you can provide any useful information.

Thanks,
Herman