def transcript(numseq):
    protseq = ''
    for l in range(0, len(numseq), 2):
        protseq += letters[ int(numseq[l : l+2]) % 20 ]
    return protseq

letters = ['V', 'A', 'R', 'N', 'D', 'C', 'Q', 'E', 'G', 'H', 'I', 'L', 'K', 'M', 'F', 'P', 'S', 'T', 'W', 'Y']

pi = '''31415926535897932384626433832795028841971693993751058209749445923078164062862
089986280348253421170679821480865132823066470938446095505822317253594081284811
174502841027019385211055596446229489549303819644288109756659334461284756482337
867831652712019091456485669234603486104543266482133936072602491412737245870066
063155881748815209209628292540917153643678925903600113305305488204665213841469
519415116094330572703657595919530921861173819326117931051185480744623799627495
673518857527248912279381830119491298336733624406566430860213949463952247371907
021798609437027705392171762931767523846748184676694051320005681271452635608277
857713427577896091736371787214684409012249534301465495853710507922796892589235
420199561121290219608640344181598136297747713099605187072113499999983729780499
510597317328160963185950244594553469083026425223082533446850352619311881710100
031378387528865875332083814206171776691473035982534904287554687311595628638823
537875937519577818577805321712268066130019278766111959092164201989380952572010
654858632788659361533818279682303019520353018529689957736225994138912497217752
834791315155748572424541506959508295331168617278558890750983817546374649393192
550604009277016711390098488240128583616035637076601047101819429555961989467678
374494482553797747268471040475346462080466842590694912933136770289891521047521
620569660240580381501935112533824300355876402474964732639141992726042699227967
823547816360093417216412199245863150302861829745557067498385054945885869269956
909272107975093029553211653449872027559602364806654991198818347977535663698074
265425278625518184175746728909777727938000816470600161452491921732172147723501
414419735685481613611573525521334757418494684385233239073941433345477624168625
189835694855620992192221842725502542568876717904946016534668049886272327917860
857843838279679766814541009538837863609506800642251252051173929848960841284886'''.split()
pi = ''.join(pi)
numseq = pi[:2080]
numseq

protseq = transcript(numseq)

with open('pi.fasta', 'w') as file:
    file.write('>pi first 2080 digits of Pi written as aminoacids \n')
    for a in range(0, 2000, 80):
        file.write(protseq[a:a+80]+'\n')