changed how we count sequences in a fastafile to allow for '>' in sequence names

[mothur.git] / seqsummarycommand.cpp

diff --git a/seqsummarycommand.cpp b/seqsummarycommand.cpp
index d613391fb5e27da1d2d933e1b4b2375b69170a95..7226f9ddac16e61ebc2239ff4d868474b990ec25 100644 (file)
--- a/seqsummarycommand.cpp
+++ b/seqsummarycommand.cpp
@@ -228,9 +228,9 @@ int SeqSummaryCommand::execute(){
                                if(processors == 1){
                                        ifstream inFASTA;
                                        openInputFile(fastafile, inFASTA);
-                                       numSeqs=count(istreambuf_iterator<char>(inFASTA),istreambuf_iterator<char>(), '>');
-                                       inFASTA.close();
-                       cout << numSeqs << endl;                
+                                       getNumSeqs(inFASTA, numSeqs);
+                                       inFASTA.close();        
+                                               
                                        lines.push_back(new linePair(0, numSeqs));
                                        
                                        driverCreateSummary(startPosition, endPosition, seqLength, ambigBases, longHomoPolymer, fastafile, summaryFile, lines[0]);
@@ -250,7 +250,7 @@ int SeqSummaryCommand::execute(){
                #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));