-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathGC.py
25 lines (20 loc) · 781 Bytes
/
GC.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
def main():
file_name = "ComputingGCContent.txt"
file = open(file_name, "r")
fasta_file = file.read()
sequencesWithName = fasta_file.split(">")
if sequencesWithName[0] == '':
sequencesWithName = sequencesWithName[1:]
for i in range(len(sequencesWithName)):
sequencesWithName[i] = sequencesWithName[i][:13] + ' ' + sequencesWithName[i][14:]
sequencesWithName[i] = sequencesWithName[i].replace("\n", '')
maxValu = 0
for i in range(len(sequencesWithName)):
percent = ((sequencesWithName[i].count('G')+sequencesWithName[i].count('C'))/len(sequencesWithName[i][14:]))*100
if percent > maxValu:
i_max = i
maxValu = percent
print(sequencesWithName[i_max][:13])
print("%0.6f" % maxValu)
if __name__ == "__main__":
main()