-// splitAtEquals(temp2, temp); //separates length=242, temp=length, temp2=242
-// convert(temp, seqLength); //converts string to int
-//
-// if (name.length() != 0) { if(name.substr(1) != seq.getName()) { m->mothurOut("sequence name mismatch btwn fasta and qual file"); m->mothurOutEndLine(); } }
-
- string rawSequence = seq.getUnaligned();
- int seqLength = seq.getNumBases();
- bool success = 0; //guilty until proven innocent
- string name;
-
- qFile >> name;
- if (name[0] == '>') { if(name.substr(1) != seq.getName()) { m->mothurOut("sequence name mismatch btwn fasta: " + seq.getName() + " and qual file: " + name); m->mothurOutEndLine(); } }
-
- while (!qFile.eof()) { char c = qFile.get(); if (c == 10 || c == 13){ break; } }
-
- int score;
- int end = seqLength;
-
- for(int i=0;i<seqLength;i++){
- qFile >> score;
-
- if(score < qThreshold){
- end = i;
- break;
- }
- }
- for(int i=end+1;i<seqLength;i++){
- qFile >> score;
- }
-
- seq.setUnaligned(rawSequence.substr(0,end));
-
- return 1;
- }
- catch(exception& e) {
- m->errorOut(e, "TrimSeqsCommand", "stripQualThreshold");
- exit(1);
- }
-}
+// seq.setUnaligned(rawSequence.substr(0,end));
+//
+// return 1;
+// }
+// catch(exception& e) {
+// m->errorOut(e, "TrimSeqsCommand", "stripQualThreshold");
+// exit(1);
+// }
+//}