[Biojava-l] Short names for Amino acid symbols

community at struck.lu community at struck.lu
Tue Aug 5 12:33:58 UTC 2008


The link should have been:

http://biojava.org/wiki/BioJava:Cookbook:Translation:OneLetterAmbi

Hope this time my webmail doesn't garble the message :-(

Greetings,
Daniel

"George Waldon" <gwaldon at geneinfinity.org> wrote:

> The link
> http://biojava.org/wiki/BioJava:Cookbook:Translation:OneLetterAmbiDaniel
does
> not seem to work.
> 
> Did you try SymbolTonenization? Something like:
> 
> Symbol s;
> SymbolTokenization  tok=
ProteinTools.getTAlphabet().getTokenization("token");
> String 	s = tok.tokenizeSymbol(s);
> 
> Should give you the short name of any given symbol.
> 
> - George
> 
> > -----Original Message-----
> > From: biojava-l-bounces at lists.open-bio.org [mailto:biojava-l-
> > bounces at lists.open-bio.org] On Behalf Of Peter Robinson
> > Sent: Sunday, July 27, 2008 8:57 AM
> > To: biojava-l at lists.open-bio.org
> > Subject: [Biojava-l] Short names for Amino acid symbols
> > 
> > Hi,
> > 
> > thanks to all on the list who helped me get started with Biojava, and by
> > the way, the online documents are quite helpful!
> > 
> > I am trying to develop some code to look for signs of positive selection
> > in human sequences by making multiple alignments of protein sequences
> > and mapping the nucleotide sequences onto this alignment and checking
> > synonymous and nonsynonymous nucleotide substitutions in several species
> > (etc).
> > 
> > A few small questions;
> > 1) I have written a class to encapsulate all I need from a given Genbank
> > mRNA sequence; the entire mRNA, the CDS and the corresponding protein
> > sequence. I have some methods such as the following:
> > 
> >  private void setCDSSequence() {
> >         Feature CDS = getCDSFeature(this.completeSequence);
> >         Location loc = CDS.getLocation();
> >         SymbolList symL = this.completeSequence.subList(loc.getMin(),
> > loc.getMax()-3); //-3 to remove stop codon
> >         this.CDS= symL;
> >     }
> > 
> > Question: Why is there (seemingly) no way in Biojava to create a
> > Sequence object instead of a SymbolList object? Or did I miss something?
> > 
> > 2)  I would then like to printout the protein alignment to check for
> > correctness, and it seems there is no way of getting from a symbol to
> > the one-letter aminoacid code. That is,
> > 
> > proteinAlignment.get(j).symbolAt(k).getName()
> > 
> > will return "Ala" instead of "A" etc. Is there a good way of getting the
> > short symbols?
> > 
> > Thanks, Peter
> > _______________________________________________
> > Biojava-l mailing list  -  Biojava-l at lists.open-bio.org
> > http://lists.open-bio.org/mailman/listinfo/biojava-l
> _______________________________________________
> Biojava-l mailing list  -  Biojava-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/biojava-l
> 

_________________________________________________________
Mail sent using root eSolutions Webmailer - www.root.lu





More information about the Biojava-l mailing list