1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
import pprint, sys
f = open('/u1/junk/kinne/GRCh38.p12.genome.fa', 'r')
count = {}
for line in f:
    for character in line:
        count.setdefault(character, 0)
        count[character] = count[character] + 1
f.close()
pprint.pprint(count)
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX