in reply to Sequence alignment
I think this is best left to blast. After blast has run, you can read the resulting output with perl. BioPerl can help with (either or both) calling the blast binary, and/or reading its output.
See the bioperl HOWTOs, esp. SearchIO
|
|---|