rsem error message: "Number of alignable reads does not match!"
0
0
Entering edit mode
7.7 years ago
ccshao ▴ 10

When I was running rsem (v1.2.12) on my paired-end RNA-seq data, it threw out this message:

Number of alignable reads does not match!

What could be problem?

The running script:

rsem-calculate-expression -p 6 --paired-end fastqP_1_fa fastqP_2_fa rsemIndex EXAMPLE

and have built the reference by rsem-prepare-reference

rsem RNA-Seq • 3.1k views
ADD COMMENT
0
Entering edit mode

What's the output of wc -l fastqP_1_fa fastqP_2_fa? My guess would be that the files are out of sync.

ADD REPLY
0
Entering edit mode

Good points. Surprisingly, when I re-run the same scripts for the same data, with only changing the number of cpu used, it works smoothly. I will submit all jobs and see the frequency of the errors.

ADD REPLY
0
Entering edit mode

That's disconcerting!

ADD REPLY
0
Entering edit mode

The error seems never appear again, not clear why it appears/disappears.

ADD REPLY

Login before adding your answer.

Traffic: 1758 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