How can I generate a pairwise alignment of two s...

To see other biotech frequently asked questions, please visit http://biotech.fyicenter.com/faq/

(Continued from previous question...)

How can I generate a pairwise alignment of two sequences?

Look at Bio::Factory::EMBOSS to see how to use the water and needle alignment programs that are part of the EMBOSS suite. Bio::Factory::EMBOSS is part of the bioperl-run package.

Or you can use the pSW module for DNA alignments or the dpAlign module for protein alignments. These are part of the bioperl-ext package; download it via Getting BioPerl.

You can also use prss34 (part of FASTA package) to assess the significance of a pairwise alignment by shuffling the sequences.

(Continued on next question...)

Other Frequently Asked Questions

• What is this FAQ?
• How do I retrieve all the features from a sequence?
• Where do I go to get the latest release?
• Who uses this toolkit?
• Cannot get an accession from GenBank when I know it is there
• I can't get sequences with Bio::DB::GenBank any more, why not?
• Why can't I easily get a list of all the methods a object can call?
• Does Bio::SearchIO parse the HTML output that BLAST creates using the -T option?
• I get the warning (old style Annotation) on new style Annotation::Collection. What is wrong?
• How do I parse the CDS join or complement statements in GenBank or EMBL files to get the sub-locations?
• More...