*
*/
-#include <iostream>
-#include <fstream>
#include "command.hpp"
-#include "utilities.hpp"
#include "fastamap.h"
#include "globaldata.hpp"
-/* The deconvolute command reads a fasta file, finds the duplicate sequences and outputs a names file
+/* The unique.seqs command reads a fasta file, finds the duplicate sequences and outputs a names file
containing 2 columns. The first being the groupname and the second the list of identical sequence names. */
-using namespace std;
class DeconvoluteCommand : public Command {
public:
- DeconvoluteCommand() {};
+ DeconvoluteCommand(string);
~DeconvoluteCommand() { delete fastamap; };
- int execute();
+ int execute();
+ void help();
private:
GlobalData* globaldata;
ifstream in;
ofstream out, outFasta;
string filename, outputFileName, outFastafile;
+ OptionParser* parser;
+ map<string, string> parameters;
+ map<string, string>::iterator it;
+ bool abort;
+
};