#!/usr/bin/python
import gzip
import csv
with gzip.open('UP000203464.swiss.gz', 'rt') as arch:
    prots = arch.read().split('//')
ecs = {}
for prot in prots:
    if 'CATALYTIC ACTIVITY'  not in prot:
        strings = prot.split('\n')
        for s in strings:
            subs = s.split(' ')
            for sub in subs:
                pos = sub.find('EC=')
                if pos != -1:
                    EC = sub[pos : ].replace(';', '')
                    if EC in ecs:
                        ecs[EC] += 1
                    else:
                        ecs[EC] = 1
with open('ECs_strange_all.csv', 'w') as out:
    wr = csv.writer(out, delimiter = '\t')
    for ec in ecs:
        wr.writerow([ec, str(ecs[ec])])