[Bioperl-l] FW: Hash table creator

Karger, Amir AKarger@CuraGen.com
Tue, 26 Jun 2001 09:02:50 -0400


This message is in MIME format. Since your mail reader does not understand
this format, some or all of this message may not be legible.

------_=_NextPart_000_01C0FE40.4DF3CFC0
Content-Type: text/plain;
	charset="iso-8859-1"

OK, Heikki. This may be useless, now that you've rewritten CodonTable, but
since I wasted so much time yesterday, I'll send it along anyway.

I decided to create
a program that creates a hash table to translate ANY three letter
codon into the correct amino acid, including codons with ambiguous
letters, as well as the ambiguous amino acid outputs.

The idea would be to put this into CodonTable->new() (with maybe some
caching
so that you only call it if you haven't yet called it for a 
codon table with this ID number yet).

If you run the attached code, you *should* get a list of 
amino acids, and the codons that can possibly create them.
(Note that the ambiguous amino acids can only be created by
ambiguous codons.)

Surprisingly, even though it has to go through 4900 or so combinations,
it only takes a second or so to run.


Amir Karger
CuraGen Corporation



------_=_NextPart_000_01C0FE40.4DF3CFC0
Content-Type: APPLICATION/octet-stream;
	name="amino_ambiguous.pl.gz"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
	filename="amino_ambiguous.pl.gz"
Content-Description: amino_ambig
Content-ID: <Pine.SOL.3.96.1010625184201.14741C@everest.curagen.com>

H4sICIu8NzsAA2FtaW5vX2FtYmlndW91cy5wbACNVtt2EzcUfc58xWHsxJNiakLe4oYSSJpmtaTQ
QEsL1EuWlbFie8aMNDhZXvPv3brMzW5C/eA90j46d106jwa5ygZjmQyWIpvTk1UQ5EqQ0pnkemi/
eZoozRJNp2cv35+fX1ye0/FzejoMgsUd7cp8yfhokjA6Dgi/iHqsZwQ+fllFxMj99j/3A/fV4zXL
t9m4ZuNtVtes3mbzB9lF0yu+yWZNNt5kV01Wb7Kzps9b7F0z3i1WNdktu1/bPsdtdtpmdZudtCPa
YMdtuxtssmlXN9nbB9h93xZKLpZzMWILmaR0THux0COeTtBKEUSMxFhoLbJagk0mGIxlnEt9F31q
aTBaO/Qmk+jCNNe0mko+JafPgKDrNPOzjCmrPxNfodZYAycYGMx27RqKZuJOtV3Yp7WNwAgxNDN1
m+zaLSyGViYRt5rktRUUXyj8ELr5Za6m9GLdhek1uKLoe4PDoAiW1vvwyoZV+n50RCc25jRXfu5T
8uRbv09J2I4KfkQqzTS5uGC/DMdbhcQRkXcTK154+8ctd2vaJR90nIklDT4yruPP68Ni0C+XNlUZ
2dtS+NH90s6Zm1SiAiGF/dLOfp9CMrkoPdwW9Eb2+6ENvkA//N9foPIxtbqrWepMXMNzNZXXugpp
t2zKXcvXoc5tS5SRrgcf/8k/Dwqf9eokHFKHcpKK8iQRXCjFsjurokNXQitKr0mvUnJGGJcT1Whf
PRWZ6CliCTGlUi6ZFpPKsUmaj21dovD0MjTbMHyJ9ISXp43B2Vs3+NsM3p6VA7OJyrBtgsUEip66
6JrdJCnykZou2mlSN23q5PXLi/P3v72/OmopmLWldmyiyzZ3xwAsh13ZvenOUE0nkeSLETav9Qmp
wnCMXCNZZlJPpd8gxJGaBZsJt2zXLTH73Bsymze0Wjv0Ck5pQWw+p2WqlDTJw9WGVkIdVlJPServ
reTFNT4r3bRAODCKYfSMUsRx2KxXH7I9Zdfp1FTSx4aZViIPkIh1t2oMfBaFS8lmzm+2RG9K0ZbK
2ZbcrJIriwvpeTPLB9COhTYlO64PWfvQqeO2ea4jLVf4I9GfhbWBotKZJ3N0OkXiVirkFguUPVS8
bzvVxOPHQzfh612OcZ5GdQ8c02G50tl3lje6qIBTvQ+9YSXXofMUJXEndB0kNmglYqmqcd3SoozC
fhRBPSigUiQT4qaRZBL/Vyu57ZsmYjOtQdlZK0Fxqsls7b4ZVOkW0ux3uxYq7KnAlGnIzWw8q9um
Q6epMDsCzzShy9PEH+hTtnF0NIr5Y1D3yCrNzO43R2wU4pxwRxiM2X1kpIwHZTHdsbO2q8qKPlCS
tnypsCAxx3PyW4t9PesCuFUm7ptcaZurqgCMedEIrbWP1VUcpY6H3GRs2CixFuq+Cuv7S1w52qGJ
GOdxDBUY4dEMZblxvJF11jz2WkRze7EqZeWmqh8aPvnVtW6X4vZEtarnue3pxhF/0MpnEPguulrO
pX1OLXPdpzGeVbYlVY6TL0G3Yhs5M2KxxJ05l4lQbq23jpvYvoJKS3t7LR9s9lsXDYw7N8x/JnSe
Jf6uxY3uq2Aua7vZ/NNvhI6eBvYOr9+Q9QVeXdXmAjBzuEVRtuiH05N3J89ttiAVuWdYn3yXKBP6
0L9n6x5pd8b9rgbBaGT0j0ZBHONRfg5gDrQDboGBOwMwB+BOAdyCBvcHgDnQDrgFDu4EwBxoB9wA
M/Z+BzAH4K4A3IKx9wuAOQB3CeAWjL3XAHAXAO2AWzD23gGYA+2AG9DG3p8AcN8BwL0CcAvGHiaZ
44y9vwDcgrH3K4A5APcTgFsw9q4AzIF2YGPgLj7u4uPGngFuwdh7C2AOwP0M4BacPe7scWPPALdg
7L0BMAfaATfwLzHSjXMBDwAA

------_=_NextPart_000_01C0FE40.4DF3CFC0--