X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=seqerrorcommand.cpp;h=0a6eae93df72e8556554bef32c5a52074a4d0a54;hb=c27369bdcc3bebc4c449b9ae2f462d2a604a5503;hp=1fe60e8954acd14dc1a1e02d37c0d8bd987cc1c0;hpb=7a18f1438113f6c52c66f97ae0044cfa365dd221;p=mothur.git

diff --git a/seqerrorcommand.cpp b/seqerrorcommand.cpp
index 1fe60e8..0a6eae9 100644
--- a/seqerrorcommand.cpp
+++ b/seqerrorcommand.cpp
@@ -298,7 +298,7 @@ SeqErrorCommand::SeqErrorCommand(string option)  {
 			}
             else{
                 if(reportFileName != ""){
-                    m->mothurOut("we are ignoring the report file if your sequences are not aligned.  we will check that the sequences in your fasta and and qual fileare the same length.");
+                    m->mothurOut("we are ignoring the report file if your sequences are not aligned.  we will check that the sequences in your fasta and and qual file are the same length.");
                     m->mothurOutEndLine();
                 }
             }
@@ -759,7 +759,10 @@ int SeqErrorCommand::driver(string filename, string qFileName, string rFileName,
             int numParentSeqs = -1;
             int closestRefIndex = -1;
                         
-            numParentSeqs = chimeraTest.analyzeQuery(query.getName(), query.getAligned(), outChimeraReport);
+            string querySeq = query.getAligned();
+            if (!aligned) {  querySeq = query.getUnaligned();  }
+            
+            numParentSeqs = chimeraTest.analyzeQuery(query.getName(), querySeq, outChimeraReport);
             
             closestRefIndex = chimeraTest.getClosestRefIndex();
             
@@ -1378,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);
 					
@@ -1438,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