This is a test version of Biostars. For the public version, visit https://www.biostars.org.
rsem error message: "Number of alignable reads does not match!"

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

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

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.

That's disconcerting!

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

0 answers

No answers yet.

Log in to answer this question.