X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=SingleReadQ.h;fp=SingleReadQ.h;h=976637fc269d79fa484cb0033421e7f7a0769c49;hb=86e650e9577999a7ba00ab454d1f6bf674b0ea70;hp=a0099d9c41054b68da45880f94a8857e8f232f47;hpb=5b867894f164d9fa3fbc1a8e659b34890b977716;p=rsem.git

diff --git a/SingleReadQ.h b/SingleReadQ.h
index a0099d9..976637f 100644
--- a/SingleReadQ.h
+++ b/SingleReadQ.h
@@ -67,6 +67,8 @@ void SingleReadQ::calc_lq() {
   int numA = 0, numT = 0, numAO = 0, numTO = 0; // numAO : number of A in overlap seed region
   int threshold_1, threshold_2;
 
+  if (len < OLEN) { low_quality = true; return; }
+
   threshold_1 = int(0.9 * len - 1.5 * sqrt(len * 1.0) + 0.5);
   threshold_2 = (OLEN - 1) / 2 + 1;
   for (int i = 0; i < len; i++) {