view test/usecase1.py @ 668:3bdcee61356d
blocks3d/www: fix log-file caused bug of building
r2w was configured to use /tmp/blocks3d-www-build.log.
It caused build error if this file has been creted by other user.
Log-file path was changed to ./output/blocks3d-www-build.log
Makefile was changed to create ./output directory before running r2w
| author |
boris (kodomo) <bnagaev@gmail.com> |
| date |
Fri, 01 Jul 2011 10:32:37 +0400 |
| parents |
07b5351b0b56 |
| children |
ddf85d0a8924 |
line source
2 from allpy import protein
3 from allpy import processors
5 # Create sequences from string representation of sequence body
6 sequence_1 = protein.Sequence.from_string("mkstf", name="E2E4")
7 sequence_2 = protein.Sequence.from_string("mstkfff", description="Longer sequence")
9 # Create alignment from sequences
10 alignment = protein.Alignment()
11 alignment.append_sequence(sequence_1)
12 alignment.append_sequence(sequence_2)
13 alignment.realign(processors.Muscle())
15 # For each sequence, print number of gaps and non-gaps in alignment
16 for row in alignment.rows():
19 for column in alignment.columns:
24 print "%s: %s gaps, %s non-gaps" % (row.sequence.name, gaps, monomers)
26 # Print number of gaps in each column
28 for column in alignment.columns:
30 for sequence in alignment.sequences:
31 if sequence not in column:
33 gaps.append(column_gaps)
34 print " ".join(map(str, gaps))
36 # Write alignment to file
37 alignment.to_file(sys.stdout)
