rsem-v1.1.14b3

[rsem.git] / EM.cpp

diff --git a/EM.cpp b/EM.cpp
index 8362f7d47dea60b762801d5877352d196bcf35f7..f9449a85cd3226a3f29356f10fa018e4ca908cba 100644 (file)
--- a/EM.cpp
+++ b/EM.cpp
@@ -102,7 +102,7 @@ void init(ReadReader<ReadType> **&readers, HitContainer<HitType> **&hitvs, doubl
                indices[i] = new ReadIndex(readFs[i]);
        }
        for (int i = 0; i < nThreads; i++) {
-               readers[i] = new ReadReader<ReadType>(s, readFs);
+               readers[i] = new ReadReader<ReadType>(s, readFs, refs.hasPolyA(), mparams.seedLen); // allow calculation of calc_lq() function
                readers[i]->setIndices(indices);
        }