]> git.donarmstrong.com Git - mothur.git/blobdiff - chimeraslayercommand.cpp
added mantel command
[mothur.git] / chimeraslayercommand.cpp
index 67b70f27cc0aa0592cd87d5b97c1471fc0edb5ff..38c097f18b8124f1dadbd657af8413fc74f895c8 100644 (file)
@@ -273,6 +273,9 @@ ChimeraSlayerCommand::ChimeraSlayerCommand(string option)  {
                        temp = validParameter.validFile(parameters, "trim", false);                             if (temp == "not found") { temp = "f"; }
                        trim = m->isTrue(temp); 
                        
+                       //temp = validParameter.validFile(parameters, "trimera", false);                                if (temp == "not found") { temp = "f"; }
+                       //trimera = m->isTrue(temp); 
+                       
                        search = validParameter.validFile(parameters, "search", false);                 if (search == "not found") { search = "distance"; }
                        
                        temp = validParameter.validFile(parameters, "iters", false);                    if (temp == "not found") { temp = "100"; }              
@@ -309,6 +312,7 @@ void ChimeraSlayerCommand::help(){
                m->mothurOut("When using MPI, the processors parameter is set to the number of MPI processes running. \n");
                #endif
                m->mothurOut("The trim parameter allows you to output a new fasta file containing your sequences with the chimeric ones trimmed to include only their longest piece, default=F. \n");
+               //m->mothurOut("The trimera parameter allows you to check both peices of a chimeric sequence for chimeras, thus looking for trimeras and quadmeras. default=F. \n");
                m->mothurOut("The window parameter allows you to specify the window size for searching for chimeras, default=50. \n");
                m->mothurOut("The increment parameter allows you to specify how far you move each window while finding chimeric sequences, default=5.\n");
                m->mothurOut("The numwanted parameter allows you to specify how many sequences you would each query sequence compared with, default=15.\n");
@@ -619,6 +623,9 @@ int ChimeraSlayerCommand::driver(linePair* filePos, string outputFName, string f
                                        Sequence* trimmed = chimera->print(out, out2);
                                        
                                        if (trim) { trimmed->printSequence(out3); delete trimmed; }
+                                       
+                                       //do you want to check both pieces for chimeras
+                                       //if (trimera) {}
                                }
                        count++;
                        }
@@ -701,6 +708,9 @@ int ChimeraSlayerCommand::driverMPI(int start, int num, MPI_File& inMPI, MPI_Fil
                                                MPI_File_write_shared(outFastaMPI, buf2, length, MPI_CHAR, &status);
                                                delete buf2;
                                        }
+                                       
+                                       //do you want to check both pieces for chimeras
+                                       //if (trimera) {}
                                                
                                }
                        }