Multiple.hs

{- |
   Multiple alignments.

-}

module Bio.Alignment.Multiple

where

import Bio.Alignment.AlignData
import Bio.Sequence
import Bio.Clustering

-- | Progressive multiple alignment.
--   Calculate a tree from agglomerative clustering, then align
--   at each branch going bottom up.  Returns a list of columns (rows?).
progressive :: (Sequence -> Sequence -> (Double,Alignment)) -> [Sequence] -> [String]
progressive = undefined

-- |  Derive alignments indirectly, i.e. calculate A|C using alignments A|B and B|C.
--    This is central for 'Coffee' evaluation of alignments, and T-Coffee construction
--    of alignments.
indirect :: Alignment -> Alignment -> Alignment
indirect (Repl x1 x2:xs) (Repl y1 y2:ys) = Repl x1 y2 : indirect xs ys -- assert x2==y1
indirect xs@(Repl _ _:_) (Ins y1:ys)     = Ins y1     : indirect xs ys
indirect (Repl x1 _:xs) (Del y1:ys)      = Del x1     : indirect xs ys

indirect (Del x1:xs) ys                  = Del x1     : indirect xs ys -- imply del+ins/=repl

indirect (Ins x1:xs) (Repl _ y2:ys)      = Ins y2     : indirect xs ys -- assert x1==y1
indirect (Ins x1:xs) (Del y1:ys)         = indirect xs ys -- assert x1 == y1
indirect xs@(Ins _:_) (Ins y1:ys)        = Ins y1 : indirect xs ys

indirect [] ys                           = ys -- assert: all Ins
indirect xs []                           = xs -- assert: all Del