if(processors == 1){
ifstream inFASTA;
openInputFile(fastafile, inFASTA);
- numSeqs=count(istreambuf_iterator<char>(inFASTA),istreambuf_iterator<char>(), '>');
- inFASTA.close();
-
+ getNumSeqs(inFASTA, numSeqs);
+ inFASTA.close();
+
lines.push_back(new linePair(0, numSeqs));
driverCreateSummary(startPosition, endPosition, seqLength, ambigBases, longHomoPolymer, fastafile, summaryFile, lines[0]);
#else
ifstream inFASTA;
openInputFile(fastafile, inFASTA);
- numSeqs=count(istreambuf_iterator<char>(inFASTA),istreambuf_iterator<char>(), '>');
+ getNumSeqs(inFASTA, numSeqs);
inFASTA.close();
lines.push_back(new linePair(0, numSeqs));
if (m->control_pressed) { in.close(); outSummary.close(); return 1; }
Sequence current(in);
+
if (current.getName() != "") {
startPosition.push_back(current.getStartPos());
endPosition.push_back(current.getEndPos());