Entering edit mode
6.3 years ago
wingsoffire72
•
0
enter code h t = open("CodonUsage.txt","r")
N=0
c=t.read(1)
while(c!='\n'):
N+=1
c=t.read(1)
#print(N)
#print("\t")
#print(N)
t.close()
print(N)
N1 =11
nucCount=([[0 for i in range(N)]for j in range(5)])
g=0
nucAlign=([['\0' for i in range(N)]for j in range(N1)])
t = open("CodonUsage.txt","r")
p=0
#print(c)
counter=[0 for i in range(10) ]
# for i in range(10):
# counter[i]=0
#
c=t.read(1)
while(c!=''):
k=0
while(c!='\n'):
nucAlign[p][k]=c
#print(p,k,c)
if c=='T' or c=='u' or c == 't' or c=='U':
nucCount[0][k] +=1
#counter[p] +=1
elif c=='c' or c=='C' :
nucCount[1][k] +=1
#counter[p] +=1
elif c=='a' or c=='A':
nucCount[2][k] +=1
#counter[p] +=1
elif c=='g' or c=='G':
nucCount[3][k] +=1
#counter[p] +=1
# else:
# nucCount[4][k] +=1
#counter[p] +=1
c=t.read(1)
k +=1
p +=1
c=t.read(1)ere