rsem-calculate-expression error
0
0
Entering edit mode
3 months ago
Bhoomika • 0

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

RNA-Seq rsem rsem-calculate-expression • 120 views
ADD COMMENT

Login before adding your answer.

Traffic: 1842 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6