/* Calculate one sequence's contribution to the alignment score */
double marginal_score(glam2_scorer *s, glam2_aln *aln,
int seq, const fasta *f) {
double score = aln->score;
unalign(aln, seq, f);
score -= aln_score(s, aln);
realign(aln, seq, f);
return score;
}
static string Unalign(const char* filename) {
char s[MAX_PATH];
strcpy(s, filename);
return string(unalign(s));
}
