[Bioperl-l] Translate sequences align and reverse translate

[Jason Stajich]

bp_pairwise_kaks does most of this for you (you can use clustalw and  
muscle).

  at the core is aa_to_dna_align() routine from Bio::Align::Utilities


-jason
On Jan 18, 2007, at 5:44 AM, Anthony Underwood wrote:

> Hi All,
>
>
>
> I would like to automate the process of taking coding DNA sequences,
> translating them to amino acids, aligning them with clustalw and then
> reverse-translating back to DNA so as to obtain the best alignment. I
> can think of how I would do this with bioperl but since it is  
> probably a
> common process  I would like to ask if anyone already has such a  
> script
> that they wouldn't mind sharing so I don't re-invent the wheel.
>
>
>
> Many thanks,
>
>
>
> Anthony
>
>
>
>
> -----------------------------------------
> *******************************************************************
> *******
> The information contained in the EMail and any attachments is
> confidential and intended solely and for the attention and use of
> the named addressee(s). It may not be disclosed to any other person
> without the express authority of the HPA, or the intended
> recipient, or both. If you are not the intended recipient, you must
> not disclose, copy, distribute or retain this message or any part
> of it. This footnote also confirms that this EMail has been swept
> for computer viruses, but please re-sweep any attachments before
> opening or saving. HTTP://www.HPA.org.uk **************************
> ************************************************
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l

--
Jason Stajich
Miller Research Fellow
University of California, Berkeley
lab: 510.642.8441
http://pmb.berkeley.edu/~taylor/people/js.html