-
Notifications
You must be signed in to change notification settings - Fork 0
/
6) Filtering-for-synonymous_tuberculosis.R
495 lines (487 loc) · 27.3 KB
/
6) Filtering-for-synonymous_tuberculosis.R
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
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
#Filtering out only synonymous rows from Snippy generated csv file
my_data <- read.csv("snps-samples_aligned_to_m.tb.csv")
my_data1 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1011C>T p.Gly337Gly")
my_data2 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1026A>G p.Glu342Glu")
my_data3 = subset(my_data, my_data$EFFECT == "synonymous_variant c.10353C>A p.Gly3451Gly")
my_data4 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1035T>C p.Gly345Gly")
my_data5 = subset(my_data, my_data$EFFECT == "synonymous_variant c.105G>A p.Ser35Ser")
my_data6 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1068G>A p.Gly356Gly")
my_data7 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1068G>A p.Leu356Leu")
my_data8 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1089A>G p.Lys363Lys")
my_data9 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1107G>A p.Gly369Gly")
my_data10 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1110G>A p.Arg370Arg")
my_data11 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1119C>T p.Gly373Gly")
my_data12 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1155C>T p.Ala385Ala")
my_data13 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1194G>C p.Gly398Gly")
my_data14 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1197C>G p.Thr399Thr")
my_data15 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1212G>C p.Leu404Leu")
my_data16 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1221G>A p.Val407Val")
my_data17 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1227A>C p.Gly409Gly")
my_data18 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1230T>G p.Gly410Gly")
my_data19 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1236C>A p.Gly412Gly")
my_data20 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1236C>G p.Arg412Arg")
my_data21 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1236G>A p.Ala412Ala")
my_data22 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1242G>T p.Val414Val")
my_data23 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1251G>T p.Pro417Pro")
my_data24 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1275C>T p.Gly425Gly")
my_data25 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1299C>T p.Gly433Gly")
my_data26 = subset(my_data, my_data$EFFECT == "synonymous_variant c.12C>T p.Asp4Asp")
my_data27 = subset(my_data, my_data$EFFECT == "synonymous_variant c.12G>A p.Lys4Lys")
my_data28 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1317T>C p.Arg439Arg")
my_data29 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1324C>T p.Leu442Leu")
my_data30 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1344C>A p.Gly448Gly")
my_data31 = subset(my_data, my_data$EFFECT == "synonymous_variant c.135G>A p.Ser45Ser")
my_data32 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1368C>T p.Asp456Asp")
my_data33 = subset(my_data, my_data$EFFECT == "synonymous_variant c.136C>T p.Leu46Leu")
my_data34 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1395G>C p.Ala465Ala")
my_data35 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1401G>A p.Pro467Pro")
my_data36 = subset(my_data, my_data$EFFECT == "synonymous_variant c.147C>G p.Ala49Ala")
my_data37 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1485G>A p.Pro495Pro")
my_data38 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1485T>C p.Arg495Arg")
my_data39 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1504T>C p.Leu502Leu")
my_data40 = subset(my_data, my_data$EFFECT == "synonymous_variant c.150A>G p.Gln50Gln")
my_data41 = subset(my_data, my_data$EFFECT == "synonymous_variant c.150T>C p.Gly50Gly")
my_data42 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1536G>A p.Pro512Pro")
my_data43 = subset(my_data, my_data$EFFECT == "synonymous_variant c.153G>T p.Leu51Leu")
my_data44 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1572A>C p.Gly524Gly")
my_data45 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1585C>T p.Leu529Leu")
my_data46 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1596A>C p.Gly532Gly")
my_data47 = subset(my_data, my_data$EFFECT == "synonymous_variant c.15T>C p.Val5Val")
my_data48 = subset(my_data, my_data$EFFECT == "synonymous_variant c.1629C>T p.Gly543Gly")
my_data49 = subset(my_data, my_data$EFFECT == "synonymous_variant c.165C>G p.Ala55Ala")
my_data50 = subset(my_data, my_data$EFFECT == "synonymous_variant c.165C>T p.Ser55Ser")
my_data51 = subset(my_data, my_data$EFFECT == "synonymous_variant c.168C>T p.Arg56Arg")
my_data52 = subset(my_data, my_data$EFFECT == "synonymous_variant c.169T>C p.Leu57Leu")
my_data53 = subset(my_data, my_data$EFFECT == "synonymous_variant c.171C>G p.Gly57Gly")
my_data54 = subset(my_data, my_data$EFFECT == "synonymous_variant c.175C>T p.Leu59Leu")
my_data55 = subset(my_data, my_data$EFFECT == "synonymous_variant c.177G>A p.Pro59Pro")
my_data56 = subset(my_data, my_data$EFFECT == "synonymous_variant c.177T>C p.Gly59Gly")
my_data57 = subset(my_data, my_data$EFFECT == "synonymous_variant c.183C>G p.Ser61Ser")
my_data58 = subset(my_data, my_data$EFFECT == "synonymous_variant c.183G>A p.Ala61Ala")
my_data59 = subset(my_data, my_data$EFFECT == "synonymous_variant c.186T>A p.Pro62Pro")
my_data60 = subset(my_data, my_data$EFFECT == "synonymous_variant c.18A>G p.Gly6Gly")
my_data61 = subset(my_data, my_data$EFFECT == "synonymous_variant c.192G>A p.Ala64Ala")
my_data62 = subset(my_data, my_data$EFFECT == "synonymous_variant c.195A>G p.Thr65Thr")
my_data63 = subset(my_data, my_data$EFFECT == "synonymous_variant c.19C>A p.Arg7Arg")
my_data64 = subset(my_data, my_data$EFFECT == "synonymous_variant c.201C>T p.Thr67Thr")
my_data65 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2034G>C p.Arg678Arg")
my_data66 = subset(my_data, my_data$EFFECT == "synonymous_variant c.204A>C p.Arg68Arg")
my_data67 = subset(my_data, my_data$EFFECT == "synonymous_variant c.204T>C p.Arg68Arg")
my_data68 = subset(my_data, my_data$EFFECT == "synonymous_variant c.204T>C p.Phe68Phe")
my_data69 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2055T>C p.Asn685Asn")
my_data70 = subset(my_data, my_data$EFFECT == "synonymous_variant c.207C>T p.Ile69Ile")
my_data71 = subset(my_data, my_data$EFFECT == "synonymous_variant c.210G>T p.Gly70Gly")
my_data72 = subset(my_data, my_data$EFFECT == "synonymous_variant c.216G>A p.Pro72Pro")
my_data73 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2187T>C p.Gly729Gly")
my_data74 = subset(my_data, my_data$EFFECT == "synonymous_variant c.219G>A p.Arg73Arg")
my_data75 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2205C>T p.Gly735Gly")
my_data76 = subset(my_data, my_data$EFFECT == "synonymous_variant c.225G>A p.Ala75Ala")
my_data77 = subset(my_data, my_data$EFFECT == "synonymous_variant c.225T>C p.Asp75Asp")
my_data78 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2340G>A p.Gln780Gln")
my_data79 = subset(my_data, my_data$EFFECT == "synonymous_variant c.237C>T p.Gly79Gly")
my_data80 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2403G>T p.Val801Val")
my_data81 = subset(my_data, my_data$EFFECT == "synonymous_variant c.243A>G p.Gly81Gly")
my_data82 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2481C>A p.Ala827Ala")
my_data83 = subset(my_data, my_data$EFFECT == "synonymous_variant c.249G>C p.Thr83Thr")
my_data84 = subset(my_data, my_data$EFFECT == "synonymous_variant c.24C>T p.Ala8Ala")
my_data85 = subset(my_data, my_data$EFFECT == "synonymous_variant c.24G>A p.Arg8Arg")
my_data86 = subset(my_data, my_data$EFFECT == "synonymous_variant c.252C>T p.Ala84Ala")
my_data87 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2646C>T p.Tyr882Tyr")
my_data88 = subset(my_data, my_data$EFFECT == "synonymous_variant c.2691T>C p.Tyr897Tyr")
my_data89 = subset(my_data, my_data$EFFECT == "synonymous_variant c.270C>T p.Ala90Ala")
my_data90 = subset(my_data, my_data$EFFECT == "synonymous_variant c.276C>T p.Ile92Ile")
my_data91 = subset(my_data, my_data$EFFECT == "synonymous_variant c.285C>T p.Arg95Arg")
my_data92 = subset(my_data, my_data$EFFECT == "synonymous_variant c.285C>T p.Leu95Leu")
my_data93 = subset(my_data, my_data$EFFECT == "synonymous_variant c.291C>T p.Val97Val")
my_data94 = subset(my_data, my_data$EFFECT == "synonymous_variant c.291G>A p.Lys97Lys")
my_data95 = subset(my_data, my_data$EFFECT == "synonymous_variant c.297T>C p.Tyr99Tyr")
my_data96 = subset(my_data, my_data$EFFECT == "synonymous_variant c.303G>A p.Ser101Ser")
my_data97 = subset(my_data, my_data$EFFECT == "synonymous_variant c.304T>C p.Leu102Leu")
my_data98 = subset(my_data, my_data$EFFECT == "synonymous_variant c.306C>T p.Arg102Arg")
my_data99 = subset(my_data, my_data$EFFECT == "synonymous_variant c.309C>T p.Pro103Pro")
my_data100 = subset(my_data, my_data$EFFECT == "synonymous_variant c.30C>G p.Ala10Ala")
my_data101 = subset(my_data, my_data$EFFECT == "synonymous_variant c.30C>G p.Thr10Thr")
my_data102 = subset(my_data, my_data$EFFECT == "synonymous_variant c.30C>T p.Leu10Leu")
my_data103 = subset(my_data, my_data$EFFECT == "synonymous_variant c.312C>T p.His104His")
my_data104 = subset(my_data, my_data$EFFECT == "synonymous_variant c.312T>C p.Asp104Asp")
my_data105 = subset(my_data, my_data$EFFECT == "synonymous_variant c.318C>G p.Ala106Ala")
my_data106 = subset(my_data, my_data$EFFECT == "synonymous_variant c.318C>T p.Ser106Ser")
my_data107 = subset(my_data, my_data$EFFECT == "synonymous_variant c.318C>T p.Val106Val")
my_data108 = subset(my_data, my_data$EFFECT == "synonymous_variant c.3192C>T p.Thr1064Thr")
my_data109 = subset(my_data, my_data$EFFECT == "synonymous_variant c.330A>G p.Gly110Gly")
my_data110 = subset(my_data, my_data$EFFECT == "synonymous_variant c.336C>T p.Ser112Ser")
my_data111 = subset(my_data, my_data$EFFECT == "synonymous_variant c.33G>C p.Val11Val")
my_data112 = subset(my_data, my_data$EFFECT == "synonymous_variant c.342G>A p.Ala114Ala")
my_data113 = subset(my_data, my_data$EFFECT == "synonymous_variant c.345G>A p.Gln115Gln")
my_data114 = subset(my_data, my_data$EFFECT == "synonymous_variant c.3513T>C p.Ser1171Ser")
my_data115 = subset(my_data, my_data$EFFECT == "synonymous_variant c.354C>G p.Thr118Thr")
my_data116 = subset(my_data, my_data$EFFECT == "synonymous_variant c.363A>C p.Thr121Thr")
my_data117 = subset(my_data, my_data$EFFECT == "synonymous_variant c.363C>T p.Pro121Pro")
my_data118 = subset(my_data, my_data$EFFECT == "synonymous_variant c.363G>A p.Ala121Ala")
my_data119 = subset(my_data, my_data$EFFECT == "synonymous_variant c.366G>A p.Thr122Thr")
my_data120 = subset(my_data, my_data$EFFECT == "synonymous_variant c.366G>C p.Arg122Arg")
my_data121 = subset(my_data, my_data$EFFECT == "synonymous_variant c.369G>C p.Leu123Leu")
my_data122 = subset(my_data, my_data$EFFECT == "synonymous_variant c.378C>T p.Ala126Ala")
my_data123 = subset(my_data, my_data$EFFECT == "synonymous_variant c.390A>G p.Glu130Glu")
my_data124 = subset(my_data, my_data$EFFECT == "synonymous_variant c.396C>A p.Gly132Gly")
my_data125 = subset(my_data, my_data$EFFECT == "synonymous_variant c.396G>A p.Leu132Leu")
my_data126 = subset(my_data, my_data$EFFECT == "synonymous_variant c.399C>T p.Ala133Ala")
my_data127 = subset(my_data, my_data$EFFECT == "synonymous_variant c.403T>C p.Leu135Leu")
my_data128 = subset(my_data, my_data$EFFECT == "synonymous_variant c.405C>T p.Gly135Gly")
my_data129 = subset(my_data, my_data$EFFECT == "synonymous_variant c.423T>C p.Arg141Arg")
my_data130 = subset(my_data, my_data$EFFECT == "synonymous_variant c.429T>G p.Gly143Gly")
my_data131 = subset(my_data, my_data$EFFECT == "synonymous_variant c.42C>T p.Arg14Arg")
my_data132 = subset(my_data, my_data$EFFECT == "synonymous_variant c.432T>C p.Asp144Asp")
my_data133 = subset(my_data, my_data$EFFECT == "synonymous_variant c.435C>T p.Gly145Gly")
my_data134 = subset(my_data, my_data$EFFECT == "synonymous_variant c.435C>T p.Val145Val")
my_data135 = subset(my_data, my_data$EFFECT == "synonymous_variant c.444A>G p.Ala148Ala")
my_data136 = subset(my_data, my_data$EFFECT == "synonymous_variant c.444G>C p.Arg148Arg")
my_data137 = subset(my_data, my_data$EFFECT == "synonymous_variant c.447G>A p.Gly149Gly")
my_data138 = subset(my_data, my_data$EFFECT == "synonymous_variant c.450A>G p.Gly150Gly")
my_data139 = subset(my_data, my_data$EFFECT == "synonymous_variant c.453C>T p.Leu151Leu")
my_data140 = subset(my_data, my_data$EFFECT == "synonymous_variant c.459C>T p.Asp153Asp")
my_data141 = subset(my_data, my_data$EFFECT == "synonymous_variant c.45G>A p.Ala15Ala")
my_data142 = subset(my_data, my_data$EFFECT == "synonymous_variant c.468G>A p.Ser156Ser")
my_data143 = subset(my_data, my_data$EFFECT == "synonymous_variant c.483T>C p.Gly161Gly")
my_data144 = subset(my_data, my_data$EFFECT == "synonymous_variant c.486C>T p.Phe162Phe")
my_data145 = subset(my_data, my_data$EFFECT == "synonymous_variant c.48C>G p.Leu16Leu")
my_data146 = subset(my_data, my_data$EFFECT == "synonymous_variant c.4926C>T p.Gly1642Gly")
my_data147 = subset(my_data, my_data$EFFECT == "synonymous_variant c.492A>G p.Ala164Ala")
my_data148 = subset(my_data, my_data$EFFECT == "synonymous_variant c.507C>G p.Ala169Ala")
my_data149 = subset(my_data, my_data$EFFECT == "synonymous_variant c.510C>T p.Asn170Asn")
my_data150 = subset(my_data, my_data$EFFECT == "synonymous_variant c.519C>G p.Val173Val")
my_data151 = subset(my_data, my_data$EFFECT == "synonymous_variant c.532C>T p.Leu178Leu")
my_data152 = subset(my_data, my_data$EFFECT == "synonymous_variant c.540A>G p.Ala180Ala")
my_data153 = subset(my_data, my_data$EFFECT == "synonymous_variant c.540A>G p.Glu180Glu")
my_data154 = subset(my_data, my_data$EFFECT == "synonymous_variant c.546C>T p.Ala182Ala")
my_data155 = subset(my_data, my_data$EFFECT == "synonymous_variant c.558G>A p.Ala186Ala")
my_data156 = subset(my_data, my_data$EFFECT == "synonymous_variant c.561G>C p.Ser187Ser")
my_data157 = subset(my_data, my_data$EFFECT == "synonymous_variant c.573C>T p.Ile191Ile")
my_data158 = subset(my_data, my_data$EFFECT == "synonymous_variant c.573T>C p.Asp191Asp")
my_data159 = subset(my_data, my_data$EFFECT == "synonymous_variant c.579C>T p.Asp193Asp")
my_data160 = subset(my_data, my_data$EFFECT == "synonymous_variant c.579T>C p.Phe193Phe")
my_data161 = subset(my_data, my_data$EFFECT == "synonymous_variant c.57C>T p.Gly19Gly")
my_data162 = subset(my_data, my_data$EFFECT == "synonymous_variant c.585C>T p.Cys195Cys")
my_data163 = subset(my_data, my_data$EFFECT == "synonymous_variant c.588G>A p.Leu196Leu")
my_data164 = subset(my_data, my_data$EFFECT == "synonymous_variant c.597C>T p.Ala199Ala")
my_data165 = subset(my_data, my_data$EFFECT == "synonymous_variant c.600T>C p.Asp200Asp")
my_data166 = subset(my_data, my_data$EFFECT == "synonymous_variant c.60C>G p.Ala20Ala")
my_data167 = subset(my_data, my_data$EFFECT == "synonymous_variant c.60C>T p.Cys20Cys")
my_data168 = subset(my_data, my_data$EFFECT == "synonymous_variant c.612G>T p.Pro204Pro")
my_data169 = subset(my_data, my_data$EFFECT == "synonymous_variant c.618C>T p.Gly206Gly")
my_data170 = subset(my_data, my_data$EFFECT == "synonymous_variant c.621C>T p.Ala207Ala")
my_data171 = subset(my_data, my_data$EFFECT == "synonymous_variant c.627C>T p.His209His")
my_data172 = subset(my_data, my_data$EFFECT == "synonymous_variant c.645A>G p.Ser215Ser")
my_data173 = subset(my_data, my_data$EFFECT == "synonymous_variant c.645G>C p.Ala215Ala")
my_data174 = subset(my_data, my_data$EFFECT == "synonymous_variant c.645T>C p.Leu215Leu")
my_data175 = subset(my_data, my_data$EFFECT == "synonymous_variant c.663G>T p.Gly221Gly")
my_data176 = subset(my_data, my_data$EFFECT == "synonymous_variant c.666C>T p.Asn222Asn")
my_data177 = subset(my_data, my_data$EFFECT == "synonymous_variant c.669G>A p.Val223Val")
my_data178 = subset(my_data, my_data$EFFECT == "synonymous_variant c.66A>G p.Ala22Ala")
my_data179 = subset(my_data, my_data$EFFECT == "synonymous_variant c.672G>A p.Ser224Ser")
my_data180 = subset(my_data, my_data$EFFECT == "synonymous_variant c.675G>A p.Ser225Ser")
my_data181 = subset(my_data, my_data$EFFECT == "synonymous_variant c.678T>C p.Tyr226Tyr")
my_data182 = subset(my_data, my_data$EFFECT == "synonymous_variant c.690C>G p.Ala230Ala")
my_data183 = subset(my_data, my_data$EFFECT == "synonymous_variant c.69T>C p.Gly23Gly")
my_data184 = subset(my_data, my_data$EFFECT == "synonymous_variant c.702C>A p.Val234Val")
my_data185 = subset(my_data, my_data$EFFECT == "synonymous_variant c.705T>C p.Phe235Phe")
my_data186 = subset(my_data, my_data$EFFECT == "synonymous_variant c.720A>G p.Glu240Glu")
my_data187 = subset(my_data, my_data$EFFECT == "synonymous_variant c.723T>C p.Arg241Arg")
my_data188 = subset(my_data, my_data$EFFECT == "synonymous_variant c.72C>G p.Pro24Pro")
my_data189 = subset(my_data, my_data$EFFECT == "synonymous_variant c.72G>C p.Gly24Gly")
my_data190 = subset(my_data, my_data$EFFECT == "synonymous_variant c.732A>G p.Glu244Glu")
my_data191 = subset(my_data, my_data$EFFECT == "synonymous_variant c.732G>A p.Ala244Ala")
my_data192 = subset(my_data, my_data$EFFECT == "synonymous_variant c.735A>G p.Glu245Glu")
my_data193 = subset(my_data, my_data$EFFECT == "synonymous_variant c.735C>T p.Thr245Thr")
my_data194 = subset(my_data, my_data$EFFECT == "synonymous_variant c.738C>A p.Arg246Arg")
my_data195 = subset(my_data, my_data$EFFECT == "synonymous_variant c.738C>T p.Ile246Ile")
my_data196 = subset(my_data, my_data$EFFECT == "synonymous_variant c.741C>T p.Gly247Gly")
my_data197 = subset(my_data, my_data$EFFECT == "synonymous_variant c.741G>A p.Ala247Ala")
my_data198 = subset(my_data, my_data$EFFECT == "synonymous_variant c.744G>A p.Leu248Leu")
my_data199 = subset(my_data, my_data$EFFECT == "synonymous_variant c.747C>T p.Ala249Ala")
my_data200 = subset(my_data, my_data$EFFECT == "synonymous_variant c.747C>T p.Thr249Thr")
my_data201 = subset(my_data, my_data$EFFECT == "synonymous_variant c.756T>A p.Thr252Thr")
my_data202 = subset(my_data, my_data$EFFECT == "synonymous_variant c.762G>A p.Leu254Leu")
my_data203 = subset(my_data, my_data$EFFECT == "synonymous_variant c.765G>A p.Thr255Thr")
my_data204 = subset(my_data, my_data$EFFECT == "synonymous_variant c.774G>A p.Thr258Thr")
my_data205 = subset(my_data, my_data$EFFECT == "synonymous_variant c.786G>A p.Pro262Pro")
my_data206 = subset(my_data, my_data$EFFECT == "synonymous_variant c.801T>C p.Leu267Leu")
my_data207 = subset(my_data, my_data$EFFECT == "synonymous_variant c.816G>A p.Arg272Arg")
my_data208 = subset(my_data, my_data$EFFECT == "synonymous_variant c.817C>A p.Arg273Arg")
my_data209 = subset(my_data, my_data$EFFECT == "synonymous_variant c.81G>T p.Arg27Arg")
my_data210 = subset(my_data, my_data$EFFECT == "synonymous_variant c.81G>T p.Leu27Leu")
my_data211 = subset(my_data, my_data$EFFECT == "synonymous_variant c.81T>C p.Cys27Cys")
my_data212 = subset(my_data, my_data$EFFECT == "synonymous_variant c.81T>C p.Ser27Ser")
my_data213 = subset(my_data, my_data$EFFECT == "synonymous_variant c.822G>A p.Leu274Leu")
my_data214 = subset(my_data, my_data$EFFECT == "synonymous_variant c.831C>T p.Ile277Ile")
my_data215 = subset(my_data, my_data$EFFECT == "synonymous_variant c.834C>T p.Ile278Ile")
my_data216 = subset(my_data, my_data$EFFECT == "synonymous_variant c.8439A>T p.Gly2813Gly")
my_data217 = subset(my_data, my_data$EFFECT == "synonymous_variant c.843G>A p.Leu281Leu")
my_data218 = subset(my_data, my_data$EFFECT == "synonymous_variant c.843G>T p.Leu281Leu")
my_data219 = subset(my_data, my_data$EFFECT == "synonymous_variant c.849C>T p.Ile283Ile")
my_data220 = subset(my_data, my_data$EFFECT == "synonymous_variant c.84C>T p.Arg28Arg")
my_data221 = subset(my_data, my_data$EFFECT == "synonymous_variant c.852G>A p.Ala284Ala")
my_data222 = subset(my_data, my_data$EFFECT == "synonymous_variant c.852G>A p.Leu284Leu")
my_data223 = subset(my_data, my_data$EFFECT == "synonymous_variant c.855G>A p.Val285Val")
my_data224 = subset(my_data, my_data$EFFECT == "synonymous_variant c.861G>T p.Ala287Ala")
my_data225 = subset(my_data, my_data$EFFECT == "synonymous_variant c.870G>A p.Ala290Ala")
my_data226 = subset(my_data, my_data$EFFECT == "synonymous_variant c.873C>G p.Ala291Ala")
my_data227 = subset(my_data, my_data$EFFECT == "synonymous_variant c.906C>G p.Thr302Thr")
my_data228 = subset(my_data, my_data$EFFECT == "synonymous_variant c.90C>T p.Gly30Gly")
my_data229 = subset(my_data, my_data$EFFECT == "synonymous_variant c.918G>A p.Leu306Leu")
my_data230 = subset(my_data, my_data$EFFECT == "synonymous_variant c.924C>T p.Ala308Ala")
my_data231 = subset(my_data, my_data$EFFECT == "synonymous_variant c.933C>T p.Gly311Gly")
my_data232 = subset(my_data, my_data$EFFECT == "synonymous_variant c.936C>G p.Ser312Ser")
my_data233 = subset(my_data, my_data$EFFECT == "synonymous_variant c.951C>G p.Leu317Leu")
my_data234 = subset(my_data, my_data$EFFECT == "synonymous_variant c.951G>C p.Pro317Pro")
my_data235 = subset(my_data, my_data$EFFECT == "synonymous_variant c.963G>A p.Pro321Pro")
my_data236 = subset(my_data, my_data$EFFECT == "synonymous_variant c.969T>C p.Ile323Ile")
my_data237 = subset(my_data, my_data$EFFECT == "synonymous_variant c.972G>C p.Ala324Ala")
my_data238 = subset(my_data, my_data$EFFECT == "synonymous_variant c.975C>T p.Gly325Gly")
my_data239 = subset(my_data, my_data$EFFECT == "synonymous_variant c.993C>T p.Thr331Thr")
my_data240 = subset(my_data, my_data$EFFECT == "synonymous_variant c.996C>A p.Gly332Gly")
my_data241 = subset(my_data, my_data$EFFECT == "synonymous_variant c.996C>T p.Asp332Asp")
my_data242 = subset(my_data, my_data$EFFECT == "synonymous_variant c.99C>T p.Ile33Ile")
x_tb <- rbind(my_data1
,my_data2
,my_data3
,my_data4
,my_data5
,my_data6
,my_data7
,my_data8
,my_data9
,my_data10
,my_data11
,my_data12
,my_data13
,my_data14
,my_data15
,my_data16
,my_data17
,my_data18
,my_data19
,my_data20
,my_data21
,my_data22
,my_data23
,my_data24
,my_data25
,my_data26
,my_data27
,my_data28
,my_data29
,my_data30
,my_data31
,my_data32
,my_data33
,my_data34
,my_data35
,my_data36
,my_data37
,my_data38
,my_data39
,my_data40
,my_data41
,my_data42
,my_data43
,my_data44
,my_data45
,my_data46
,my_data47
,my_data48
,my_data49
,my_data50
,my_data51
,my_data52
,my_data53
,my_data54
,my_data55
,my_data56
,my_data57
,my_data58
,my_data59
,my_data60
,my_data61
,my_data62
,my_data63
,my_data64
,my_data65
,my_data66
,my_data67
,my_data68
,my_data69
,my_data70
,my_data71
,my_data72
,my_data73
,my_data74
,my_data75
,my_data76
,my_data77
,my_data78
,my_data79
,my_data80
,my_data81
,my_data82
,my_data83
,my_data84
,my_data85
,my_data86
,my_data87
,my_data88
,my_data89
,my_data90
,my_data91
,my_data92
,my_data93
,my_data94
,my_data95
,my_data96
,my_data97
,my_data98
,my_data99
,my_data100
,my_data101
,my_data102
,my_data103
,my_data104
,my_data105
,my_data106
,my_data107
,my_data108
,my_data109
,my_data110
,my_data111
,my_data112
,my_data113
,my_data114
,my_data115
,my_data116
,my_data117
,my_data118
,my_data119
,my_data120
,my_data121
,my_data122
,my_data123
,my_data124
,my_data125
,my_data126
,my_data127
,my_data128
,my_data129
,my_data130
,my_data131
,my_data132
,my_data133
,my_data134
,my_data135
,my_data136
,my_data137
,my_data138
,my_data139
,my_data140
,my_data141
,my_data142
,my_data143
,my_data144
,my_data145
,my_data146
,my_data147
,my_data148
,my_data149
,my_data150
,my_data151
,my_data152
,my_data153
,my_data154
,my_data155
,my_data156
,my_data157
,my_data158
,my_data159
,my_data160
,my_data161
,my_data162
,my_data163
,my_data164
,my_data165
,my_data166
,my_data167
,my_data168
,my_data169
,my_data170
,my_data171
,my_data172
,my_data173
,my_data174
,my_data175
,my_data176
,my_data177
,my_data178
,my_data179
,my_data180
,my_data181
,my_data182
,my_data183
,my_data184
,my_data185
,my_data186
,my_data187
,my_data188
,my_data189
,my_data190
,my_data191
,my_data192
,my_data193
,my_data194
,my_data195
,my_data196
,my_data197
,my_data198
,my_data199
,my_data200
,my_data201
,my_data202
,my_data203
,my_data204
,my_data205
,my_data206
,my_data207
,my_data208
,my_data209
,my_data210
,my_data211
,my_data212
,my_data213
,my_data214
,my_data215
,my_data216
,my_data217
,my_data218
,my_data219
,my_data220
,my_data221
,my_data222
,my_data223
,my_data224
,my_data225
,my_data226
,my_data227
,my_data228
,my_data229
,my_data230
,my_data231
,my_data232
,my_data233
,my_data234
,my_data235
,my_data236
,my_data237
,my_data238
,my_data239
,my_data240
,my_data241
,my_data242)
write.csv(x_tb,'x_tb.csv')