{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
47{0.0984401, -0.109199, -0.0350913, 0.0339016, 0.321584,
480.216223, 0.216223, 0.159065, 0.544727, 0.0339016,
491.176, -0.882389, 1.70788, 0.216223, -0.882389},
50{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
51{-0.405465, 0.287682, 0.693147, 0.441833, 0.693147,
521.13498, 0.287682, 0.575364, 0.105361, 0.287682,
531.44036, -0.405465, 0.693147, 0.575364, -0.405465},
54{-2.18605, -2.18605, -2.18605, -2.18605, -2.18605,
55-2.18605, -2.18605, -2.18605, -0.576613, -1.08744,
56-25.2119, -0.576613, -25.2119, 0.116534, 0.211844},
57{-2.30259, -2.30259, -2.30259, -2.30259, -2.30259,
58-2.30259, -2.30259, -2.30259, -1.20397, -0.356675,
59-25.3284, -0.356675, -25.3284, 0.262364, 0.336472},
60{0.842183, 0.474458, 0.675129, 0.675129, 0.0689929,
610.223144, 1.16761, 0.842183, -0.336472, -0.113329,
62-24.7486, 0.842183, -24.7486, 0.0689929, -0.336472},
63{-1.10691, -1.10691, -1.39459, -0.701446, -1.39459,
640.0717439, -1.39459, -1.80006, 0.451234, 1.03316,
65-0.883768, -0.547295, 1.17036, -0.19062, -0.547295},
66{-1.22378, -1.22378, -1.22378, -1.22378, -1.22378,
67-1.22378, -1.22378, -1.22378, 0.385662, -1.22378,
68-24.2496, 0.567984, -24.2496, 0.162519, -0.530628},
69{0.70657, 0.70657, 0.0779615, -0.209721, 0.396415,
70-0.392042, -0.615186, 0.0779615, -0.392042, -2.00148,
710.301105, -0.392042, -25.0273, 0.0779615, -0.0555699},
72{-2.4248, -2.4248, -2.4248, -2.4248, -2.4248,
73-2.4248, -2.4248, -2.4248, -2.4248, -1.03851,
74-25.4507, -1.73166, -25.4507, -0.0269075, -0.227578},
75{1.76947, 1.7263, 1.78346, 1.87624, 1.8635,
761.31346, 1.66568, 1.39861, -0.19062, 0.642289,
77-0.413764, 0.502527, -2.49321, -0.413764, -1.10691},
78{-0.993252, 0.105361, 0.952658, 0.393043, -0.993252,
790.798508, -0.300105, -0.300105, -0.993252, -0.993252,
80-24.0191, -0.993252, -24.0191, -0.993252, -0.300105},
81{-1.96009, -1.96009, -1.96009, -1.96009, -1.96009,
82-1.96009, -1.96009, -1.96009, -0.861482, -0.861482,
83-24.9859, 0.34249, -24.9859, -0.5738, -0.0141846},
84{-1.30833, -2.00148, -1.30833, -2.00148, -2.00148,
85-0.615186, -2.00148, 0.0779615, 0.994252, 0.637577,
86-25.0273, -2.00148, -0.902868, -2.00148, 1.08956},
87{-1.84055, -1.84055, -1.84055, -1.84055, -1.84055,
88-0.0487902, -1.84055, -1.84055, 0.462035, 0.238892,
89-24.8664, 1.04982, -0.741937, 1.10389, 0.462035},
90{-1.335, -2.02815, -2.02815, -2.02815, -2.02815,
91-2.02815, -2.02815, -2.02815, -0.0822381, -0.641854,
92-25.054, 0.679902, -25.054, 0.456758, 0.169076},
93{-0.236389, -1.335, -0.354172, -0.641854, 0.131336,
94-0.131028, 0.274437, 0.338975, 0.824483, -0.0357181,
950.701881, 0.3996, 0.562119, 0.274437, -0.131028},
96{-1.57898, 0.0304592, -0.662688, -0.885832, -0.662688,
970.292823, -0.326216, -0.326216, 0.212781, -0.480366,
980.561087, -0.192684, -0.480366, -1.17351, 0.0304592},
99{0.588787, 0.811931, 0.301105, 0.483427, 0.158004,
1000.301105, -0.00904984, 0.888892, -2.40695, 0.0779615,
1011.05879, -1.30833, -25.4328, -0.327504, 0.426268},
102{0.798508, 0.510826, 0.510826, -0.587787, -0.587787,
1030.105361, 1.20397, 0.510826, -0.587787, 0.510826,
104-23.6136, 1.60944, -23.6136, 0.105361, -0.587787},
105{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
106{-1.72277, -1.72277, -0.336472, -1.72277, -1.72277,
107-1.72277, -0.624154, -1.72277, -1.72277, -1.02962,
108-24.7486, -0.113329, -24.7486, -1.72277, 0.223144},
109{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
110{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
111{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
117{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
118{1.13943, 0.916291, 0.916291, 1.03407, 0.628609,
1190.782759, 0.446287, 0.628609, 0.782759, 0.782759,
1202.0149, -0.470004, 2.27084, 1.72722, 0.223144},
121{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
124-23.0259, 0, -23.0259, 0, 0},
125{-0.693147, -0.693147, -0.693147, -0.693147, -0.693147,
126-0.693147, -0.693147, -0.693147, -0.693147, -0.693147,
127-23.719, -0.693147, -23.719, 0, 1.38629},
128{-0.788457, -0.788457, -0.788457, -0.788457, -0.788457,
129-0.788457, -0.788457, -0.788457, -0.788457, -0.788457,
130-23.8143, -0.788457, -23.8143, 0.597837, 1.29098},
131{0.430783, 1.12393, 0.836248, 1.12393, -0.262364,
132-0.262364, 1.81708, -0.262364, 1.12393, -0.262364,
133-23.2882, 1.68355, -23.2882, -0.262364, -0.262364},
134{0.393043, -0.300105, -0.300105, -0.300105, 0.105361,
1350.616186, -0.300105, 0.393043, -0.300105, -0.300105,
136-24.0191, -0.300105, -0.300105, -0.993252, -0.993252},
137{0.223144, 0.223144, 0.223144, 0.223144, 0.223144,
1380.223144, 0.223144, 0.223144, 0.223144, 0.223144,
139-22.8027, 2.16905, -22.8027, 0.223144, 0.223144},
140{0.567984, -0.530628, 1.07881, -0.530628, 1.07881,
141-0.530628, -0.530628, 0.567984, -0.530628, -0.530628,
142-23.5565, -0.530628, -23.5565, -0.530628, 0.162519},
143{-0.916291, -0.916291, -0.916291, -0.916291, -0.916291,
144-0.916291, -0.916291, -0.916291, -0.916291, -0.916291,
145-23.9421, -0.223144, -23.9421, 0.182322, -0.916291},
146{1.08619, 1.40464, 1.20397, 1.08619, 1.20397,
1471.5717, -0.993252, -0.993252, -0.300105, -0.300105,
148-0.993252, -0.300105, -24.0191, -0.993252, -0.993252},
149{0.510826, 1.20397, 0.510826, 0.510826, 1.60944,
1501.20397, 1.60944, 0.510826, 0.510826, 1.20397,
151-22.515, 1.89712, -22.515, 0.510826, 0.510826},
152{-0.470004, -0.470004, -0.470004, -0.470004, -0.470004,
153-0.470004, -0.470004, -0.470004, -0.470004, -0.470004,
154-23.4959, 0.628609, -23.4959, -0.470004, 0.916291},
155{-0.530628, -0.530628, -0.530628, -0.530628, -0.530628,
156-0.530628, 0.162519, 0.567984, 1.07881, 0.162519,
157-23.5565, -0.530628, -23.5565, -0.530628, 1.07881},
158{-0.336472, -0.336472, -0.336472, -0.336472, -0.336472,
159-0.336472, -0.336472, -0.336472, 0.356675, 0.356675,
160-23.3623, 0.76214, -23.3623, -0.336472, -0.336472},
161{-0.530628, -0.530628, -0.530628, -0.530628, -0.530628,
162-0.530628, -0.530628, -0.530628, -0.530628, -0.530628,
163-23.5565, -0.530628, -23.5565, -0.530628, -0.530628},
164{-0.955511, -0.955511, -0.955511, 0.430783, 0.430783,
165-0.955511, 0.653926, 1.75254, 0.653926, 1.12393,
1660.653926, -0.262364, -0.262364, -0.955511, -0.955511},
167{-0.0953102, -0.788457, 0.597837, -0.0953102, -0.0953102,
168-0.0953102, -0.0953102, -0.0953102, 0.820981, -0.788457,
1690.310155, -0.788457, -0.788457, -0.788457, -0.0953102},
170{0.693147, 1.02962, -0.916291, 0.182322, -0.916291,
1710.470004, 1.02962, -0.916291, -0.916291, 0.470004,
1720.182322, -0.916291, -23.9421, -0.223144, -0.916291},
173{0.916291, 0.916291, 0.916291, 0.916291, 0.916291,
1740.916291, 0.916291, 0.916291, 0.916291, 0.916291,
175-22.1096, 0.916291, -22.1096, 0.916291, 0.916291},
176{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
177{-0.262364, -0.262364, -0.262364, -0.262364, -0.262364,
178-0.262364, -0.262364, -0.262364, -0.262364, 0.836248,
179-23.2882, -0.262364, -23.2882, -0.262364, -0.262364},
180{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
181{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
182{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}
186 template<
classSeq>
191 if(seq.size() < 15) {
192 throwruntime_error(
"Sequence length must be at least 15 " 193 "to predict signal sequence");
196 throwruntime_error(
"max_pos must be at least 12 " 197 "to predict signal sequence");
200 constdouble (*Mat)[15];
209 doublemax_score = -1e6;
211 for(
unsigned int i= 0;
i<= max_index;
i++) {
213 for(
unsigned intj = 0; j < 15; j++) {
214 intsub = (
unsignedchar)seq[
i+ j];
217 if(sum > max_score) {
244 stringseq_ncbistdaa;
247vec.GetSeqData(0, vec.size(), seq_ncbistdaa);
friend void x_PredictSignalSeq(const Seq &seq, CSignalSeq::EDomain domain, TSeqPos max_pos, TSeqPos &pos, double &score)
static void Predict(const string &seq, EDomain domain, TSeqPos max_pos, TSeqPos &pos, double &score)
Find the most likely predicted signal sequence cleavage site (pos) and the associated score (>3....
unsigned int TSeqPos
Type for sequence locations and lengths.
#define END_NCBI_SCOPE
End previously defined NCBI scope.
#define BEGIN_NCBI_SCOPE
Define ncbi namespace.
void x_PredictSignalSeq(const Seq &seq, CSignalSeq::EDomain domain, TSeqPos max_pos, TSeqPos &pos, double &score)
const double const_BacMat[26][15]
const double const_EukMat[26][15]
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4