A test to ensure Bio::PrimarySeqI->trunc() doesn't use clone() for a Bio::Seq::RichSe...
[bioperl-live.git] / t / data / hmmpfam_cs.out
blob28e8a728190cbe84a3657f79594707e2a3f10daa
1 hmmpfam - search one or more sequences against HMM database
2 HMMER 2.3.2 (Oct 2003)
3 Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
4 Freely distributed under the GNU General Public License (GPL)
5 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
6 HMM file:                 ../Shared/Pfam_fs
7 Sequence file:            single_porphyra_AA.fa
8 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
10 Query sequence: gi|90819130|dbj|BAE92499.1|
11 Accession:      [none]
12 Description:    glutamate synthase [Porphyra yezoensis]
14 Scores for sequence family classification (score includes all domains):
15 Model           Description                             Score    E-value  N 
16 --------        -----------                             -----    ------- ---
17 Glu_synthase    Conserved region in glutamate synthas   858.6   3.6e-255   2
18 GATase_2        Glutamine amidotransferases class-II    731.8   3.9e-226   1
19 Glu_syn_central Glutamate synthase central domain       649.1   7.9e-213   1
20 GXGXG           GXGXG motif                             367.3   2.7e-107   1
21 HdeA            hns-dependent expression protein A (H     9.6      0.015   1
22 GDC-P           Glycine cleavage system P-protein         7.1      0.086   1
23 Cache_1         Cache domain                              7.0       0.14   1
24 IBN_N           Importin-beta N-terminal domain           8.2       0.17   1
25 DUF1200         Protein of unknown function (DUF1200)     6.7       0.42   1
26 cobW            CobW/HypB/UreG, nucleotide-binding do     5.1       0.45   1
27 PUF             Pumilio-family RNA binding repeat         6.5       0.47   1
28 Arch_flagellin  Archaebacterial flagellin                 4.1       0.66   1
29 FMN_dh          FMN-dependent dehydrogenase               3.2       0.89   1
30 RNA_pol_Rpb2_4  RNA polymerase Rpb2, domain 4             4.6        1.4   1
31 DUF477          Domain of unknown function (DUF477)       3.8        1.7   1
32 FRG1            FRG1-like family                          0.2        1.7   1
33 DUF1393         Protein of unknown function (DUF1393)     3.1          2   1
34 tRNA_anti       OB-fold nucleic acid binding domain       4.9          2   1
35 SelT            Selenoprotein T                           3.1        2.2   1
36 RNase_PH_C      3' exoribonuclease family, domain 2       4.2        2.3   1
37 Pencillinase_R  Penicillinase repressor                   3.9        2.5   1
38 Hormone_4       Neurohypophysial hormones, N-terminal     4.4        2.5   1
39 DSRB            Dextransucrase DSRB                       2.7        2.7   1
40 FtsK_SpoIIIE    FtsK/SpoIIIE family                       2.6        3.1   1
41 UBA             UBA/TS-N domain                           4.2        3.1   1
42 DUF1981         Domain of unknown function (DUF1981)      3.6        3.3   1
43 Gla             Vitamin K-dependent carboxylation/gam     4.0        3.5   1
44 Scm3            Centromere protein Scm3                   2.2        3.5   1
45 Ribosomal_S6    Ribosomal protein S6                      3.3        3.7   1
46 Cystatin        Cystatin domain                           2.4        3.9   1
47 Phage_prot_Gp6  Phage portal protein, SPP1 Gp6-like       1.0          4   1
48 DUF1976         Domain of unknown function (DUF1976)     -1.5        4.3   1
49 DUF37           Domain of unknown function DUF37          3.0        4.5   1
50 Flavodoxin_NdrI NrdI Flavodoxin like                      2.1        4.6   1
51 Bac_rhodopsin   Bacteriorhodopsin                         0.9        4.9   1
52 Nitro_FeMo-Co   Dinitrogenase iron-molybdenum cofacto     2.1        5.3   1
53 MoCF_biosynth   Probable molybdopterin binding domain     1.3        5.6   1
54 PaaA_PaaC       Phenylacetic acid catabolic protein       0.4        5.6   1
55 Albicidin_res   Albicidin resistance domain               1.7        5.7   1
56 DUF1514         Protein of unknown function (DUF1514)     3.5        5.7   1
57 T5orf172        T5orf172 domain                           2.0        6.1   1
58 Nup133_N        Nup133 N terminal like                   -0.6        6.5   1
59 BicD            Microtubule-associated protein Bicaud    -1.6        6.8   1
60 Sel1            Sel1 repeat                               2.5          7   1
61 CAP_C           DE   Adenylate cyclase associated (CA     1.3        7.4   1
62 Colicin         Colicin pore forming domain               1.4        7.5   1
63 MADF_DNA_bdg    Alcohol dehydrogenase transcription f     1.8        8.2   1
64 DUF258          Protein of unknown function, DUF258       0.3        8.3   1
65 PspB            Phage shock protein B                     0.4        8.4   1
66 GspM            General secretion pathway, M protein      1.0        8.6   1
67 Coq4            Coenzyme Q (ubiquinone) biosynthesis     -0.3        9.1   1
68 P22_AR_N        P22_AR N-terminal domain                 -0.2        9.5   1
69 C1_2            C1 domain                                 1.1        9.6   1
70 Phage_Mu_P      Bacteriophage Mu P protein               -0.4         10   1
72 Parsed for domains:
73 Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
74 --------        ------- ----- -----    ----- -----      -----  -------
75 GATase_2          1/1      34   404 ..     1   385 []   731.8 3.9e-226
76 FRG1              1/1      88   107 ..   151   173 ..     0.2      1.7
77 C1_2              1/1     191   210 ..     9    27 ..     1.1      9.6
78 MADF_DNA_bdg      1/1     235   261 ..    57    95 .]     1.8      8.2
79 PaaA_PaaC         1/1     258   269 ..     1    13 [.     0.4      5.6
80 Albicidin_res     1/1     274   289 ..    50    65 ..     1.7      5.7
81 UBA               1/1     311   331 ..    18    38 .]     4.2      3.1
82 Gla               1/1     342   357 ..    27    42 .]     4.0      3.5
83 RNA_pol_Rpb2_4    1/1     369   381 ..     1    13 [.     4.6      1.4
84 MoCF_biosynth     1/1     371   396 ..    23    49 ..     1.3      5.6
85 DUF1200           1/1     389   401 ..     1    13 [.     6.7     0.42
86 Nup133_N          1/1     397   419 ..   475   498 .]    -0.6      6.5
87 DUF1976           1/1     428   448 ..  1296  1319 .]    -1.5      4.3
88 Bac_rhodopsin     1/1     445   472 ..   219   250 .]     0.9      4.9
89 Coq4              1/1     459   481 ..    60    82 ..    -0.3      9.1
90 Glu_syn_central   1/1     478   773 ..     1   301 []   649.1 7.9e-213
91 Flavodoxin_NdrI   1/1     488   497 ..   122   131 .]     2.1      4.6
92 P22_AR_N          1/1     524   541 ..   110   126 .]    -0.2      9.5
93 Cache_1           1/1     537   557 ..     1    23 [.     7.0     0.14
94 Glu_synthase      1/2     650   676 ..   297   323 ..     1.3        3
95 HdeA              1/1     727   749 ..    58    79 .]     9.6    0.015
96 Sel1              1/1     729   745 ..    32    49 .]     2.5        7
97 DUF1981           1/1     765   787 ..    62    88 .]     3.6      3.3
98 tRNA_anti         1/1     818   839 ..    54    85 .]     4.9        2
99 Cystatin          1/1     826   859 ..     1    38 [.     2.4      3.9
100 RNase_PH_C        1/1     827   846 ..    64    84 .]     4.2      2.3
101 Glu_synthase      2/2     830  1216 ..     1   412 []   857.3   9e-255
102 DUF258            1/1     839   860 ..   282   305 .]     0.3      8.3
103 Pencillinase_R    1/1     856   894 ..    84   118 .]     3.9      2.5
104 SelT              1/1     872   885 ..    96   111 .]     3.1      2.2
105 Nitro_FeMo-Co     1/1     879   897 ..    87   105 .]     2.1      5.3
106 DUF37             1/1     927   934 ..    61    68 .]     3.0      4.5
107 Scm3              1/1     953   963 ..   103   113 .]     2.2      3.5
108 cobW              1/1    1038  1058 ..   202   222 .]     5.1     0.45
109 Arch_flagellin    1/1    1050  1072 ..   197   219 .]     4.1     0.66
110 DUF1393           1/1    1055  1068 ..     1    14 [.     3.1        2
111 FtsK_SpoIIIE      1/1    1107  1143 ..   163   198 ..     2.6      3.1
112 FMN_dh            1/1    1109  1148 ..   291   330 ..     3.2     0.89
113 DSRB              1/1    1120  1134 ..     1    16 [.     2.7      2.7
114 Phage_Mu_P        1/1    1122  1131 ..     1    10 [.    -0.4       10
115 Hormone_4         1/1    1168  1176 ..     1     9 []     4.4      2.5
116 GDC-P             1/1    1205  1225 ..    10    30 ..     7.1    0.086
117 PspB              1/1    1268  1276 ..     1     9 [.     0.4      8.4
118 T5orf172          1/1    1271  1293 ..    35    58 ..     2.0      6.1
119 CAP_C             1/1    1283  1292 ..   161   170 .]     1.3      7.4
120 GXGXG             1/1    1290  1485 ..     1   228 []   367.3 2.7e-107
121 DUF1514           1/1    1453  1469 ..    50    66 .]     3.5      5.7
122 Colicin           1/1    1456  1467 ..   192   203 .]     1.4      7.5
123 Ribosomal_S6      1/1    1461  1481 ..    16    36 ..     3.3      3.7
124 BicD              1/1    1465  1481 ..     1    17 [.    -1.6      6.8
125 PUF               1/1    1470  1486 ..    19    35 .]     6.5     0.47
126 DUF477            1/1    1472  1495 ..     1    24 [.     3.8      1.7
127 Phage_prot_Gp6    1/1    1479  1492 ..     1    14 [.     1.0        4
128 IBN_N             1/1    1498  1516 ..     1    20 [.     8.2     0.17
129 GspM              1/1    1506  1520 ..     1    15 [.     1.0      8.6
131 Alignments of top-scoring domains:
132 GATase_2: domain 1 of 1, from 34 to 404: score 731.8, E = 3.9e-226
133                 CS    EEEEEEEEETSSHSBHHHHHHHHHHHHHGGGGSSCSTTSSCECEEEE
134                    *->CGvlGfiAhikgkpshkivedaleaLerLeHRGavgADgktGDGAGI
135                       CGv GfiA+ ++ ++hkiv +aleaL+++eHRGa++AD ++GDGAGI
136   gi|9081913    34    CGV-GFIADVNNVANHKIVVQALEALTCMEHRGACSADRDSGDGAGI 79   
138                 CS EEECTCCCHHHHHHHCT----S GC-EEEEEEE-SSHHHHHHHHHHHHHH
139                    ltqiPdgFFrevakelGieLpe.gqYAVGmvFLPqdelaraearkifEki
140                     t+iP+++F++  ++++i++ ++   +VGm+FLP   l+    + i+E +
141   gi|9081913    80 TTAIPWNLFQKSLQNQNIKFEQnDSVGVGMLFLPAHKLKES--KLIIETV 127  
143                 CS HHHTT-EEEEEEE--B-GGGS-HHHHHC--EEEEEEEE-TT--HHHHHHC
144                    aeeeGLeVLGWReVPvnnsvLGetAlatePvIeQvFvgapsgdgedfErr
145                    ++ee+Le++GWR VP+  +vLG++A  + P++eQvF+ +++ +++ +E++
146   gi|9081913   128 LKEENLEIIGWRLVPTVQEVLGKQAYLNKPHVEQVFCKSSNLSKDRLEQQ 177  
148                 CS EEEEECHSCHHHHTHHH.    BEEEEEESSEEEEEECC-GGGHHHHBHG
149                    LyviRkrieksivaenvn....fYiCSLSsrTIVYKGMLtseQLgqFYpD
150                    L+++Rk+iek+i+  + +  ++fYiCSLS++TIVYKGM++s++LgqFY+D
151   gi|9081913   178 LFLVRKKIEKYIGINGKDwaheFYICSLSCYTIVYKGMMRSAVLGQFYQD 227  
153                 CS GGSTTEEBSEEEEEECESSSSSCTGGGSSCEEECCCTTCEEEEEEEEETT
154                    LqderfeSalAivHsRFSTNTfPsWplAQPfRVnslwgggivlAHNGEIN
155                    L++++++S++Ai+H+RFSTNT+P+WplAQP+R         ++ HNGEIN
156   gi|9081913   228 LYHSEYTSSFAIYHRRFSTNTMPKWPLAQPMR---------FVSHNGEIN 268  
158                 CS THHHHHHHHHHTSCCCSSTTCGHHHHCC-SSS-TTSCHHHHHHHHHHHHH
159                    TlrgNrnwMraRegvlksplFgddldkLkPIvneggSDSaalDnvlEllv
160                    Tl gN nwM++Re +l+s++++d++++LkPI n+++SDSa+lD ++Ell+
161   gi|9081913   269 TLLGNLNWMQSREPLLQSKVWKDRIHELKPITNKDNSDSANLDAAVELLI 318  
163                 CS HTT--HHHHHHHHS----TT-GGGTST-HHHHHHHHHHHHHHCCHCCEEE
164                    raGRslpeAlMMlIPEAWqnnpdmdkdrpekraFYeylsglmEPWDGPAa
165                    ++GRs++eAlM+l+PEA+qn+pd   +++e+ +FYey+sgl+EPWDGPA+
166   gi|9081913   319 ASGRSPEEALMILVPEAFQNQPDFA-NNTEISDFYEYYSGLQEPWDGPAL 367  
168                 CS EEEETSSEEEEEEETTTSCESEEEEEEEEEE.TTEEEEEESSC   
169                    lvftDGryavgAtLDRNGLTRPaRygiTrdldkDglvvvaSEa<-*
170                    +vft+G++ +gAtLDRNGL RPaRy+iT    kD+lv+v+SE+   
171   gi|9081913   368 VVFTNGKV-IGATLDRNGL-RPARYVIT----KDNLVIVSSES    404  
173 FRG1: domain 1 of 1, from 88 to 107: score 0.2, E = 1.7
174                    *->FQkfKvDLqdrklrinekDkkel<-*
175                       FQk+   Lq+  +  +++D+ ++   
176   gi|9081913    88    FQKS---LQNQNIKFEQNDSVGV    107  
178 C1_2: domain 1 of 1, from 191 to 210: score 1.1, E = 9.6
179                    *->idgfyg...fYsCkkccddftl<-*
180                       i+g+++ ++fY C+  c  +t+   
181   gi|9081913   191    INGKDWaheFYICSLSC--YTI    210  
183 MADF_DNA_bdg: domain 1 of 1, from 235 to 261: score 1.8, E = 8.2
184                    *->drYrrelrkirqgnsegsstgsgesykskWryyeelsFL<-*
185                       +++  ++r+               ++ +kW+++  ++F    
186   gi|9081913   235    SSFAIYHRRFS------------TNTMPKWPLAQPMRFV    261  
188 PaaA_PaaC: domain 1 of 1, from 258 to 269: score 0.4, E = 5.6
189                 CS    X............   
190                    *->MYnFvEHGGvint<-*
191                       M  Fv H G int   
192   gi|9081913   258    M-RFVSHNGEINT    269  
194 Albicidin_res: domain 1 of 1, from 274 to 289: score 1.7, E = 5.7
195                    *->LrlmharEPsLrkgtG<-*
196                       L+ m+ rEP L+ +++   
197   gi|9081913   274    LNWMQSREPLLQSKVW    289  
199 UBA: domain 1 of 1, from 311 to 331: score 4.2, E = 3.1
200                 CS    HHHHHHHHHTTT-HHHHHHHH   
201                    *->eeakkALeatngnverAvewL<-*
202                       ++a++ L a++ ++e+A+++L   
203   gi|9081913   311    DAAVELLIASGRSPEEALMIL    331  
205 Gla: domain 1 of 1, from 342 to 357: score 4.0, E = 3.5
206                 CS    CSSHHHHHHHHHHCTC   
207                    *->fednegtkefwrkYfg<-*
208                       f++n+++  f++ Y g   
209   gi|9081913   342    FANNTEISDFYEYYSG    357  
211 RNA_pol_Rpb2_4: domain 1 of 1, from 369 to 381: score 4.6, E = 1.4
212                 CS    EEETTEEEEEESS   
213                    *->VYvNGklvGthrn<-*
214                       V+ NGk++G + +   
215   gi|9081913   369    VFTNGKVIGATLD    381  
217 MoCF_biosynth: domain 1 of 1, from 371 to 396: score 1.3, E = 5.6
218                 CS    CHHHHHHHHHHHTTTCEEEEEEEE-SS   
219                    *->tNgpmLaalLresaGaevirygiVpDd<-*
220                       tNg+ + a L +  G  ++ry+i +D+   
221   gi|9081913   371    TNGKVIGATLDR-NGLRPARYVITKDN    396  
223 DUF1200: domain 1 of 1, from 389 to 401: score 6.7, E = 0.42
224                    *->kYvltedtLlIks<-*
225                       +Yv+t+d L+I+s   
226   gi|9081913   389    RYVITKDNLVIVS    401  
228 Nup133_N: domain 1 of 1, from 397 to 419: score -0.6, E = 6.5
229                    *->lylltrnsGvvrIeHaleedstne<-*
230                       l++ + +sGvv++e +  + s  +   
231   gi|9081913   397    LVIVSSESGVVQVE-PGNVKSKGR    419  
233 DUF1976: domain 1 of 1, from 428 to 448: score -1.5, E = 4.3
234                    *->VsvYiyFkevtdnksLsEysVtyk<-*
235                       V++++   ++++nk ++  sVt k   
236   gi|9081913   428    VDIFS--HKILNNKEIK-TSVTTK    448  
238 Bac_rhodopsin: domain 1 of 1, from 445 to 472: score 0.9, E = 4.9
239                 CS    HHHHHHHHHHHHHHHHHCHHHTC---------   
240                    *->vvAKVgFgfilLrsravlertvavgsalaage<-*
241                       v++K+++g +l ++r++le  +   + l+++    
242   gi|9081913   445    VTTKIPYGELLTDARQILE--HK--PFLSDQQ    472  
244 Coq4: domain 1 of 1, from 459 to 481: score -0.3, E = 9.1
245                    *->rrILkEkPRissetldlkkLrkL<-*
246                       r+IL  kP  s  ++d kkL +L   
247   gi|9081913   459    RQILEHKPFLSDQQVDIKKLMQL    481  
249 Glu_syn_central: domain 1 of 1, from 478 to 773: score 649.1, E = 7.9e-213
250                 CS    HHHHHHCTT--HHHHHCTCHHHHHHSS--EE-S---S--CCC-SS--
251                    *->llrrQkAFGYTyEdvelvllPMAetGkEalGSMGdDtPLAVLSekpr
252                       l+++Q+AFGYT+Edvelv+++MA+++kE++++MGdD+PL +LSek++
253   gi|9081913   478    LMQLQTAFGYTNEDVELVIEHMASQAKEPTFCMGDDIPLSILSEKSH 524  
255                 CS -GGGCEEE----SSS----TTTTGGG-B--EEES--S-TTS-SGGGC-CE
256                    lLYdYFKQlFAQVTNPPIDPIREelVMSLetylGpegNlLeptpeqarrl
257                    +LYdYFKQ+FAQVTNP+IDP+RE+lVMSL+ ++G+++NlL+  p+ a+++
258   gi|9081913   525 ILYDYFKQRFAQVTNPAIDPLRESLVMSLAIQIGHKSNLLDDQPTLAKHI 574  
260                 CS EESSSB--HHHHHH.HHHH....CCCCEEEEESEEESTTSTTCHHHHHHH
261                    kLesPILsnselekmlknidairegfkaatIditFdveeGvdgLeaaLdr
262                    kLesP+++++el++ + +     +++++  I+++F  e+G++ ++  + +
263   gi|9081913   575 KLESPVINEGELNA-IFE-----SKLSCIRINTLFQLEDGPKNFKQQIQQ 618  
265                 CS HHHHHHHHHHCT-SEEEEESTCG--CTTEEE--HHHHHHHHHHHHHCTT-
266                    lceeAeeAirsGaniivLSDRndildeervaIPaLLAvGAVHhHLIrkgL
267                    lce A++Ai +G ni+vLSD+n+ ld+e+v+IP+LLAvGAVHhHLI kgL
268   gi|9081913   619 LCENASQAILDGNNILVLSDKNNSLDSEKVSIPPLLAVGAVHHHLINKGL 668  
270                 CS CCC-EEEEEESS--SHHHHHHHHCTT-SEEEEHCCHHHHHHHHCCCCCCC
271                    RtkvslvVETGEaREvHHFAvLiGYGAsAInPYLAyETirdWWlirrGll
272                    R+ +s+ VET++++++HHFA+LiGYGAsAI+PYLA+ET r+WW + ++++
273   gi|9081913   669 RQEASILVETAQCWSTHHFACLIGYGASAICPYLAFETARHWWSNPKTKM 718  
275                 CS CHTTTS- T--HHHHHHHHHHHHHHHHHHHHHCTT--BHHHHCCS--EEE
276                    lmskGkl.elsleeavkNYrkAiekGlLKIMSKMGISTlqSYrGAQIFEA
277                    lmskG+l++++++ea++NY+kA+e+GlLKI+SKMGIS+l+SY+GAQIFE+
278   gi|9081913   719 LMSKGRLpACNIQEAQANYKKAVEAGLLKILSKMGISLLSSYHGAQIFEI 768  
280                 CS SSB-H   
281                    vGLsk<-*
282                    +GL++   
283   gi|9081913   769 LGLGS    773  
285 Flavodoxin_NdrI: domain 1 of 1, from 488 to 497: score 2.1, E = 4.6
286                 CS    -HHHHHHHHH   
287                    *->TneDVerVrk<-*
288                       TneDVe V +   
289   gi|9081913   488    TNEDVELVIE    497  
291 P22_AR_N: domain 1 of 1, from 524 to 541: score -0.2, E = 9.5
292                    *->dVLydYWtrkGkAv..NPR<-*
293                       ++LydY+  + +A  +NP+   
294   gi|9081913   524    HILYDYFK-QRFAQvtNPA    541  
296 Cache_1: domain 1 of 1, from 537 to 557: score 7.0, E = 0.14
297                    *->wTePYvdaalktgdlViTiaqPv<-*
298                       +T+P++d +  +++lV ++a+++   
299   gi|9081913   537    VTNPAIDPL--RESLVMSLAIQI    557  
301 Glu_synthase: domain 1 of 2, from 650 to 676: score 1.3, E = 3
302                 CS    --HHHHHHHHHHHHHCTT-CCCSEEEE   
303                    *->lPwelgLaevhqtLvengLRdrVsLia<-*
304                       +P  l++ +vh  L++ gLR + s+ +   
305   gi|9081913   650    IPPLLAVGAVHHHLINKGLRQEASILV    676  
307 HdeA: domain 1 of 1, from 727 to 749: score 9.6, E = 0.015
308                    *->ACk.QdkkAsFkdKvkaEldKvk<-*
309                       AC  Q+ +A++k+ v+a l K+    
310   gi|9081913   727    ACNiQEAQANYKKAVEAGLLKIL    749  
312 Sel1: domain 1 of 1, from 729 to 745: score 2.5, E = 7
313                 CS    .HHH.HHHHHHHHHHTT-   
314                    *->DyekeAlkwyekAAeqGn<-*
315                       ++++ A + y+kA e+G    
316   gi|9081913   729    NIQE-AQANYKKAVEAGL    745  
318 DUF1981: domain 1 of 1, from 765 to 787: score 3.6, E = 3.3
319                    *->iFgvltlaakeesesivklAfqiid.qi<-*
320                       iF++l+l++       v+lAf+ +++qi   
321   gi|9081913   765    IFEILGLGSEV-----VNLAFKGTTsQI    787  
323 tRNA_anti: domain 1 of 1, from 818 to 839: score 4.9, E = 2
324                 CS    EEEEEEETTSSTSTCTCTT..EEEEEEEEEEE   
325                    *->tGkvkkrpggeqNnlkTGeKAlelvveeievl<-*
326                       +G v+ rpgge          ++++ +e+      
327   gi|9081913   818    YGFVQYRPGGE----------YHINNPEMSKA    839  
329 Cystatin: domain 1 of 1, from 826 to 859: score 2.4, E = 3.9
330                 CS    ECEEEEET.STSHHHHHHHHHHHHHHHHHSSSSEEEEE   
331                    *->GglspvdpNendpevqealdfAlakyNeksndnylfel<-*
332                       Gg   +++    pe  +al+ A+  yN +  +ny++ l   
333   gi|9081913   826    GGEYHINN----PEMSKALHQAVRGYNPEYYNNYQSLL    859  
335 RNase_PH_C: domain 1 of 1, from 827 to 846: score 4.2, E = 2.3
336                 CS    SSSS.B.HHHHHHHHHHHHHH   
337                    *->GkgnglteelleealelAkeg<-*
338                       G +++++ +++ +al++A+ g   
339   gi|9081913   827    G-EYHINNPEMSKALHQAVRG    846  
341 Glu_synthase: domain 2 of 2, from 830 to 1216: score 857.3, E = 9e-255
342                 CS    -SS-HHHHHHHHHHHHC--T-HHHHHHHHHHHHTS.-S-SGGGGEEE
343                    *->hrnepeviktlqkavqvpveskpsydkYreplnertpigalrdlLef
344                       h n+pe++k l++av+    +   y +Y+ +l +r p++alrdlL++
345   gi|9081913   830    HINNPEMSKALHQAVRG--YNPEYYNNYQSLLQNR-PPTALRDLLKL 873  
347                 CS --SS--......--GGGS--HHHHHTTEEEEEB-CTTC-HHHHHHHHHHH
348                    kyaeepldtdkiipieevepaleikkrfctgaMSyGALSeeAheALAiAm
349                     ++++p      i+i+eve+++ i + fctg+MS+GALS+e+he+LAiAm
350   gi|9081913   874 QSNRAP------ISIDEVESIEDILQKFCTGGMSLGALSRETHETLAIAM 917  
352                 CS HHCT-EEEETTT---GGGCSB-TTS-T S BTTSTT--S--TT-B---SE
353                    nriGtksNtGEGGedperlkpaadlds.G.SpTlpHLkGLqnednarSAI
354                    nriG+ksN+GEGGedp r+k + d++s+G+Sp lpHLkGL+n+d+a+SAI
355   gi|9081913   918 NRIGGKSNSGEGGEDPVRFKILNDVNSsGtSPLLPHLKGLKNGDTASSAI 967  
357                 CS EEE-TT-TT--............HHHHCC-SEEEEE---TTSTTT--EE-
358                    kQvASGRFGVtkRnGefWeefkRseYLvnAdalEIKiAQGAKPGeGGhLP
359                    kQ+ASGRFGVt            +eYL+nA++lEIKiAQGAKPGeGG+LP
360   gi|9081913   968 KQIASGRFGVT------------PEYLMNAKQLEIKIAQGAKPGEGGQLP 1005 
362                 CS GGG--HHHHHHHTS-TT--EE--SS-TT-SSHHHHHHHHHHHHHH-.TTS
363                    GeKVspeIAriRnstPGvgliSPpPHHDIysiEDLaqLIydLkeindpkA
364                    G+K+sp+IA +R ++PGv liSPpPHHDIysiEDL+qLI+dL++in pkA
365   gi|9081913  1006 GKKISPYIATLRKCKPGVPLISPPPHHDIYSIEDLSQLIFDLHQIN-PKA 1054 
367                 CS EEEEEEE-STTHHHHHHH...HHHTT-SEEEEE-TT---SSEECCHHHHC
368                    pisVKLVsehgvgtiaaGhmqvakAnADiIlIdGhdGGTGASpktsikha
369                    +isVKLVse g+gtiaaG   vak+nADiI+I+GhdGGTGASp++sikha
370   gi|9081913  1055 KISVKLVSEIGIGTIAAG---VAKGNADIIQISGHDGGTGASPLSSIKHA 1101 
372                 CS ---HHHHHHHHHHHHHCTT-CCCSEEEEESS--SHHHHHHHHHCT-SEEE
373                    GlPwelgLaevhqtLvengLRdrVsLiadGGLrTGaDVakAaaLGAdavg
374                    G PwelgL+evhq+L en+LRdrV+L++dGGLrTG D+++Aa++GA+++g
375   gi|9081913  1102 GSPWELGLSEVHQLLAENQLRDRVTLRVDGGLRTGSDIVLAAIMGAEEFG 1151 
377                 CS -SHHHHHHCT--S---CCCT--TTSSS---CCHH..CT----HHHHHHHH
378                    iGTaaLiAlGCimaRvCHtntCPvGvATQDPeLrKrlkfegaperVvNyf
379                    +GT+a+iA+GCimaR+CHtn+CPvGvATQ++eLr   +f g+pe +vN+f
380   gi|9081913  1152 FGTVAMIATGCIMARICHTNKCPVGVATQREELR--ARFSGVPEALVNFF 1199 
382                 CS HHHHHHHHHHHHHHT-S   
383                    iflaeEvrellaqlGfr<-*
384                    +f+  Evre+la+lG++   
385   gi|9081913  1200 LFIGNEVREILASLGYK    1216 
387 DUF258: domain 1 of 1, from 839 to 860: score 0.3, E = 8.3
388                 CS    HHHHHHHCTSS-HHHHHHHHHHHH   
389                    *->AVkaAveeGeIseeRYesYlklle<-*
390                       A+ +Av    +++e Y++Y+ ll+   
391   gi|9081913   839    ALHQAVR--GYNPEYYNNYQSLLQ    860  
393 Pencillinase_R: domain 1 of 1, from 856 to 894: score 3.9, E = 2.5
394                 CS    XXXXXXXXXXXXXXXXXXX    XXXXXXXXXXXXXXXX   
395                    *->drlfggsvgalvanfleee....klSeddieeLrelLde<-*
396                       + l++++++ ++ ++l+ ++++ ++S d++e ++++L++   
397   gi|9081913   856    QSLLQNRPPTALRDLLKLQsnraPISIDEVESIEDILQK    894  
399 SelT: domain 1 of 1, from 872 to 885: score 3.1, E = 2.2
400                    *->KLqtGrvYAPPtpqEL<-*
401                       KLq++r   P++++E+   
402   gi|9081913   872    KLQSNRA--PISIDEV    885  
404 Nitro_FeMo-Co: domain 1 of 1, from 879 to 897: score 2.1, E = 5.3
405                 CS    EEE-TTSSBHHHHHHHHHC   
406                    *->pikagegetieeaiealqe<-*
407                       pi   e e+ie+ + ++ +   
408   gi|9081913   879    PISIDEVESIEDILQKFCT    897  
410 DUF37: domain 1 of 1, from 927 to 934: score 3.0, E = 4.5
411                    *->hpGGyDPV<-*
412                       ++GG DPV   
413   gi|9081913   927    GEGGEDPV    934  
415 Scm3: domain 1 of 1, from 953 to 963: score 2.2, E = 3.5
416                    *->HLraLeteddi<-*
417                       HL++L+++d++   
418   gi|9081913   953    HLKGLKNGDTA    963  
420 cobW: domain 1 of 1, from 1038 to 1058: score 5.1, E = 0.45
421                 CS    ...HHHHHHHHHH-SSS-EEE   
422                    *->adlekleadlrrlnpeapiip<-*
423                       +dl++l+ dl+++np+a+i     
424   gi|9081913  1038    EDLSQLIFDLHQINPKAKISV    1058 
426 Arch_flagellin: domain 1 of 1, from 1050 to 1072: score 4.1, E = 0.66
427                    *->inpstkvrgeVvpenGapgtief<-*
428                       inp  k+++++v+e+G+ ++      
429   gi|9081913  1050    INPKAKISVKLVSEIGIGTIAAG    1072 
431 DUF1393: domain 1 of 1, from 1055 to 1068: score 3.1, E = 2
432                    *->klSvKtVVAiGIGA<-*
433                       k+SvK V  iGIG+   
434   gi|9081913  1055    KISVKLVSEIGIGT    1068 
436 FtsK_SpoIIIE: domain 1 of 1, from 1107 to 1143: score 2.6, E = 3.1
437                    *->lviDnydeLaeenlL.ervtsLknqGlsygvhvmata<-*
438                       l++ + ++L +en+L++rvt+ + +Gl +g +++++a   
439   gi|9081913  1107    LGLSEVHQLLAENQLrDRVTLRVDGGLRTGSDIVLAA    1143 
441 FMN_dh: domain 1 of 1, from 1109 to 1148: score 3.2, E = 0.89
442                 CS    HHHHHHHHHCHHTTTSSEEEEESS-SSHHHHHHHHHHTSS   
443                    *->LpeVvPIlkeaAvkgdieVllDgGvRRGtDVlKALALGAr<-*
444                       L eV  +l e  + +++   +DgG R+G+D++ A  +GA+   
445   gi|9081913  1109    LSEVHQLLAENQLRDRVTLRVDGGLRTGSDIVLAAIMGAE    1148 
447 DSRB: domain 1 of 1, from 1120 to 1134: score 2.7, E = 2.7
448                    *->mKvndrvtvKtDGgpR<-*
449                        ++ drvt + DGg R   
450   gi|9081913  1120    -QLRDRVTLRVDGGLR    1134 
452 Phage_Mu_P: domain 1 of 1, from 1122 to 1131: score -0.4, E = 10
453                    *->sntVtLrvgG<-*
454                        ++VtLrv+G   
455   gi|9081913  1122    RDRVTLRVDG    1131 
457 Hormone_4: domain 1 of 1, from 1168 to 1176: score 4.4, E = 2.5
458                 CS    X-TT--TT-   
459                    *->CyirnCPrG<-*
460                       C  + CP+G   
461   gi|9081913  1168    CHTNKCPVG    1176 
463 GDC-P: domain 1 of 1, from 1205 to 1225: score 7.1, E = 0.086
464                    *->eqqeMLstiGlssLddLidat<-*
465                       e++e+L+++G++sLdd ++++   
466   gi|9081913  1205    EVREILASLGYKSLDDITGQN    1225 
468 PspB: domain 1 of 1, from 1268 to 1276: score 0.4, E = 8.4
469                    *->MsaffLagP<-*
470                       M+ ++La+P   
471   gi|9081913  1268    MDDDILAIP    1276 
473 T5orf172: domain 1 of 1, from 1271 to 1293: score 2.0, E = 6.1
474                    *->dvvalievedaraklEklLHkrFk<-*
475                       d+ a+ ev++a  klE+++ k+Fk   
476   gi|9081913  1271    DILAIPEVSNAI-KLETEITKHFK    1293 
478 CAP_C: domain 1 of 1, from 1283 to 1292: score 1.3, E = 7.4
479                 CS    EEEEEE----   
480                    *->KLvTevveha<-*
481                       KL+Te++ h    
482   gi|9081913  1283    KLETEITKHF    1292 
484 GXGXG: domain 1 of 1, from 1290 to 1485: score 367.3, E = 2.7e-107
485                 CS    EEEEE-TT--STTHHHHHHHHHHCTTTS.S-TTCEEEEEEEEE-TTT
486                    *->keeaiiNtdrlvgtrlsgeiakkygeegalpkdtgkivfnGsAGqsf
487                       k+++i Nt+r+vgtrlsg iak yg+ g + k+ +k++f+GsAGqsf
488   gi|9081913  1290    KHFKIANTNRTVGTRLSGIIAKNYGNTG-F-KGLIKLNFYGSAGQSF 1334 
490                 CS TTT-BTTEEEEEEEEE-S.TTTTT-ECCEEEEE--TT-.......SS-GG
491                    GafmagGvtLeleGdAnddyvGkgmsGGeIvikgnagdpvGnnMdageyv
492                    Gaf+a+G++L l+G+And yvGkgm+GG+Ivi+++ag         +e +
493   gi|9081913  1335 GAFLASGINLKLMGEAND-YVGKGMNGGSIVIVPPAGT-------IYEDN 1376 
495                 CS GSEEC-SSTTTT--CEEEEESSEE-TTTTTT-.....CCEEEEESEB.-S
496                    gnviaGNtclyGatGGkifiaGdAGerfgvrnkayKdsgatiVveGvaGd
497                    ++vi+GNtclyGatGG++f++G+AGerf+vrn     s a+ VveGv Gd
498   gi|9081913  1377 NQVIIGNTCLYGATGGYLFAQGQAGERFAVRN-----SLAESVVEGV-GD 1420 
500                 CS STTTT-EEEEEEESS-B-SSBTTT--CCEEEEE-TTS.......THHHHB
501                    hggEYMtGGtivVlGdaGrnvGagMtGGiaYvlgeiedfsyMiatlpgkv
502                    h++EYMtGG+ivVlG+aGrnvGagMtGG+aY+l+e+e        + ++v
503   gi|9081913  1421 HACEYMTGGVIVVLGKAGRNVGAGMTGGLAYFLDEDE-------NFIDRV 1463 
505                 CS -CCCEEEE...ES-S......CCHHHHHHHH   
506                    nleiVeledlkrievkrkklLpegekqlkel<-*
507                    n+eiV+ +   r+ +      ++ge+qlk+l   
508   gi|9081913  1464 NSEIVKIQ---RVIT------KAGEEQLKNL    1485 
510 DUF1514: domain 1 of 1, from 1453 to 1469: score 3.5, E = 5.7
511                    *->LeeyrieveRikkevkk<-*
512                       L e+++ ++R++ e+ k   
513   gi|9081913  1453    LDEDENFIDRVNSEIVK    1469 
515 Colicin: domain 1 of 1, from 1456 to 1467: score 1.4, E = 7.5
516                 CS    SHHHHHHHHHCH   
517                    *->DdkfveklNkli<-*
518                       D++f++ +N +i   
519   gi|9081913  1456    DENFIDRVNSEI    1467 
521 Ribosomal_S6: domain 1 of 1, from 1461 to 1481: score 3.3, E = 3.7
522                 CS    CCHHHHHHHHHHHHHCTT-EE   
523                    *->EqvkqeiekYqkvLtnngAei<-*
524                       ++v++ei k+q+v+t++g+e+   
525   gi|9081913  1461    DRVNSEIVKIQRVITKAGEEQ    1481 
527 BicD: domain 1 of 1, from 1465 to 1481: score -1.6, E = 6.8
528                    *->gqaysnqrkvAkdGeer<-*
529                        + +++qr+ +k Gee+   
530   gi|9081913  1465    SEIVKIQRVITKAGEEQ    1481 
532 PUF: domain 1 of 1, from 1470 to 1486: score 6.5, E = 0.47
533                    *->lQkllevateeqkqlil<-*
534                       +Q+++++a+eeq ++++   
535   gi|9081913  1470    IQRVITKAGEEQLKNLI    1486 
537 DUF477: domain 1 of 1, from 1472 to 1495: score 3.8, E = 1.7
538                    *->gtLspserarLeqalaalEqktga<-*
539                       ++++++  ++L   ++  ++ktg+   
540   gi|9081913  1472    RVITKAGEEQLKNLIENHAAKTGS    1495 
542 Phage_prot_Gp6: domain 1 of 1, from 1479 to 1492: score 1.0, E = 4
543                    *->eEmikkFidkHklr<-*
544                       eE +k++i+ H+++   
545   gi|9081913  1479    EEQLKNLIENHAAK    1492 
547 IBN_N: domain 1 of 1, from 1498 to 1516: score 8.2, E = 0.17
548                 CS    HHHHHHHHHCCTHHCHHHHH   
549                    *->AEkqLeqlekqklPgfllaL<-*
550                       A++ Le+++++ lP+f++ +   
551   gi|9081913  1498    AHTILEKWNSY-LPQFWQVV    1516 
553 GspM: domain 1 of 1, from 1506 to 1520: score 1.0, E = 8.6
554                 CS    XXXXXXXXXXXXXXX   
555                    *->mneLqawWqgrspRE<-*
556                       ++ L ++Wq ++p+E   
557   gi|9081913  1506    NSYLPQFWQVVPPSE    1520