-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathGenescan
340 lines (139 loc) · 8.51 KB
/
Genescan
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
GENSCAN 1.0 Date run: 19-May-118 Time: 06:40:04
Sequence /tmp/05_19_18-06:40:03.fasta : 72278 bp : 31.69% C+G : Isochore 1 ( 0 - 43 C+G%)
Parameter matrix: Arabidopsis.smat
Predicted genes/exons:
Gn.Ex Type S .Begin ...End .Len Fr Ph I/Ac Do/T CodRg P.... Tscr..
----- ---- - ------ ------ ---- -- -- ---- ---- ----- ----- ------
1.01 Intr + 77 267 191 2 2 0 18 254 0.963 12.91
1.02 Term + 345 742 398 1 2 10 45 304 0.810 17.55
1.03 PlyA + 1175 1180 6 -4.04
2.03 PlyA - 1274 1269 6 1.05
2.02 Term - 1492 1344 149 2 2 16 44 180 0.250 8.48
2.01 Init - 2174 2096 79 2 1 54 5 109 0.194 5.47
2.00 Prom - 2899 2860 40 -8.15
3.04 PlyA - 3073 3068 6 1.05
3.03 Term - 6698 6559 140 0 2 38 48 140 0.167 7.14
3.02 Intr - 11864 11800 65 1 2 48 58 57 0.234 1.14
3.01 Init - 12596 12151 446 2 2 49 70 365 0.437 31.33
3.00 Prom - 12958 12919 40 -3.95
4.08 PlyA - 14494 14489 6 1.05
4.07 Term - 17516 16323 1194 2 0 38 41 667 0.863 52.72
4.06 Intr - 18766 18668 99 1 0 15 31 141 0.947 5.49
4.05 Intr - 20516 20011 506 0 2 73 24 276 0.887 16.47
4.04 Intr - 20933 20571 363 0 0 32 21 163 0.404 3.23
4.03 Intr - 21577 21423 155 0 2 48 42 137 0.578 8.89
4.02 Intr - 21883 21701 183 0 0 47 19 232 0.541 15.18
4.01 Init - 24807 24734 74 0 2 55 115 20 0.540 6.99
4.00 Prom - 25501 25462 40 -10.15
5.04 PlyA - 25796 25791 6 -3.94
5.03 Term - 26312 25798 515 0 2 69 41 238 0.600 15.72
5.02 Intr - 26588 26499 90 0 0 -24 98 104 0.439 4.25
5.01 Init - 27789 26955 835 0 1 68 76 295 0.551 26.20
5.00 Prom - 28638 28599 40 -12.33
6.00 Prom + 28948 28987 40 -10.65
6.01 Init + 29212 29325 114 0 0 18 76 65 0.278 3.56
6.02 Intr + 30112 30304 193 0 1 8 95 194 0.942 15.24
6.03 Intr + 30352 30467 116 2 2 -13 30 119 0.549 0.25
6.04 Intr + 30752 30916 165 1 0 46 95 87 0.872 9.44
6.05 Intr + 31056 31135 80 2 2 64 21 61 0.811 -0.57
6.06 Intr + 31239 31364 126 0 0 18 61 105 0.324 4.67
6.07 Term + 41795 41951 157 2 1 41 42 185 0.163 10.62
6.08 PlyA + 42095 42100 6 1.05
7.05 PlyA - 42123 42118 6 -1.75
7.04 Term - 42975 42891 85 2 1 72 39 89 0.069 3.55
7.03 Intr - 49020 48731 290 0 2 64 103 53 0.696 4.52
7.02 Intr - 49238 49102 137 0 2 65 94 -3 0.555 2.37
7.01 Init - 50289 50283 7 0 1 57 113 0 0.837 5.64
7.00 Prom - 56523 56484 40 -5.45
8.03 PlyA - 56609 56604 6 1.05
8.02 Term - 57244 57070 175 0 1 17 49 218 0.884 12.05
8.01 Init - 59245 59238 8 1 2 69 97 0 0.936 4.38
8.00 Prom - 61774 61735 40 -6.55
9.00 Prom + 64275 64314 40 -6.75
9.01 Init + 64827 64829 3 2 0 71 101 0 0.640 4.65
9.02 Intr + 65253 65306 54 2 0 107 78 11 0.662 5.26
9.03 Intr + 66159 66261 103 2 1 73 53 126 0.727 11.43
9.04 Intr + 66369 66714 346 1 1 119 56 241 0.646 22.53
9.05 Intr + 66790 67155 366 1 0 70 106 225 0.858 20.84
9.06 Intr + 67400 67478 79 2 1 6 68 47 0.982 -1.87
9.07 Term + 67554 67766 213 2 0 55 49 150 0.982 9.05
9.08 PlyA + 68276 68281 6 1.05
10.00 Prom + 68329 68368 40 -7.95
10.01 Init + 69119 69229 111 1 0 57 80 51 0.345 6.46
10.02 Intr + 69914 70081 168 1 0 -2 65 96 0.589 2.52
10.03 Intr + 70167 70269 103 2 1 84 36 146 0.997 12.83
10.04 Intr + 70364 70709 346 0 1 36 57 215 0.584 11.73
10.05 Intr + 70900 71247 348 1 0 29 84 381 0.579 30.35
10.06 Intr + 71329 71407 79 1 1 59 93 53 0.797 6.53
Suboptimal exons with probability > 1.000
Exnum Type S .Begin ...End .Len Fr Ph B/Ac Do/T CodRg P.... Tscr..
----- ---- - ------ ------ ---- -- -- ---- ---- ----- ----- ------
NO EXONS FOUND AT GIVEN PROBABILITY CUTOFF
Predicted peptide sequence(s):
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_1|196_aa
XFGVPESIVTDNGANLNSQLMREICEQFKITHRNSTAYRPQMNGTVEAANKNIKKILRKI
IDNRQTVIPAEVKIPSLRIIQETELSDIEWVRKRIEQLTLIDEKRMVVVYHGQLYRQRMI
RAFHKKVRVRIFEIGQLVLKCIFLHQAEYKGKFAPNWQGPYMVRKVLSRGALVLSEMDGT
EWPKPMNSHVVKRYYV
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_2|75_aa
MHESLDEVIPVPNEDGASPDNICDLLGKYNVWKRSLPSVMSRQTLANIFEIDASQSIRLF
VGKYSKLFQTADEPY
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_3|216_aa
MDFGKNIIEKGKTIIVNSFTNLSTSKSKKNKKIGSTSKSKTNKTSQLRINTDDYTDVNET
VFNIDSNSGLDPYHELQRRFGNFDEELPEDDEHDDVNDYIDSFNNNFDEYDDDEIEPVTA
TTPTPSPSSPAPAPFRCLAPVPPRILGLSIQRSYSREKDLEELAKMIVVIDWKAAFGRDQ
KFSDAKIETDKPKEKPTCLLIGVWFDCGEVALQLAS
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_4|857_aa
MDCDIVDVYDFLYLLINVEDGFGGCCNRSEQIQCEEISDQVVASPDHSDKYDARELSTDD
NIQATPISEIQEVVVQEDVAPTEHRSVATSSLKEHATYSEAVKDYRWVEAMQAEVQDLEN
NKTWVITDLLHGKKPIGCIEFTRSSKGILMNQRKYALELITEVGMSATKPAGTPMDVNVK
LTSRQYDEQVEDRKVSEDPLIDQAAYKKLIGKLLYLNMTRPNISFSTQTLSQFLQQPKKS
HMEAALMVVKYLEREPGQDWASCPLTRRSVTGYMVKIGDSLVSWKAKKQNTVSRSSAEAE
YMSLASTVSELVWLLGMLKEVGVEVQLPVQVYSDSKTAIQIATNPVYHERTKHIEIDCQF
IRKRIQQGLIKVDYIPTQEQPADVLRKGLSRLQHEYLLSKLGVLNIFVPRSLKGSNKIFP
VDKLVKEIVGVDVRFFGSVLSNFDSILRFSVLEKRAPYTKNQFLHEKNRDPLGRRRDFRS
WVIPGFDFGSNFESTQSVLEAVGVLTAIIVVHESGHFLAAYLQGIHVSKFAVGFGPILAK
FNAKNVEYSLRAFPLGGFVGFPDNDPDSDIPPDDKNLLKNRPIFDRVIVISAGVIANIIF
AYVIIFTQVLLVGLPVQESFPGVLVPDVRPFSAASRDGLLPGDVILGVNGIDLGKNGPSL
VTEVVDVIKKSPKRNVLLKIGRGGGSVDVRVTPDENSDGTGRIGVQLSPNFKISKVQPKN
ILEAFSFSGREFWGLTYNVLDSLKQTFMNFSQTASKVSGPVAIIAVGAEVAKSNVDGLYQ
FAAVLNINLAVINLLPLPALDGGTLALILVEAARGGKKLPLEVEQGIMSSGIMFVIIVGL
FLLVRDTLNLDFIRDLL
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_5|479_aa
MLTDRAADDLSQNFNNGKNVGVSSPLQASSSIGQGIDYNHPLFLNPTDITGISIISFQFL
GAENYTLWSRSITLALLGRNKIGLVNGSCSKEVYNEELWGQWERVNAIVLSWLMNSVSKN
LLIGIVFASTAAQVWSDLKERFDRVDGYRTYSLHKEIVSLQQGTNSVSMYYTKLKSLWDE
AEVLVHAPCCKCEKSRGFVVHLNRGKLYQFLMGLNETYHQARGQILMMDPLPTINHAYAM
IVGDESQKVVVSYIGSMGLNSVSMDSVAMYSKTGSSSGDQYDNILKGYQQKTNPTATDCS
ATHAAYTADSSVISNVFYLPEFKHNLMSVSKVTKELGCSVTFFPNCCVFQELYSGKVKEV
GEKSGGLYTHKTVTSETTTALAVTQVFSDMEVWHQRLGHVSSAVLARIFDMNKESQCKVA
KCLVCPYAKQTRLSFPSSRIKSSTCFDLIHVDLCGPYNTSTFDGNKYFLTIVDDLVVFT
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_6|316_aa
MPLYLFLKFLYFIFRSIFIENVAGKVVVITGASSGIGEHLAYEYAKKGARLVLAARRRKS
LEQVADMAYWLGSPHVISVHADVSKVEDCQRLIEETIRNFGRLKMIIANFRVVDHLVSNA
AVTPLYMFEDLVEVTNAAPAMASKAALISFFETLRVELGTQIGITIVTPGLIESEMTKGK
FLTTEGKLEVDQVMRDVEMSVTPILPVEKCARSIVKSACRGDKETSNRDPQQDTFGRDRI
TEIRVSRVSPISTYKSGLIEHQTQFQLPTKLPVKEIVLHPDTKTLIPRNRLKDFAVSTGR
GRCGQEVTIEELVLSF
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_7|172_aa
MNDDWTNDSRFQIINSKSSSNLNSMYHAKLIGQPWLLLTLLTRSLCQLRKLEKIEESEKK
TLMKSRNKNGKQSPLIDIMNDSPIVGLAMGEFRDPIFRNIQKRITCQAKYLVTPESCEDL
LRGQIKILLQKVEKDILFCHFDVNIDCEMFVAVFVEFLNDEINIPSYDFQAD
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_8|60_aa
MTNNSGEIQKSLKRGREGSWTLQVRRTTTFLRICRIDATGRWISGILGVKKRRDGVGLVS
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_9|387_aa
MTFDADIIHKFWRCLDDYFDLKAILDGGNDLPFPDGLLVNGRGSNGLTFTVDQGRTYRFR
ISNVGLTTAINFRIQGHKMVLVEAEGTHTLQNTYESLDIHLGQSYSVLVTMDQPGQDYYI
VASTRFTSPVLTATSILHYSNSAGGVSGPPPGGPTIEIDWSLNQARSVRQNLTASGPRPN
PQGSYHYGLVNTTRTIRLANSAPMINGKKRYAVNSVSFIPADTPLKLADYFKIPGVFNLG
SIQDYPTGGGGYLQTSVLAADFRAYVEIVFENPEDTVQSWHIDGHIFFVVGMDGGQWSAA
SRLNYNIRDGISRCTIQVYPRSWTALYMPLDNVGMWNIRSENWARQYLGQQFYLRVYSPV
NSWRDEYPIPIGALLCGRASGRKTRPL
>/tmp/05_19_18-06:40:03.fasta|GENSCAN_predicted_peptide_10|385_aa
MVMMFLHVNFIHGEDPYRFYTWNVTYGDIYPLGVKQQDQIGSFYYFPSLAFHKAAGGFGS
INIASRSVIPVPFPPPAGEFSILTGDWFKQNHSDLQAILDGGHDLPFPDGLLINGRGSNG
YTFTVDQGKIYRFRISNVGLTTSVNFRIQGHKMMVVEVEGTHTVQNTYDSLDIHLGQSYS
VLLTADQPAQDYYIAVSTRFTSQVLTATSTLRYSNSVGSVTGPPPGGPTIEIDWSFNQAR
SLSGPRPNPQGSYHYGLINTTRTIRLANSAPIINGKQRYAINSVSFVPADTPLKLADHFN
IPGVFTLGSIPDSPTGSGAYLQTSVMAADFRAYTEVIFENLEDSVQSYHIDGHHFFVVGM
GRGEWTPASRLTYNLRDTISRSTVQ