next up previous contents
Next: Refine Shake Up: Dynamic Programming Previous: GlobalAlign Match

Refine Log

Calling Sequences:
LogDelLocalRefine(
m, DM)
Parameters:
m : Match
DM : DayMatrix

Returns: Match

Synposis: This function computes the alignment with maximum similarity score for the Match m at the PAM distance defined by the DayMatrix DM. It computes a local alignment by applying dynamic programming in a forwards and backwards manner until the score stabilizes. The LogDelLocalRefine function uses a logarithmic deletion cost.

Examples:

> DB := ReadDb('~cbrg/DB/SwissProt'):   CreateDayMatrices():
> m := Match(op(Sequence(Entry(8588))), op(Sequence(Entry(8577)))):
> m := LogDelLocalRefine(m, DM);           # compute best local alignment at 250 PAM
> print(m);




Gaston Gonnet
1998-09-15