Query         012032
Match_columns 472
No_of_seqs    122 out of 266
Neff          5.6 
Searched_HMMs 46136
Date          Fri Mar 29 07:54:08 2013
Command       hhsearch -i /work/01045/syshi/csienesis_hhblits_a3m/012032.a3m -d /work/01045/syshi/HHdatabase/Cdd.hhm -o /work/01045/syshi/hhsearch_cdd/012032hhsearch_cdd -cpu 12 -v 0 

 No Hit                             Prob E-value P-value  Score    SS Cols Query HMM  Template HMM
  1 KOG2233 Alpha-N-acetylglucosam 100.0  3E-149  6E-154 1145.2  39.6  439    1-472   134-574 (666)
  2 PF05089 NAGLU:  Alpha-N-acetyl 100.0  3E-113  7E-118  862.8  16.9  258    1-265    75-332 (333)
  3 PF12972 NAGLU_C:  Alpha-N-acet 100.0 1.1E-50 2.3E-55  403.3  13.1  178  271-472     1-178 (267)
  4 PF02449 Glyco_hydro_42:  Beta-  94.4   0.058 1.2E-06   56.4   5.2   93   21-119    47-140 (374)
  5 PF00150 Cellulase:  Cellulase   93.7    0.36 7.8E-06   47.0   9.0  133   16-184    57-200 (281)
  6 cd06593 GH31_xylosidase_YicI Y  81.2      21 0.00045   36.3  11.9  113   24-149    69-204 (308)
  7 PF12891 Glyco_hydro_44:  Glyco  80.7     1.7 3.7E-05   43.4   3.6  120   24-148    27-176 (239)
  8 PF13812 PPR_3:  Pentatricopept  71.1     2.6 5.6E-05   27.7   1.4   20   18-37     15-34  (34)
  9 PF01229 Glyco_hydro_39:  Glyco  67.7     8.3 0.00018   42.0   5.2  103   25-146    88-198 (486)
 10 cd06592 GH31_glucosidase_KIAA1  58.5      65  0.0014   32.8   9.5   82   24-117    73-166 (303)
 11 PF02065 Melibiase:  Melibiase;  58.1      29 0.00064   37.1   7.1  122   24-153   107-236 (394)
 12 PRK00035 hemH ferrochelatase;   53.4 2.4E+02  0.0051   29.0  12.7  152    8-180    69-268 (333)
 13 PF05119 Terminase_4:  Phage te  47.7      21 0.00046   29.8   3.3   31    8-38     55-85  (100)
 14 PRK15014 6-phospho-beta-glucos  44.9      47   0.001   36.4   6.3   68   23-117   112-179 (477)
 15 cd03198 GST_C_CLIC GST_C famil  44.1      65  0.0014   29.3   6.1   43  397-439     5-47  (134)
 16 cd06594 GH31_glucosidase_YihQ   43.5   2E+02  0.0043   29.6  10.3  105   24-143    74-204 (317)
 17 COG1422 Predicted membrane pro  41.3      85  0.0018   30.7   6.6   63  386-448    66-129 (201)
 18 cd04257 AAK_AK-HSDH AAK_AK-HSD  40.7      18 0.00039   37.0   2.2   68   23-114   169-236 (294)
 19 cd06388 PBP1_iGluR_AMPA_GluR4   40.3      89  0.0019   32.7   7.3   92   22-140   190-287 (371)
 20 cd04247 AAK_AK-Hom3 AAK_AK-Hom  40.2      14  0.0003   38.3   1.2   70   21-114   176-246 (306)
 21 PF03705 CheR_N:  CheR methyltr  39.6      24 0.00051   26.5   2.1   20   18-37     24-43  (57)
 22 TIGR00676 fadh2 5,10-methylene  36.5      39 0.00084   34.0   3.7   25   25-49    174-199 (272)
 23 PF10629 DUF2475:  Protein of u  36.3      16 0.00034   29.9   0.7   12   38-49      4-15  (71)
 24 cd06389 PBP1_iGluR_AMPA_GluR2   35.2 1.3E+02  0.0028   31.3   7.5   59   21-94    187-250 (370)
 25 PF00232 Glyco_hydro_1:  Glycos  34.8 1.3E+02  0.0028   32.6   7.6   68   21-116    99-166 (455)
 26 COG5005 Mu-like prophage prote  33.8      11 0.00025   34.2  -0.5   42  175-216    62-105 (140)
 27 PF13410 GST_C_2:  Glutathione   33.7      63  0.0014   24.6   3.8   28  418-446     3-30  (69)
 28 cd06394 PBP1_iGluR_Kainate_KA1  33.6      92   0.002   32.3   6.1   67   21-102   197-268 (333)
 29 PF02836 Glyco_hydro_2_C:  Glyc  33.6 1.9E+02   0.004   29.1   8.1  106   22-145    37-152 (298)
 30 cd06565 GH20_GcnA-like Glycosy  32.9 5.2E+02   0.011   26.3  11.4  121    8-152    54-184 (301)
 31 TIGR00825 EIIBC-GUT PTS system  31.7     9.3  0.0002   39.3  -1.6   42    5-47     10-58  (331)
 32 PF12876 Cellulase-like:  Sugar  31.4      72  0.0016   26.3   4.0   53   99-152     2-64  (88)
 33 PF08858 IDEAL:  IDEAL domain;   29.1 2.1E+02  0.0046   20.3   5.8   33  394-426     4-37  (37)
 34 TIGR03356 BGL beta-galactosida  29.0 3.4E+02  0.0073   29.3   9.6   99   19-145    92-210 (427)
 35 cd03309 CmuC_like CmuC_like. P  28.4 2.9E+02  0.0063   28.7   8.6  115    5-135   177-320 (321)
 36 TIGR00756 PPR pentatricopeptid  26.3      44 0.00096   21.2   1.5   19   19-37     15-33  (35)
 37 COG0276 HemH Protoheme ferro-l  25.9 7.9E+02   0.017   25.8  13.8  158    9-181    69-264 (320)
 38 PF14433 SUKH-3:  SUKH-3 immuni  25.5 1.4E+02  0.0031   26.7   5.1   46  252-298     3-49  (142)
 39 cd03201 GST_C_DHAR GST_C famil  25.3 2.2E+02  0.0048   24.8   6.2    8  445-452    91-98  (121)
 40 cd06599 GH31_glycosidase_Aec37  25.2 5.1E+02   0.011   26.5   9.7  109   23-143    75-211 (317)
 41 PRK10203 hypothetical protein;  25.1      64  0.0014   29.1   2.7   38    1-51     21-61  (122)
 42 cd03190 GST_C_ECM4_like GST_C   24.1 2.8E+02  0.0061   24.6   6.8   37  409-446    25-61  (142)
 43 PF07328 VirD1:  T-DNA border e  23.5   3E+02  0.0065   25.4   6.6   53  385-439    73-128 (147)
 44 cd06380 PBP1_iGluR_AMPA N-term  22.6 3.7E+02   0.008   27.5   8.2   17   21-37    194-210 (382)
 45 COG4334 Uncharacterized protei  22.6      46   0.001   29.7   1.3   39  159-208     6-45  (131)
 46 cd06597 GH31_transferase_CtsY   21.9 7.2E+02   0.016   25.8  10.2   85   24-119    88-190 (340)
 47 cd04243 AAK_AK-HSDH-like AAK_A  21.6      58  0.0013   33.3   2.0   67   24-114   169-235 (293)
 48 cd03186 GST_C_SspA GST_N famil  21.5 4.4E+02  0.0094   21.5   7.1   29  417-446    31-59  (107)
 49 TIGR01558 sm_term_P27 phage te  21.2      98  0.0021   27.0   3.1   31    8-38     64-94  (116)
 50 PRK13511 6-phospho-beta-galact  20.7 2.3E+02   0.005   31.0   6.5   71   18-117    89-161 (469)
 51 PF13041 PPR_2:  PPR repeat fam  20.5      74  0.0016   22.9   1.9   22   18-39     17-38  (50)
 52 PF09350 DUF1992:  Domain of un  20.2      74  0.0016   25.8   1.9   38    1-51     15-56  (71)

No 1  
>KOG2233 consensus Alpha-N-acetylglucosaminidase [Intracellular trafficking, secretion, and vesicular transport]
Probab=100.00  E-value=2.8e-149  Score=1145.18  Aligned_cols=439  Identities=46%  Similarity=0.870  Sum_probs=427.5

Q ss_pred             CCCCCCCCCCCCHHHHHHHHHHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCcccccccc
Q 012032            1 MSNLHGWGGPLPQSWLDQQLVLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLL   80 (472)
Q Consensus         1 MgNi~gwgGPLp~~wi~~q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~L   80 (472)
                      ||||++||||||++|+..|.-||||||+|||+||||||||+|+||||++++++||++++++++.|+.|+  ++|||.++|
T Consensus       134 MGNl~awgGpLs~aw~~~ql~LqkrIidrm~~lGmTpvLPaFaG~VP~al~rlfPesnf~rl~rWn~f~--s~~~C~l~v  211 (666)
T KOG2233|consen  134 MGNLHAWGGPLSPAWMLNQLLLQKRIIDRMLELGMTPVLPAFAGHVPDALERLFPESNFTRLPRWNNFT--SRYSCMLLV  211 (666)
T ss_pred             hcCccccCCCCCHHHHHHHHHHHHHHHHHHHHcCCCccchhhccccHHHHHHhCchhceeeccccCCCC--cceeeeEEc
Confidence            999999999999999999999999999999999999999999999999999999999999999999996  699999999


Q ss_pred             CCCChHHHHHHHHHHHHHHHHhCCCCcccccCCCCCCCCCCCChHHHHHHHHHHHHHHhccCCCceEEEecccCCCCCCC
Q 012032           81 DATDPLFIEIGRAFIEQQLKEYGRTSHIYNCDTFDENTPPVDSPEYISSLGAAIYSGMQSGDSDAVWLMQGWLFSYDPFW  160 (472)
Q Consensus        81 dP~DplF~~I~~~F~~eq~~~fG~~~h~Y~~D~FnE~~pp~~dp~~L~~~~~~iy~am~~~dP~AvWvmQgW~F~~~~fW  160 (472)
                      .|+||||.+||.+|+++|+|+||.++|+|+||||||+.||+++|+|+++.+.+||+||+++||+||||||||+|.+|+||
T Consensus       212 ~P~dplF~eIgs~Flr~~~kefG~~tniy~~DpFNE~~Pp~sepey~~staaAiyesm~kvdknaVWllQgWlF~~d~FW  291 (666)
T KOG2233|consen  212 SPFDPLFQEIGSTFLRHQIKEFGGVTNIYSADPFNEILPPESEPEYVKSTAAAIYESMKKVDKNAVWLLQGWLFTYDPFW  291 (666)
T ss_pred             cCCcchHHHHHHHHHHHHHHHhCCcccccccCcccccCCCCCChHHHHHHHHHHHHHHhccCcceEEeeecceeecCCCC
Confidence            99999999999999999999999889999999999999999999999999999999999999999999999999999999


Q ss_pred             CchhHHHhHhCCCCCCEEEEecCCCcccccccccCcCCCCceeeeccCCCCccccccchhhhhcChHHhhhCCCCceEEe
Q 012032          161 RPPQMKALLNSVPLGKLVVLDLFAEVKPIWSTSKQFYGVPYIWCMLHNFAGNIEMYGILDSIAFGPVEARTSENTTMVGV  240 (472)
Q Consensus       161 ~~~~~~a~L~~Vp~~~mliLDL~~E~~p~W~~t~~f~G~pwIWc~LhNFGGn~gl~G~l~~i~~~~~~a~~~~~~~m~Gi  240 (472)
                      +++++||||++||.||||||||++|..|+|++|.+|+|+|||||||||||||.+|+|.++.|++||.+|+..+||+|||+
T Consensus       292 ~~~~ikafL~avP~GrllVLDLyaEv~P~~~~Tasf~GqpfIWCmLHNFGGn~~lfGal~~InsG~e~Ar~~~nStlVGt  371 (666)
T KOG2233|consen  292 PPWQIKAFLSAVPRGRLLVLDLYAEVVPAWQMTASFQGQPFIWCMLHNFGGNRELFGALQKINSGPELARMKPNSTLVGT  371 (666)
T ss_pred             ChHHHHHHHhcCCcCcEEEEehhhhhhhHHHhhhhhcCCceeEEEeeccCCcHhhhhhHHHhccCHHHHhcCCCCceeec
Confidence            99999999999999999999999999999999999999999999999999999999999999999999999999999999


Q ss_pred             eeCccccccChhHHHHHHhhcCCCCCCCHHHHHHHHHhcccCCCChhHHHHHHHHHhcccCCCCCCCCCC-CcceecccC
Q 012032          241 GMSMEGIEQNPVVYDLMSEMAFQHEKVDVKAWINQYSVRRYGRSVPAIQDAWNVLYHTVYNCTDGATDKN-RDVIVAFPD  319 (472)
Q Consensus       241 G~tpEGie~NpvvYeL~~d~aW~~~~id~~~W~~~Ya~rRYG~~~~~~~~AW~iL~~tvY~~~~~~~~~~-~~~~~~~P~  319 (472)
                      ||+||||+||||||+||.|++|+++++|++.|+++|++||||+.++.+++||.+|.+|||||+++...++ ...+..||+
T Consensus       372 Gm~~EgI~QN~VvYsf~~e~~wsde~ldl~~Wlksys~sRY~~~~~~~eaaW~lL~~tvYn~~~~~~~~~~~~~l~rRp~  451 (666)
T KOG2233|consen  372 GMSMEGIDQNYVVYSFMIERGWSDEPLDLNNWLKSYSESRYGVDFKVAEAAWTLLAGTVYNCPGKWATRGFSYFLYRRPS  451 (666)
T ss_pred             ccCccccccchhhHHHHHHhccCCCCCCHHHHHHHHHHhhhccchHHHHHHHHHHhhhhcCCCcccccCCCceeeEeccc
Confidence            9999999999999999999999999999999999999999999999999999999999999998777666 666777787


Q ss_pred             CCCccccccccccccCCCCcccchhccccCCCCCCCCccCCHHHHHHHHHHHHhccccCCCC-CcccchHHHHHHHHHHH
Q 012032          320 VDPSIISVTEGKYQNYGKPVSKEAVLKSETSSYDHPHLWYSTSEVIRALELFIASGNELSAS-NTYRYDLIDLTRQALAK  398 (472)
Q Consensus       320 l~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~Yd~~~l~~A~~lll~~~~~l~~~-~~y~yDLvDvtRQvL~n  398 (472)
                      ++                               .++.+||+.+++..||++||.+.+.+..+ ++|||||||||||+|+.
T Consensus       452 f~-------------------------------~k~~~Wyn~sev~~a~rllL~a~~~l~~e~~~fr~DlvDltRq~lqe  500 (666)
T KOG2233|consen  452 FQ-------------------------------RKTEYWYNVSEVFSAWRLLLTALVHLLGEHPLFRYDLVDLTRQMLQE  500 (666)
T ss_pred             cc-------------------------------CCceeeecHHHHHHHHHHHHhhhhhhccCchhhHhHHHHHHHHHHHH
Confidence            65                               23679999999999999999998887666 99999999999999999


Q ss_pred             HHHHHHHHHHHHHHhCChHHHHHHHHHHHHHHHHHHHHhccCCCCChhHHHHHHHhhCCCHHHHHhhhCcccCC
Q 012032          399 YANELFLNIIEAYQLNDAHGVFQLSRRFLELVEDMDGLLACHDGFLLGPWLESAKQLAQNEEQEKQVRCPYVSQ  472 (472)
Q Consensus       399 ~~~~~~~~~~~Ay~~~d~~~~~~~~~~~l~li~dlD~LL~t~~~FlLg~Wl~~Ar~~a~~~~ek~~yE~NAr~q  472 (472)
                      .++++|.++.+||..||...+.++|.++++|++|||.+|+++.+||||+||++||+.|.+.+|++.||+|||+|
T Consensus       501 lanq~y~e~~sAflkkd~~sl~~~~~~llelf~dle~~lasd~nfLlg~WleqAk~~A~n~~er~~fe~nar~q  574 (666)
T KOG2233|consen  501 LANQAYLEARSAFLKKDKQSLGALSEKLLELFQDLESYLASDDNFLLGQWLEQAKSAAPNSEERQVFEVNARDQ  574 (666)
T ss_pred             HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHccCcchhHHHHHHHHHhhCCChhhhhhccccccce
Confidence            99999999999999999999999999999999999999999999999999999999999999999999999998


No 2  
>PF05089 NAGLU:  Alpha-N-acetylglucosaminidase (NAGLU) tim-barrel domain;  InterPro: IPR024733 Alpha-N-acetylglucosaminidase is a lysosomal enzyme that is required for the stepwise degradation of heparan sulphate []. Mutations on the alpha-N-acetylglucosaminidase gene can lead to mucopolysaccharidosis type IIIB (MPS IIIB; or Sanfilippo syndrome type B), characterised by neurological dysfunction but relatively mild somatic manifestations [].  Alpha-N-acetylglucosaminidase is composed of three domains. This entry represents the central domain, which has a tim barrel fold [].; PDB: 4A4A_A 2VC9_A 2VCC_A 2VCB_A 2VCA_A.
Probab=100.00  E-value=3.3e-113  Score=862.79  Aligned_cols=258  Identities=59%  Similarity=1.153  Sum_probs=188.8

Q ss_pred             CCCCCCCCCCCCHHHHHHHHHHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCcccccccc
Q 012032            1 MSNLHGWGGPLPQSWLDQQLVLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLL   80 (472)
Q Consensus         1 MgNi~gwgGPLp~~wi~~q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~L   80 (472)
                      ||||+|||||||++||++|++||||||+|||||||+||||||+||||++||+++|+|+|+++++|++|      |||++|
T Consensus        75 MgNl~gwgGPLp~~w~~~q~~Lq~kIl~RmreLGm~PVLPaF~G~VP~~~~~~~P~a~i~~~~~W~~f------~~~~~L  148 (333)
T PF05089_consen   75 MGNLQGWGGPLPQSWIDQQAELQKKILDRMRELGMTPVLPAFAGHVPRAFKRKYPNANITRQGNWNGF------CRPYFL  148 (333)
T ss_dssp             TTS--STT----TTHHHHHHHHHHHHHHHHHHHT-EEEEE--S-EE-TTHHHHSTT--EE---EETTE------E--EEE
T ss_pred             hCCcccCCCCCCHHHHHHHHHHHHHHHHHHHHcCCcccCCCcCCCCChHHHhcCCCCEEeeCCCcCCC------CCCcee
Confidence            99999999999999999999999999999999999999999999999999999999999999999986      489999


Q ss_pred             CCCChHHHHHHHHHHHHHHHHhCCCCcccccCCCCCCCCCCCChHHHHHHHHHHHHHHhccCCCceEEEecccCCCCCCC
Q 012032           81 DATDPLFIEIGRAFIEQQLKEYGRTSHIYNCDTFDENTPPVDSPEYISSLGAAIYSGMQSGDSDAVWLMQGWLFSYDPFW  160 (472)
Q Consensus        81 dP~DplF~~I~~~F~~eq~~~fG~~~h~Y~~D~FnE~~pp~~dp~~L~~~~~~iy~am~~~dP~AvWvmQgW~F~~~~fW  160 (472)
                      ||+||+|++||++||+||+|+|| ++|||++|||||++||+++|+||+++|++||++|+++||+||||||||+|++++||
T Consensus       149 ~P~dplF~~i~~~F~~~q~~~yG-~~~~Y~~D~FnE~~p~~~~~~~l~~~s~~v~~am~~~dp~AvWvmQgWlF~~~~fW  227 (333)
T PF05089_consen  149 DPTDPLFAEIAKLFYEEQIKLYG-TDHIYAADPFNEGGPPSGDPEYLANVSKAVYKAMQAADPDAVWVMQGWLFYYDPFW  227 (333)
T ss_dssp             -SS--HHHHHHHHHHHHHHHHH----SEEE--TTTTS---TTS---HHHHHHHHHHHHHHH-TT-EEEEEE--------B
T ss_pred             CCCCchHHHHHHHHHHHHHHhcC-CCceeCCCccCCCCCCCCchHHHHHHHHHHHHHHHhhCCCcEEEEccccccccccc
Confidence            99999999999999999999999 78999999999999999999999999999999999999999999999999999999


Q ss_pred             CchhHHHhHhCCCCCCEEEEecCCCcccccccccCcCCCCceeeeccCCCCccccccchhhhhcChHHhhhCCCCceEEe
Q 012032          161 RPPQMKALLNSVPLGKLVVLDLFAEVKPIWSTSKQFYGVPYIWCMLHNFAGNIEMYGILDSIAFGPVEARTSENTTMVGV  240 (472)
Q Consensus       161 ~~~~~~a~L~~Vp~~~mliLDL~~E~~p~W~~t~~f~G~pwIWc~LhNFGGn~gl~G~l~~i~~~~~~a~~~~~~~m~Gi  240 (472)
                      ++++++|||++||++|||||||+||..|+|+++++|+|||||||||||||||+||+|+++.|++++.+|+++++++||||
T Consensus       228 ~~~~~~a~L~~Vp~~~mliLDL~se~~p~w~~~~~f~G~pwIwc~L~NFGG~~gl~G~~~~i~~~~~~a~~~~~~~m~G~  307 (333)
T PF05089_consen  228 TPNPIKALLSGVPKGRMLILDLFSERFPQWKRTESFYGKPWIWCMLHNFGGNTGLYGNLENIASGPIEARASPNSNMVGI  307 (333)
T ss_dssp             TTBS-HHHHTT-SGGGEEEEETTTTTS---HHHHCTT---EEEEE---STT--SS---HHHHHHHHHHHHHT--S-EEEE
T ss_pred             CcchHHHHHcCCCCCCeEEEEccccccchhccccccccchhhhhcccCCCCCCCCcccHHHHHhhHHHHHHhcCCCceEE
Confidence            99999999999999999999999999999999999999999999999999999999999999999999999989999999


Q ss_pred             eeCccccccChhHHHHHHhhcCCCC
Q 012032          241 GMSMEGIEQNPVVYDLMSEMAFQHE  265 (472)
Q Consensus       241 G~tpEGie~NpvvYeL~~d~aW~~~  265 (472)
                      |+|||||||||||||||+|++|+++
T Consensus       308 G~tpEgi~~NpvvYeL~~e~aW~~~  332 (333)
T PF05089_consen  308 GLTPEGIEQNPVVYELMLEMAWRKD  332 (333)
T ss_dssp             EE--S-S-S-HHHHHHHHHGGG-SS
T ss_pred             EecccccccCHHHHHHHHHhcccCC
Confidence            9999999999999999999999975


No 3  
>PF12972 NAGLU_C:  Alpha-N-acetylglucosaminidase (NAGLU) C-terminal domain;  InterPro: IPR024732 Alpha-N-acetylglucosaminidase is a lysosomal enzyme required for the stepwise degradation of heparan sulphate []. Mutations on the alpha-N-acetylglucosaminidase (NAGLU) gene can lead to Mucopolysaccharidosis type IIIB (MPS IIIB, or Sanfilippo syndrome type B) characterised by neurological dysfunction but relatively mild somatic manifestations []. The structure shows that the enzyme is composed of three domains. This C-terminal domain has an all alpha helical fold [].; PDB: 2VC9_A 2VCC_A 2VCB_A 2VCA_A 4A4A_A.
Probab=100.00  E-value=1.1e-50  Score=403.32  Aligned_cols=178  Identities=38%  Similarity=0.753  Sum_probs=146.2

Q ss_pred             HHHHHHHhcccCCCChhHHHHHHHHHhcccCCCCCCCCCCCcceecccCCCCccccccccccccCCCCcccchhccccCC
Q 012032          271 AWINQYSVRRYGRSVPAIQDAWNVLYHTVYNCTDGATDKNRDVIVAFPDVDPSIISVTEGKYQNYGKPVSKEAVLKSETS  350 (472)
Q Consensus       271 ~W~~~Ya~rRYG~~~~~~~~AW~iL~~tvY~~~~~~~~~~~~~~~~~P~l~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  350 (472)
                      +||++|++||||+.++++.+||++|++|||++++...+++++++|+||+++..                        ..+
T Consensus         1 ~Wi~~Y~~rRYG~~~~~~~~AW~iL~~tvY~~~~~~~~~~~si~~~rPsl~~~------------------------~~~   56 (267)
T PF12972_consen    1 EWIKDYATRRYGKYDPEAREAWQILLRTVYNNTGGQYGRGESIFCARPSLNGN------------------------SAS   56 (267)
T ss_dssp             HHHHHHHHHHHT---HHHHHHHHHHHHTTTS---S---SS--GGGS-S-SS---------------------------SS
T ss_pred             ChHHHHHHccCCCCCHHHHHHHHHHHhhhCCCCCCCCCCCcCceeeCCCCCcc------------------------ccc
Confidence            69999999999999999999999999999999887777789999999998743                        235


Q ss_pred             CCCCCCccCCHHHHHHHHHHHHhccccCCCCCcccchHHHHHHHHHHHHHHHHHHHHHHHHHhCChHHHHHHHHHHHHHH
Q 012032          351 SYDHPHLWYSTSEVIRALELFIASGNELSASNTYRYDLIDLTRQALAKYANELFLNIIEAYQLNDAHGVFQLSRRFLELV  430 (472)
Q Consensus       351 ~~~~~~~~Yd~~~l~~A~~lll~~~~~l~~~~~y~yDLvDvtRQvL~n~~~~~~~~~~~Ay~~~d~~~~~~~~~~~l~li  430 (472)
                      +|+++.++|||.+|++||++|++++++|+++++|||||||||||||+|+++.+|.++++||+++|.+.|++.+++||+||
T Consensus        57 ~~~~~~~~Yd~~~l~~A~~~ll~a~~~l~~~~~yryDlvDvtRQvL~n~~~~~~~~~~~ay~~~d~~~~~~~~~~~l~ll  136 (267)
T PF12972_consen   57 TWGPPGIWYDPADLEKAWRLLLKAADELKDSETYRYDLVDVTRQVLSNYADELYQQLVDAYNSKDIEAFKALSARFLELL  136 (267)
T ss_dssp             TT------S-HHHHHHHHHHHHHCHHHHTT-HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHTT-HHHHHHHHHHHHHHH
T ss_pred             CCCCCCCCCCHHHHHHHHHHHHHhHHhcCCCchHHhHHHHHHHHHHHHHHHHHHHHHHHHHHhCCHHHHHHHHHHHHHHH
Confidence            67889999999999999999999999999999999999999999999999999999999999999999999999999999


Q ss_pred             HHHHHHhccCCCCChhHHHHHHHhhCCCHHHHHhhhCcccCC
Q 012032          431 EDMDGLLACHDGFLLGPWLESAKQLAQNEEQEKQVRCPYVSQ  472 (472)
Q Consensus       431 ~dlD~LL~t~~~FlLg~Wl~~Ar~~a~~~~ek~~yE~NAr~q  472 (472)
                      +|||+||+||++|+||+||++||++|.|++||++||+|||+|
T Consensus       137 ~dlD~lL~t~~~f~Lg~Wi~~Ar~~g~~~~e~~~yE~NAR~q  178 (267)
T PF12972_consen  137 DDLDRLLATNPEFLLGKWIEDARAWGTTPEEKDLYEYNARNQ  178 (267)
T ss_dssp             HHHHHHHTT-GGGBHHHHHHHHHHSSTT--HHHHHHHHHHHH
T ss_pred             HHHHHHHCcCCCCCHHHHHHHHHHHCCCHHHHHHHHHHHHHH
Confidence            999999999999999999999999999999999999999987


No 4  
>PF02449 Glyco_hydro_42:  Beta-galactosidase;  InterPro: IPR013529 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. This group of beta-galactosidase enzymes (3.2.1.23 from EC) belong to the glycosyl hydrolase 42 family GH42 from CAZY. The enzyme catalyses the hydrolysis of terminal, non-reducing terminal beta-D-galactosidase residues.; GO: 0004565 beta-galactosidase activity, 0005975 carbohydrate metabolic process, 0009341 beta-galactosidase complex; PDB: 1KWK_A 1KWG_A 3U7V_A.
Probab=94.36  E-value=0.058  Score=56.38  Aligned_cols=93  Identities=18%  Similarity=0.211  Sum_probs=70.5

Q ss_pred             HHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceec-cCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHH
Q 012032           21 VLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQ-LGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQL   99 (472)
Q Consensus        21 ~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~-~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~   99 (472)
                      ..-.++|+++++.||..||--+.+..|.-+.+++|++..+. .|...++..      ....++.+|.|.+..+.|+++..
T Consensus        47 ~~lD~~l~~a~~~Gi~viL~~~~~~~P~Wl~~~~Pe~~~~~~~g~~~~~g~------~~~~~~~~p~yr~~~~~~~~~l~  120 (374)
T PF02449_consen   47 SWLDRVLDLAAKHGIKVILGTPTAAPPAWLYDKYPEILPVDADGRRRGFGS------RQHYCPNSPAYREYARRFIRALA  120 (374)
T ss_dssp             HHHHHHHHHHHCTT-EEEEEECTTTS-HHHHCCSGCCC-B-TTTSBEECCC------STT-HCCHHHHHHHHHHHHHHHH
T ss_pred             HHHHHHHHHHHhccCeEEEEecccccccchhhhcccccccCCCCCcCccCC------ccccchhHHHHHHHHHHHHHHHH
Confidence            34567899999999999998778889999999999986643 333322211      12357899999999999999999


Q ss_pred             HHhCCCCcccccCCCCCCCC
Q 012032          100 KEYGRTSHIYNCDTFDENTP  119 (472)
Q Consensus       100 ~~fG~~~h~Y~~D~FnE~~p  119 (472)
                      +.|++..++.+.+..||..-
T Consensus       121 ~~y~~~p~vi~~~i~NE~~~  140 (374)
T PF02449_consen  121 ERYGDHPAVIGWQIDNEPGY  140 (374)
T ss_dssp             HHHTTTTTEEEEEECCSTTC
T ss_pred             hhccccceEEEEEeccccCc
Confidence            99997777889999999765


No 5  
>PF00150 Cellulase:  Cellulase (glycosyl hydrolase family 5);  InterPro: IPR001547 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. Glycoside hydrolase family 5 GH5 from CAZY comprises enzymes with several known activities; endoglucanase (3.2.1.4 from EC); beta-mannanase (3.2.1.78 from EC); exo-1,3-glucanase (3.2.1.58 from EC); endo-1,6-glucanase (3.2.1.75 from EC); xylanase (3.2.1.8 from EC); endoglycoceramidase (3.2.1.123 from EC). The microbial degradation of cellulose and xylans requires several types of enzymes. Fungi and bacteria produces a spectrum of cellulolytic enzymes (cellulases) and xylanases which, on the basis of sequence similarities, can be classified into families. One of these families is known as the cellulase family A [] or as the glycosyl hydrolases family 5 []. One of the conserved regions in this family contains a conserved glutamic acid residue which is potentially involved [] in the catalytic mechanism.; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 3NDY_A 3NDZ_B 1LF1_A 1TVP_B 1TVN_A 3AYR_A 3AYS_A 1QI0_A 1W3K_A 1OCQ_A ....
Probab=93.67  E-value=0.36  Score=47.03  Aligned_cols=133  Identities=13%  Similarity=0.140  Sum_probs=81.9

Q ss_pred             HHHHHHHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHH
Q 012032           16 LDQQLVLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFI   95 (472)
Q Consensus        16 i~~q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~   95 (472)
                      -+...+-=++||+..++.||..||--...                  +.|..-.         -..+..+...+--+.|.
T Consensus        57 ~~~~~~~ld~~v~~a~~~gi~vild~h~~------------------~~w~~~~---------~~~~~~~~~~~~~~~~~  109 (281)
T PF00150_consen   57 DETYLARLDRIVDAAQAYGIYVILDLHNA------------------PGWANGG---------DGYGNNDTAQAWFKSFW  109 (281)
T ss_dssp             THHHHHHHHHHHHHHHHTT-EEEEEEEES------------------TTCSSST---------STTTTHHHHHHHHHHHH
T ss_pred             cHHHHHHHHHHHHHHHhCCCeEEEEeccC------------------ccccccc---------cccccchhhHHHHHhhh
Confidence            34555555899999999999998732222                  4463211         01222222333334477


Q ss_pred             HHHHHHhCCCCcccccCCCCCCCCCCC-------ChHHHHHHHHHHHHHHhccCCCceEEEec--ccCCCCCCCCchhHH
Q 012032           96 EQQLKEYGRTSHIYNCDTFDENTPPVD-------SPEYISSLGAAIYSGMQSGDSDAVWLMQG--WLFSYDPFWRPPQMK  166 (472)
Q Consensus        96 ~eq~~~fG~~~h~Y~~D~FnE~~pp~~-------dp~~L~~~~~~iy~am~~~dP~AvWvmQg--W~F~~~~fW~~~~~~  166 (472)
                      +...+.|.+...+.+.|++||-.....       .+..+....++++++|+++||+..=++.|  |.        .....
T Consensus       110 ~~la~~y~~~~~v~~~el~NEP~~~~~~~~w~~~~~~~~~~~~~~~~~~Ir~~~~~~~i~~~~~~~~--------~~~~~  181 (281)
T PF00150_consen  110 RALAKRYKDNPPVVGWELWNEPNGGNDDANWNAQNPADWQDWYQRAIDAIRAADPNHLIIVGGGGWG--------ADPDG  181 (281)
T ss_dssp             HHHHHHHTTTTTTEEEESSSSGCSTTSTTTTSHHHTHHHHHHHHHHHHHHHHTTSSSEEEEEEHHHH--------TBHHH
T ss_pred             hhhccccCCCCcEEEEEecCCccccCCccccccccchhhhhHHHHHHHHHHhcCCcceeecCCCccc--------cccch
Confidence            788888943456789999999776533       23567889999999999999995444555  52        22222


Q ss_pred             HhHhCCCC--CCEEEEecCC
Q 012032          167 ALLNSVPL--GKLVVLDLFA  184 (472)
Q Consensus       167 a~L~~Vp~--~~mliLDL~~  184 (472)
                      +... .|.  +..+|++...
T Consensus       182 ~~~~-~P~~~~~~~~~~~H~  200 (281)
T PF00150_consen  182 AAAD-NPNDADNNDVYSFHF  200 (281)
T ss_dssp             HHHH-STTTTTTSEEEEEEE
T ss_pred             hhhc-CcccccCceeEEeeE
Confidence            2223 663  4667777643


No 6  
>cd06593 GH31_xylosidase_YicI YicI alpha-xylosidase is a glycosyl hydrolase family 31 (GH31) enzyme that catalyzes the release of an alpha-xylosyl residue from the non-reducing end of alpha-xyloside substrates such as alpha-xylosyl fluoride and isoprimeverose. YicI forms a homohexamer (a trimer of dimers). All GH31 enzymes cleave a terminal carbohydrate moiety from a substrate that varies considerably in size, depending on the enzyme, and may be either a starch or a glycoprotein. The YicI family corresponds to subgroup 4 in the Ernst et al classification of GH31 enzymes.
Probab=81.19  E-value=21  Score=36.33  Aligned_cols=113  Identities=19%  Similarity=0.290  Sum_probs=63.0

Q ss_pred             HHHHHHHHHcCCeecc---CCCCCCCchhhHhh--------CCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHH
Q 012032           24 KKILVRIYELGMNPVL---PAFSGNVPAALQNV--------FPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGR   92 (472)
Q Consensus        24 kkIl~RmrelGM~PVL---PgF~G~VP~~~k~~--------~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~   92 (472)
                      +++++++++.||+.|+   |...-- .+.+++.        .++-+......|.+-        ..++|.+.|.-.+-=+
T Consensus        69 ~~~i~~l~~~G~~~~~~~~P~i~~~-~~~~~e~~~~g~~v~~~~g~~~~~~~w~g~--------~~~~Dftnp~a~~w~~  139 (308)
T cd06593          69 EGMLSRLKEKGFKVCLWINPYIAQK-SPLFKEAAEKGYLVKKPDGSVWQWDLWQPG--------MGIIDFTNPDACKWYK  139 (308)
T ss_pred             HHHHHHHHHCCCeEEEEecCCCCCC-chhHHHHHHCCeEEECCCCCeeeecccCCC--------cccccCCCHHHHHHHH
Confidence            6899999999999976   543221 1123332        222223333345431        3578888887654332


Q ss_pred             HHHHHHHHHhCCCCcccccCCCCCCCCCCCC-----------hHHHHHHHHHHHHHHhccCCCc-eEEE
Q 012032           93 AFIEQQLKEYGRTSHIYNCDTFDENTPPVDS-----------PEYISSLGAAIYSGMQSGDSDA-VWLM  149 (472)
Q Consensus        93 ~F~~eq~~~fG~~~h~Y~~D~FnE~~pp~~d-----------p~~L~~~~~~iy~am~~~dP~A-vWvm  149 (472)
                      .-+++..+ +| . ..+-+| |+|..|....           ..|--..++++|+++++..++- .+++
T Consensus       140 ~~~~~~~~-~G-i-d~~~~D-~~e~~p~~~~~~~g~~~~~~hn~y~~~~~~~~~~~~~~~~~~~r~~~~  204 (308)
T cd06593         140 DKLKPLLD-MG-V-DCFKTD-FGERIPTDVVYYDGSDGEKMHNYYALLYNKAVYEATKEVKGEGEAVVW  204 (308)
T ss_pred             HHHHHHHH-hC-C-cEEecC-CCCCCCccccccCCCCcceeeeHHHHHHHHHHHHHHHHhcCCCCeEEE
Confidence            22332222 46 3 456677 6675554321           1244567889999998887752 4443


No 7  
>PF12891 Glyco_hydro_44:  Glycoside hydrolase family 44;  InterPro: IPR024745 This is a family of putative bacterial glycoside hydrolases.; PDB: 3IK2_A 3ZQ9_A 2YJQ_B 2YKK_A 2YIH_A 2EEX_A 2EQD_A 2E0P_A 2E4T_A 2EO7_A ....
Probab=80.66  E-value=1.7  Score=43.38  Aligned_cols=120  Identities=18%  Similarity=0.236  Sum_probs=68.6

Q ss_pred             HHHHHHHHHcCCeec--cCCCCCCCchhh-----HhhCCCCceeccCCCC--------CCCCCCccccccccCCCChHHH
Q 012032           24 KKILVRIYELGMNPV--LPAFSGNVPAAL-----QNVFPSAKITQLGNWF--------SVKSDPRWCCTYLLDATDPLFI   88 (472)
Q Consensus        24 kkIl~RmrelGM~PV--LPgF~G~VP~~~-----k~~~P~a~i~~~~~W~--------gf~~~~~~~~~~~LdP~DplF~   88 (472)
                      .+.+++-++-|+.++  || -.|.|++.=     ...||.++.-+|..|.        |......  .+.-++|..+.=.
T Consensus        27 ~~f~~~~~~~ga~~m~T~p-m~G~Vakd~~~~~~~~~fp~~~y~~Q~~~d~~~~~~Gng~~~~~~--~~~~~~P~~~~~~  103 (239)
T PF12891_consen   27 DTFIDQNLAAGAYSMMTLP-MIGYVAKDANSVSESESFPSWRYGPQQWFDPWNPDCGNGVKPDKT--ALTSNDPDTPDNP  103 (239)
T ss_dssp             HHHHHHHHHTT-EEEEEE---SSEEES-BSEGBGGGTSSSTTEEEBS-EETTEEEEE-SEESTSS--S--SSSGGSSSSE
T ss_pred             HHHHHHhhhcCcceeEeec-ccceEecCCCCcccccCCChhhcccccccCcCcCCCCccccCCCC--CCCCCCCCCCccH
Confidence            456677778898888  44 378888876     6789999888866553        1110000  1222344422122


Q ss_pred             HHHHHHHHHHHHHhCCC-----CcccccCCCCCCC----------CCCCChHHHHHHHHHHHHHHhccCCCceEE
Q 012032           89 EIGRAFIEQQLKEYGRT-----SHIYNCDTFDENT----------PPVDSPEYISSLGAAIYSGMQSGDSDAVWL  148 (472)
Q Consensus        89 ~I~~~F~~eq~~~fG~~-----~h~Y~~D~FnE~~----------pp~~dp~~L~~~~~~iy~am~~~dP~AvWv  148 (472)
                      .....|+..+++.||+.     -+||.+|  ||..          |-.-..+.|.+-.-+.=++++++||.|.=+
T Consensus       104 ~y~~ewV~~l~~~~g~a~~~~gvk~y~lD--NEP~LW~~TH~dVHP~~~t~~El~~r~i~~AkaiK~~DP~a~v~  176 (239)
T PF12891_consen  104 VYMDEWVNYLVNKYGNASTNGGVKYYSLD--NEPDLWHSTHRDVHPEPVTYDELRDRSIEYAKAIKAADPDAKVF  176 (239)
T ss_dssp             EEHHHHHHHHHHHH--TTSTTS--EEEES--S-GGGHHHHTTTT--S---HHHHHHHHHHHHHHHHHH-TTSEEE
T ss_pred             hHHHHHHHHHHHHHhccccCCCceEEEec--CchHhhcccccccCCCCCCHHHHHHHHHHHHHHHHhhCCCCeEe
Confidence            23577888899999964     2799999  8853          222233456677777778999999999754


No 8  
>PF13812 PPR_3:  Pentatricopeptide repeat domain
Probab=71.07  E-value=2.6  Score=27.67  Aligned_cols=20  Identities=15%  Similarity=0.235  Sum_probs=17.3

Q ss_pred             HHHHHHHHHHHHHHHcCCee
Q 012032           18 QQLVLQKKILVRIYELGMNP   37 (472)
Q Consensus        18 ~q~~LQkkIl~RmrelGM~P   37 (472)
                      .+.+...++++.|++.|++|
T Consensus        15 g~~~~a~~~~~~M~~~gv~P   34 (34)
T PF13812_consen   15 GDPDAALQLFDEMKEQGVKP   34 (34)
T ss_pred             CCHHHHHHHHHHHHHhCCCC
Confidence            35677899999999999998


No 9  
>PF01229 Glyco_hydro_39:  Glycosyl hydrolases family 39;  InterPro: IPR000514 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. Glycoside hydrolase family 39 GH39 from CAZY comprises enzymes with several known activities; alpha-L-iduronidase (3.2.1.76 from EC); beta-xylosidase (3.2.1.37 from EC). The most highly conserved regions in these enzymes are located in their N-terminal sections. These contain a glutamic acid residue which, on the basis of similarities with other families of glycosyl hydrolases [], probably acts as the proton donor in their catalytic mechanism.; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 2BS9_D 2BFG_E 1W91_B 1UHV_D 1PX8_A.
Probab=67.73  E-value=8.3  Score=42.03  Aligned_cols=103  Identities=20%  Similarity=0.326  Sum_probs=63.4

Q ss_pred             HHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCC-hHHHHHHHHHHHHHHHHhC
Q 012032           25 KILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATD-PLFIEIGRAFIEQQLKEYG  103 (472)
Q Consensus        25 kIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~D-plF~~I~~~F~~eq~~~fG  103 (472)
                      +|++-+.+.||+|.+-  -|+-|.++.....     ..-.|.+          .+-.|.| ..+.++-+.|.++.++.||
T Consensus        88 ~i~D~l~~~g~~P~ve--l~f~p~~~~~~~~-----~~~~~~~----------~~~pp~~~~~W~~lv~~~~~h~~~RYG  150 (486)
T PF01229_consen   88 QILDFLLENGLKPFVE--LGFMPMALASGYQ-----TVFWYKG----------NISPPKDYEKWRDLVRAFARHYIDRYG  150 (486)
T ss_dssp             HHHHHHHHCT-EEEEE--E-SB-GGGBSS-------EETTTTE----------E-S-BS-HHHHHHHHHHHHHHHHHHHH
T ss_pred             HHHHHHHHcCCEEEEE--EEechhhhcCCCC-----ccccccC----------CcCCcccHHHHHHHHHHHHHHHHhhcC
Confidence            6999999999999654  3778887754211     1111222          1244666 6789999999999999999


Q ss_pred             CCCcc--cccCCCCCCCCC-CC---C-hHHHHHHHHHHHHHHhccCCCce
Q 012032          104 RTSHI--YNCDTFDENTPP-VD---S-PEYISSLGAAIYSGMQSGDSDAV  146 (472)
Q Consensus       104 ~~~h~--Y~~D~FnE~~pp-~~---d-p~~L~~~~~~iy~am~~~dP~Av  146 (472)
                      . ..+  ......||-.-. .+   . .+|+ +.-++++++++++||.+.
T Consensus       151 ~-~ev~~W~fEiWNEPd~~~f~~~~~~~ey~-~ly~~~~~~iK~~~p~~~  198 (486)
T PF01229_consen  151 I-EEVSTWYFEIWNEPDLKDFWWDGTPEEYF-ELYDATARAIKAVDPELK  198 (486)
T ss_dssp             H-HHHTTSEEEESS-TTSTTTSGGG-HHHHH-HHHHHHHHHHHHH-TTSE
T ss_pred             C-ccccceeEEeCcCCCcccccCCCCHHHHH-HHHHHHHHHHHHhCCCCc
Confidence            3 222  234567884432 11   1 2455 577889999999999986


No 10 
>cd06592 GH31_glucosidase_KIAA1161 KIAA1161 is an uncharacterized Homo sapiens protein with a glycosyl hydrolase family 31 (GH31) domain that is homologous to the Escherichia coli YihQ glucosidase. Orthologs of KIA1161 are found in eukaryotes and prokaryotes. In bacteria, YihQ (along with YihO) is important for bacterial O-antigen capsule assembly and translocation. Enzymes of the GH31 family possess a wide range of different hydrolytic activities including alpha-glucosidase (glucoamylase and sucrase-isomaltase), alpha-xylosidase, 6-alpha-glucosyltransferase, 3-alpha-isomaltosyltransferase and alpha-1,4-glucan lyase. All GH31 enzymes cleave a terminal carbohydrate moiety from a substrate that varies considerably in size, depending on the enzyme, and may be either a starch or a glycoprotein.
Probab=58.45  E-value=65  Score=32.85  Aligned_cols=82  Identities=18%  Similarity=0.337  Sum_probs=49.4

Q ss_pred             HHHHHHHHHcCCeecc---CCCCCCCchhhHhhC-CCCcee-------ccCC-CCCCCCCCccccccccCCCChHHHHHH
Q 012032           24 KKILVRIYELGMNPVL---PAFSGNVPAALQNVF-PSAKIT-------QLGN-WFSVKSDPRWCCTYLLDATDPLFIEIG   91 (472)
Q Consensus        24 kkIl~RmrelGM~PVL---PgF~G~VP~~~k~~~-P~a~i~-------~~~~-W~gf~~~~~~~~~~~LdP~DplF~~I~   91 (472)
                      +++++++++.||++|+   |.++ .-.+.+++.. .+.=|.       -.+. |.+        .+.++|.+.|.-.+.=
T Consensus        73 ~~mi~~l~~~G~k~~l~i~P~i~-~~s~~~~e~~~~g~~vk~~~g~~~~~~~~w~g--------~~~~~Dftnp~a~~w~  143 (303)
T cd06592          73 KGMIDQLHDLGFRVTLWVHPFIN-TDSENFREAVEKGYLVSEPSGDIPALTRWWNG--------TAAVLDFTNPEAVDWF  143 (303)
T ss_pred             HHHHHHHHHCCCeEEEEECCeeC-CCCHHHHhhhhCCeEEECCCCCCCcccceecC--------CcceEeCCCHHHHHHH
Confidence            6788999999999998   7554 2234455411 111111       1122 332        2568999999866555


Q ss_pred             HHHHHHHHHHhCCCCcccccCCCCCC
Q 012032           92 RAFIEQQLKEYGRTSHIYNCDTFDEN  117 (472)
Q Consensus        92 ~~F~~eq~~~fG~~~h~Y~~D~FnE~  117 (472)
                      ..-+++..+..| . ..+-+| |+|.
T Consensus       144 ~~~~~~~~~~~G-v-dg~w~D-~~E~  166 (303)
T cd06592         144 LSRLKSLQEKYG-I-DSFKFD-AGEA  166 (303)
T ss_pred             HHHHHHHHHHhC-C-cEEEeC-CCCc
Confidence            555555555668 4 567778 6885


No 11 
>PF02065 Melibiase:  Melibiase;  InterPro: IPR000111 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. Glycosyl hydrolase family 27, family 31 and family 36 alpha-galactosidases form the glycosyl hydrolase clan GH-D (acc_GH from CAZY), a superfamily of alpha-galactosidases, alpha-N-acetylgalactosaminidases, and isomaltodextranases which are likely to share a common catalytic mechanism and structural topology. Alpha-galactosidase (3.2.1.22 from EC) (melibiase) [] catalyzes the hydrolysis of melibiose into galactose and glucose. In man, the deficiency of this enzyme is the cause of Fabry's disease (X-linked sphingolipidosis). Alpha-galactosidase is present in a variety of organisms. There is a considerable degree of similarity in the sequence of alpha-galactosidase from various eukaryotic species. Escherichia coli alpha-galactosidase (gene melA), which requires NAD and magnesium as cofactors, is not structurally related to the eukaryotic enzymes; by contrast, an Escherichia coli plasmid encoded alpha-galactosidase (gene rafA P16551 from SWISSPROT) [] contains a region of about 50 amino acids which is similar to a domain of the eukaryotic alpha-galactosidases. Alpha-N-acetylgalactosaminidase (3.2.1.49 from EC) [] catalyzes the hydrolysis of terminal non-reducing N-acetyl-D-galactosamine residues in N-acetyl-alpha-D- galactosaminides. In man, the deficiency of this enzyme is the cause of Schindler and Kanzaki diseases. The sequence of this enzyme is highly related to that of the eukaryotic alpha-galactosidases.; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 1KTC_A 1KTB_A 1UAS_A 3H55_A 3H53_A 3IGU_B 3H54_A 3LRM_A 3LRL_A 3LRK_A ....
Probab=58.12  E-value=29  Score=37.11  Aligned_cols=122  Identities=18%  Similarity=0.110  Sum_probs=68.6

Q ss_pred             HHHHHHHHHcCCeeccCCCCCCCch--hhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHH
Q 012032           24 KKILVRIYELGMNPVLPAFSGNVPA--ALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKE  101 (472)
Q Consensus        24 kkIl~RmrelGM~PVLPgF~G~VP~--~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~  101 (472)
                      +-+.++.+++||+|=|=-=..+|-+  .+.+.+|+--+. .+.......    ..-++||++.|.-.+--...+....+.
T Consensus       107 ~~l~~~i~~~Gmk~GlW~ePe~v~~~S~l~~~hPdw~l~-~~~~~~~~~----r~~~vLD~~~pev~~~l~~~i~~ll~~  181 (394)
T PF02065_consen  107 KPLADYIHSLGMKFGLWFEPEMVSPDSDLYREHPDWVLR-DPGRPPTLG----RNQYVLDLSNPEVRDYLFEVIDRLLRE  181 (394)
T ss_dssp             HHHHHHHHHTT-EEEEEEETTEEESSSCHCCSSBGGBTC-CTTSE-ECB----TTBEEB-TTSHHHHHHHHHHHHHHHHH
T ss_pred             HHHHHHHHHCCCeEEEEeccccccchhHHHHhCccceee-cCCCCCcCc----ccceEEcCCCHHHHHHHHHHHHHHHHh
Confidence            4577788999999976432233321  244556644222 111111110    124689999998776666667777778


Q ss_pred             hCCCCcccccCCCCCCCCCCCCh---H---HHHHHHHHHHHHHhccCCCceEEEeccc
Q 012032          102 YGRTSHIYNCDTFDENTPPVDSP---E---YISSLGAAIYSGMQSGDSDAVWLMQGWL  153 (472)
Q Consensus       102 fG~~~h~Y~~D~FnE~~pp~~dp---~---~L~~~~~~iy~am~~~dP~AvWvmQgW~  153 (472)
                      +| . .|...|--.....+.+..   .   |...+ -.+++.|++..|+-++-.=+|=
T Consensus       182 ~g-i-dYiK~D~n~~~~~~~~~~~~~~~~~~~~~~-y~l~~~L~~~~P~v~iE~CssG  236 (394)
T PF02065_consen  182 WG-I-DYIKWDFNRDITEAGSPSLPEGYHRYVLGL-YRLLDRLRARFPDVLIENCSSG  236 (394)
T ss_dssp             TT---SEEEEE-TS-TTS-SSTTS-GHHHHHHHHH-HHHHHHHHHHTTTSEEEE-BTT
T ss_pred             cC-C-CEEEeccccCCCCCCCCCchHHHHHHHHHH-HHHHHHHHHhCCCcEEEeccCC
Confidence            99 4 788888654444332221   1   22222 3578899999999999887763


No 12 
>PRK00035 hemH ferrochelatase; Reviewed
Probab=53.37  E-value=2.4e+02  Score=29.03  Aligned_cols=152  Identities=20%  Similarity=0.259  Sum_probs=79.7

Q ss_pred             CCCCCHHHHHHHHHHHHHH-----------------------HHHHHHcCC-----eeccCCCCCC----CchhhHhh--
Q 012032            8 GGPLPQSWLDQQLVLQKKI-----------------------LVRIYELGM-----NPVLPAFSGN----VPAALQNV--   53 (472)
Q Consensus         8 gGPLp~~wi~~q~~LQkkI-----------------------l~RmrelGM-----~PVLPgF~G~----VP~~~k~~--   53 (472)
                      |.||...+.++...||++.                       |++|++-|.     .|..|-|+.+    ++..+++.  
T Consensus        69 gSPl~~~t~~q~~~L~~~l~~~~~~~~V~~am~y~~P~i~eal~~l~~~G~~~IivlPL~p~~s~~t~gs~~~~i~~~~~  148 (333)
T PRK00035         69 GSPLNVITRRQAEALQAELAARGPDLPVYLGMRYWNPSIEEALEALKADGVDRIVVLPLYPQYSYSTTASYFEDLARALA  148 (333)
T ss_pred             CChhHHHHHHHHHHHHHHHhccCCCceEEEeecCCCCCHHHHHHHHHhcCCCEEEEEECCCccccccHHHHHHHHHHHHH
Confidence            6799888888888888764                       345655554     4666766644    44444322  


Q ss_pred             -CC-CCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHHhCCC-CcccccCCCCCCCC----CCCChHH
Q 012032           54 -FP-SAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKEYGRT-SHIYNCDTFDENTP----PVDSPEY  126 (472)
Q Consensus        54 -~P-~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~fG~~-~h~Y~~D~FnE~~p----p~~dp~~  126 (472)
                       ++ ..++.-.             +++.   .+|.|.+.-..-+++..+.++.. ...--+=.+|= .|    ..+|| |
T Consensus       149 ~~~~~~~i~~i-------------~~~~---~~p~~i~~l~~~I~~~~~~~~~~~~~~~llfs~HG-~P~~~~~~gd~-Y  210 (333)
T PRK00035        149 KLRLQPEIRFI-------------RSYY---DHPGYIEALAESIREALAKHGEDPEPDRLLFSAHG-LPQRYIDKGDP-Y  210 (333)
T ss_pred             hcCCCCcEEEe-------------CCcc---CCHHHHHHHHHHHHHHHHhcCcccCCcEEEEecCC-CchHHhhcCCC-h
Confidence             22 1111111             1211   46667666666666665555521 01112222331 11    12243 4


Q ss_pred             ---HHHHHHHHHHHHhccCCCceEEEecccCCCCC-CCCchhHHHhHhCCCC---CCEEEE
Q 012032          127 ---ISSLGAAIYSGMQSGDSDAVWLMQGWLFSYDP-FWRPPQMKALLNSVPL---GKLVVL  180 (472)
Q Consensus       127 ---L~~~~~~iy~am~~~dP~AvWvmQgW~F~~~~-fW~~~~~~a~L~~Vp~---~~mliL  180 (472)
                         +...++.|.+.+.-  +..- +..+|+-.-.+ =|..|.+...|..+.+   .+++|+
T Consensus       211 ~~~~~~t~~~l~~~l~~--~~~~-~~~~fqs~~g~~~Wl~P~~~~~l~~l~~~g~k~V~v~  268 (333)
T PRK00035        211 QQQCEETARLLAEALGL--PDED-YDLTYQSRFGPEPWLEPYTDDTLEELAEKGVKKVVVV  268 (333)
T ss_pred             HHHHHHHHHHHHHHhCC--CCCC-eEEEeeCCCCCCccCCCCHHHHHHHHHHcCCCeEEEE
Confidence               44456666666542  1222 34577743333 5998999999888765   356664


No 13 
>PF05119 Terminase_4:  Phage terminase, small subunit;  InterPro: IPR006448 This entry is represented by Streptococcus phage 7201, Orf21. The characteristics of the protein distribution suggest prophage matches in addition to the phage matches. This group of sequences describe the distinct family of phage (and integrated prophage) putative terminase small subunit sequnces. Members tend to be encoded by the gene adjacent to the phage terminase large subunit gene.
Probab=47.68  E-value=21  Score=29.78  Aligned_cols=31  Identities=23%  Similarity=0.277  Sum_probs=28.3

Q ss_pred             CCCCCHHHHHHHHHHHHHHHHHHHHcCCeec
Q 012032            8 GGPLPQSWLDQQLVLQKKILVRIYELGMNPV   38 (472)
Q Consensus         8 gGPLp~~wi~~q~~LQkkIl~RmrelGM~PV   38 (472)
                      |++.+.-++..+....++|.+=.++||++|.
T Consensus        55 G~~~~nP~~~~~~~~~~~~~~l~~~lGLtP~   85 (100)
T PF05119_consen   55 GNPKKNPAVSILNKAMKQMRSLASELGLTPA   85 (100)
T ss_pred             CCcccCHHHHHHHHHHHHHHHHHHHcCCCHH
Confidence            4688888999999999999999999999994


No 14 
>PRK15014 6-phospho-beta-glucosidase BglA; Provisional
Probab=44.94  E-value=47  Score=36.39  Aligned_cols=68  Identities=19%  Similarity=0.457  Sum_probs=45.8

Q ss_pred             HHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHHh
Q 012032           23 QKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKEY  102 (472)
Q Consensus        23 QkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~f  102 (472)
                      =+++++.+++.||+|++==|-=-.|..|.++|        |.|.+         +.+.    ..|.+    |-+..-+.|
T Consensus       112 Y~~lid~l~~~GI~P~vTL~H~dlP~~L~~~y--------GGW~n---------~~~~----~~F~~----Ya~~~f~~f  166 (477)
T PRK15014        112 YDDMFDELLKYNIEPVITLSHFEMPLHLVQQY--------GSWTN---------RKVV----DFFVR----FAEVVFERY  166 (477)
T ss_pred             HHHHHHHHHHcCCEEEEEeeCCCCCHHHHHhc--------CCCCC---------hHHH----HHHHH----HHHHHHHHh
Confidence            47999999999999998776667888887655        55632         1111    34544    444555778


Q ss_pred             CCCCcccccCCCCCC
Q 012032          103 GRTSHIYNCDTFDEN  117 (472)
Q Consensus       103 G~~~h~Y~~D~FnE~  117 (472)
                      |+--+++.  +|||.
T Consensus       167 gdrVk~Wi--T~NEp  179 (477)
T PRK15014        167 KHKVKYWM--TFNEI  179 (477)
T ss_pred             cCcCCEEE--EecCc
Confidence            86445554  88884


No 15 
>cd03198 GST_C_CLIC GST_C family, Chloride Intracellular Channel (CLIC) subfamily; composed of CLIC1-5, p64, parchorin, and similar proteins. They are auto-inserting, self-assembling intracellular anion channels involved in a wide variety of functions including regulated secretion, cell division, and apoptosis. They can exist in both water-soluble and membrane-bound states and are found in various vesicles and membranes. Biochemical studies of the C. elegans homolog, EXC-4, show that the membrane localization domain is present in the N-terminal part of the protein. The structure of soluble human CLIC1 reveals that it is monomeric and adopts a fold similar to GSTs, containing an N-terminal domain with a thioredoxin fold and a C-terminal alpha helical domain. Upon oxidation, the N-terminal domain of CLIC1 undergoes a structural change to form a non-covalent dimer stabilized by the formation of an intramolecular disulfide bond between two cysteines that are far apart in the reduced form. T
Probab=44.12  E-value=65  Score=29.30  Aligned_cols=43  Identities=12%  Similarity=0.200  Sum_probs=37.6

Q ss_pred             HHHHHHHHHHHHHHHHhCChHHHHHHHHHHHHHHHHHHHHhcc
Q 012032          397 AKYANELFLNIIEAYQLNDAHGVFQLSRRFLELVEDMDGLLAC  439 (472)
Q Consensus       397 ~n~~~~~~~~~~~Ay~~~d~~~~~~~~~~~l~li~dlD~LL~t  439 (472)
                      ...+..+|.++..-.++++.+.-+...+++++-++-||..|+.
T Consensus         5 ~~~~~~~f~~~~~~~~~~~~~~~e~~~~~l~~~L~~ld~~L~~   47 (134)
T cd03198           5 NTAGEDIFAKFSAYIKNSNPALNENLEKGLLKALKKLDDYLNS   47 (134)
T ss_pred             hhhHHHHHHHHHHHHcCCChhhhHHHHHHHHHHHHHHHHHHcc
Confidence            3467788999998888998888888889999999999999986


No 16 
>cd06594 GH31_glucosidase_YihQ YihQ is a bacterial alpha-glucosidase with a conserved glycosyl hydrolase family 31 (GH31) domain that catalyzes the release of an alpha-glucosyl residue from the non-reducing end of alpha-glucoside substrates such as alpha-glucosyl fluoride. Orthologs of YihQ that have not yet been functionally characterized are present in plants and fungi. YihQ has sequence similarity to other GH31 enzymes such as CtsZ, a 6-alpha-glucosyltransferase from Bacillus globisporus, and YicI, an alpha-xylosidase from Echerichia coli. In bacteria, YihQ (along with YihO) is important for bacterial O-antigen capsule assembly and translocation.
Probab=43.50  E-value=2e+02  Score=29.64  Aligned_cols=105  Identities=15%  Similarity=0.179  Sum_probs=59.2

Q ss_pred             HHHHHHHHHcCCeecc---CCCCCCCchhhHhh--------CCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHH
Q 012032           24 KKILVRIYELGMNPVL---PAFSGNVPAALQNV--------FPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGR   92 (472)
Q Consensus        24 kkIl~RmrelGM~PVL---PgF~G~VP~~~k~~--------~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~   92 (472)
                      +++++.+++.|++.|+   |++.--.+.-+++.        .++-+......|.+.        +.++|-+.|.-    +
T Consensus        74 ~~mi~~Lh~~G~~~~~~i~P~v~~~~~~~y~~~~~~g~~vk~~~g~~~~~~~w~g~--------~~~~Dftnp~a----~  141 (317)
T cd06594          74 DELIEELKARGIRVLTYINPYLADDGPLYYEEAKDAGYLVKDADGSPYLVDFGEFD--------CGVLDLTNPAA----R  141 (317)
T ss_pred             HHHHHHHHHCCCEEEEEecCceecCCchhHHHHHHCCeEEECCCCCeeeeccCCCC--------ceeeecCCHHH----H
Confidence            6889999999999875   55442222212332        122112222334442        35788888874    4


Q ss_pred             HHHHHHHH----HhCCCCcccccCCCCCCCCCCC------Ch-----HHHHHHHHHHHHHHhccCC
Q 012032           93 AFIEQQLK----EYGRTSHIYNCDTFDENTPPVD------SP-----EYISSLGAAIYSGMQSGDS  143 (472)
Q Consensus        93 ~F~~eq~~----~fG~~~h~Y~~D~FnE~~pp~~------dp-----~~L~~~~~~iy~am~~~dP  143 (472)
                      .++.++.+    ..| . ..+-.| |+|..|+..      +.     .|--..++++|++++++.|
T Consensus       142 ~ww~~~~~~~~~~~G-v-dg~w~D-~~E~~p~d~~~~~g~~~~~~hN~y~~~~~~~~~~~~~~~~~  204 (317)
T cd06594         142 DWFKQVIKEMLLDLG-L-SGWMAD-FGEYLPFDAVLHSGEDAATMHNRYPELWAKLNREAVEEAGK  204 (317)
T ss_pred             HHHHHHHHHHhhhcC-C-cEEEec-CCCCCCCcceecCCCCHHHHhhHHHHHHHHHHHHHHHHhcc
Confidence            44444433    246 4 345666 788766522      11     2555678888888888765


No 17 
>COG1422 Predicted membrane protein [Function unknown]
Probab=41.29  E-value=85  Score=30.72  Aligned_cols=63  Identities=22%  Similarity=0.187  Sum_probs=49.5

Q ss_pred             chHHHH-HHHHHHHHHHHHHHHHHHHHHhCChHHHHHHHHHHHHHHHHHHHHhccCCCCChhHH
Q 012032          386 YDLIDL-TRQALAKYANELFLNIIEAYQLNDAHGVFQLSRRFLELVEDMDGLLACHDGFLLGPW  448 (472)
Q Consensus       386 yDLvDv-tRQvL~n~~~~~~~~~~~Ay~~~d~~~~~~~~~~~l~li~dlD~LL~t~~~FlLg~W  448 (472)
                      +=+||- .-+-+.+.++++.+++.+|++++|..+++++-++=.++++|+-++..-+-.=++=.|
T Consensus        66 ~~liD~ekm~~~qk~m~efq~e~~eA~~~~d~~~lkkLq~~qmem~~~Q~elmk~qfkPM~~~~  129 (201)
T COG1422          66 KLLIDQEKMKELQKMMKEFQKEFREAQESGDMKKLKKLQEKQMEMMDDQRELMKMQFKPMLYIS  129 (201)
T ss_pred             HHhccHHHHHHHHHHHHHHHHHHHHHHHhCCHHHHHHHHHHHHHHHHHHHHHHHHhhhhHHHHH
Confidence            344553 235688999999999999999999999999999999999999988765544444333


No 18 
>cd04257 AAK_AK-HSDH AAK_AK-HSDH: Amino Acid Kinase Superfamily (AAK), AK-HSDH; this CD includes the N-terminal catalytic domain of aspartokinase (AK) of the bifunctional enzyme AK - homoserine dehydrogenase (HSDH). These aspartokinases are found in bacteria (E. coli AKI-HSDHI, ThrA  and E. coli AKII-HSDHII, MetL) and higher plants (Z. mays AK-HSDH). AK and HSDH are the first and third enzymes in the biosynthetic pathway of the aspartate family of amino acids. AK catalyzes the phosphorylation of Asp to P-aspartyl phosphate. HSDH catalyzes the NADPH-dependent conversion of Asp 3-semialdehyde to homoserine. ThrA and MetL are involved in threonine and methionine biosynthesis, respectively. In E. coli, ThrA is subject to allosteric regulation by the end product L-threonine and the native enzyme is reported to be tetrameric. As with bacteria, plant AK and HSDH are feedback inhibited by pathway end products. Maize AK-HSDH is a Thr-sensitive 180-kD enzyme. Arabidopsis AK-HSDH is an alanine-act
Probab=40.70  E-value=18  Score=37.04  Aligned_cols=68  Identities=18%  Similarity=0.222  Sum_probs=43.2

Q ss_pred             HHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHHh
Q 012032           23 QKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKEY  102 (472)
Q Consensus        23 QkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~f  102 (472)
                      .++|-+.+.+.|..||.|||.|-=+        +-+++.+| -+|               .|-+...+|..+=-+....+
T Consensus       169 ~~~l~~~~~~~~~v~Vv~Gfig~~~--------~G~~ttlG-RGG---------------SD~~A~~lA~~l~a~~l~i~  224 (294)
T cd04257         169 KERIKAWFSSNGKVIVVTGFIASNP--------QGETTTLG-RNG---------------SDYSAAILAALLDADQVEIW  224 (294)
T ss_pred             HHHHHHHHhcCCCEEEecCcccCCC--------CCCEEECC-CCc---------------hHHHHHHHHHHhCCCEEEEE
Confidence            3445555555699999999988322        33345544 222               34455666766554555567


Q ss_pred             CCCCcccccCCC
Q 012032          103 GRTSHIYNCDTF  114 (472)
Q Consensus       103 G~~~h~Y~~D~F  114 (472)
                      -+++.+|.+||=
T Consensus       225 tdVdGvyt~DP~  236 (294)
T cd04257         225 TDVDGVYSADPR  236 (294)
T ss_pred             eCCCccCCCCCC
Confidence            778889999994


No 19 
>cd06388 PBP1_iGluR_AMPA_GluR4 N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the GluR4 subunit of the AMPA receptor. N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the GluR4 subunit of the AMPA (alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor. The AMPA receptor is a member of the glutamate-receptor ion channels (iGluRs) which are the major mediators of excitatory synaptic transmission in the central nervous system. AMPA receptors are composed of four types of subunits (GluR1, GluR2, GluR3, and GluR4) which combine to form a tetramer and play an important role in mediating the rapid excitatory synaptic current. Furthermore, this N-terminal domain of the iGluRs has homology with LIVBP, a bacterial periplasmic binding protein, as well as with the structurally related glutamate-binding domain of the G-protein-coupled metabotropic receptors (mGluRs).
Probab=40.32  E-value=89  Score=32.65  Aligned_cols=92  Identities=16%  Similarity=0.292  Sum_probs=52.8

Q ss_pred             HHHHHHHHHHHcCCee-----ccCCCCCCCchhhHh-hCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHH
Q 012032           22 LQKKILVRIYELGMNP-----VLPAFSGNVPAALQN-VFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFI   95 (472)
Q Consensus        22 LQkkIl~RmrelGM~P-----VLPgF~G~VP~~~k~-~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~   95 (472)
                      .-..||+..+++||..     ||-+. |+--..+.+ ++-.++|+.      |         .+++|.+|...++-+.|-
T Consensus       190 ~~~~il~qa~~~gm~~~~y~~il~~~-~~~~~~l~~~~~g~~nitg------~---------~~~~~~~~~v~~~~~~~~  253 (371)
T cd06388         190 RLQNILEQIVSVGKHVKGYHYIIANL-GFKDISLERFMHGGANVTG------F---------QLVDFNTPMVTKLMQRWK  253 (371)
T ss_pred             HHHHHHHHHHhcCccccceEEEEccC-ccccccHHHHhccCCceEE------E---------EeecCCChhHHHHHHHHH
Confidence            4478999999999986     77443 443334433 244444432      2         246888888877666664


Q ss_pred             HHHHHHhCCCCcccccCCCCCCCCCCCChHHHHHHHHHHHHHHhc
Q 012032           96 EQQLKEYGRTSHIYNCDTFDENTPPVDSPEYISSLGAAIYSGMQS  140 (472)
Q Consensus        96 ~eq~~~fG~~~h~Y~~D~FnE~~pp~~dp~~L~~~~~~iy~am~~  140 (472)
                      +..++.|      ++.+     .+|.....+.-++...+.+||++
T Consensus       254 ~~~~~~~------~~~~-----~~~~~~aAl~YDaV~l~a~A~~~  287 (371)
T cd06388         254 KLDQREY------PGSE-----SPPKYTSALTYDGVLVMAEAFRN  287 (371)
T ss_pred             hcCcccc------CCCC-----CCccchHHHHHHHHHHHHHHHHH
Confidence            4444343      1222     13444455555666666666664


No 20 
>cd04247 AAK_AK-Hom3 AAK_AK-Hom3: Amino Acid Kinase Superfamily (AAK), AK-Hom3; this CD includes the N-terminal catalytic domain of the aspartokinase HOM3, a monofunctional class enzyme found in Saccharomyces cerevisiae and other related AK domains. Aspartokinase, the first enzyme in the aspartate metabolic pathway, catalyzes the conversion of aspartate and ATP to aspartylphosphate and ADP, and in fungi, is responsible for the production of threonine, isoleucine and methionine. S. cerevisiae has a single aspartokinase isoenzyme type, which is regulated by feedback, allosteric inhibition by L-threonine. Recent studies show that the allosteric transition triggered by binding of threonine to AK involves a large change in the conformation of the native hexameric enzyme that is converted to an inactive one of different shape and substantially smaller hydrodynamic size.
Probab=40.24  E-value=14  Score=38.27  Aligned_cols=70  Identities=23%  Similarity=0.284  Sum_probs=45.1

Q ss_pred             HHHHHHHHHHHHc-CCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHH
Q 012032           21 VLQKKILVRIYEL-GMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQL   99 (472)
Q Consensus        21 ~LQkkIl~Rmrel-GM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~   99 (472)
                      ++++++.+++.++ |-.||.|||.|..|..        +++.+|. +               -+|-....+|+..=-+..
T Consensus       176 ~~~~~~~~~~~~~~~~v~Vv~GFig~~~~G--------~~ttLGR-g---------------GsD~~A~~la~~l~a~~v  231 (306)
T cd04247         176 ELAQVLGEKITACENRVPVVTGFFGNVPGG--------LLSQIGR-G---------------YTDLCAALCAVGLNADEL  231 (306)
T ss_pred             HHHHHHHHHhhccCCceEEeeccEecCCCC--------CeEEeCC-C---------------chHHHHHHHHHHcCCCEE
Confidence            4677777777765 6789999999976652        3444441 1               233344555554433445


Q ss_pred             HHhCCCCcccccCCC
Q 012032          100 KEYGRTSHIYNCDTF  114 (472)
Q Consensus       100 ~~fG~~~h~Y~~D~F  114 (472)
                      ..+.+++.+|++||=
T Consensus       232 ~i~tdVdGvyt~DP~  246 (306)
T cd04247         232 QIWKEVDGIFTADPR  246 (306)
T ss_pred             EEeecCCeeECCCCC
Confidence            567778889999994


No 21 
>PF03705 CheR_N:  CheR methyltransferase, all-alpha domain;  InterPro: IPR022641  CheR proteins are part of the chemotaxis signaling mechanism which methylates the chemotaxis receptor at specific glutamate residues. This entry refers to the N-terminal domain of the CherR-type MCP methyltransferases, which are found in bacteria, archaea and green plants. This entry is found in association with PF01739 from PFAM.  Methyl transfer from the ubiquitous S-adenosyl-L-methionine (AdoMet) to either nitrogen, oxygen or carbon atoms is frequently employed in diverse organisms ranging from bacteria to plants and mammals. The reaction is catalysed by methyltransferases (Mtases) and modifies DNA, RNA, proteins and small molecules, such as catechol for regulatory purposes. The various aspects of the role of DNA methylation in prokaryotic restriction-modification systems and in a number of cellular processes in eukaryotes including gene regulation and differentiation is well documented. Three classes of DNA Mtases transfer the methyl group from AdoMet to the target base to form either N-6-methyladenine, or N-4-methylcytosine, or C-5- methylcytosine. In C-5-cytosine Mtases, ten conserved motifs are arranged in the same order []. Motif I (a glycine-rich or closely related consensus sequence; FAGxGG in M.HhaI []), shared by other AdoMet-Mtases [], is part of the cofactor binding site and motif IV (PCQ) is part of the catalytic site. In contrast, sequence comparison among N-6-adenine and N-4-cytosine Mtases indicated two of the conserved segments [], although more conserved segments may be present. One of them corresponds to motif I in C-5-cytosine Mtases, and the other is named (D/N/S)PP(Y/F). Crystal structures are known for a number of Mtases [, , , ]. The cofactor binding sites are almost identical and the essential catalytic amino acids coincide. The comparable protein folding and the existence of equivalent amino acids in similar secondary and tertiary positions indicate that many (if not all) AdoMet-Mtases have a common catalytic domain structure. This permits tertiary structure prediction of other DNA, RNA, protein, and small-molecule AdoMet-Mtases from their amino acid sequences []. Flagellated bacteria swim towards favourable chemicals and away from deleterious ones. Sensing of chemoeffector gradients involves chemotaxis receptors, transmembrane (TM) proteins that detect stimuli through their periplasmic domains and transduce the signals via their cytoplasmic domains []. Signalling outputs from these receptors are influenced both by the binding of the chemoeffector ligand to their periplasmic domains and by methylation of specific glutamate residues on their cytoplasmic domains. Methylation is catalysed by CheR, an S-adenosylmethionine-dependent methyltransferase [], which reversibly methylates specific glutamate residues within a coiled coil region, to form gamma-glutamyl methyl ester residues [, ]. The structure of the Salmonella typhimurium chemotaxis receptor methyltransferase CheR, bound to S-adenosylhomocysteine, has been determined to a resolution of 2.0 A []. The structure reveals CheR to be a two-domain protein, with a smaller N-terminal helical domain linked via a single polypeptide connection to a larger C-terminal alpha/beta domain. The C-terminal domain has the characteristics of a nucleotide-binding fold, with an insertion of a small anti-parallel beta-sheet subdomain. The S-adenosylhomocysteine-binding site is formed mainly by the large domain, with contributions from residues within the N-terminal domain and the linker region [].; PDB: 1AF7_A 1BC5_A.
Probab=39.64  E-value=24  Score=26.51  Aligned_cols=20  Identities=20%  Similarity=0.529  Sum_probs=15.4

Q ss_pred             HHHHHHHHHHHHHHHcCCee
Q 012032           18 QQLVLQKKILVRIYELGMNP   37 (472)
Q Consensus        18 ~q~~LQkkIl~RmrelGM~P   37 (472)
                      ++.-|++||..||+++|+.-
T Consensus        24 K~~~l~rRl~~rm~~~~~~~   43 (57)
T PF03705_consen   24 KRSLLERRLARRMRALGLPS   43 (57)
T ss_dssp             GHHHHHHHHHHHHHHHT---
T ss_pred             hHHHHHHHHHHHHHHcCCCC
Confidence            46779999999999999864


No 22 
>TIGR00676 fadh2 5,10-methylenetetrahydrofolate reductase, prokaryotic form. This protein is an FAD-containing flavoprotein.
Probab=36.54  E-value=39  Score=34.03  Aligned_cols=25  Identities=16%  Similarity=0.371  Sum_probs=19.4

Q ss_pred             HHHHHHHHcCCe-eccCCCCCCCchh
Q 012032           25 KILVRIYELGMN-PVLPAFSGNVPAA   49 (472)
Q Consensus        25 kIl~RmrelGM~-PVLPgF~G~VP~~   49 (472)
                      +.++++|+.|++ ||+||+.+..-..
T Consensus       174 ~~~~~~~~~gi~~PIi~Gi~p~~s~k  199 (272)
T TIGR00676       174 RFVDRCRAAGIDVPIIPGIMPITNFK  199 (272)
T ss_pred             HHHHHHHHcCCCCCEecccCCcCCHH
Confidence            457889999876 8898888776655


No 23 
>PF10629 DUF2475:  Protein of unknown function (DUF2475);  InterPro: IPR018902  This entry represents both UPF0573 and UPF0605 families. Both these families of proteins have no known function. 
Probab=36.29  E-value=16  Score=29.87  Aligned_cols=12  Identities=33%  Similarity=0.933  Sum_probs=10.1

Q ss_pred             ccCCCCCCCchh
Q 012032           38 VLPAFSGNVPAA   49 (472)
Q Consensus        38 VLPgF~G~VP~~   49 (472)
                      .+|||+||||..
T Consensus         4 ~iPgY~G~vP~~   15 (71)
T PF10629_consen    4 MIPGYTGYVPGY   15 (71)
T ss_pred             CCCCccCcCCcc
Confidence            379999999984


No 24 
>cd06389 PBP1_iGluR_AMPA_GluR2 N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the GluR2 subunit of the AMPA receptor. N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the GluR2 subunit of the AMPA (alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor. The AMPA receptor is a member of the glutamate-receptor ion channels (iGluRs) which are the major mediators of excitatory synaptic transmission in the central nervous system. AMPA receptors are composed of four types of subunits (GluR1, GluR2, GluR3, and GluR4) which combine to form a tetramer and play an important role in mediating the rapid excitatory synaptic current. Furthermore, this N-terminal domain of the iGluRs has homology with LIVBP, a bacterial periplasmic binding protein, as well as with the structurally related glutamate-binding domain of the G-protein-coupled metabotropic receptors (mGluRs).
Probab=35.22  E-value=1.3e+02  Score=31.31  Aligned_cols=59  Identities=19%  Similarity=0.269  Sum_probs=35.7

Q ss_pred             HHHHHHHHHHHHcCCeec----cCCCCCCCchhhHh-hCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHH
Q 012032           21 VLQKKILVRIYELGMNPV----LPAFSGNVPAALQN-VFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAF   94 (472)
Q Consensus        21 ~LQkkIl~RmrelGM~PV----LPgF~G~VP~~~k~-~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F   94 (472)
                      +.-..|++.++++||++-    +-+-.|+--..|.+ ++..++|+      ||         .+.+|.+|...++-+.|
T Consensus       187 ~~~~~il~qa~~~gm~~~~y~~il~~~~~~~~~l~~~~~~~~nit------g~---------~~~~~~~~~v~~f~~~~  250 (370)
T cd06389         187 DKVNDIVDQVITIGKHVKGYHYIIANLGFTDGDLSKIQFGGANVS------GF---------QIVDYDDPLVSKFIQRW  250 (370)
T ss_pred             HHHHHHHHHHHHhCccccceEEEEccCCccccchhhhccCCcceE------EE---------EEecCCCchHHHHHHHH
Confidence            344689999999999876    22333554445532 25555553      22         24678888766555554


No 25 
>PF00232 Glyco_hydro_1:  Glycosyl hydrolase family 1;  InterPro: IPR001360 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. Glycoside hydrolase family 1 GH1 from CAZY comprises enzymes with a number of known activities; beta-glucosidase (3.2.1.21 from EC); beta-galactosidase (3.2.1.23 from EC); 6-phospho-beta-galactosidase (3.2.1.85 from EC); 6-phospho-beta-glucosidase (3.2.1.86 from EC); lactase-phlorizin hydrolase (3.2.1.62 from EC), (3.2.1.108 from EC); beta-mannosidase (3.2.1.25 from EC); myrosinase (3.2.1.147 from EC). ; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 1QVB_A 3AHY_D 2E9L_A 2ZOX_A 2JFE_X 2E9M_A 3FIZ_A 3FIY_A 3CMJ_A 3FJ0_A ....
Probab=34.79  E-value=1.3e+02  Score=32.61  Aligned_cols=68  Identities=19%  Similarity=0.535  Sum_probs=43.6

Q ss_pred             HHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHH
Q 012032           21 VLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLK  100 (472)
Q Consensus        21 ~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~  100 (472)
                      +-=+++++.+++-||+|++-=|---+|..|.+ +        |.|.+         +.++    ..|.+    |-+..-+
T Consensus        99 ~~Y~~~i~~l~~~gi~P~vtL~H~~~P~~l~~-~--------ggw~~---------~~~~----~~F~~----Ya~~~~~  152 (455)
T PF00232_consen   99 DFYRDLIDELLENGIEPIVTLYHFDLPLWLED-Y--------GGWLN---------RETV----DWFAR----YAEFVFE  152 (455)
T ss_dssp             HHHHHHHHHHHHTT-EEEEEEESS--BHHHHH-H--------TGGGS---------THHH----HHHHH----HHHHHHH
T ss_pred             hhhHHHHHHHHhhccceeeeeeecccccceee-c--------ccccC---------HHHH----HHHHH----HHHHHHH
Confidence            33478999999999999999999999999986 3        44532         1111    35544    4455556


Q ss_pred             HhCCCCcccccCCCCC
Q 012032          101 EYGRTSHIYNCDTFDE  116 (472)
Q Consensus       101 ~fG~~~h~Y~~D~FnE  116 (472)
                      .||+--+++.  +|||
T Consensus       153 ~~gd~V~~w~--T~NE  166 (455)
T PF00232_consen  153 RFGDRVKYWI--TFNE  166 (455)
T ss_dssp             HHTTTBSEEE--EEET
T ss_pred             HhCCCcceEE--eccc
Confidence            7886334443  6777


No 26 
>COG5005 Mu-like prophage protein gpG [General function prediction only]
Probab=33.76  E-value=11  Score=34.21  Aligned_cols=42  Identities=17%  Similarity=0.305  Sum_probs=29.5

Q ss_pred             CCEEEE--ecCCCcccccccccCcCCCCceeeeccCCCCccccc
Q 012032          175 GKLVVL--DLFAEVKPIWSTSKQFYGVPYIWCMLHNFAGNIEMY  216 (472)
Q Consensus       175 ~~mliL--DL~~E~~p~W~~t~~f~G~pwIWc~LhNFGGn~gl~  216 (472)
                      |+.|.+  +|....++.+..+..--|+.=-+..+|||||.+|+-
T Consensus        62 Gk~L~~~GrL~~sltt~y~n~~AlvGtne~YaaiHqfGG~~gr~  105 (140)
T COG5005          62 GKILQDSGRLAGSLTTDYGNNTALVGTNEEYAAIHQFGGKTGRP  105 (140)
T ss_pred             CCceeecchhhhcccccCCCceeeecccchhHHHHHhcCcCCCC
Confidence            444444  233334666667777778888899999999999963


No 27 
>PF13410 GST_C_2:  Glutathione S-transferase, C-terminal domain; PDB: 4DEJ_H 3IC8_A 2JL4_A 2V6K_B 3CBU_B 1JLW_B 3F6D_B 3G7I_A 3F63_A 3G7J_B ....
Probab=33.71  E-value=63  Score=24.58  Aligned_cols=28  Identities=21%  Similarity=0.330  Sum_probs=23.2

Q ss_pred             HHHHHHHHHHHHHHHHHHHhccCCCCChh
Q 012032          418 GVFQLSRRFLELVEDMDGLLACHDGFLLG  446 (472)
Q Consensus       418 ~~~~~~~~~l~li~dlD~LL~t~~~FlLg  446 (472)
                      ..+..-+++.+.++.+|..|+.++ |++|
T Consensus         3 ~~~~~~~~~~~~l~~le~~L~~~~-fl~G   30 (69)
T PF13410_consen    3 AVERARAQLEAALDALEDHLADGP-FLFG   30 (69)
T ss_dssp             HHHHHHHHHHHHHHHHHHHHTTSS-BTTB
T ss_pred             HHHHHHHHHHHHHHHHHHHHhhCC-CCCC
Confidence            456677789999999999999998 7766


No 28 
>cd06394 PBP1_iGluR_Kainate_KA1_2 N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the KA1 and KA2 subunits of Kainate receptor. N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the KA1 and KA2 subunits of Kainate receptor. While this N-terminal domain belongs to the periplasmic-binding fold type I superfamily, the glutamate-binding domain of the iGluR is structurally homologous to the periplasmic-binding fold type II. The LIVBP-like domain of iGluRs is thought to play a role in the initial assembly of iGluR subunits, but it is not well understood how this domain is arranged and functions in intact iGluR. There are five types of kainate receptors, GluR5, GluR6, GluR7, KA1, and KA2, which are structurally similar to AMPA and NMDA subunits of ionotropic glutamate receptors. KA1 and KA2 subunits can only form functional receptors with one of the GluR5-7 subunits. Moreover, GluR5-7 can also form functional homomeric receptor channels act
Probab=33.61  E-value=92  Score=32.29  Aligned_cols=67  Identities=16%  Similarity=0.245  Sum_probs=46.1

Q ss_pred             HHHHHHHHHHHHcCCeeccCCCC--CCCch--hhHh-hCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHH
Q 012032           21 VLQKKILVRIYELGMNPVLPAFS--GNVPA--ALQN-VFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFI   95 (472)
Q Consensus        21 ~LQkkIl~RmrelGM~PVLPgF~--G~VP~--~~k~-~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~   95 (472)
                      +...+|++.++++||..-.-+|-  |.-..  .|.+ .+|.++|+.      |.         +.+|++|...++-+.|-
T Consensus       197 ~~a~~il~qa~~lGm~~~~y~~i~T~l~~~~~~L~~~~~~~~niTg------F~---------l~d~~~~~v~~f~~~~~  261 (333)
T cd06394         197 SMSHTILLKASELGMTSAFYKYILTTMDFPLLRLDSIVDDRSNILG------FS---------MFNQSHAFYQEFIRSLN  261 (333)
T ss_pred             HHHHHHHHHHHHcCCCCCceEEEEecCCcccccHHHhhcCCcceEE------EE---------eecCCcHHHHHHHHHHH
Confidence            45678999999999987766666  66666  4433 466666543      32         58899999877777766


Q ss_pred             HHHHHHh
Q 012032           96 EQQLKEY  102 (472)
Q Consensus        96 ~eq~~~f  102 (472)
                      +++.+.+
T Consensus       262 ~~~~~~~  268 (333)
T cd06394         262 QSWRENC  268 (333)
T ss_pred             Hhhhhhc
Confidence            6554433


No 29 
>PF02836 Glyco_hydro_2_C:  Glycosyl hydrolases family 2, TIM barrel domain;  InterPro: IPR006103 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. Glycoside hydrolase family 2 GH2 from CAZY comprises enzymes with several known activities; beta-galactosidase (3.2.1.23 from EC); beta-mannosidase (3.2.1.25 from EC); beta-glucuronidase (3.2.1.31 from EC). These enzymes contain a conserved glutamic acid residue which has been shown [], in Escherichia coli lacZ (P00722 from SWISSPROT), to be the general acid/base catalyst in the active site of the enzyme. Beta-galactosidase from E. coli has a TIM-barrel-like core surrounded by four other largely beta domains [].; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 3CMG_A 3FN9_C 1YQ2_A 3K4D_B 3LPG_B 3LPF_A 3K4A_B 3K46_B 3GM8_A 3DEC_A ....
Probab=33.57  E-value=1.9e+02  Score=29.10  Aligned_cols=106  Identities=15%  Similarity=0.185  Sum_probs=64.4

Q ss_pred             HHHHHHHHHHHcCCeeccCCCCCCCch--hh---HhhCCCCceeccC-----CCCCCCCCCccccccccCCCChHHHHHH
Q 012032           22 LQKKILVRIYELGMNPVLPAFSGNVPA--AL---QNVFPSAKITQLG-----NWFSVKSDPRWCCTYLLDATDPLFIEIG   91 (472)
Q Consensus        22 LQkkIl~RmrelGM~PVLPgF~G~VP~--~~---k~~~P~a~i~~~~-----~W~gf~~~~~~~~~~~LdP~DplF~~I~   91 (472)
                      ..++=+..|+++|+.-|=-   .|.|.  .|   .+++-=.-+.+..     .|..+.       ..-....||.|.+..
T Consensus        37 ~~~~d~~l~k~~G~N~iR~---~h~p~~~~~~~~cD~~GilV~~e~~~~~~~~~~~~~-------~~~~~~~~~~~~~~~  106 (298)
T PF02836_consen   37 AMERDLELMKEMGFNAIRT---HHYPPSPRFYDLCDELGILVWQEIPLEGHGSWQDFG-------NCNYDADDPEFRENA  106 (298)
T ss_dssp             HHHHHHHHHHHTT-SEEEE---TTS--SHHHHHHHHHHT-EEEEE-S-BSCTSSSSTS-------CTSCTTTSGGHHHHH
T ss_pred             HHHHHHHHHHhcCcceEEc---ccccCcHHHHHHHhhcCCEEEEeccccccCccccCC-------ccccCCCCHHHHHHH
Confidence            4445567899999998864   33443  23   3443211122222     233221       002456899999999


Q ss_pred             HHHHHHHHHHhCCCCcccccCCCCCCCCCCCChHHHHHHHHHHHHHHhccCCCc
Q 012032           92 RAFIEQQLKEYGRTSHIYNCDTFDENTPPVDSPEYISSLGAAIYSGMQSGDSDA  145 (472)
Q Consensus        92 ~~F~~eq~~~fG~~~h~Y~~D~FnE~~pp~~dp~~L~~~~~~iy~am~~~dP~A  145 (472)
                      ..=++++.+.+-+--.+-.=-.+||+        .-....+.+++.+++.||..
T Consensus       107 ~~~~~~~v~~~~NHPSIi~W~~gNE~--------~~~~~~~~l~~~~k~~DptR  152 (298)
T PF02836_consen  107 EQELREMVRRDRNHPSIIMWSLGNES--------DYREFLKELYDLVKKLDPTR  152 (298)
T ss_dssp             HHHHHHHHHHHTT-TTEEEEEEEESS--------HHHHHHHHHHHHHHHH-TTS
T ss_pred             HHHHHHHHHcCcCcCchheeecCccC--------ccccchhHHHHHHHhcCCCC
Confidence            99999999988865456555677999        23445578889999999986


No 30 
>cd06565 GH20_GcnA-like Glycosyl hydrolase family 20 (GH20) catalytic domain of N-acetyl-beta-D-glucosaminidase (GcnA, also known as BhsA) and related proteins. GcnA  is an exoglucosidase which cleaves N-acetyl-beta-D-galactosamine (NAG) and N-acetyl-beta-D-galactosamine residues from 4-methylumbelliferylated (4MU) substrates, as well as cleaving NAG from chito-oligosaccharides (i.e. NAG polymers).  In contrast, sulfated forms of the substrate are unable to be cleaved and act instead as mild competitive inhibitors. Additionally, the enzyme is known to be poisoned by several first-row transition metals as well as by mercury.  GcnA forms a homodimer with subunits comprised of three domains, an N-terminal zincin-like domain, this central catalytic GH20 domain, and a C-terminal alpha helical domain.  The GH20 hexosaminidases are thought to act via a catalytic mechanism in which the catalytic nucleophile is not provided by solvent or the enzyme, but by the substrate itself.
Probab=32.90  E-value=5.2e+02  Score=26.33  Aligned_cols=121  Identities=16%  Similarity=0.099  Sum_probs=73.7

Q ss_pred             CCCCCHHHHHHHHHHHHHHHHHHHHcCCeec-cCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChH
Q 012032            8 GGPLPQSWLDQQLVLQKKILVRIYELGMNPV-LPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPL   86 (472)
Q Consensus         8 gGPLp~~wi~~q~~LQkkIl~RmrelGM~PV-LPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~Dpl   86 (472)
                      +|+++++=+       |.|++--+++||+-| .--..||.=..++  +|.  ......+..        .+..|+|.+|.
T Consensus        54 ~~~yT~~ei-------~ei~~yA~~~gI~vIPeid~pGH~~~~l~--~~~--~~~l~~~~~--------~~~~l~~~~~~  114 (301)
T cd06565          54 RGAYTKEEI-------REIDDYAAELGIEVIPLIQTLGHLEFILK--HPE--FRHLREVDD--------PPQTLCPGEPK  114 (301)
T ss_pred             CCCcCHHHH-------HHHHHHHHHcCCEEEecCCCHHHHHHHHh--Ccc--cccccccCC--------CCCccCCCChh
Confidence            566654433       568899999999866 3334688766554  342  222222221        13469999999


Q ss_pred             HHHHHHHHHHHHHHHhCCCCcc-cccCCCCCCCCCC--------CChHHHHHHHHHHHHHHhccCCCceEEEecc
Q 012032           87 FIEIGRAFIEQQLKEYGRTSHI-YNCDTFDENTPPV--------DSPEYISSLGAAIYSGMQSGDSDAVWLMQGW  152 (472)
Q Consensus        87 F~~I~~~F~~eq~~~fG~~~h~-Y~~D~FnE~~pp~--------~dp~~L~~~~~~iy~am~~~dP~AvWvmQgW  152 (472)
                      =-++-+..++|..+.|.. ..| -++|=+++.+...        +..+-....-+.|.+-+++..+    -++.|
T Consensus       115 t~~fi~~li~ev~~~f~s-~~~HIG~DE~~~~g~~~~~~~~~~~~~~~l~~~~~~~v~~~v~~~g~----~~~~W  184 (301)
T cd06565         115 TYDFIEEMIRQVLELHPS-KYIHIGMDEAYDLGRGRSLRKHGNLGRGELYLEHLKKVLKIIKKRGP----KPMMW  184 (301)
T ss_pred             HHHHHHHHHHHHHHhCCC-CeEEECCCcccccCCCHHHHHhcCCCHHHHHHHHHHHHHHHHHHcCC----EEEEE
Confidence            999999999999999973 233 6788777643211        1111122333456666666666    35567


No 31 
>TIGR00825 EIIBC-GUT PTS system, glucitol/sorbitol-specific, IIBC component. Bacterial PTS transporters transport and concomitantly phosphorylate their sugar substrates, and typically consist of multiple subunits or protein domains. The Gut family consists only of glucitol-specific permeases, but these occur both in Gram-negative and Gram-positive bacteria.E. coli consists of IIA protein, a IIC protein and a IIBC protein. This family is specific for the IIBC component.
Probab=31.69  E-value=9.3  Score=39.26  Aligned_cols=42  Identities=29%  Similarity=0.616  Sum_probs=29.8

Q ss_pred             CCCCCCCCHHHHHHHHHH------HHHHHHHHHHc-CCeeccCCCCCCCc
Q 012032            5 HGWGGPLPQSWLDQQLVL------QKKILVRIYEL-GMNPVLPAFSGNVP   47 (472)
Q Consensus         5 ~gwgGPLp~~wi~~q~~L------QkkIl~Rmrel-GM~PVLPgF~G~VP   47 (472)
                      .||||||--.=-+.++.+      -+-|.+|+.|| ||++| -||.-.||
T Consensus        10 gGwGGPL~i~pt~~KKivyiTGG~~p~i~~kiaeLTG~eaV-dGFk~~~p   58 (331)
T TIGR00825        10 GGWGGPLTVKPTEGKKIVYITAGTEPAIVDKLVNLTGWKAV-DGFKTGEP   58 (331)
T ss_pred             CCcCCCEEEecCCCeEEEEEcCCCCCHHHHHHHHhhCCeec-ccccCCCC
Confidence            469999854322333111      45688999999 99987 79999998


No 32 
>PF12876 Cellulase-like:  Sugar-binding cellulase-like;  InterPro: IPR024778 O-Glycosyl hydrolases 3.2.1. from EC are a widespread group of enzymes that hydrolyse the glycosidic bond between two or more carbohydrates, or between a carbohydrate and a non-carbohydrate moiety. A classification system for glycosyl hydrolases, based on sequence similarity, has led to the definition of 85 different families [, ]. This classification is available on the CAZy (CArbohydrate-Active EnZymes) web site. This entry represents a family of putative cellulase enzymes.; PDB: 3GYC_B.
Probab=31.40  E-value=72  Score=26.32  Aligned_cols=53  Identities=13%  Similarity=0.232  Sum_probs=34.6

Q ss_pred             HHHhCCCCcccccCCCCC-CCCC-----C--C--ChHHHHHHHHHHHHHHhccCCCceEEEecc
Q 012032           99 LKEYGRTSHIYNCDTFDE-NTPP-----V--D--SPEYISSLGAAIYSGMQSGDSDAVWLMQGW  152 (472)
Q Consensus        99 ~~~fG~~~h~Y~~D~FnE-~~pp-----~--~--dp~~L~~~~~~iy~am~~~dP~AvWvmQgW  152 (472)
                      .+.||..+.+=.-|.+|| -.+.     .  .  ..+.+...-+.+.+.|+++||..- |.-|.
T Consensus         2 v~~~~~~~~Il~Wdl~NE~p~~~~~~~~~~~~~~~~~~~~~~l~~~~~~iR~~dP~~p-vt~g~   64 (88)
T PF12876_consen    2 VTRFGYDPRILAWDLWNEPPNNWADGYPAEWGDPKAEAYAEWLKEAFRWIRAVDPSQP-VTSGF   64 (88)
T ss_dssp             HHHTT-GGGEEEEESSTTTT-TT-TT-TT-TT-TTSHHHHHHHHHHHHHHHTT-TTS--EE--B
T ss_pred             chhhcCCCCEEEEEeecCCCCcccccccccccchhHHHHHHHHHHHHHHHHHhCCCCc-EEeec
Confidence            356888788999999999 3311     1  1  235677778899999999999887 54453


No 33 
>PF08858 IDEAL:  IDEAL domain;  InterPro: IPR014957 This entry represents the C-terminal domain of Bacteriophage SPP1, p90. The characteristics of the protein distribution suggest prophage matches in addition to the phage matches. his domain may also be referred to as the IDEAL domain, after the sequence of the most conserved region of the domain.; PDB: 3DO9_A.
Probab=29.10  E-value=2.1e+02  Score=20.27  Aligned_cols=33  Identities=24%  Similarity=0.249  Sum_probs=25.2

Q ss_pred             HHHH-HHHHHHHHHHHHHHHhCChHHHHHHHHHH
Q 012032          394 QALA-KYANELFLNIIEAYQLNDAHGVFQLSRRF  426 (472)
Q Consensus       394 QvL~-n~~~~~~~~~~~Ay~~~d~~~~~~~~~~~  426 (472)
                      +++. ..-..++.+|..|-..+|.+.|..++.++
T Consensus         4 ~~~~~~~~~~L~~~ID~ALd~~D~e~F~~Ls~eL   37 (37)
T PF08858_consen    4 ESLREFRKEQLLELIDEALDNRDKEWFYELSEEL   37 (37)
T ss_dssp             HHHHHHHHHHHHHHHHHHHHTT-HHHHHHHHHHH
T ss_pred             HHHHHHHHHHHHHHHHHHHHcCCHHHHHHHHhhC
Confidence            3444 34467899999999999999999998764


No 34 
>TIGR03356 BGL beta-galactosidase.
Probab=29.03  E-value=3.4e+02  Score=29.26  Aligned_cols=99  Identities=16%  Similarity=0.305  Sum_probs=60.1

Q ss_pred             HHHHHHHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHH
Q 012032           19 QLVLQKKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQ   98 (472)
Q Consensus        19 q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq   98 (472)
                      ..+-=+++++.+++.||+||+==+-=-+|..|.+        . |.|.+         +.+.    ..|.+.++...   
T Consensus        92 ~~~~y~~~i~~l~~~gi~pivtL~Hfd~P~~l~~--------~-gGw~~---------~~~~----~~f~~ya~~~~---  146 (427)
T TIGR03356        92 GLDFYDRLVDELLEAGIEPFVTLYHWDLPQALED--------R-GGWLN---------RDTA----EWFAEYAAVVA---  146 (427)
T ss_pred             HHHHHHHHHHHHHHcCCeeEEeeccCCccHHHHh--------c-CCCCC---------hHHH----HHHHHHHHHHH---
Confidence            3445568999999999999966655456776653        1 33532         1111    46666665554   


Q ss_pred             HHHhCCCCcccccCCCCCCC------------CCCC-Ch--H-----HHHHHHHHHHHHHhccCCCc
Q 012032           99 LKEYGRTSHIYNCDTFDENT------------PPVD-SP--E-----YISSLGAAIYSGMQSGDSDA  145 (472)
Q Consensus        99 ~~~fG~~~h~Y~~D~FnE~~------------pp~~-dp--~-----~L~~~~~~iy~am~~~dP~A  145 (472)
                       +.||+--.++  =+|||..            ||.. +.  .     .+-.+..++++.+++.+|++
T Consensus       147 -~~~~d~v~~w--~t~NEp~~~~~~~y~~G~~~P~~~~~~~~~~~~hnll~Aha~A~~~~~~~~~~~  210 (427)
T TIGR03356       147 -ERLGDRVKHW--ITLNEPWCSAFLGYGLGVHAPGLRDLRAALQAAHHLLLAHGLAVQALRANGPGA  210 (427)
T ss_pred             -HHhCCcCCEE--EEecCcceecccchhhccCCCCCccHHHHHHHHHHHHHHHHHHHHHHHHhCCCC
Confidence             7788644555  4999953            3321 11  1     23345566778888888864


No 35 
>cd03309 CmuC_like CmuC_like. Proteins similar to the putative corrinoid methyltransferase CmuC. Its function has been inferred from sequence similarity to the methyltransferases CmuA and MtaA. Mutants of Methylobacterium sp. disrupted in cmuC and purU appear deficient in some step of chloromethane metabolism.
Probab=28.35  E-value=2.9e+02  Score=28.65  Aligned_cols=115  Identities=15%  Similarity=0.125  Sum_probs=58.2

Q ss_pred             CCCCCC----CCHHHHHHH-HHHHHHHHHHHHHcCCeeccCCCCCCC-----------------------chhhHhhCCC
Q 012032            5 HGWGGP----LPQSWLDQQ-LVLQKKILVRIYELGMNPVLPAFSGNV-----------------------PAALQNVFPS   56 (472)
Q Consensus         5 ~gwgGP----Lp~~wi~~q-~~LQkkIl~RmrelGM~PVLPgF~G~V-----------------------P~~~k~~~P~   56 (472)
                      +.|+|=    ||.+.+++= .---|||++..++.|-.||.-=.+|..                       +..+++++++
T Consensus       177 Ddwa~~~~~~LSpe~f~efv~P~~krIi~~ik~~~g~piilH~cG~~~~~l~~~~e~g~dvl~~d~~~~dl~eak~~~g~  256 (321)
T cd03309         177 DDLGSQKGSFISPATFREFILPRMQRIFDFLRSNTSALIVHHSCGAAASLVPSMAEMGVDSWNVVMTANNTAELRRLLGD  256 (321)
T ss_pred             CCCccccCCccCHHHHHHHHHHHHHHHHHHHHhccCCceEEEeCCCcHHHHHHHHHcCCCEEEecCCCCCHHHHHHHhCC
Confidence            457665    777777543 234577777777774345444445543                       1223334432


Q ss_pred             CceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHHhCC-CCcccccCCCCCCCCCCCChHHHHHHHHHHH
Q 012032           57 AKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKEYGR-TSHIYNCDTFDENTPPVDSPEYISSLGAAIY  135 (472)
Q Consensus        57 a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~fG~-~~h~Y~~D~FnE~~pp~~dp~~L~~~~~~iy  135 (472)
                       +++=+|   +++       |.+|.....  .+.-+...++..+.+|. ..|++.-+  .+. |.+..|+.+..+++.++
T Consensus       257 -k~~l~G---NlD-------p~~L~~~~t--~E~i~~~v~~~l~~~g~~~~fIf~~~--~~~-~~~~~~~~~~~~~~~~~  320 (321)
T cd03309         257 -KVVLAG---AID-------DVALDTATW--PEEDARGVAKAAAECAPIHPFISAPT--AGL-PFSIFPEVLRRVSAFLD  320 (321)
T ss_pred             -CeEEEc---CCC-------hHHhcCCCC--HHHHHHHHHHHHHHhCCCCCEEeCcc--CCC-CcccCHHHHHHHHHhhc
Confidence             122222   222       223332221  35556777777788874 33444433  443 22334667777777654


No 36 
>TIGR00756 PPR pentatricopeptide repeat domain (PPR motif). This family has a similar consensus to the TPR domain (tetratricopeptide), pfam pfam00515, a 33-residue repeat. It is predicted to form a pair of antiparallel helices similar to that of TPR.
Probab=26.32  E-value=44  Score=21.20  Aligned_cols=19  Identities=16%  Similarity=0.396  Sum_probs=15.8

Q ss_pred             HHHHHHHHHHHHHHcCCee
Q 012032           19 QLVLQKKILVRIYELGMNP   37 (472)
Q Consensus        19 q~~LQkkIl~RmrelGM~P   37 (472)
                      +.+-..+++++|++.|+.|
T Consensus        15 ~~~~a~~~~~~M~~~g~~p   33 (35)
T TIGR00756        15 RVEEALELFKEMLERGIEP   33 (35)
T ss_pred             CHHHHHHHHHHHHHcCCCC
Confidence            3456789999999999887


No 37 
>COG0276 HemH Protoheme ferro-lyase (ferrochelatase) [Coenzyme metabolism]
Probab=25.93  E-value=7.9e+02  Score=25.81  Aligned_cols=158  Identities=17%  Similarity=0.226  Sum_probs=91.1

Q ss_pred             CCCCHHHHHHHHHHHHH---------------------HHHHHHHcCC-----eeccCCCCCCCchh-hHhhCCCCceec
Q 012032            9 GPLPQSWLDQQLVLQKK---------------------ILVRIYELGM-----NPVLPAFSGNVPAA-LQNVFPSAKITQ   61 (472)
Q Consensus         9 GPLp~~wi~~q~~LQkk---------------------Il~RmrelGM-----~PVLPgF~G~VP~~-~k~~~P~a~i~~   61 (472)
                      .||-.---.++.+|+++                     .++.|++-|.     -|..|=|+-.+=.+ +++.. .| +.+
T Consensus        69 sPL~~~T~~q~~~L~~~L~~~~~~V~~amry~~P~i~~~v~~l~~~gv~~iv~~pLyPqyS~sTt~s~~~~~~-~a-l~~  146 (320)
T COG0276          69 SPLNVITRAQAAALEERLDLPDFKVYLAMRYGPPFIEEAVEELKKDGVERIVVLPLYPQYSSSTTGSYVDELA-RA-LKE  146 (320)
T ss_pred             CccHHHHHHHHHHHHHHhCCCCccEEEeecCCCCcHHHHHHHHHHcCCCeEEEEECCcccccccHHHHHHHHH-HH-HHh
Confidence            48866556666777776                     5788999998     46678887777644 22111 11 111


Q ss_pred             cCCCCCCCCCCccccccccC--CCChHHHHHHHHHHHHHHHHhCCCCcccccCCCCCCC---CCCCChHHHHHHHHH---
Q 012032           62 LGNWFSVKSDPRWCCTYLLD--ATDPLFIEIGRAFIEQQLKEYGRTSHIYNCDTFDENT---PPVDSPEYISSLGAA---  133 (472)
Q Consensus        62 ~~~W~gf~~~~~~~~~~~Ld--P~DplF~~I~~~F~~eq~~~fG~~~h~Y~~D~FnE~~---pp~~dp~~L~~~~~~---  133 (472)
                      .+.+-.+         .+++  +++|.|-+.=..=+++-.+.+| -++..-+=.||=.-   ---+|| |...+-++   
T Consensus       147 ~~~~~~i---------~~I~~~~~~p~yI~a~a~~I~~~~~~~~-~~~~~llfSaHglP~~~~~~GDp-Y~~q~~~t~~l  215 (320)
T COG0276         147 LRGQPKI---------STIPDYYDEPLYIEALADSIREKLAKHP-RDDDVLLFSAHGLPKRYIDEGDP-YPQQCQETTRL  215 (320)
T ss_pred             cCCCCce---------EEecCccCChHHHHHHHHHHHHHHHhcC-CCCeEEEEecCCCchhhhhcCCc-hHHHHHHHHHH
Confidence            1111111         1222  4788998888888888888888 34555555566221   011343 55555444   


Q ss_pred             HHHHHhccCCCceEEEecccCCCCCCCCchhHHHhHhCCCCC---CEEEEe
Q 012032          134 IYSGMQSGDSDAVWLMQGWLFSYDPFWRPPQMKALLNSVPLG---KLVVLD  181 (472)
Q Consensus       134 iy~am~~~dP~AvWvmQgW~F~~~~fW~~~~~~a~L~~Vp~~---~mliLD  181 (472)
                      |-+.+.--..+-+..-|+ -|..- =|-.|.+...|..+++.   +++|.=
T Consensus       216 i~e~lg~~~~~~~~~~QS-~~G~~-~WL~P~t~~~l~~L~~~g~k~iiv~p  264 (320)
T COG0276         216 IAEALGLPEEEYDLTFQS-RFGPE-PWLQPYTDDLLEELGEKGVKKIIVVP  264 (320)
T ss_pred             HHHHcCCCchheeEEeec-CCCCC-CCCCCCHHHHHHHHHhcCCCeEEEEC
Confidence            444444333444444566 44433 38888888888887763   777764


No 38 
>PF14433 SUKH-3:  SUKH-3 immunity protein
Probab=25.50  E-value=1.4e+02  Score=26.74  Aligned_cols=46  Identities=17%  Similarity=0.282  Sum_probs=37.6

Q ss_pred             hHHHHHHhhcCCC-CCCCHHHHHHHHHhcccCCCChhHHHHHHHHHhc
Q 012032          252 VVYDLMSEMAFQH-EKVDVKAWINQYSVRRYGRSVPAIQDAWNVLYHT  298 (472)
Q Consensus       252 vvYeL~~d~aW~~-~~id~~~W~~~Ya~rRYG~~~~~~~~AW~iL~~t  298 (472)
                      .|-+++...+|+. ..||++.|.+.|...+|-. .|.+.++|+-+..=
T Consensus         3 ~v~~~L~~aGW~~~R~idi~~~~~~~~~~g~~~-~paa~~fL~efGgL   49 (142)
T PF14433_consen    3 KVIELLRAAGWYEGRKIDISLWEKILEEEGYPV-FPAAVEFLAEFGGL   49 (142)
T ss_pred             HHHHHHHHcCCCCCcccCHHHHHHHHHhcCCCC-CHHHHHHHHHcCCe
Confidence            5778999999985 4679999999999997766 67888899865444


No 39 
>cd03201 GST_C_DHAR GST_C family, Dehydroascorbate Reductase (DHAR) subfamily; composed of plant-specific DHARs, monomeric enzymes catalyzing the reduction of DHA into ascorbic acid (AsA) using glutathione as the reductant. DHAR allows plants to recycle oxidized AsA before it is lost. AsA serves as a cofactor of violaxanthin de-epoxidase in the xanthophyll cycle and as an antioxidant in the detoxification of reactive oxygen species. Because AsA is the major reductant in plants, DHAR serves to regulate their redox state. It has been suggested that a significant portion of DHAR activity is plastidic, acting to reduce the large amounts of ascorbate oxidized during hydrogen peroxide scavenging by ascorbate peroxidase. DHAR contains a conserved cysteine in its active site and in addition to its reductase activity, shows thiol transferase activity similar to glutaredoxins.
Probab=25.30  E-value=2.2e+02  Score=24.78  Aligned_cols=8  Identities=0%  Similarity=0.380  Sum_probs=3.8

Q ss_pred             hhHHHHHH
Q 012032          445 LGPWLESA  452 (472)
Q Consensus       445 Lg~Wl~~A  452 (472)
                      |..|++..
T Consensus        91 l~~w~~rl   98 (121)
T cd03201          91 VKSYMKAL   98 (121)
T ss_pred             HHHHHHHH
Confidence            44444444


No 40 
>cd06599 GH31_glycosidase_Aec37 Glycosyl hydrolase family 31 (GH31) domain of a bacterial protein family represented by Escherichia coli protein Aec37. The gene encoding Aec37 (aec-37) is located within a genomic island (AGI-3) isolated from the extraintestinal avian pathogenic Escherichia coli strain BEN2908. The function of Aec37 and its orthologs is unknown; however, deletion of a region of the genome that includes aec-37 affects the assimilation of seven carbohydrates, decreases growth rate of the strain in minimal medium containing galacturonate or trehalose, and attenuates the virulence of E. coli BEN2908 in chickens.  All GH31 enzymes cleave a terminal carbohydrate moiety from a substrate that varies considerably in size, depending on the enzyme, and may be either a starch or a glycoprotein.
Probab=25.16  E-value=5.1e+02  Score=26.53  Aligned_cols=109  Identities=17%  Similarity=0.159  Sum_probs=55.7

Q ss_pred             HHHHHHHHHHcCCeecc---CCCCCCCchhhHhhCCC--------CceeccC-CCCCCCCCCccccccccCCCChHHHHH
Q 012032           23 QKKILVRIYELGMNPVL---PAFSGNVPAALQNVFPS--------AKITQLG-NWFSVKSDPRWCCTYLLDATDPLFIEI   90 (472)
Q Consensus        23 QkkIl~RmrelGM~PVL---PgF~G~VP~~~k~~~P~--------a~i~~~~-~W~gf~~~~~~~~~~~LdP~DplF~~I   90 (472)
                      =+++++++++.|++-|+   |.....-| .+++.-..        -+....+ -|.|.        ..++|-+.|...+.
T Consensus        75 p~~mi~~L~~~g~k~~~~i~P~i~~~~~-~y~e~~~~g~~v~~~~g~~~~~~~~w~g~--------~~~~Dftnp~a~~w  145 (317)
T cd06599          75 PAAFVAKFHERGIRLAPNIKPGLLQDHP-RYKELKEAGAFIKPPDGREPSIGQFWGGV--------GSFVDFTNPEGREW  145 (317)
T ss_pred             HHHHHHHHHHCCCEEEEEeCCcccCCCH-HHHHHHHCCcEEEcCCCCCcceecccCCC--------eEeecCCChHHHHH
Confidence            45889999999998876   55443222 24332211        1111112 23332        34788888876543


Q ss_pred             HHHHHHHHHHHhCCCCcccccCCCCCCCCCCCC----------------hHHHHHHHHHHHHHHhccCC
Q 012032           91 GRAFIEQQLKEYGRTSHIYNCDTFDENTPPVDS----------------PEYISSLGAAIYSGMQSGDS  143 (472)
Q Consensus        91 ~~~F~~eq~~~fG~~~h~Y~~D~FnE~~pp~~d----------------p~~L~~~~~~iy~am~~~dP  143 (472)
                      =+.-++++....| . ..+-+| ++|......+                +.|--..++++|+++++..|
T Consensus       146 w~~~~~~~~~~~G-v-dg~w~D-~~E~~~~~~~~~~~~~g~~~~~~~~~n~y~~l~~~a~~~~~~~~~~  211 (317)
T cd06599         146 WKEGVKEALLDLG-I-DSTWND-NNEYEIWDDDAVCDGFGKPGTIGELRPVQPNLMARASHEAQAEHYP  211 (317)
T ss_pred             HHHHHHHHHhcCC-C-cEEEec-CCCCccCCCcceecCCCCccchhhcccchHHHHHHHHHHHHHHhCC
Confidence            3333333333346 3 456666 5664211000                12333467788888776655


No 41 
>PRK10203 hypothetical protein; Provisional
Probab=25.08  E-value=64  Score=29.14  Aligned_cols=38  Identities=32%  Similarity=0.496  Sum_probs=27.2

Q ss_pred             CCCCCCCCCCCCH---HHHHHHHHHHHHHHHHHHHcCCeeccCCCCCCCchhhH
Q 012032            1 MSNLHGWGGPLPQ---SWLDQQLVLQKKILVRIYELGMNPVLPAFSGNVPAALQ   51 (472)
Q Consensus         1 MgNi~gwgGPLp~---~wi~~q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k   51 (472)
                      +-||.|.|=|||.   +.+..-..+..||++-             +|+||+.+.
T Consensus        21 fdnLpG~GKPL~~~d~~~~p~e~r~~~rilkn-------------ag~lP~el~   61 (122)
T PRK10203         21 FDNLPGSGEPLILDDDSHVPPELRAGYRLLKN-------------AGCLPPELE   61 (122)
T ss_pred             ccCCCCCCCCCCCccCCCCCHHHHHHHHHHhh-------------CCCCCHHHH
Confidence            3599999999973   3344556677777754             589999874


No 42 
>cd03190 GST_C_ECM4_like GST_C family, ECM4-like subfamily; composed of predominantly uncharacterized and taxonomically diverse proteins with similarity to the translation product of the Saccharomyces cerevisiae gene ECM4.  ECM4, a gene of unknown function, is involved in cell surface biosynthesis and architecture. S. cerevisiae ECM4 mutants show increased amounts of the cell wall hexose, N-acetylglucosamine. More recently, global gene expression analysis shows that ECM4 is upregulated during genotoxic conditions and together with the expression profiles of 18 other genes could potentially differentiate between genotoxic and cytotoxic insults in yeast.
Probab=24.10  E-value=2.8e+02  Score=24.62  Aligned_cols=37  Identities=19%  Similarity=0.285  Sum_probs=26.9

Q ss_pred             HHHHhCChHHHHHHHHHHHHHHHHHHHHhccCCCCChh
Q 012032          409 EAYQLNDAHGVFQLSRRFLELVEDMDGLLACHDGFLLG  446 (472)
Q Consensus       409 ~Ay~~~d~~~~~~~~~~~l~li~dlD~LL~t~~~FlLg  446 (472)
                      ..+..++.+.++....++.+.++.+|..|+.+ .|+.|
T Consensus        25 ~~~~~~~~~~~~~~~~~l~~~l~~LE~~L~~~-~yl~G   61 (142)
T cd03190          25 KAGFATTQEAYDEAVDELFEALDRLEELLSDR-RYLLG   61 (142)
T ss_pred             HHhhccCHHHHHHHHHHHHHHHHHHHHHHccC-CeeeC
Confidence            34445667777888888888999999988765 45555


No 43 
>PF07328 VirD1:  T-DNA border endonuclease VirD1;  InterPro: IPR009933 This family consists of several T-DNA border endonuclease VirD1 proteins, which appear to be found exclusively in Agrobacterium species. Agrobacterium, a plant pathogen, is capable to stably transform the plant cell with a segment of its own DNA called T-DNA (transferred DNA). This process depends, among others, on the specialised bacterial virulence proteins VirD1 and VirD2 that excise the T-DNA from its adjacent sequences. VirD1 is thought to interact with VirD2 in this process [].; GO: 0004519 endonuclease activity
Probab=23.50  E-value=3e+02  Score=25.40  Aligned_cols=53  Identities=25%  Similarity=0.239  Sum_probs=38.6

Q ss_pred             cchHHHHHHHHHHHHHHHHHHHHHHHHHhC---ChHHHHHHHHHHHHHHHHHHHHhcc
Q 012032          385 RYDLIDLTRQALAKYANELFLNIIEAYQLN---DAHGVFQLSRRFLELVEDMDGLLAC  439 (472)
Q Consensus       385 ~yDLvDvtRQvL~n~~~~~~~~~~~Ay~~~---d~~~~~~~~~~~l~li~dlD~LL~t  439 (472)
                      |-+|-|+.|++=+-..+  ..++..||+..   |.+.|...-..|=+.+.++|.||.+
T Consensus        73 r~~l~~il~sIg~la~N--in~i~~Aa~~~~~pd~e~f~aER~~fGk~fA~ld~lLr~  128 (147)
T PF07328_consen   73 RQKLEDILRSIGGLATN--INQILKAANRTPRPDYEAFRAERKAFGKEFADLDALLRS  128 (147)
T ss_pred             HHHHHHHHHHHHHHHHH--HHHHHHHHhcCCCccHHHHHHHHHHHHHHHHHHHHHHHH
Confidence            56677777766443322  34577788744   6788999999999999999998864


No 44 
>cd06380 PBP1_iGluR_AMPA N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the AMPA receptor. N-terminal leucine/isoleucine/valine-binding protein (LIVBP)-like domain of the AMPA (alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor, a member of the glutamate-receptor ion channels (iGluRs). AMPA receptors are the major mediators of excitatory synaptic transmission in the central nervous system.  While this N-terminal domain belongs to the periplasmic-binding fold type I superfamily, the glutamate-binding domain of the iGluR is structurally homologous to the periplasmic-binding fold type II. The LIVBP-like domain of iGluRs is thought to play a role in the initial assembly of iGluR subunits, but it is not well understood how this domain is arranged and functions in intact iGluR.  AMPA receptors consist of four types of subunits (GluR1, GluR2, GluR3, and GluR4) which combine to form a tetramer and play an important roles in mediating the rapid excita
Probab=22.63  E-value=3.7e+02  Score=27.51  Aligned_cols=17  Identities=35%  Similarity=0.421  Sum_probs=14.1

Q ss_pred             HHHHHHHHHHHHcCCee
Q 012032           21 VLQKKILVRIYELGMNP   37 (472)
Q Consensus        21 ~LQkkIl~RmrelGM~P   37 (472)
                      +-..+|++.++++||.+
T Consensus       194 ~~~~~i~~qa~~~gm~~  210 (382)
T cd06380         194 ERLNKILEQIVDVGKNR  210 (382)
T ss_pred             HHHHHHHHHHHHhhhcc
Confidence            44578999999999995


No 45 
>COG4334 Uncharacterized protein conserved in bacteria [Function unknown]
Probab=22.58  E-value=46  Score=29.65  Aligned_cols=39  Identities=23%  Similarity=0.596  Sum_probs=21.0

Q ss_pred             CCCchhHHHhHhCCCCCCEEEEecCCCcccccccccCcCCCC-ceeeeccC
Q 012032          159 FWRPPQMKALLNSVPLGKLVVLDLFAEVKPIWSTSKQFYGVP-YIWCMLHN  208 (472)
Q Consensus       159 fW~~~~~~a~L~~Vp~~~mliLDL~~E~~p~W~~t~~f~G~p-wIWc~LhN  208 (472)
                      .|+-++++.+..        +-||+-.-.|+=.||   ||+| ||||..-+
T Consensus         6 ~w~~D~l~ki~~--------~dDl~IsPfre~grt---ygtptWIW~v~vD   45 (131)
T COG4334           6 IWKIDELKKISK--------IDDLYISPFREEGRT---YGTPTWIWFVYVD   45 (131)
T ss_pred             ccchHHHHhhhh--------ccceEecCccccCcc---cCCccEEEEEEEC
Confidence            355555554332        234544444443443   5665 99998654


No 46 
>cd06597 GH31_transferase_CtsY CtsY (cyclic tetrasaccharide-synthesizing enzyme Y) is a bacterial 3-alpha-isomaltosyltransferase, first identified in  Arthrobacter globiformis, that produces cyclic tetrasaccharides together with a closely related enzyme CtsZ. CtsY and CtsZ both have a glycosyl hydrolase family 31 (GH31) catalytic domain.  All GH31 enzymes cleave a terminal carbohydrate moiety from a substrate that varies considerably in size, depending on the enzyme, and may be either a starch or a glycoprotein.
Probab=21.91  E-value=7.2e+02  Score=25.83  Aligned_cols=85  Identities=22%  Similarity=0.282  Sum_probs=44.3

Q ss_pred             HHHHHHHHHcCCeecc---CCCC-CCCc---------hhhHhh----CCCCceecc-CCCCCCCCCCccccccccCCCCh
Q 012032           24 KKILVRIYELGMNPVL---PAFS-GNVP---------AALQNV----FPSAKITQL-GNWFSVKSDPRWCCTYLLDATDP   85 (472)
Q Consensus        24 kkIl~RmrelGM~PVL---PgF~-G~VP---------~~~k~~----~P~a~i~~~-~~W~gf~~~~~~~~~~~LdP~Dp   85 (472)
                      ++.++++.+.|++-+|   |... +..|         ++.++-    .++-+.... +.|.|-        +.++|-+.|
T Consensus        88 ~~mi~~Lh~~G~kv~l~v~P~i~~~~~~~~~~~~~~~~~~~~g~~vk~~~G~~~~~~~~W~g~--------~~~~Dftnp  159 (340)
T cd06597          88 KGMIDELHEQGVKVLLWQIPIIKLRPHPHGQADNDEDYAVAQNYLVQRGVGKPYRIPGQWFPD--------SLMLDFTNP  159 (340)
T ss_pred             HHHHHHHHHCCCEEEEEecCccccccccccccchhHHHHHHCCEEEEcCCCCccccccccCCC--------ceeecCCCH
Confidence            6888999999999875   6554 2111         111111    111111111 245442        468999998


Q ss_pred             HHHHHHHHHHHHHHHHhCCCCcccccCCCCCCCC
Q 012032           86 LFIEIGRAFIEQQLKEYGRTSHIYNCDTFDENTP  119 (472)
Q Consensus        86 lF~~I~~~F~~eq~~~fG~~~h~Y~~D~FnE~~p  119 (472)
                      .-.+-=+.-+++..+.+| . .-+-+| ++|..+
T Consensus       160 ~a~~Ww~~~~~~~~~~~G-i-dg~w~D-~~E~~~  190 (340)
T cd06597         160 EAAQWWMEKRRYLVDELG-I-DGFKTD-GGEHVW  190 (340)
T ss_pred             HHHHHHHHHHHHHHHhcC-C-cEEEec-CCCccC
Confidence            865532222333333467 4 456677 778543


No 47 
>cd04243 AAK_AK-HSDH-like AAK_AK-HSDH-like: Amino Acid Kinase Superfamily (AAK), AK-HSDH-like; this family includes the N-terminal catalytic domain of aspartokinase (AK) of the bifunctional enzyme AK- homoserine dehydrogenase (HSDH). These aspartokinases are found in such bacteria as E. coli (AKI-HSDHI, ThrA  and  AKII-HSDHII, MetL) and in higher plants (Z. mays AK-HSDH). AK and HSDH are the first and third enzymes in the biosynthetic pathway of the aspartate family of amino acids. AK catalyzes the phosphorylation of Asp to P-aspartyl phosphate. HSDH catalyzes the NADPH-dependent conversion of Asp 3-semialdehyde to homoserine. ThrA and MetL are involved in threonine and methionine biosynthesis, respectively. In E. coli, ThrA is subject to allosteric regulation by the end product L-threonine and the native enzyme is reported to be tetrameric. As with bacteria, plant AK and HSDH are feedback inhibited by pathway end products. Maize AK-HSDH is a Thr-sensitive 180-kD enzyme. Arabidopsis AK-
Probab=21.64  E-value=58  Score=33.34  Aligned_cols=67  Identities=15%  Similarity=0.170  Sum_probs=40.2

Q ss_pred             HHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHHHHHHHHhC
Q 012032           24 KKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFIEQQLKEYG  103 (472)
Q Consensus        24 kkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~~eq~~~fG  103 (472)
                      ++|.+.+.+.|..||.|||-|.=+        .-+++.+| -+|               +|-.-..||+.+=-+....+-
T Consensus       169 ~~~~~~~~~~~~v~Vv~Gfig~~~--------~G~~ttLG-Rgg---------------sD~~A~~~a~~l~a~~~~i~t  224 (293)
T cd04243         169 ERLAQLLAEHGKVVVTQGFIASNE--------DGETTTLG-RGG---------------SDYSAALLAALLDAEEVEIWT  224 (293)
T ss_pred             HHHHHHHhcCCCEEEecCccccCC--------CCCEEEeC-CCC---------------cHHHHHHHHHHcCCCEEEEEe
Confidence            344444444499999999998422        22344444 111               244455566555444445667


Q ss_pred             CCCcccccCCC
Q 012032          104 RTSHIYNCDTF  114 (472)
Q Consensus       104 ~~~h~Y~~D~F  114 (472)
                      +++.+|.+||-
T Consensus       225 dvdGiyt~dP~  235 (293)
T cd04243         225 DVDGVYTADPR  235 (293)
T ss_pred             CCCccCCCCCC
Confidence            77889999993


No 48 
>cd03186 GST_C_SspA GST_N family, Stringent starvation protein A (SspA) subfamily; SspA is a RNA polymerase (RNAP)-associated protein required for the lytic development of phage P1 and for stationary phase-induced acid tolerance of E. coli. It is implicated in survival during nutrient starvation. SspA adopts the GST fold with an N-terminal thioredoxin-fold domain and a C-terminal alpha helical domain, but it does not bind glutathione (GSH) and lacks GST activity. SspA is highly conserved among gram-negative bacteria. Related proteins found in Neisseria (called RegF), Francisella and Vibrio regulate the expression of virulence factors necessary for pathogenesis.
Probab=21.53  E-value=4.4e+02  Score=21.54  Aligned_cols=29  Identities=7%  Similarity=0.177  Sum_probs=19.4

Q ss_pred             HHHHHHHHHHHHHHHHHHHHhccCCCCChh
Q 012032          417 HGVFQLSRRFLELVEDMDGLLACHDGFLLG  446 (472)
Q Consensus       417 ~~~~~~~~~~l~li~dlD~LL~t~~~FlLg  446 (472)
                      +.....-.++.+.+..+|..|+.++ |++|
T Consensus        31 ~~~~~~~~~~~~~l~~le~~L~~~~-~l~G   59 (107)
T cd03186          31 KEAEKARKELRESLLALAPVFAHKP-YFMS   59 (107)
T ss_pred             HHHHHHHHHHHHHHHHHHHHHcCCC-cccC
Confidence            3344555668888888999887543 5554


No 49 
>TIGR01558 sm_term_P27 phage terminase, small subunit, putative, P27 family. Members tend to be adjacent to the phage terminase large subunit gene.
Probab=21.18  E-value=98  Score=27.00  Aligned_cols=31  Identities=23%  Similarity=0.193  Sum_probs=25.1

Q ss_pred             CCCCCHHHHHHHHHHHHHHHHHHHHcCCeec
Q 012032            8 GGPLPQSWLDQQLVLQKKILVRIYELGMNPV   38 (472)
Q Consensus         8 gGPLp~~wi~~q~~LQkkIl~RmrelGM~PV   38 (472)
                      |++.+.-.+.-+.+..++|.+=..+|||+|.
T Consensus        64 g~~k~nPa~~i~~~a~~~~~~l~~elGLtP~   94 (116)
T TIGR01558        64 GSPKANPALTVVEDAFKQLRSIGSALGLTPS   94 (116)
T ss_pred             CCeecChHHHHHHHHHHHHHHHHHHcCCCHH
Confidence            4455556777788889999999999999986


No 50 
>PRK13511 6-phospho-beta-galactosidase; Provisional
Probab=20.68  E-value=2.3e+02  Score=30.96  Aligned_cols=71  Identities=15%  Similarity=0.361  Sum_probs=47.4

Q ss_pred             HHHHHH--HHHHHHHHHcCCeeccCCCCCCCchhhHhhCCCCceeccCCCCCCCCCCccccccccCCCChHHHHHHHHHH
Q 012032           18 QQLVLQ--KKILVRIYELGMNPVLPAFSGNVPAALQNVFPSAKITQLGNWFSVKSDPRWCCTYLLDATDPLFIEIGRAFI   95 (472)
Q Consensus        18 ~q~~LQ--kkIl~RmrelGM~PVLPgF~G~VP~~~k~~~P~a~i~~~~~W~gf~~~~~~~~~~~LdP~DplF~~I~~~F~   95 (472)
                      .++.|+  +++++.|++.||+||+==|-=-+|..|.+         .|.|.+         +.++    ..|.+.|+.- 
T Consensus        89 N~~gl~~Y~~lid~l~~~GI~P~VTL~H~dlP~~L~~---------~GGW~n---------~~~v----~~F~~YA~~~-  145 (469)
T PRK13511         89 NPKGVEYYHRLFAECHKRHVEPFVTLHHFDTPEALHS---------NGDWLN---------RENI----DHFVRYAEFC-  145 (469)
T ss_pred             CHHHHHHHHHHHHHHHHcCCEEEEEecCCCCcHHHHH---------cCCCCC---------HHHH----HHHHHHHHHH-
Confidence            344554  69999999999999987776678999875         245643         1111    4666655554 


Q ss_pred             HHHHHHhCCCCcccccCCCCCC
Q 012032           96 EQQLKEYGRTSHIYNCDTFDEN  117 (472)
Q Consensus        96 ~eq~~~fG~~~h~Y~~D~FnE~  117 (472)
                         -+.||+ -.++.  +|||-
T Consensus       146 ---~~~fgd-Vk~W~--T~NEP  161 (469)
T PRK13511        146 ---FEEFPE-VKYWT--TFNEI  161 (469)
T ss_pred             ---HHHhCC-CCEEE--Eccch
Confidence               567896 45443  78883


No 51 
>PF13041 PPR_2:  PPR repeat family 
Probab=20.50  E-value=74  Score=22.93  Aligned_cols=22  Identities=14%  Similarity=0.215  Sum_probs=18.3

Q ss_pred             HHHHHHHHHHHHHHHcCCeecc
Q 012032           18 QQLVLQKKILVRIYELGMNPVL   39 (472)
Q Consensus        18 ~q~~LQkkIl~RmrelGM~PVL   39 (472)
                      ++.+-..+++++|++.|+.|=.
T Consensus        17 ~~~~~a~~l~~~M~~~g~~P~~   38 (50)
T PF13041_consen   17 GKFEEALKLFKEMKKRGIKPDS   38 (50)
T ss_pred             cCHHHHHHHHHHHHHcCCCCCH
Confidence            4567788999999999999844


No 52 
>PF09350 DUF1992:  Domain of unknown function (DUF1992);  InterPro: IPR018961  This entry represents a family of proteins that may have a role in protein folding or as a chaperone. DnaJ is a member of the J-protein family, which are defined by the presence of a J domain that can regulate the activity of 70kDa heat-shock proteins []. Some of the proteins in this entry contain a J domain.
Probab=20.18  E-value=74  Score=25.84  Aligned_cols=38  Identities=29%  Similarity=0.565  Sum_probs=25.2

Q ss_pred             CCCCCCCCCCCCH----HHHHHHHHHHHHHHHHHHHcCCeeccCCCCCCCchhhH
Q 012032            1 MSNLHGWGGPLPQ----SWLDQQLVLQKKILVRIYELGMNPVLPAFSGNVPAALQ   51 (472)
Q Consensus         1 MgNi~gwgGPLp~----~wi~~q~~LQkkIl~RmrelGM~PVLPgF~G~VP~~~k   51 (472)
                      +-||.|.|=||+.    ..+..-..+..|||++             +|++|+.+.
T Consensus        15 FdnLpG~GKPL~~~~~~~~~~~~~~~~~~iLk~-------------~g~lPp~i~   56 (71)
T PF09350_consen   15 FDNLPGAGKPLPLDDDNPYWPAEERMANRILKN-------------AGYLPPWIE   56 (71)
T ss_pred             ccCCCCCCCCCCCCCCCcCCCHHHHHHHHhhcc-------------cCCCCHHHH
Confidence            3589999999884    2233445566777665             677887763


Done!