

CodonLookupDictionary = {"GAC":.20
                   ,"GCU":.21
                   ,"CGU":.22
                   ,"AAA":.23
                   ,"UCU":.24
                   ,"GAC":.25
                   ,"CAG":.26
                   ,"AAC":.27
                   ,"GUA":.28
                   ,"CGU":.29
                    }
CodonProduct = 1

CodeString = open("GenCode.txt", "r").read()
CodeString = CodeString.replace("\n", "")
CodeString = CodeString.replace(" ", "")
CodeString = CodeString.replace("t", "u").upper()

print
print "Source", len(CodeString), "chars:", CodeString
print

TmpString = CodeString[:150]
print "First", len(TmpString), "chars:", TmpString
print

TmpString = CodeString[150:-60]
print "Target", len(TmpString), "chars:"

CodonCount = len(TmpString)/3

while TmpString:
    Codon = TmpString[:3]
    CodonValue = CodonLookupDictionary[Codon]
    CodonProduct = CodonProduct * CodonValue
    print " ", Codon, "%5.4f" % CodonValue
    TmpString = TmpString[3:]
    
print
print "Product = %10.9f" % CodonProduct
print "MagicNumber = %10.9f" % (CodonProduct/CodonCount)
print

TmpString = CodeString[-60:]
print "Final", len(TmpString), "chars:", TmpString
print

