[BioPython] Re: translating genes with >1 exon

Cayte katel@worldpath.net
Wed, 12 Apr 2000 23:35:01 -0700


----- Original Message ----- 
From: Andrew Dalke <dalke@acm.org>
To: <biopython@biopython.org>
Sent: Monday, April 10, 2000 9:41 AM
Subject: [BioPython] Re: translating genes with >1 exon


> On the bioperl list, Keith James asked about translating
> an exon:
> > Get the PrimarySeq object of each subsequence
> > Get the sequence string of each PrimarySeq object
> > Join the sequence strings together
> > Make a new PrimarySeq object from the string
> > Translate that, making yet another PrimarySeq object
> 
> In the proposal code I've been working on, I added
> __add__ (and I should add __radd__), so you can do
> 
> s1 = Seq("ATGCATCACAATCG", Alphabet.IUPAC.unambiguous_dna)
> s2 = Seq("W", Alphabet.IUPAC.ambiguous_dna)
> t = s1 + s2
> 
> and get that t is
> 
> Seq("ATGCATCACAATCGW", IUPACAmbiguousDNA())
> 
> Thus, translation of subsets is:
>   translate(seq[5:20] + seq[29:65] + seq[100:200])
> 
   Should there be a program to guess where the exon start and end?

                                              Cayte