int main() { int evalcost; char *align = "R7I2R2D1R3I1R3"; char *seq1 = "acgtagatatatagat"; char *seq2 = "agaaagaggtaagaggga"; alignment alg = alignment_new(seq1, seq2, align); alignment_show(alg); evalcost = alignment_evalcost(alg); printf("\tCosts: %d\n\n", evalcost); printf("\tadding one deletion...\n"); alignment_add_operation(&alg, 1, 'D'); alignment_show(alg); printf("\tadding one insertion...\n"); alignment_add_operation(&alg, 1, 'I'); alignment_show(alg); printf("\tadding one replacement...\n"); alignment_add_operation(&alg, 1, 'R'); alignment_show(alg); return(EXIT_SUCCESS); }
int alignment_add_operations (alignment *a, char operation, int length) { int i; int s = 0; for (i = 0; i < length && s == 0; i++) { s += alignment_add_operation (a, operation); } return s; }
int alignment_add_replacement (alignment* a) { return alignment_add_operation (a, 'R'); }
int alignment_add_insertion (alignment* a) { return alignment_add_operation (a, 'I'); }
int alignment_add_deletion (alignment* a) { return alignment_add_operation (a, 'D'); }