struct Compare {
int AA, AT, AG, AC, TA, TT, TG, TC, GA, GT, GG, GC, CA, CT, CG, CC, NA, NT, NG, NC, Ai, Ti, Gi, Ci, Ni, dA, dT, dG, dC;
- string refName, queryName;
+ string refName, queryName, sequence;
double errorRate;
int weight, matches, mismatches, total;
mismatches = 0;
total = 0;
errorRate = 1.0000;
+ sequence = "";
}
};
class SeqErrorCommand : public Command {
public:
SeqErrorCommand(string);
+ SeqErrorCommand();
~SeqErrorCommand();
+ vector<string> getRequiredParameters();
+ vector<string> getValidParameters();
+ vector<string> getRequiredFiles();
+ map<string, vector<string> > getOutputFiles() { return outputTypes; }
int execute();
void help();
map<string,int> getWeights();
Compare getErrors(Sequence, Sequence);
void printErrorHeader();
- void printErrorData(Compare);
-
- string queryFileName, referenceFileName, namesFileName, errorFileName, outputDir;
+ void printErrorData(Compare, int);
+ void printSubMatrix();
+ void printErrorFRFile(map<char, vector<int> >, map<char, vector<int> >);
+ void printErrorQuality(map<char, vector<int> >);
+ void printQualityFR(vector<vector<int> >, vector<vector<int> >);
+
+ string queryFileName, referenceFileName, qualFileName, reportFileName, namesFileName, outputDir;
double threshold;
+ bool ignoreChimeras;
int numRefs;
- ofstream errorFile;
+ ofstream errorSummaryFile, errorSeqFile;
+ vector<string> outputNames;
+ map<string, vector<string> > outputTypes;
vector<Sequence> referenceSeqs;
+ vector<vector<int> > substitutionMatrix;
};
#endif