rsem v1.1.13, speed up EM by only updating model parameters for first 10 iterations...

[rsem.git] / PairedEndReadQ.h

diff --git a/PairedEndReadQ.h b/PairedEndReadQ.h
index f04afb7089e66a853aa798afbc5415b678ee8ca5..36b703b059ed4bd4fad7a3eb873c216f49f22e19 100644 (file)
--- a/PairedEndReadQ.h
+++ b/PairedEndReadQ.h
@@ -61,6 +61,7 @@ void PairedEndReadQ::write(int argc, std::ostream* argv[]) {
 
 void PairedEndReadQ::calc_lq() {
        low_quality = mate1.isLowQuality() && mate2.isLowQuality();
+       if (mate1.getReadLength() < OLEN || mate2.getReadLength() < OLEN) low_quality = true;
 }
 
 #endif /* PAIREDENDREADQ_H_ */