我怎样才能缩短这个节目??
选择 | 换行 | 行号
- 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)