myTestRepository/codon5.m at master · dschslava/myTestRepository · GitHub

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
function [NoC, lastcodon, nuc] = codon5(dna, x, NoC)


% TCT = 0;
% TCG = 0;
% TCC = 0;
% TCA = 0;
% CCC = 0;
% CCA = 0;
% CCT = 0;
% ACC = 0;
% GCC = 0;
% GCT = 0;

% AAA = 0;
% AAC = 0;
% AAG = 0;
% AAT = 0;
% ACA = 0;
% ACC = 0;
% ACG = 0;
% ACT = 0;
% AGA = 0;
% AGC = 0;
% AGG = 0;
% AGT = 0;
% AGA = 0;
% AGC = 0;
% AGG = 0;
% AGT = 0;
% ATA = 0;
% ATC = 0;
% ATG = 0;
% ATT = 0;
%
% CAA = 0;
% CAC = 0;
% CAG = 0;
% CAT = 0;
% CCA = 0;
% CCC = 0;
% CCG = 0;
% CCT = 0;
% CGA = 0;
% CGC = 0;
% CGG = 0;
% CGT = 0;
% CGA = 0;
% CGC = 0;
% CGG = 0;
% CGT = 0;
% CTA = 0;
% CTC = 0;
% CTG = 0;
% CTT = 0;
%
% GAA = 0;
% GAC = 0;
% GAG = 0;
% GAT = 0;
% GCA = 0;
% GCC = 0;
% GCG = 0;
% GCT = 0;
% GGA = 0;
% GGC = 0;
% GGG = 0;
% GGT = 0;
% GGA = 0;
% GGC = 0;
% GGG = 0;
% GGT = 0;
% GTA = 0;
% GTC = 0;
% GTG = 0;
% GTT = 0;
%
% TAA = 0;
% TAC = 0;
% TAG = 0;
% TAT = 0;
% TCA = 0;
% TCC = 0;
% TCG = 0;
% TCT = 0;
% TGA = 0;
% TGC = 0;
% TGG = 0;
% TGT = 0;
% TGA = 0;
% TGC = 0;
% TGG = 0;
% TGT = 0;
% TTA = 0;
% TTC = 0;
% TTG = 0;
% TTT = 0;
% nucleotide
nuc = {'AAA', 'AAC', 'AAG', 'AAT', 'ACA', 'ACC', 'ACG', 'ACT', 'AGA',...
    'AGC', 'AGG', 'AGT', 'ATA', 'ATC',...
    'ATG', 'ATT', 'CAA', 'CAC', 'CAG', 'CAT', 'CCA', 'CCC', 'CCG',...
    'CCT', 'CGA', 'CGC', 'CGG', 'CGT',...
    'CTA', 'CTC', 'CTG', 'CTT', 'GAA', 'GAC', 'GAG', 'GAT', 'GCA',...
    'GCC', 'GCG', 'GCT', 'GGA', 'GGC',...
    'GGG', 'GGT', 'GTA', 'GTC', 'GTG', 'GTT', 'TAA', 'TAC', 'TAG',...
    'TAT', 'TCA', 'TCC', 'TCG', 'TCT',...
    'TGA', 'TGC', 'TGG', 'TGT', 'TTA', 'TTC', 'TTG', 'TTT'};

place = find(strcmp([dna(x) dna(x+1) dna(x+2)], nuc));
lastcodon = [dna(x) dna(x+1) dna(x+2)];
NoC(place) = NoC(place) + 1;
end