[BioPython] help for local alignment

Peter biopython at maubp.freeserve.co.uk
Mon Dec 22 17:47:16 UTC 2008


On Mon, Dec 22, 2008 at 4:38 PM, Chandan Kumar <biopython.chen at gmail.com> wrote:
> Dear all,
>                 can any one provide me simple code for local alignment
> python code which can be applied for protein or nucleotide sequence. Please
> provide me the simplest code as I am new to python and from biology
> background.
>
> Thanking you.
>
> Kind regards
> Chen

Hi Chen,

Are you wanting to do pairwise alignments (aligning two sequences to
each other), or multiple sequence alignments?

For multiple sequence alignments, you might want to use a 3rd party
tool like ClustalW, or MUSCLE. Biopython can parse several alignment
formats including ClustalW format.  See our tutorial for examples
using ClustalW.

Biopython's Bio.pairwise2 can do pairwise alignments, although we only
have the built in documentation for this at the moment (nothing in our
tutorial).  This documentation is also available online:
http://biopython.org/DIST/docs/api/Bio.pairwise2-module.html

For pairwise sequence alignments I personally use the EMBOSS tools
"water" (Smith-Waterman algorithm for local alignment) or "needle"
(Needleman-Wunsch for global alignment).  Biopython's Bio.AlignIO
module can parse their output.

Peter



More information about the Biopython mailing list