Question: rsem error message: "Number of alignable reads does not match!"
0
gravatar for chunxuan
5.0 years ago by
chunxuan10
Germany
chunxuan10 wrote:

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 • 1.9k views
ADD COMMENTlink modified 4.9 years ago by Biostar ♦♦ 20 • written 5.0 years ago by chunxuan10

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 REPLYlink modified 5.0 years ago • written 5.0 years ago by Devon Ryan89k

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 REPLYlink written 5.0 years ago by chunxuan10

That's disconcerting!

ADD REPLYlink written 5.0 years ago by Devon Ryan89k

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

ADD REPLYlink written 5.0 years ago by chunxuan10
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: 1520 users visited in the last hour