rsem v1.1.11, allow spaces appear in field seqname and attributes gene_id, transcript...

[rsem.git] / Refs.h

diff --git a/Refs.h b/Refs.h
index 41fcf061112bc4b6ca1cc96cb4395ca4a819e944..63b67e16e0d705d8c66d3946085f76b79496724c 100644 (file)
--- a/Refs.h
+++ b/Refs.h
@@ -52,7 +52,7 @@ class Refs {
   }
 
   bool isValid(int sid, int dir, int pos, const std::string& readseq, int LEN, int C) {
-    if (sid <= 0 || sid > M || dir != 0 && dir != 1 ||  pos < 0 || pos + LEN > seqs[sid].getTotLen() || LEN > (int)readseq.length()) return false;
+    if (sid <= 0 || sid > M || (dir != 0 && dir != 1) ||  pos < 0 || pos + LEN > seqs[sid].getTotLen() || LEN > (int)readseq.length()) return false;
     const std::string& seq = seqs[sid].getSeq(dir);
     return countMismatch(seq, pos, readseq, LEN, C) <= C;
   }
@@ -107,7 +107,7 @@ void Refs::makeRefs(char *inpF, RefSeqPolicy& policy, PolyARules& rules) {
 }
 
 //inpF in fasta format, with sequence all in one line together
-//option 0 read all, 1 do not read sequences and names
+//option 0 read all, 1 do not read sequences
 void Refs::loadRefs(char *inpF, int option) {
   std::ifstream fin;
   RefSeq seq;