




Description 
Qualityaware alignments
Generally, quality data are ignored for alignment/pattern searching
like SmithWaterman, NeedlemanWunsch, or BLAST(pnx). I believe
that accounting for quality will at the very least affect things like
BLAST statistics, and e.g. is crucial for good EST annotation using Blastx.
This module performs sequences alignments, takes quality values into
account.
See also http://bioinformatics.oxfordjournals.org/cgi/content/abstract/btn052v1.


Synopsis 

local_score :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > Double   local_align :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > (Double, EditList)   global_score :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > Double   global_align :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > (Double, EditList)   overlap_score :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > Double   overlap_align :: QualMx t Double > (Double, Double) > Sequence t > Sequence t > (Double, EditList)   qualMx :: Qual > Qual > (Chr, Chr) > Double   test :: IO () 



SmithWaterman


Locally optimal alignment with affine gaps, i.e. best infix match.



Calculate local edit distance (SmithWaterman alignment score)



Calculate local alignment (SmithWaterman)
(can we replace uncurry max' with fst  a local alignment must always end on a subst, no?)


NeedlemanWunsch


Globally optimal alignment with affine gaps, the whole sequences are matched.



Calculate global edit distance (NeedlemanWunsch alignment score)



Calculate global alignment (NeedlemanWunsch)


Overlapping alignment.


The suffix of one sequence matches a prefix of another.



Calucalte best overlap score, where gaps at the edges are free
The starting point is like for local score (0 cost for initial indels),
the result is the maximum anywhere in the last column or bottom row of the matrix.



Calucalte best overlap score, where gaps at the edges are free
The starting point is like for local score (0 cost for initial indels),
the result is the maximum anywhere in the last column or bottom row of the matrix.


Matrix construction




Interactive testing of alignments




Produced by Haddock version 2.6.1 