-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathadd-names.py
100 lines (71 loc) · 2.16 KB
/
add-names.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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
#this script adds species names to fasta headers
#add species namaes to Aedes aegypti file
fout1 = open("./named-families/Aaeg-named.fa","w")
for line in open("./TE-Families/ae-aegypti-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Aaeg_" + oldline
fout1.write(newline)
else:
fout1.write(line)
fout1.close()
#add species names to Aedes albopictus file
fout2 = open("./named-families/Aalb-named.fa","w")
for line in open("./TE-Families/ae-albopictus-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Aalb_" + oldline
fout2.write(newline)
else:
fout2.write(line)
fout2.close()
#add species names to Anopheles coluzzii file
fout3 = open("./named-families/Acol-named.fa","w")
for line in open("./TE-Families/an-coluzzii-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Acol_" + oldline
fout3.write(newline)
else:
fout3.write(line)
fout3.close()
#add species names to Anopheles gambiae file
fout4 = open("./named-families/Agam-named.fa","w")
for line in open("./TE-Families/an-gambiae-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Agam_" + oldline
fout4.write(newline)
else:
fout4.write(line)
fout4.close()
#add species names to Culex quinquefasciatus file
fout5 = open("./named-families/Cqui-named.fa","w")
for line in open("./TE-Families/cu-qu-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Cqui_" + oldline
fout5.write(newline)
else:
fout5.write(line)
fout5.close()
#add species names to Culex tarsalis file
fout6 = open("./named-families/Ctar-named.fa","w")
for line in open("./TE-Families/cu-tarsalis-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Ctar_" + oldline
fout6.write(newline)
else:
fout6.write(line)
fout6.close()
#add species names to Wyeomyia smithii file
fout7 = open("./named-families/Wsmi-named.fa","w")
for line in open("./TE-Families/wyeomyia-families.fa"):
if line.startswith('>'):
oldline = line.split('>')[1]
newline = ">Wsmi_" + oldline
fout7.write(newline)
else:
fout7.write(line)
fout7.close()