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

[rsem.git] / PairedEndRead.h

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