try {
int seqLength = errorSeq.size();
+ cout << start << '\t' << stop << '\t' << seqLength << endl;
+
int qIndex = start - 1;
for(int i=0;i<seqLength;i++){
else if(errorSeq[i] == 'd') { /* there are no qScores for deletions */ }
if(errorSeq[i] != 'd') { qIndex++; }
-
+ if(qIndex > stop){ break; }
}
}
catch(exception& e) {