2 #include "textfile.h"
\r
4 #include "objscore.h"
\r
10 TextFile f(g_pstrSPFileName);
\r
15 ALPHA Alpha = ALPHA_Undefined;
\r
19 Alpha = a.GuessAlpha();
\r
22 case SEQTYPE_Protein:
\r
23 Alpha = ALPHA_Amino;
\r
35 Quit("Invalid SeqType");
\r
42 const unsigned uSeqCount = a.GetSeqCount();
\r
44 Quit("No sequences in input file %s", g_pstrSPFileName);
\r
46 MSA::SetIdCount(uSeqCount);
\r
47 for (unsigned uSeqIndex = 0; uSeqIndex < uSeqCount; ++uSeqIndex)
\r
48 a.SetSeqId(uSeqIndex, uSeqIndex);
\r
50 SetSeqWeightMethod(g_SeqWeight1);
\r
52 TreeFromMSA(a, tree, g_Cluster2, g_Distance2, g_Root2);
\r
53 SetMuscleTree(tree);
\r
54 SetMSAWeightsMuscle((MSA &) a);
\r
56 SCORE SP = ObjScoreSP(a);
\r
58 Log("File=%s;SP=%.4g\n", g_pstrSPFileName, SP);
\r
59 fprintf(stderr, "File=%s;SP=%.4g\n", g_pstrSPFileName, SP);
\r