[Biopython] slicing an alignment

Alan Bergland Alan_Bergland at brown.edu
Fri Feb 12 18:35:59 UTC 2010


Hi all,

	I am a newbie, so my sincere apologies if this is a really naive  
question.  Is there an easy way to slice an alignment.  For instance,  
if I have imported a fasta alignment like this:

	>seq1
	ACCCGT
	>seq2
	ACCGGT

	Is there a single commmand that would slice out site 2 from the whole  
alignment and return "C, C"?  Or, do I have to call each record and  
slice individually?

Thanks!
Alan





More information about the Biopython mailing list