Code:
map = {}
dna = [ "A", "C", "G", "T"]
with open("matrix.txt", "r") as myfile:
line = myfile.readline()
mylist = line.split()
map{"AA"} = mylist[0]
map{"AC"} = mylist[1]
map{"AG"} = mylist[2]
map{"AT"} = mylist[3]
line = myfile.readline()
mylist = line.split()
map{"CA"} = mylist[0]
map{"CC"} = mylist[1]
map{"CG"} = mylist[2]
map{"CT"} = mylist[3]
line = myfile.readline()
mylist = line.split()
map{"GA"} = mylist[0]
map{"GC"} = mylist[1]
map{"GG"} = mylist[2]
map{"GT"} = mylist[3]
line = myfile.readline()
mylist = line.split()
map{"TA"} = mylist[0]
map{"TC"} = mylist[1]
map{"TG"} = mylist[2]
map{"TT"} = mylist[3]
gap = -2
map{"A-"} = gap
map{"C-"} = gap
map{"T-"} = gap
map{"G-"} = gap
map{"-A"} = gap
map{"-C"} = gap
map{"-T"} = gap
map{"-G"} = gap
with open("dna.fasta.txt","r") as myfile:
data = myfile.readlines()
myfile.close()
for i in range(len(data)):
data[i] = data[i].rstrip("\n")
seq1 = data[1]
seq2 = data[3]
score = 0
for i in range(0,len(seq1),1):
key = seq1[i] + seq2[i]
score = score + int(map[key])
print(score)
Comment