X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=refchimeratest.h;h=5ea0801f17f437af48efbc148df9ef04875e3f31;hp=41df4ebba55842eae2c36ff9f02af303d20b64c7;hb=615301e57c25e241356a9c2380648d117709458d;hpb=3c5a709812a966bf3cabc8bb15f86f0c36cf1e34 diff --git a/refchimeratest.h b/refchimeratest.h index 41df4eb..5ea0801 100644 --- a/refchimeratest.h +++ b/refchimeratest.h @@ -16,13 +16,22 @@ class RefChimeraTest { public: - RefChimeraTest(vector&, bool); + RefChimeraTest(){}; + RefChimeraTest(vector&, bool); int printHeader(ofstream&); - int analyzeQuery(string, string, ofstream&); + int analyzeQuery(string, string, ofstream&); int getClosestRefIndex(); + string getClosestRefAlignment(); + string getQueryAlignment(); + private: - int getMismatches(string&, vector >&, vector >&, int&); - int getChimera(vector >&, vector >&, int&, int&, int&, vector&, vector&, vector&, vector&); + int getAlignedMismatches(string&, vector >&, vector >&, int&); + int analyzeAlignedQuery(string, string, ofstream&); + int analyzeUnalignedQuery(string, string, ofstream&); + double alignQueryToReferences(string, string, string&, string&, double&); + int getUnalignedDiffs(string, string, vector&, vector&, vector&, vector&); + + int getChimera(vector >&, vector >&, int&, int&, int&, vector&, vector&, vector&, vector&); int getTrimera(vector >&, vector >&, int&, int&, int&, int&, int&, vector&, vector&, vector&, vector&); string stitchBimera(int, int, int); string stitchTrimera(int, int, int, int, int); @@ -33,6 +42,8 @@ private: int numRefSeqs; int alignLength; int bestMatch; + string bestRefAlignment; + string bestQueryAlignment; //ofstream chimeraReportFile; bool aligned;