Revision as of 09:51, 20 January 2007 by Jblucks
from Bio.Seq import Seq,translate #create a sequence object of some DNA my_seq = Seq('CATGTAGATAG') #print out some details about it print 'seq is %i bases long' % len(my_seq) print 'reverse complement is %s' % my_seq.reverse_complement().tostring() #or see the whole record print 'sequence record:', my_seq #translate the sequence into a protein my_protein = translate(my_seq) print 'protein translation is %s' % my_protein.tostring() print 'protein record:', my_protein
seq is 11 bases long reverse complement is CTATCTACATG sequence record: Seq('CATGTAGATAG', Alphabet()) protein translation is HVD protein record: Seq('HVD', HasStopCodon(IUPACProtein(), '*'))
- Learn how to program in Python
- Browse the Biopython Tutorial
- Examine the Class Diagram if you'd like to know more about the relationships between the modules.
- Use the Wiki Search tools to find more information on specific topics.