Now, for
something completely different in this node -- this time somewhat more useless. Next is a
Python script to calculate a
DNA sequence for any stream given in
standard input.
#!/usr/bin/python
import sys
g = ('a','c','g','t')
gencode=lambda x:[ sys.stdout.write("%s%s%s%s "%(g[int(ord(y))/64],\
g[(int(ord(y))%64)/16],\
g[((int(ord(y))%64)%16)/4],\
g[((int(ord(y))%64)%16)%4]))\
for y in tuple(x)]
if __name__ == "__main__":
while 1:
try: gencode(raw_input())
except: break
print
Example session:
$ ./dna.py < dna.py
agat agac agtt ctcc ctat ctag agtt cgag cggc cgtg agtt ctaa
ctgc ctca cgga cgtt cgtg cggc cgtc ctaa cgtt ctag ctca agaa
ctat ctgc ctat cgct agaa attc agaa agga ...