X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=seqerrorcommand.cpp;h=e48d50794d15b4b77dae6a15134d9bec54c8a911;hb=81ad9ebabe7f951eb2d1c922b1342e9d6e521c4e;hp=67e43aa0d4b7f6a6ff251bc45a8254ca49b8b162;hpb=20071b183e619c122bf9f63b4bb42722507c4e4a;p=mothur.git

diff --git a/seqerrorcommand.cpp b/seqerrorcommand.cpp
index 67e43aa..e48d507 100644
--- a/seqerrorcommand.cpp
+++ b/seqerrorcommand.cpp
@@ -673,10 +673,10 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName,
 		int numSeqs = 0;
 		
 		map<string, int>::iterator it;
-		qScoreErrorMap['m'].assign(41, 0);
-		qScoreErrorMap['s'].assign(41, 0);
-		qScoreErrorMap['i'].assign(41, 0);
-		qScoreErrorMap['a'].assign(41, 0);
+		qScoreErrorMap['m'].assign(101, 0);
+		qScoreErrorMap['s'].assign(101, 0);
+		qScoreErrorMap['i'].assign(101, 0);
+		qScoreErrorMap['a'].assign(101, 0);
 		
 		errorForward['m'].assign(maxLength,0);
 		errorForward['s'].assign(maxLength,0);
@@ -712,8 +712,8 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName,
 			qualForwardMap.resize(maxLength);
 			qualReverseMap.resize(maxLength);
 			for(int i=0;i<maxLength;i++){
-				qualForwardMap[i].assign(41,0);
-				qualReverseMap[i].assign(41,0);
+				qualForwardMap[i].assign(101,0);
+				qualReverseMap[i].assign(101,0);
 			}	
 		}
 		else if(qFileName != "" && !aligned){
@@ -724,8 +724,8 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName,
 			qualForwardMap.resize(maxLength);
 			qualReverseMap.resize(maxLength);
 			for(int i=0;i<maxLength;i++){
-				qualForwardMap[i].assign(41,0);
-				qualReverseMap[i].assign(41,0);
+				qualForwardMap[i].assign(101,0);
+				qualReverseMap[i].assign(101,0);
 			}	
         }
         
@@ -1271,7 +1271,7 @@ void SeqErrorCommand::printErrorQuality(map<char, vector<int> > qScoreErrorMap){
 		outputNames.push_back(errorQualityFileName);  outputTypes["errorquality"].push_back(errorQualityFileName);
 
 		errorQualityFile << "qscore\tmatches\tsubstitutions\tinsertions\tambiguous" << endl;
-		for(int i=0;i<41;i++){
+		for(int i=0;i<101;i++){
 			errorQualityFile << i << '\t' << qScoreErrorMap['m'][i] << '\t' << qScoreErrorMap['s'][i] << '\t' << qScoreErrorMap['i'][i] << '\t'<< qScoreErrorMap['a'][i] << endl;
 		}
 		errorQualityFile.close();
@@ -1381,6 +1381,10 @@ int SeqErrorCommand::setLines(string filename, string qfilename, string rfilenam
 					
 					string sname = "";  nameStream >> sname;
 					sname = sname.substr(1);
+                    
+                    for (int i = 0; i < sname.length(); i++) {
+                        if (sname[i] == ':') { sname[i] = '_'; m->changedSeqNames = true; }
+                    }
 					
 					map<string, int>::iterator it = firstSeqNames.find(sname);
 					
@@ -1441,6 +1445,10 @@ int SeqErrorCommand::setLines(string filename, string qfilename, string rfilenam
                     istringstream nameStream(input);
                     string sname = "";  nameStream >> sname;
                     
+                    for (int i = 0; i < sname.length(); i++) {
+                        if (sname[i] == ':') { sname[i] = '_'; m->changedSeqNames = true; }
+                    }
+                    
                     map<string, int>::iterator it = firstSeqNamesReport.find(sname);
                 
                     if(it != firstSeqNamesReport.end()) { //this is the start of a new chunk