[Bioperl-l] IUPAC support for DNA alignment
aaron.j.mackey at gsk.com
aaron.j.mackey at gsk.com
Fri Jun 27 14:09:57 UTC 2008
You could replicate what they do here with EST_GENOME (re-engineered to
accept ambiguity codes):
http://www.genome.org/cgi/content/short/17/2/212
But I think the answer is user-dependent -- some might want the "full
score" (as in the above case), others might want the "(probabilistically)
averaged score", etc. So, let the scoring matrix be subclass-able (or
mix-able), so that users can specify the exact desired behavior via a
handful of predefined (and useful) behaviors.
-Aaron
More information about the Bioperl-l
mailing list