openInputFile(fastafile, inFASTA);
numSeqs=count(istreambuf_iterator<char>(inFASTA),istreambuf_iterator<char>(), '>');
inFASTA.close();
-
+ cout << numSeqs << endl;
lines.push_back(new linePair(0, numSeqs));
driverCreateSummary(startPosition, endPosition, seqLength, ambigBases, longHomoPolymer, fastafile, summaryFile, lines[0]);
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());