Query         018216
Match_columns 359
No_of_seqs    348 out of 1543
Neff          7.1 
Searched_HMMs 46136
Date          Fri Mar 29 07:07:14 2013
Command       hhsearch -i /work/01045/syshi/csienesis_hhblits_a3m/018216.a3m -d /work/01045/syshi/HHdatabase/Cdd.hhm -o /work/01045/syshi/hhsearch_cdd/018216hhsearch_cdd -cpu 12 -v 0 

 No Hit                             Prob E-value P-value  Score    SS Cols Query HMM  Template HMM
  1 KOG1141 Predicted histone meth 100.0 8.5E-49 1.8E-53  397.1   6.4  179   14-231   615-826 (1262)
  2 KOG1082 Histone H3 (Lys9) meth 100.0 1.1E-44 2.4E-49  357.1  16.4  274   55-359    50-358 (364)
  3 KOG4442 Clathrin coat binding  100.0 2.6E-43 5.6E-48  355.6  13.4  165  178-357    92-262 (729)
  4 KOG1080 Histone H3 (Lys4) meth 100.0 3.1E-29 6.8E-34  267.9  10.3  136  204-355   866-1005(1005)
  5 KOG1079 Transcriptional repres 100.0 3.5E-29 7.5E-34  252.1   8.7  174  132-334   521-711 (739)
  6 smart00317 SET SET (Su(var)3-9  99.9 1.6E-21 3.5E-26  159.5  11.9  112  205-331     1-116 (116)
  7 KOG1083 Putative transcription  99.8 2.8E-20   6E-25  194.7   3.9  128  193-334  1166-1296(1306)
  8 KOG1085 Predicted methyltransf  99.8 1.1E-18 2.5E-23  161.2   8.4  124  199-334   251-379 (392)
  9 PF05033 Pre-SET:  Pre-SET moti  99.7 1.2E-17 2.7E-22  136.0   5.0  100   69-196     1-103 (103)
 10 smart00468 PreSET N-terminal t  99.6 4.2E-16 9.1E-21  126.0   6.8   39   67-105     1-40  (98)
 11 COG2940 Proteins containing SE  99.6 2.3E-16 5.1E-21  161.3   3.1  157  184-355   312-480 (480)
 12 PF00856 SET:  SET domain;  Int  99.4 5.2E-13 1.1E-17  113.9   7.4   54  276-332   109-162 (162)
 13 KOG1081 Transcription factor N  99.0 9.8E-11 2.1E-15  118.8   2.5  150  179-356   288-438 (463)
 14 KOG1141 Predicted histone meth  98.9 2.7E-09 5.8E-14  110.8   7.4  298   19-354   823-1261(1262)
 15 KOG2589 Histone tail methylase  98.6 6.3E-08 1.4E-12   93.2   5.2   59  283-347   194-252 (453)
 16 cd01395 HMT_MBD Methyl-CpG bin  98.0 1.2E-06 2.6E-11   64.0   0.4   30   14-43     20-49  (60)
 17 KOG2461 Transcription factor B  97.9 8.9E-06 1.9E-10   81.2   3.6  116  202-334    26-145 (396)
 18 smart00570 AWS associated with  97.1 0.00029 6.3E-09   49.8   1.7   31  168-202    20-50  (51)
 19 smart00508 PostSET Cysteine-ri  97.0 0.00043 9.3E-09   41.9   1.5   17  340-356     2-18  (26)
 20 smart00391 MBD Methyl-CpG bind  96.1  0.0017 3.7E-08   50.1   0.4   31   14-44     24-54  (77)
 21 cd00122 MBD MeCP2, MBD1, MBD2,  95.7  0.0035 7.7E-08   46.2   0.6   29   15-43     22-50  (62)
 22 cd01397 HAT_MBD Methyl-CpG bin  94.9   0.011 2.3E-07   45.1   1.0   29   14-42     21-49  (73)
 23 PF01429 MBD:  Methyl-CpG bindi  94.7   0.011 2.3E-07   45.6   0.5   31   15-45     28-58  (77)
 24 KOG2084 Predicted histone tail  91.3    0.36 7.9E-06   48.7   5.7   40  288-334   208-248 (482)
 25 cd01396 MeCP2_MBD MeCP2, MBD1,  90.1    0.11 2.5E-06   40.0   0.6   29   15-43     23-51  (77)
 26 KOG1337 N-methyltransferase [G  82.1     1.1 2.3E-05   46.2   2.9   41  288-334   239-279 (472)
 27 PF03638 TCR:  Tesmin/TSO1-like  79.0    0.96 2.1E-05   30.7   0.9   37  149-197     4-40  (42)
 28 PF12937 F-box-like:  F-box-lik  61.4     5.4 0.00012   27.0   1.5   31   31-61      7-37  (47)
 29 PF08666 SAF:  SAF domain;  Int  57.5     6.4 0.00014   28.1   1.5   15  314-328     3-17  (63)
 30 PF00646 F-box:  F-box domain;   49.5     9.1  0.0002   25.7   1.1   33   30-62      8-40  (48)
 31 PF11403 Yeast_MT:  Yeast metal  45.0      22 0.00048   22.8   2.2   18  149-171    22-39  (40)
 32 smart00256 FBOX A Receptor for  42.8      18  0.0004   22.9   1.8   31   31-61      4-34  (41)
 33 KOG2155 Tubulin-tyrosine ligas  39.9      17 0.00036   37.1   1.7   49  285-334   204-252 (631)
 34 smart00858 SAF This domain fam  30.5      34 0.00073   24.2   1.6   16  314-329     3-18  (64)
 35 PF01847 VHL:  von Hippel-Linda  24.4      84  0.0018   27.6   3.2   38  285-330    16-53  (156)
 36 smart00317 SET SET (Su(var)3-9  24.0      86  0.0019   24.3   3.1   16  216-231    98-113 (116)
 37 cd05468 pVHL von Hippel-Landau  23.2 1.1E+02  0.0024   26.2   3.7   37  286-330    11-47  (141)
 38 KOG1081 Transcription factor N  22.4      26 0.00056   36.1  -0.3  138  211-358   120-267 (463)
 39 KOG4289 Cadherin EGF LAG seven  21.8 1.5E+02  0.0032   35.1   5.2   78  148-230  1789-1887(2531)
 40 TIGR03569 NeuB_NnaB N-acetylne  20.3      55  0.0012   32.2   1.5   20  312-331   277-296 (329)

No 1  
>KOG1141 consensus Predicted histone methyl transferase [Chromatin structure and dynamics]
Probab=100.00  E-value=8.5e-49  Score=397.07  Aligned_cols=179  Identities=21%  Similarity=0.306  Sum_probs=148.1

Q ss_pred             cccceeeeeeCCchhhhhhHHHHhccCCcchhhhHhhh-------hhhhhhhcccccccccccccCCcCCCCeeeEeCCC
Q 018216           14 HKRTIEVDEEQKQNQFLQWARLILPWLTPGELANVSLT-------CRTLSQISKSITLSRSLDASRSVENFPIPFHNAAD   86 (359)
Q Consensus        14 ~~~~~~~~~~~~~~~~~~~~~~v~~~l~~~~~~~~~~~-------~~~~~~~~~~~~~~~~~DiS~G~E~~pI~~vN~vD   86 (359)
                      +++.+||.|++|||+.||.+++|.+||.++.+.++...       ....+...+.....++.||+.|+|.+||+++|++|
T Consensus       615 ~~~~fhv~yktpcg~~lr~~~el~ryL~et~c~flf~~~f~~~~yV~~~r~~~p~kp~~~~~Di~~g~e~vpis~~neid  694 (1262)
T KOG1141|consen  615 YIRDFHVEYKTPCGMPLRMRIELYRYLVETRCKFLFVIGFDRAFYVVRHRAPNPLKPGNRCTDIPCGREHVPISEKNEID  694 (1262)
T ss_pred             hhhcceeeccCCCccchHHHHHHHHHHHHhcCcEEEEeecccchheeecccCCCcCCcceeccccCCccccccceeeccc
Confidence            78899999999999999999999999999988776654       22334455556678899999999999999999999


Q ss_pred             CCCCCCeEEcccccCCCCCCCCCCCCCccccccccccccccccccccCcccccccccCCCCCCCCCcCCCccCCCCCCCc
Q 018216           87 KTPYAYFIYTPSQIIPPPCPAQFPPRQFWASTNAAADAESNSSLSRLGFDSVSLVCESDESESGCDCEECFEVGLGDGVF  166 (359)
Q Consensus        87 ~~~~~~f~Y~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~c~C~~~~~~~gC~C~~~~~~~~~~~~~  166 (359)
                      +.++|.|.|....|....            ...+  +++++                    ..+|+|..+|.+.     .
T Consensus       695 s~~lpq~ay~K~~ip~~~------------nl~n--~~~~f--------------------l~scdc~~gcid~-----~  735 (1262)
T KOG1141|consen  695 SHRLPQAAYKKHMIPTNN------------NLSN--RRKDF--------------------LQSCDCPTGCIDS-----M  735 (1262)
T ss_pred             CcCCccchhheeeccCCC------------cccc--cChhh--------------------hhcCCCCcchhhh-----h
Confidence            999999999997652111            1111  12222                    2789999999877     8


Q ss_pred             CCCCCCC-------------------------cccccceeecCCCCCCCC-CCCCccccCCCcccEEEEEeCCCCcccee
Q 018216          167 GCPCFSG-------------------------LEDVGIVSECGPSCGCGS-ECGNRLTQRGISVRLKIVRSVNKGWGLYA  220 (359)
Q Consensus       167 ~c~C~~~-------------------------~~~~~~i~EC~~~C~C~~-~C~Nrv~Q~g~~~~leV~~t~~kGwGlrA  220 (359)
                      +|+|.++                         ...+..+|||+..|+|.+ .|.||++|+|++++|++|++..||||+|+
T Consensus       736 kcachQltvk~~~t~p~~~v~~t~gykyKRl~e~~ptg~yEc~k~ckc~~~~C~nrmvqhg~qvRlq~fkt~~kGWg~rc  815 (1262)
T KOG1141|consen  736 KCACHQLTVKKKTTGPNQNVASTNGYKYKRLIEIRPTGPYECLKACKCCGPDCLNRMVQHGYQVRLQRFKTIHKGWGRRC  815 (1262)
T ss_pred             hhhHHHHHHHhhccCCCcccccCcchhhHHHHHhcCCCHHHHHHhhccCcHHHHHHHhhcCceeEeeeccccccccceEe
Confidence            8999864                         012466899999999876 79999999999999999999999999999


Q ss_pred             ccccCCCcEEE
Q 018216          221 DQFIKQGQFIC  231 (359)
Q Consensus       221 ~~~I~~G~fI~  231 (359)
                      +++|.+|.|||
T Consensus       816 lddi~~g~fVc  826 (1262)
T KOG1141|consen  816 LDDITGGNFVC  826 (1262)
T ss_pred             eeecCCceEEE
Confidence            99999999999


No 2  
>KOG1082 consensus Histone H3 (Lys9) methyltransferase SUV39H1/Clr4, required for transcriptional silencing [Chromatin structure and dynamics; Transcription]
Probab=100.00  E-value=1.1e-44  Score=357.07  Aligned_cols=274  Identities=32%  Similarity=0.498  Sum_probs=190.6

Q ss_pred             hhhhcccccccccccccCCcCCCCeeeEeCCCCCCCCCeEEcccccCCCCCCCCCCCCCccccccccccccccccccccC
Q 018216           55 LSQISKSITLSRSLDASRSVENFPIPFHNAADKTPYAYFIYTPSQIIPPPCPAQFPPRQFWASTNAAADAESNSSLSRLG  134 (359)
Q Consensus        55 ~~~~~~~~~~~~~~DiS~G~E~~pI~~vN~vD~~~~~~f~Y~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  134 (359)
                      .....+..+.....||+.|.|++||+++|+||...+++|.|++..+++.+...     .+           .        
T Consensus        50 ~~~~~~~~~~~~~~d~~~~~e~~~v~~~n~id~~~~~~f~y~~~~~~~~~~~~-----~~-----------~--------  105 (364)
T KOG1082|consen   50 ANDKDKLEAKSELEDIALGSENLPVPLVNRIDEDAPLYFQYIATEIVDPGELS-----DC-----------E--------  105 (364)
T ss_pred             cccccccccccccccccCccccCceeeeeeccCCccccceeccccccCccccc-----cC-----------c--------
Confidence            33344444566789999999999999999999876689999999876553100     00           0        


Q ss_pred             cccccccccCCCCCC---CCCcCCCccCCCCCCCcCCCCC---CCcccccceeecCCCCCCCCCCCCccccCCCcccEEE
Q 018216          135 FDSVSLVCESDESES---GCDCEECFEVGLGDGVFGCPCF---SGLEDVGIVSECGPSCGCGSECGNRLTQRGISVRLKI  208 (359)
Q Consensus       135 ~~~~~c~C~~~~~~~---gC~C~~~~~~~~~~~~~~c~C~---~~~~~~~~i~EC~~~C~C~~~C~Nrv~Q~g~~~~leV  208 (359)
                       ....|.|...+...   .|.|.....    ....+..+.   ........+|||++.|+|++.|.||++|+|++.+|+|
T Consensus       106 -~~~~c~C~~~~~~~~~~~C~C~~~n~----~~~~~~~~~~~~~~~~~~~~i~EC~~~C~C~~~C~nRv~q~g~~~~leI  180 (364)
T KOG1082|consen  106 -NSTGCRCCSSCSSVLPLTCLCERHNG----GLVAYTCDGDCGTLGKFKEPVFECSVACGCHPDCANRVVQKGLQFHLEV  180 (364)
T ss_pred             -cccCCCccCCCCCCCCccccChHhhC----CccccccCCccccccccCccccccccCCCCCCcCcchhhccccccceEE
Confidence             00113332222111   266655431    111222221   2233456899999999999999999999999999999


Q ss_pred             EEeCCCCccceeccccCCCcEEE----EecCHHHHHHHHHHhhcccCCCCCCchhhhhhhcCC-------------CCCc
Q 018216          209 VRSVNKGWGLYADQFIKQGQFIC----ELLTTKEARRRQQIYDGLASSPRNSSALLVIREHLP-------------SGKA  271 (359)
Q Consensus       209 ~~t~~kGwGlrA~~~I~~G~fI~----E~i~~~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~-------------~~~~  271 (359)
                      |+++.+|||||+++.|++|+||+    |+++..+++.|...++.......  .+.....+...             ....
T Consensus       181 frt~~kGwgvRs~~~I~~G~fvcEyaGe~~t~~e~~~~~~~~~~~~~~~~--~~~~~~~~~~~~~~~~~~~~~~~~~~~~  258 (364)
T KOG1082|consen  181 FRTPEKGWGVRTLDPIPAGEFVCEYAGEVLTSEEAQRRTHLREYLDDDCD--AYSIADREWVDESPVGNTFVAPSLPGGP  258 (364)
T ss_pred             EecCCceeeecccccccCCCeeEEEeeEecChHHhhhccccccccccccc--cchhhhccccccccccccccccccccCC
Confidence            99999999999999999999999    56677888877543332111111  00000000000             0111


Q ss_pred             ceeEEEeccccCCcccccccCCCCCce-EEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCCCC-----------CCC
Q 018216          272 CLRMNIDATRIGNIARFINHSCDGGNL-STTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEIRA-----------RPR  339 (359)
Q Consensus       272 ~~~~~IDa~~~GN~aRFiNHSC~PNN~-~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~~~-----------~~~  339 (359)
                      ...+.|||+..||++|||||||.||+. +.++..+.+..++|++|||+++|+||||||||||....           ...
T Consensus       259 ~~~~~ida~~~GNv~RfinHSC~PN~~~~~v~~~~~~~~~~~i~ffa~~~I~p~~ELT~dYg~~~~~~~~~~~~~~~~~~  338 (364)
T KOG1082|consen  259 GRELLIDAKPHGNVARFINHSCSPNLLYQAVFQDEFVLLYLRIGFFALRDISPGEELTLDYGKAYKLLVQDGANIYTPVM  338 (364)
T ss_pred             CcceEEchhhcccccccccCCCCccceeeeeeecCCccchheeeeeeccccCCCcccchhhccccccccccccccccccc
Confidence            246899999999999999999999543 45666655678899999999999999999999996521           346


Q ss_pred             CceeecCCCCCccccCCCCC
Q 018216          340 GLPCYCGSTSCFGILPSENT  359 (359)
Q Consensus       340 ~~~C~CGs~~Crg~l~~~~~  359 (359)
                      ...|.||+.+||+.+...++
T Consensus       339 ~~~c~c~~~~cr~~~~~~~~  358 (364)
T KOG1082|consen  339 KKNCNCGLEKCRGLLGSAPC  358 (364)
T ss_pred             chhhcCCCHHhCcccCCCcc
Confidence            78899999999999987653


No 3  
>KOG4442 consensus Clathrin coat binding protein/Huntingtin interacting protein HIP1, involved in regulation of endocytosis [Intracellular trafficking, secretion, and vesicular transport]
Probab=100.00  E-value=2.6e-43  Score=355.63  Aligned_cols=165  Identities=35%  Similarity=0.599  Sum_probs=149.9

Q ss_pred             cceeecCC-CCC-CCCCCCCccccCCCcccEEEEEeCCCCccceeccccCCCcEEE----EecCHHHHHHHHHHhhcccC
Q 018216          178 GIVSECGP-SCG-CGSECGNRLTQRGISVRLKIVRSVNKGWGLYADQFIKQGQFIC----ELLTTKEARRRQQIYDGLAS  251 (359)
Q Consensus       178 ~~i~EC~~-~C~-C~~~C~Nrv~Q~g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~----E~i~~~e~~~r~~~~~~~~~  251 (359)
                      ....||++ .|. |+..|+|+.+|+....+++||+|++|||||||.++|++|+||+    |+|..+|.++|...|+.   
T Consensus        92 ~t~iECs~~~C~~cg~~C~NQRFQkkqyA~vevF~Te~KG~GLRA~~dI~~g~FI~EY~GEVI~~~Ef~kR~~~Y~~---  168 (729)
T KOG4442|consen   92 MTSIECSDRECPRCGVYCKNQRFQKKQYAKVEVFLTEKKGCGLRAEEDIPKGQFILEYIGEVIEEKEFEKRVKRYAK---  168 (729)
T ss_pred             hhhcccCCccCCCccccccchhhhhhccCceeEEEecCcccceeeccccCCCcEEeeeccccccHHHHHHHHHHHHh---
Confidence            36789999 698 9999999999999999999999999999999999999999999    55677899999999997   


Q ss_pred             CCCCCchhhhhhhcCCCCCcceeEEEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeC
Q 018216          252 SPRNSSALLVIREHLPSGKACLRMNIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSY  331 (359)
Q Consensus       252 ~~~~~~~~~~~~~~~~~~~~~~~~~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dY  331 (359)
                      .+..+.|++.+..         ..+|||+.+||.||||||||+| |+....|.+.|..  ||+|||.|.|++||||||||
T Consensus       169 d~~kh~Yfm~L~~---------~e~IDAT~KGnlaRFiNHSC~P-Na~~~KWtV~~~l--RvGiFakk~I~~GEEITFDY  236 (729)
T KOG4442|consen  169 DGIKHYYFMALQG---------GEYIDATKKGNLARFINHSCDP-NAEVQKWTVPDEL--RVGIFAKKVIKPGEEITFDY  236 (729)
T ss_pred             cCCceEEEEEecC---------CceecccccCcHHHhhcCCCCC-CceeeeeeeCCee--EEEEeEecccCCCceeeEec
Confidence            3455666666554         4689999999999999999999 7888888888864  99999999999999999999


Q ss_pred             CCCCCCCCCceeecCCCCCccccCCC
Q 018216          332 GEIRARPRGLPCYCGSTSCFGILPSE  357 (359)
Q Consensus       332 g~~~~~~~~~~C~CGs~~Crg~l~~~  357 (359)
                      ++.++.....+|+||+++||||||++
T Consensus       237 qf~rYGr~AQ~CyCgeanC~G~IGgk  262 (729)
T KOG4442|consen  237 QFDRYGRDAQPCYCGEANCRGWIGGK  262 (729)
T ss_pred             ccccccccccccccCCcccccccCCC
Confidence            99999899999999999999999987


No 4  
>KOG1080 consensus Histone H3 (Lys4) methyltransferase complex, subunit SET1 and related methyltransferases [Chromatin structure and dynamics; Transcription]
Probab=99.96  E-value=3.1e-29  Score=267.92  Aligned_cols=136  Identities=35%  Similarity=0.626  Sum_probs=118.8

Q ss_pred             ccEEEEEeCCCCccceeccccCCCcEEEEecCH----HHHHHHHHHhhcccCCCCCCchhhhhhhcCCCCCcceeEEEec
Q 018216          204 VRLKIVRSVNKGWGLYADQFIKQGQFICELLTT----KEARRRQQIYDGLASSPRNSSALLVIREHLPSGKACLRMNIDA  279 (359)
Q Consensus       204 ~~leV~~t~~kGwGlrA~~~I~~G~fI~E~i~~----~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~IDa  279 (359)
                      ..|...++..+||||||.+.|.+|++|+||+++    ..|+.|+..|..   .+...+|+|.+++         ..+|||
T Consensus       866 k~~~F~~s~iH~wglfa~~~i~~~dmViEY~Ge~vR~~iad~RE~~Y~~---~gi~~sYlfrid~---------~~ViDA  933 (1005)
T KOG1080|consen  866 KYVKFGRSGIHGWGLFAMENIAAGDMVIEYRGELVRSSIADLREARYER---MGIGDSYLFRIDD---------EVVVDA  933 (1005)
T ss_pred             hhhccccccccccceeeccCccccceEEEeeceehhhhHHHHHHHHHhc---cCcccceeeeccc---------ceEEec
Confidence            347777899999999999999999999988765    467888888886   4557789988876         479999


Q ss_pred             cccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCCCCCCCCceeecCCCCCccccC
Q 018216          280 TRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEIRARPRGLPCYCGSTSCFGILP  355 (359)
Q Consensus       280 ~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~~~~~~~~~C~CGs~~Crg~l~  355 (359)
                      ++.||+||||||||+| ||....+.++|..  +|+|||.|+|.+||||||||.+.... .+.+|+|||++|||+|.
T Consensus       934 tk~gniAr~InHsC~P-NCyakvi~V~g~~--~IvIyakr~I~~~EElTYDYkF~~e~-~kipClCgap~Crg~~n 1005 (1005)
T KOG1080|consen  934 TKKGNIARFINHSCNP-NCYAKVITVEGDK--RIVIYSKRDIAAGEELTYDYKFPTED-DKIPCLCGAPNCRGFLN 1005 (1005)
T ss_pred             cccCchhheeecccCC-CceeeEEEecCee--EEEEEEecccccCceeeeeccccccc-cccccccCCCccccccC
Confidence            9999999999999999 8888888888875  99999999999999999999986333 38999999999999984


No 5  
>KOG1079 consensus Transcriptional repressor EZH1 [Transcription]
Probab=99.95  E-value=3.5e-29  Score=252.11  Aligned_cols=174  Identities=36%  Similarity=0.631  Sum_probs=146.4

Q ss_pred             ccCcccccccccCCC--CCCCCCcCCCccCCCCCCCcCCCCCCCcccccceeecCCC-CC-CC---------CCCCCccc
Q 018216          132 RLGFDSVSLVCESDE--SESGCDCEECFEVGLGDGVFGCPCFSGLEDVGIVSECGPS-CG-CG---------SECGNRLT  198 (359)
Q Consensus       132 ~~~~~~~~c~C~~~~--~~~gC~C~~~~~~~~~~~~~~c~C~~~~~~~~~i~EC~~~-C~-C~---------~~C~Nrv~  198 (359)
                      ...++..+|.|.-+|  ...||.|...|. +     ..|+|..      ...||+|. |. |+         .+|.|.-+
T Consensus       521 n~~~CEk~C~C~~dC~nrF~GC~Ck~QC~-t-----kqCpC~~------A~rECdPd~Cl~cg~~~~~d~~~~~C~N~~l  588 (739)
T KOG1079|consen  521 NETFCEKFCYCSPDCRNRFPGCRCKAQCN-T-----KQCPCYL------AVRECDPDVCLMCGNVDHFDSSKISCKNTNL  588 (739)
T ss_pred             cCcchhhcccCCHHHHhcCCCCCcccccc-c-----CcCchhh------hccccCchHHhccCcccccccCccccccchh
Confidence            345777888886444  348999976662 2     4699986      57899987 74 33         28999999


Q ss_pred             cCCCcccEEEEEeCCCCccceeccccCCCcEEE----EecCHHHHHHHHHHhhcccCCCCCCchhhhhhhcCCCCCccee
Q 018216          199 QRGISVRLKIVRSVNKGWGLYADQFIKQGQFIC----ELLTTKEARRRQQIYDGLASSPRNSSALLVIREHLPSGKACLR  274 (359)
Q Consensus       199 Q~g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~----E~i~~~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  274 (359)
                      |+|.+.++.|..|...|||||+++.+.|++||.    |+|+++||++|+.+|+.     ...+|+|.+..         +
T Consensus       589 ~~~~qkr~llapSdVaGwGlFlKe~v~KnefisEY~GE~IS~dEADrRGkiYDr-----~~cSflFnln~---------d  654 (739)
T KOG1079|consen  589 QRGEQKRVLLAPSDVAGWGLFLKESVSKNEFISEYTGEIISHDEADRRGKIYDR-----YMCSFLFNLNN---------D  654 (739)
T ss_pred             hhhhhcceeechhhccccceeeccccCCCceeeeecceeccchhhhhccccccc-----ccceeeeeccc---------c
Confidence            999999999999999999999999999999999    77999999999999996     33456666554         4


Q ss_pred             EEEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          275 MNIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       275 ~~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                      |+|||++.||.+||+|||-+| ||....+.+.|.+  ||+|||.|.|.+||||||||+|.
T Consensus       655 yviDs~rkGnk~rFANHS~nP-NCYAkvm~V~Gdh--RIGifAkRaIeagEELffDYrYs  711 (739)
T KOG1079|consen  655 YVIDSTRKGNKIRFANHSFNP-NCYAKVMMVAGDH--RIGIFAKRAIEAGEELFFDYRYS  711 (739)
T ss_pred             ceEeeeeecchhhhccCCCCC-CcEEEEEEecCCc--ceeeeehhhcccCceeeeeeccC
Confidence            899999999999999999999 6766666677764  99999999999999999999987


No 6  
>smart00317 SET SET (Su(var)3-9, Enhancer-of-zeste, Trithorax) domain. Putative methyl transferase, based on outlier plant homologues
Probab=99.86  E-value=1.6e-21  Score=159.54  Aligned_cols=112  Identities=45%  Similarity=0.737  Sum_probs=85.0

Q ss_pred             cEEEEEeCCCCccceeccccCCCcEEEEe----cCHHHHHHHHHHhhcccCCCCCCchhhhhhhcCCCCCcceeEEEecc
Q 018216          205 RLKIVRSVNKGWGLYADQFIKQGQFICEL----LTTKEARRRQQIYDGLASSPRNSSALLVIREHLPSGKACLRMNIDAT  280 (359)
Q Consensus       205 ~leV~~t~~kGwGlrA~~~I~~G~fI~E~----i~~~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~IDa~  280 (359)
                      ++++++++++|+||+|+++|++|++|+++    +...+...+...+....   ....+++....         .+.||++
T Consensus         1 ~~~~~~~~~~G~gl~a~~~i~~g~~i~~~~g~~~~~~~~~~~~~~~~~~~---~~~~~~~~~~~---------~~~id~~   68 (116)
T smart00317        1 KLEVFKSPGKGWGVRATEDIPKGEFIGEYVGEIITSEEAEERSKAYDTDG---ADSFYLFEIDS---------DLCIDAR   68 (116)
T ss_pred             CcEEEecCCCcEEEEECCccCCCCEEEEEEeEEECHHHHHHHHHHHHhcC---CCCEEEEECCC---------CEEEeCC
Confidence            36788899999999999999999999954    44555555544344311   11223332221         3789999


Q ss_pred             ccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeC
Q 018216          281 RIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSY  331 (359)
Q Consensus       281 ~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dY  331 (359)
                      ..||++|||||||.| |++...+...+.  .++.|+|+|||++|||||+||
T Consensus        69 ~~~~~~~~iNHsc~p-N~~~~~~~~~~~--~~~~~~a~r~I~~GeEi~i~Y  116 (116)
T smart00317       69 RKGNIARFINHSCEP-NCELLFVEVNGD--SRIVIFALRDIKPGEELTIDY  116 (116)
T ss_pred             ccCcHHHeeCCCCCC-CEEEEEEEECCC--cEEEEEECCCcCCCCEEeecC
Confidence            999999999999999 787777666554  389999999999999999999


No 7  
>KOG1083 consensus Putative transcription factor ASH1/LIN-59 [Transcription]
Probab=99.79  E-value=2.8e-20  Score=194.67  Aligned_cols=128  Identities=34%  Similarity=0.512  Sum_probs=100.9

Q ss_pred             CCCccccC-CCcccEEEEEeCCCCccceeccccCCCcEEEEecCHHHHHHHHHH--hhcccCCCCCCchhhhhhhcCCCC
Q 018216          193 CGNRLTQR-GISVRLKIVRSVNKGWGLYADQFIKQGQFICELLTTKEARRRQQI--YDGLASSPRNSSALLVIREHLPSG  269 (359)
Q Consensus       193 C~Nrv~Q~-g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~E~i~~~e~~~r~~~--~~~~~~~~~~~~~~~~~~~~~~~~  269 (359)
                      |.|+.+|+ +.-..|+||+.+.+||||++.++|++|+||+||+++....++.+.  ...+.  ..-..|.+.+..     
T Consensus      1166 c~nqrm~r~e~cp~L~v~~gp~~G~~v~tk~PikagtfI~EYvGeVit~ke~e~~mmtl~~--~d~~~~cL~I~p----- 1238 (1306)
T KOG1083|consen 1166 CSNQRMQRHEECPPLEVFRGPKKGWGVRTKEPIKAGTFIMEYVGEVITEKEFEPRMMTLYH--NDDDHYCLVIDP----- 1238 (1306)
T ss_pred             hhhHHhhhhccCCCcceeccCCCCccccccccccccchHHHHHHHHHHHHhhcccccccCC--CCCcccccccCc-----
Confidence            77776665 567789999999999999999999999999988887654433221  11100  111223333332     


Q ss_pred             CcceeEEEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          270 KACLRMNIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       270 ~~~~~~~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                          +.+||+.++||.+||+||||.| |+..+.|.++|.  .||++||+|||.+||||||||++.
T Consensus      1239 ----~l~id~~R~~n~~RfinhscKP-Nc~~qkwSVNG~--~Rv~L~A~rDi~kGEELtYDYN~k 1296 (1306)
T KOG1083|consen 1239 ----GLFIDIPRMGNGARFINHSCKP-NCEMQKWSVNGE--YRVGLFALRDLPKGEELTYDYNFK 1296 (1306)
T ss_pred             ----cccCChhhccccccccccccCC-CCccccccccce--eeeeeeecCCCCCCceEEEecccc
Confidence                4689999999999999999999 899999999986  599999999999999999999875


No 8  
>KOG1085 consensus Predicted methyltransferase (contains a SET domain) [General function prediction only]
Probab=99.76  E-value=1.1e-18  Score=161.24  Aligned_cols=124  Identities=33%  Similarity=0.444  Sum_probs=101.2

Q ss_pred             cCCCcccEEEEEeCCCCccceeccccCCCcEEEEecC----HHHHHHHHHHhhcccCCCCCCchhhhhhhcCCCCCccee
Q 018216          199 QRGISVRLKIVRSVNKGWGLYADQFIKQGQFICELLT----TKEARRRQQIYDGLASSPRNSSALLVIREHLPSGKACLR  274 (359)
Q Consensus       199 Q~g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~E~i~----~~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  274 (359)
                      .+|....|.+..-.+||.||+|...+.+|+||.||.+    -.|+..|+..|......+.++.|+-...         ..
T Consensus       251 l~g~~egl~~~~~dgKGRGv~a~~~F~rgdFVVEY~Gdliei~eAk~rE~~Ya~De~~GcYMYyF~h~s---------k~  321 (392)
T KOG1085|consen  251 LKGTNEGLLEVYKDGKGRGVRAKVNFERGDFVVEYRGDLIEISEAKVREEQYANDEEIGCYMYYFEHNS---------KK  321 (392)
T ss_pred             HhccccceeEEeeccccceeEeecccccCceEEEEecceeeechHHHHHHHhccCcccceEEEeeeccC---------ee
Confidence            4566677888888889999999999999999997754    4788999999987554444443332211         26


Q ss_pred             EEEeccccC-CcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          275 MNIDATRIG-NIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       275 ~~IDa~~~G-N~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                      |+|||++.- -++|.||||=.+ |+++.++.++|.  ||+.+.|.|||.+||||+||||+.
T Consensus       322 yCiDAT~et~~lGRLINHS~~g-Nl~TKvv~Idg~--pHLiLvA~rdIa~GEELlYDYGDR  379 (392)
T KOG1085|consen  322 YCIDATKETPWLGRLINHSVRG-NLKTKVVEIDGS--PHLILVARRDIAQGEELLYDYGDR  379 (392)
T ss_pred             eeeecccccccchhhhcccccC-cceeeEEEecCC--ceEEEEeccccccchhhhhhcccc
Confidence            899999854 579999999998 899999999885  799999999999999999999975


No 9  
>PF05033 Pre-SET:  Pre-SET motif;  InterPro: IPR007728 This region is found in a number of histone lysine methyltransferases (HMTase), N-terminal to the SET domain; it is generally described as the pre-SET domain. Histone lysine methylation is part of the histone code that regulated chromatin function and epigenetic control of gene function. Histone lysine methyltransferases (HMTase) differ both in their substrate specificity for the various acceptor lysines as well as in their product specificity for the number of methyl groups (one, two, or three) they transfer. With just one exception [], the HMTases belong to SET family that can be classified according to the sequences surrounding the SET domain [, ]. Structural studies on the human SET7/9, a mono-methylase, have revealed the molecular basis for the specificity of the enzyme for the histone-target and the roles of the invariant residues in the SET domain in determining the methylation specificities [].  The pre-SET domain, as found in the SUV39 SET family, contains nine invariant cysteine residues that are grouped into two segments separated by a region of variable length. These 9 cysteines coordinate 3 zinc ions to form a triangular cluster, where each of the zinc ions is coordinated by 4 four cysteines to give a tetrahedral configuration. The function of this domain is structural, holding together 2 long segments of random coils and stabilising the SET domain. The C-terminal region including the post-SET domain is disordered when not interacting with a histone tail and in the absence of zinc. The three conserved cysteines in the post-SET domain form a zinc-binding site [] when coupled to a fourth conserved cysteine in the knot-like structure close to the SET domain active site []. The structured post-SET region brings in the C-terminal residues that participate in S-adenosylmethine-binding and histone tail interactions. The three conserved cysteine residues are essential for HMTase activity, as replacement with serine abolishes HMTase activity []. ; GO: 0008270 zinc ion binding, 0018024 histone-lysine N-methyltransferase activity, 0034968 histone lysine methylation, 0005634 nucleus; PDB: 3K5K_A 2O8J_D 3RJW_B 1ML9_A 1PEG_B 1MVH_A 1MVX_A 3BO5_A 2RFI_B 3MO5_B ....
Probab=99.70  E-value=1.2e-17  Score=135.98  Aligned_cols=100  Identities=31%  Similarity=0.498  Sum_probs=52.5

Q ss_pred             cccCCcCCCCeeeEeCCCCCCC-CCeEEcccccCCCCCCCCCCCCCccccccccccccccccccccCcccccccccCCC-
Q 018216           69 DASRSVENFPIPFHNAADKTPY-AYFIYTPSQIIPPPCPAQFPPRQFWASTNAAADAESNSSLSRLGFDSVSLVCESDE-  146 (359)
Q Consensus        69 DiS~G~E~~pI~~vN~vD~~~~-~~f~Y~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~c~C~~~~-  146 (359)
                      |||+|+|++||+++|+||++++ +.|+|++++++++++.. +++...                       ..|.|...| 
T Consensus         1 Dis~g~e~~pI~~~N~vd~~~~p~~F~Yi~~~~~~~~~~~-~~~~~~-----------------------~~C~C~~~C~   56 (103)
T PF05033_consen    1 DISRGKENVPIPVVNDVDDEPPPPNFEYIPENIYGEGVPD-IDPEFL-----------------------QGCDCSGDCS   56 (103)
T ss_dssp             -TTCTSSSS-EEEEESSSS--SSTSSEE-SS-EESTTSS--TBGGGT-----------------------S----SSSST
T ss_pred             CCCCCccCCCEEEEeCCCCCCCCCCeEEeeeEEcCCCccc-cccccC-----------------------ccCccCCCCC
Confidence            8999999999999999999876 69999999987665431 111100                       123343333 


Q ss_pred             CCCCCCcCCCccCCCCCCCcCCCCCCC-cccccceeecCCCCCCCCCCCCc
Q 018216          147 SESGCDCEECFEVGLGDGVFGCPCFSG-LEDVGIVSECGPSCGCGSECGNR  196 (359)
Q Consensus       147 ~~~gC~C~~~~~~~~~~~~~~c~C~~~-~~~~~~i~EC~~~C~C~~~C~Nr  196 (359)
                      ....|.|......    ...|..=..+ .....+||||++.|+|+.+|.||
T Consensus        57 ~~~~C~C~~~~~~----~~~Y~~~g~l~~~~~~~i~EC~~~C~C~~~C~NR  103 (103)
T PF05033_consen   57 NPSNCECLQRNGG----IFAYDSNGRLRIPDKPPIFECNDNCGCSPSCRNR  103 (103)
T ss_dssp             CTTTSHHHCCTSS----S-SB-TTSSBSSSSTSEEE---TTSSS-TTSTT-
T ss_pred             CCCCCcCccccCc----cccccCCCcCccCCCCeEEeCCCCCCCCCCCCCC
Confidence            3345555433311    1122221111 13567899999999999999998


No 10 
>smart00468 PreSET N-terminal to some SET domains. A Cys-rich putative Zn2+-binding domain that occurs N-terminal to some SET domains. Function is unknown. Unpublished.
Probab=99.63  E-value=4.2e-16  Score=126.02  Aligned_cols=39  Identities=33%  Similarity=0.489  Sum_probs=34.3

Q ss_pred             cccccCCcCCCCeeeEeCCCCCCC-CCeEEcccccCCCCC
Q 018216           67 SLDASRSVENFPIPFHNAADKTPY-AYFIYTPSQIIPPPC  105 (359)
Q Consensus        67 ~~DiS~G~E~~pI~~vN~vD~~~~-~~f~Y~~~~~~~~~~  105 (359)
                      ..|||+|+|++||++||+||++++ +.|+|+++++++.+.
T Consensus         1 ~~Dis~G~E~~pI~~vN~vD~~~~p~~F~Yi~~~~~~~gv   40 (98)
T smart00468        1 CLDISNGKENVPVPLVNEVDEDPPPPDFEYISEYIYGQGV   40 (98)
T ss_pred             CccccCCccCCCcceEecCCCCCCCCCcEECcceEcCCCc
Confidence            379999999999999999999766 689999999876654


No 11 
>COG2940 Proteins containing SET domain [General function prediction only]
Probab=99.60  E-value=2.3e-16  Score=161.26  Aligned_cols=157  Identities=31%  Similarity=0.497  Sum_probs=115.0

Q ss_pred             CCCCCCCCCCCCccccCCCcccEEEEEeCCCCccceeccccCCCcEEEEec----CHHHHHHHHHHhhcccCCCCCCchh
Q 018216          184 GPSCGCGSECGNRLTQRGISVRLKIVRSVNKGWGLYADQFIKQGQFICELL----TTKEARRRQQIYDGLASSPRNSSAL  259 (359)
Q Consensus       184 ~~~C~C~~~C~Nrv~Q~g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~E~i----~~~e~~~r~~~~~~~~~~~~~~~~~  259 (359)
                      .........+.|...+........+..+..+||||||++.|++|++|.++.    ...++..+...++.+.    ....+
T Consensus       312 ~s~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~g~fa~~~i~~~e~i~~~~~~~~~~~~~~~~~~~~~~~~----~~~~~  387 (480)
T COG2940         312 KSNVSKLKELLNSNGCKKRREPNVVQESEIKGYGVFALESIKKGEFIIEYHGEIIRRKEAREREENYDLLG----NEFSF  387 (480)
T ss_pred             cccCccccchhhhcccccccchhhhhhhcccccceeehhhccchHHHHHhcCcccchHHHHhhhccccccc----cccch
Confidence            333333445666666677777788888999999999999999999988554    4455555555443211    11111


Q ss_pred             hhhhhcCCCCCcceeEEEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCCCCC--
Q 018216          260 LVIREHLPSGKACLRMNIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEIRAR--  337 (359)
Q Consensus       260 ~~~~~~~~~~~~~~~~~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~~~~--  337 (359)
                      ..+...        ..++|+...|+++|||||||.| |+........|  .-++.++|+|||.+||||++||+.....  
T Consensus       388 ~~~~~~--------~~~~d~~~~g~~~r~~nHS~~p-N~~~~~~~~~g--~~~~~~~~~rDI~~geEl~~dy~~~~~~~~  456 (480)
T COG2940         388 GLLEDK--------DKVRDSQKAGDVARFINHSCTP-NCEASPIEVNG--IFKISIYAIRDIKAGEELTYDYGPSLEDNR  456 (480)
T ss_pred             hhcccc--------chhhhhhhcccccceeecCCCC-Ccceecccccc--cceeeecccccchhhhhhccccccccccch
Confidence            111111        3679999999999999999999 78777666655  3489999999999999999999976322  


Q ss_pred             ------CCCceeecCCCCCccccC
Q 018216          338 ------PRGLPCYCGSTSCFGILP  355 (359)
Q Consensus       338 ------~~~~~C~CGs~~Crg~l~  355 (359)
                            .....|.||+..|+++|.
T Consensus       457 ~~~~~~~~~~~~~~~~~~~~~~~~  480 (480)
T COG2940         457 ELKKLLEKRWGCACGEDRCSHTMS  480 (480)
T ss_pred             hhhhhhhhhhccccCCCccCCCCC
Confidence                  246789999999999874


No 12 
>PF00856 SET:  SET domain;  InterPro: IPR001214 The SET domain appears generally as one part of a larger multidomain protein, and recently there were described three structures of very different proteins with distinct domain compositions: Neurospora crassa DIM-5, a member of the Su(var) family of HKMTs which methylate histone H3 on lysine 9,human SET7 (also called SET9), which methylates H3 on lysine 4 and garden pea Rubisco LSMT, an enzyme that does not modify histones, but instead methylates lysine 14 in the flexible tail of the large subunit of the enzyme Rubisco. The SET domain itself turned out to be an uncommon structure. Although in all three studies, electron density maps revealed the location of the AdoMet or AdoHcy cofactor, the SET domain bears no similarity at all to the canonical/AdoMet-dependent methyltransferase fold. Strictly conserved in the C-terminal motif of the SET domain tyrosine could be involved in abstracting a proton from the protonated amino group of the substrate lysine, promoting its nucleophilic attack on the sulphonium methyl group of the AdoMet cofactor. In contrast to the AdoMet-dependent protein methyltranferases of the classical type, which tend to bind their polypeptide substrates on top of the cofactor, it is noted from the Rubisco LSMT structure that the AdoMet seems to bind in a separate cleft, suggesting how a polypeptide substrate could be subjected to multiple rounds of methylation without having to be released from the enzyme. In contrast, SET7/9 is able to add only a single methyl group to its substrate. It has been demonstrated that association of SET domain and myotubularin-related proteins modulates growth control []. The SET domain-containing Drosophila melanogaster (Fruit fly) protein, enhancer of zeste, has a function in segment determination and the mammalian homologue may be involved in the regulation of gene transcription and chromatin structure. Histone lysine methylation is part of the histone code that regulated chromatin function and epigenetic control of gene function. Histone lysine methyltransferases (HMTase) differ both in their substrate specificity for the various acceptor lysines as well as in their product specificity for the number of methyl groups (one, two, or three) they transfer. With just one exception [], the HMTases belong to SET family that can be classified according to the sequences surrounding the SET domain [, ]. Structural studies on the human SET7/9, a mono-methylase, have revealed the molecular basis for the specificity of the enzyme for the histone-target and the roles of the invariant residues in the SET domain in determining the methylation specificities [].  The pre-SET domain, as found in the SUV39 SET family, contains nine invariant cysteine residues that are grouped into two segments separated by a region of variable length. These 9 cysteines coordinate 3 zinc ions to form to form a triangular cluster, where each of the zinc ions is coordinated by 4 four cysteines to give a tetrahedral configuration. The function of this domain is structural, holding together 2 long segments of random coils. The C-terminal region including the post-SET domain is disordered when not interacting with a histone tail and in the absence of zinc. The three conserved cysteines in the post-SET domain form a zinc-binding site when coupled to a fourth conserved cysteine in the knot-like structure close to the SET domain active site []. The structured post-SET region brings in the C-terminal residues that participate in S-adenosylmethine-binding and histone tail interactions. The three conserved cysteine residues are essential for HMTase activity, as replacement with serine abolishes HMTase activity [], []. ; GO: 0005515 protein binding; PDB: 3TG5_A 3S7F_A 3RIB_B 3TG4_A 3S7J_A 3S7D_A 3S7B_A 3H6L_A 3SMT_A 3K5K_A ....
Probab=99.40  E-value=5.2e-13  Score=113.94  Aligned_cols=54  Identities=35%  Similarity=0.530  Sum_probs=40.9

Q ss_pred             EEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCC
Q 018216          276 NIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYG  332 (359)
Q Consensus       276 ~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg  332 (359)
                      ..++.....++.|+||||.| |+...+....+.  ..+.|.|.|+|++|||||++||
T Consensus       109 ~~~~~~l~p~~d~~NHsc~p-n~~~~~~~~~~~--~~~~~~a~r~I~~GeEi~isYG  162 (162)
T PF00856_consen  109 DRDGIALYPFADMLNHSCDP-NCEVSFDFDGDG--GCLVVRATRDIKKGEEIFISYG  162 (162)
T ss_dssp             EEEEEEEETGGGGSEEESST-SEEEEEEEETTT--TEEEEEESS-B-TTSBEEEEST
T ss_pred             cccccccCcHhHhecccccc-ccceeeEeeccc--ceEEEEECCccCCCCEEEEEEC
Confidence            45667778899999999999 665544422222  3899999999999999999998


No 13 
>KOG1081 consensus Transcription factor NSD1 and related SET domain proteins [Transcription]
Probab=99.02  E-value=9.8e-11  Score=118.78  Aligned_cols=150  Identities=30%  Similarity=0.476  Sum_probs=105.4

Q ss_pred             ceeecCC-CCCCCCCCCCccccCCCcccEEEEEeCCCCccceeccccCCCcEEEEecCHHHHHHHHHHhhcccCCCCCCc
Q 018216          179 IVSECGP-SCGCGSECGNRLTQRGISVRLKIVRSVNKGWGLYADQFIKQGQFICELLTTKEARRRQQIYDGLASSPRNSS  257 (359)
Q Consensus       179 ~i~EC~~-~C~C~~~C~Nrv~Q~g~~~~leV~~t~~kGwGlrA~~~I~~G~fI~E~i~~~e~~~r~~~~~~~~~~~~~~~  257 (359)
                      ..+||-+ .|.+...|.|+..-+.....      +.+    +|..+|.+|  |.+++...+...|...-..   ......
T Consensus       288 ~~~~~~p~~~~~~~~~~~~~~sk~~~~e------~~~----~~~~~~~k~--vg~~i~~~e~~~~~~~~~~---~~~~~~  352 (463)
T KOG1081|consen  288 LAYEVHPKVCSAEERCHNQQFSKESYPE------PQK----TAKADIRKG--VGEVIDDKECKARLQRVKE---SDLVDF  352 (463)
T ss_pred             hhhhhcccccccccccccchhhhhcccc------cch----hhHHhhhcc--cCcccchhhheeehhhhhc---cchhhh
Confidence            4566654 49988899998774443322      222    888899998  7888888776655433222   011111


Q ss_pred             hhhhhhhcCCCCCcceeEEEeccccCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCCCCC
Q 018216          258 ALLVIREHLPSGKACLRMNIDATRIGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEIRAR  337 (359)
Q Consensus       258 ~~~~~~~~~~~~~~~~~~~IDa~~~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~~~~  337 (359)
                      ++..+..+         ..||+...||.+||+||||.| ++....+...+.  .++++||.++|++|||||++|..... 
T Consensus       353 ~~~~~e~~---------~~id~~~~~n~sr~~nh~~~~-~v~~~k~~~~~~--t~~~~~a~~~i~~g~e~t~~~n~~~~-  419 (463)
T KOG1081|consen  353 YMVFIQKD---------RIIDAGPKGNYSRFLNHSCQP-NVETEKWQVIGD--TRVGLFAPRQIEAGEELTFNYNGNCE-  419 (463)
T ss_pred             hhhhhhcc---------cccccccccchhhhhcccCCC-ceeechhheecc--cccccccccccccchhhhheeecccc-
Confidence            22222221         279999999999999999999 665555544443  48999999999999999999986533 


Q ss_pred             CCCceeecCCCCCccccCC
Q 018216          338 PRGLPCYCGSTSCFGILPS  356 (359)
Q Consensus       338 ~~~~~C~CGs~~Crg~l~~  356 (359)
                      ...+.|.||+.+|.++++.
T Consensus       420 ~~~~~~~~~~e~~~~~~~k  438 (463)
T KOG1081|consen  420 GNEKRCCCGSENCTETKGK  438 (463)
T ss_pred             CCcceEeecccccccCCcc
Confidence            4678999999999998864


No 14 
>KOG1141 consensus Predicted histone methyl transferase [Chromatin structure and dynamics]
Probab=98.89  E-value=2.7e-09  Score=110.78  Aligned_cols=298  Identities=24%  Similarity=0.323  Sum_probs=203.9

Q ss_pred             eeeeeCCchhhhhhHHHHh------ccCCcchhhhHhhhhhhhhhhcccccccccccccCCcCCCCeeeEeCCCCCCCCC
Q 018216           19 EVDEEQKQNQFLQWARLIL------PWLTPGELANVSLTCRTLSQISKSITLSRSLDASRSVENFPIPFHNAADKTPYAY   92 (359)
Q Consensus        19 ~~~~~~~~~~~~~~~~~v~------~~l~~~~~~~~~~~~~~~~~~~~~~~~~~~~DiS~G~E~~pI~~vN~vD~~~~~~   92 (359)
                      .++|-.|+|..++-..++-      +||...|++.|.-...+-.-+...-..+..+|.+.|.+.+|||.+|.+|+..+|.
T Consensus       823 ~fVciy~g~~l~~~~sdks~~~~~~~~~~~id~~~f~~~~dt~~~~tvD~~g~d~~d~~~g~sg~~~p~~~~~d~~~~~~  902 (1262)
T KOG1141|consen  823 NFVCIYPGGALLHQISDKSEYIHVTRSLLTIDCFSFDARIDTATYITVDDKGLDVADFSLGTSGIPIPLVNSVDNDEPPS  902 (1262)
T ss_pred             eEEEEecchhhhhhhchhhhhcccchhhhcccccchhccccccceeeccccccchhhhhccccCCCCccccccccCCCcc
Confidence            4677789999988765544      6666667766655544443343333456679999999999999999999988887


Q ss_pred             eEEcccccCCCCCCCCCCCCCccccccccccccccccccccCcccccccccCCCCCCCCCcCCCccCCCCCCCcCCCCCC
Q 018216           93 FIYTPSQIIPPPCPAQFPPRQFWASTNAAADAESNSSLSRLGFDSVSLVCESDESESGCDCEECFEVGLGDGVFGCPCFS  172 (359)
Q Consensus        93 f~Y~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~c~C~~~~~~~gC~C~~~~~~~~~~~~~~c~C~~  172 (359)
                      .+|......+-           |.+..         +....+|+            +||+|...+.|.     ..|.|.+
T Consensus       903 c~d~~~~~~~~-----------~~~~~---------s~~~~~~~------------~~~s~d~hp~d~-----~~~~~~~  945 (1262)
T KOG1141|consen  903 CEDSKRRFQYN-----------DQVDI---------SSVSRDFC------------SGCSCDGHPSDA-----SKCECQQ  945 (1262)
T ss_pred             ccccceeeccc-----------ccchh---------hhhccccc------------cccccCCCCccc-----CcccCCC
Confidence            77776543111           22111         11223443            567777666554     4555553


Q ss_pred             C--------------------cc--------cccceeecCCCCCCCCCCCCccccCCCcccE--------EEEEeCCCCc
Q 018216          173 G--------------------LE--------DVGIVSECGPSCGCGSECGNRLTQRGISVRL--------KIVRSVNKGW  216 (359)
Q Consensus       173 ~--------------------~~--------~~~~i~EC~~~C~C~~~C~Nrv~Q~g~~~~l--------eV~~t~~kGw  216 (359)
                      +                    .+        -....+||+..|.|...|.||++|.+.+.+.        +||++...||
T Consensus       946 ~~~~~~~~cpp~~s~d~~~~~~eS~~~~ns~~~~~f~e~~~hss~~~~e~~~~v~~~~~~~me~~s~~~l~i~~~~~~~~ 1025 (1262)
T KOG1141|consen  946 LSIEAMKRCPPNLSFDGHDELYESSEKQNSFLKLFFFECNDHSSCHRKEYNRVVQNNIKYPMEVSSFNDLQIFKTAQSGW 1025 (1262)
T ss_pred             CChhhhcCCCCccccCchhhhhhhhhhcchhhhccceeccccchhcccccchhhhcCCccceeeeecccccccccccccc
Confidence            2                    00        0235689999999999999999999887764        5566778899


Q ss_pred             cceeccccCCCcEEEEecCH----HHHH-HHH---HHhhc--------------------ccC-----------CC----
Q 018216          217 GLYADQFIKQGQFICELLTT----KEAR-RRQ---QIYDG--------------------LAS-----------SP----  253 (359)
Q Consensus       217 GlrA~~~I~~G~fI~E~i~~----~e~~-~r~---~~~~~--------------------~~~-----------~~----  253 (359)
                      |+++..+|+.-+||+++++.    ..+. -|.   ..|..                    ...           .+    
T Consensus      1026 ~~~edtD~~~~~~~~~~~~~ppt~~l~~~~r~aqad~~sn~~D~~~~~~l~es~~~~~T~~r~~t~~~~~~~~~d~dd~q 1105 (1262)
T KOG1141|consen 1026 GVREDTDIPQSTFICTYVGAPPTDDLADELRNAQADQYSNDLDLKDTVELEESREDHETDFRGDTSDYDDEEGSDGDDGQ 1105 (1262)
T ss_pred             cccccccCCCCcccccccCCCCchhhHHHHhhhhhccccCccchhhhhhhhhcccccccccCCCCCCCcccccccCccHH
Confidence            99999999999999977643    1111 111   00000                    000           00    


Q ss_pred             --------------CC----------------------------CchhhhhhhcCCCCC-----------cceeEEEecc
Q 018216          254 --------------RN----------------------------SSALLVIREHLPSGK-----------ACLRMNIDAT  280 (359)
Q Consensus       254 --------------~~----------------------------~~~~~~~~~~~~~~~-----------~~~~~~IDa~  280 (359)
                                    ..                            ...+...+ ...+..           ...-|+|||+
T Consensus      1106 ~I~k~ve~qd~~~~~~~T~~~~RQ~~~~s~k~~~~~s~~~~~~ts~~~~~~d-kges~~~~~~~~~~y~~~~~~yvIDAk 1184 (1262)
T KOG1141|consen 1106 DIMKMVERQDSSESGEETKRLTRQKRKQSKKSGKGGSVEKDDTTSRDSMEKD-KGESKDEPVFNWDKYFEPFPLYVIDAK 1184 (1262)
T ss_pred             HHHHHhhcccccccccccchhhhhhhhhhhhcccCccccccccCccchhhhc-cCccCcccccchhhccCCCceEEEecc
Confidence                          00                            00000000 000000           0013799999


Q ss_pred             ccCCcccccccCCCCCc-eEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC--CCCCCCceeecCCCCCcccc
Q 018216          281 RIGNIARFINHSCDGGN-LSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI--RARPRGLPCYCGSTSCFGIL  354 (359)
Q Consensus       281 ~~GN~aRFiNHSC~PNN-~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~--~~~~~~~~C~CGs~~Crg~l  354 (359)
                      ..||++||+||||.||. +|.|++++|+.++|+++|||.|-|+||+||||||+|.  .+..+...|+||+.+|||+|
T Consensus      1185 ~eGNlGRfLNHSC~PNl~VQnVfvdTHdlrfPwVAFFt~kyVkAgtELTWDY~Ye~g~v~~keL~C~CGa~~CrgrL 1261 (1262)
T KOG1141|consen 1185 QEGNLGRFLNHSCDPNLHVQNVFVDTHDLRFPWVAFFTRKYVKAGTELTWDYQYEQGQVATKELTCHCGAENCRGRL 1261 (1262)
T ss_pred             cccchhhhhccCCCccceeeeeeeeccccCCchhhhhhhhhhccCceeeeeccccccccccceEEEecChhhhhccc
Confidence            99999999999999953 3889999999999999999999999999999999997  67889999999999999998


No 15 
>KOG2589 consensus Histone tail methylase [Chromatin structure and dynamics]
Probab=98.57  E-value=6.3e-08  Score=93.17  Aligned_cols=59  Identities=29%  Similarity=0.437  Sum_probs=48.6

Q ss_pred             CCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCCCCCCCCceeecCC
Q 018216          283 GNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEIRARPRGLPCYCGS  347 (359)
Q Consensus       283 GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~~~~~~~~~C~CGs  347 (359)
                      -.-|+||||.|.| ||.++  .+++.   ++.+.++|||+||||||--||...+..+...|.|-+
T Consensus       194 LGPaafINHDCrp-nCkFv--s~g~~---tacvkvlRDIePGeEITcFYgs~fFG~~N~~CeC~T  252 (453)
T KOG2589|consen  194 LGPAAFINHDCRP-NCKFV--STGRD---TACVKVLRDIEPGEEITCFYGSGFFGENNEECECVT  252 (453)
T ss_pred             eccHHhhcCCCCC-Cceee--cCCCc---eeeeehhhcCCCCceeEEeecccccCCCCceeEEee
Confidence            3458999999999 66543  34333   899999999999999999999988888888999854


No 16 
>cd01395 HMT_MBD Methyl-CpG binding domains (MBD) present in putative histone methyltransferases (HMT) such as CLLD8 and SETDB1 proteins; CLLD8 contains a MBD, a PreSET and a bifurcated SET domain, suggesting that CLLD8 might be associated with methylation-mediated transcriptional repression. SETDB1 and other proteins in this group have a similar domain architecture. SETDB1 is a novel KAP-1-associated histone H3, lysine 9-specific methyltransferase that contributes to HP1-mediated silencing of euchromatic genes by KRAB zinc-finger proteins.
Probab=98.03  E-value=1.2e-06  Score=64.02  Aligned_cols=30  Identities=13%  Similarity=-0.031  Sum_probs=27.5

Q ss_pred             cccceeeeeeCCchhhhhhHHHHhccCCcc
Q 018216           14 HKRTIEVDEEQKQNQFLQWARLILPWLTPG   43 (359)
Q Consensus        14 ~~~~~~~~~~~~~~~~~~~~~~v~~~l~~~   43 (359)
                      ++...+|.|+||||++||++.||+.||..+
T Consensus        20 ~~~k~~V~Y~aPCGr~Lr~~~EV~~YL~~t   49 (60)
T cd01395          20 GKVKKHVIYKAPCGRSLRNMSEVHRYLRET   49 (60)
T ss_pred             CCcccceEEECCcchhhhcHHHHHHHHHhc
Confidence            356788999999999999999999999988


No 17 
>KOG2461 consensus Transcription factor BLIMP-1/PRDI-BF1, contains C2H2-type Zn-finger and SET domains [Transcription]
Probab=97.88  E-value=8.9e-06  Score=81.16  Aligned_cols=116  Identities=21%  Similarity=0.219  Sum_probs=75.7

Q ss_pred             CcccEEEEEe--CCCCccceeccccCCCcEEEEecCHHHHHHHHHHhhcccCCCCCCchhhhhhhcCCCCCcceeEEEec
Q 018216          202 ISVRLKIVRS--VNKGWGLYADQFIKQGQFICELLTTKEARRRQQIYDGLASSPRNSSALLVIREHLPSGKACLRMNIDA  279 (359)
Q Consensus       202 ~~~~leV~~t--~~kGwGlrA~~~I~~G~fI~E~i~~~e~~~r~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~IDa  279 (359)
                      +...|.++.+  ...|.||++...|++|+--+-+.++.-.....+..     +..+.-.++..+        ..-++||+
T Consensus        26 LP~~l~i~~Ssv~~~~lgV~s~~~i~~G~~FGP~~G~~~~~~~~~~~-----n~~y~W~I~~~d--------~~~~~iDg   92 (396)
T KOG2461|consen   26 LPPELRIKPSSVPVTGLGVWSNASILPGTSFGPFEGEIIASIDSKSA-----NNRYMWEIFSSD--------NGYEYIDG   92 (396)
T ss_pred             CCCceEeeccccCCccccccccccccCcccccCccCccccccccccc-----cCcceEEEEeCC--------CceEEecc
Confidence            6677888876  47889999999999988777444331011111100     111111222111        02378998


Q ss_pred             cc--cCCcccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          280 TR--IGNIARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       280 ~~--~GN~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                      +.  ..|+.||+|=+++. -.|+.+....+.   .|.++|.|+|++||||.+.|+.+
T Consensus        93 ~d~~~sNWmRYV~~Ar~~-eeQNL~A~Q~~~---~Ifyrt~r~I~p~eELlVWY~~e  145 (396)
T KOG2461|consen   93 TDEEHSNWMRYVNSARSE-EEQNLLAFQIGE---NIFYRTIRDIRPNEELLVWYGSE  145 (396)
T ss_pred             CChhhcceeeeecccCCh-hhhhHHHHhccC---ceEEEecccCCCCCeEEEEeccc
Confidence            76  78999999988886 334433333343   79999999999999999999976


No 18 
>smart00570 AWS associated with SET domains. subdomain of PRESET
Probab=97.06  E-value=0.00029  Score=49.83  Aligned_cols=31  Identities=32%  Similarity=0.767  Sum_probs=25.6

Q ss_pred             CCCCCCcccccceeecCCCCCCCCCCCCccccCCC
Q 018216          168 CPCFSGLEDVGIVSECGPSCGCGSECGNRLTQRGI  202 (359)
Q Consensus       168 c~C~~~~~~~~~i~EC~~~C~C~~~C~Nrv~Q~g~  202 (359)
                      ..|+++    .+.+||+..|+|++.|+||.+|+..
T Consensus        20 sdClNR----~l~~EC~~~C~~G~~C~NqrFqk~~   50 (51)
T smart00570       20 SDCLNR----MLLIECSSDCPCGSYCSNQRFQKRQ   50 (51)
T ss_pred             hHHHHH----HHhhhcCCCCCCCcCccCcccccCc
Confidence            347775    5789998889999999999999863


No 19 
>smart00508 PostSET Cysteine-rich motif following a subset of SET domains.
Probab=96.96  E-value=0.00043  Score=41.92  Aligned_cols=17  Identities=41%  Similarity=1.087  Sum_probs=15.1

Q ss_pred             CceeecCCCCCccccCC
Q 018216          340 GLPCYCGSTSCFGILPS  356 (359)
Q Consensus       340 ~~~C~CGs~~Crg~l~~  356 (359)
                      .+.|+|||.+|||+|+.
T Consensus         2 ~~~C~CGs~~CRG~l~~   18 (26)
T smart00508        2 KQPCLCGAPNCRGFLGX   18 (26)
T ss_pred             CeeeeCCCccccceecc
Confidence            47999999999999954


No 20 
>smart00391 MBD Methyl-CpG binding domain. Methyl-CpG binding domain, also known as the TAM (TTF-IIP5, ARBP, MeCP1) domain
Probab=96.12  E-value=0.0017  Score=50.09  Aligned_cols=31  Identities=10%  Similarity=-0.003  Sum_probs=27.7

Q ss_pred             cccceeeeeeCCchhhhhhHHHHhccCCcch
Q 018216           14 HKRTIEVDEEQKQNQFLQWARLILPWLTPGE   44 (359)
Q Consensus        14 ~~~~~~~~~~~~~~~~~~~~~~v~~~l~~~~   44 (359)
                      .....+|.|.+|||+.||++.+|..||..+.
T Consensus        24 ~~~~~dV~Y~sP~GkklRs~~ev~~YL~~~~   54 (77)
T smart00391       24 SAGKFDVYYISPCGKKLRSKSELARYLHKNG   54 (77)
T ss_pred             CCCcccEEEECCCCCeeeCHHHHHHHHHhCC
Confidence            4567999999999999999999999997653


No 21 
>cd00122 MBD MeCP2, MBD1, MBD2, MBD3, MBD4, CLLD8-like, and BAZ2A-like proteins constitute a family of proteins that share the methyl-CpG-binding domain (MBD). The MBD consists of about 70 residues and is defined as the minimal region required for binding to methylated DNA by a methyl-CpG-binding protein which binds specifically to methylated DNA. The MBD can recognize a single symmetrically methylated CpG either as naked DNA or within chromatin.  MeCP2, MBD1 and MBD2 (and likely MBD3) form complexes with histone deacetylase and are involved in histone deacetylase-dependent repression of transcription. MBD4 is an endonuclease that forms a complex with the DNA mismatch-repair protein MLH1. The MBDs present in putative chromatin remodelling subunit, BAZ2A, and putative histone methyltransferase, CLLD8, represent two phylogenetically distinct groups within the MBD protein family.
Probab=95.71  E-value=0.0035  Score=46.21  Aligned_cols=29  Identities=10%  Similarity=-0.074  Sum_probs=26.9

Q ss_pred             ccceeeeeeCCchhhhhhHHHHhccCCcc
Q 018216           15 KRTIEVDEEQKQNQFLQWARLILPWLTPG   43 (359)
Q Consensus        15 ~~~~~~~~~~~~~~~~~~~~~v~~~l~~~   43 (359)
                      ....+|.|.+|||+.||++.+|..||..+
T Consensus        22 ~~k~dv~Y~sP~Gk~~Rs~~ev~~yL~~~   50 (62)
T cd00122          22 AGKGDVYYYSPCGKKLRSKPEVARYLEKT   50 (62)
T ss_pred             CCcceEEEECCCCceecCHHHHHHHHHhC
Confidence            56889999999999999999999999875


No 22 
>cd01397 HAT_MBD Methyl-CpG binding domains (MBD) present in putative chromatin remodelling factor such as BAZ2A; BAZ2A contains a MBD, DDT, PHD-type zinc finger and Bromo domain suggesting that BAZ2A might be associated with histone acetyltransferase (HAT) activity. The Drosophila melanogaster toutatis protein, a putative subunit of the chromatin-remodeling complex, and other such proteins in this group share a similar domain architecture with BAZ2A, as does the Caenorhabditis elegans flectin homolog.
Probab=94.91  E-value=0.011  Score=45.09  Aligned_cols=29  Identities=14%  Similarity=0.029  Sum_probs=26.1

Q ss_pred             cccceeeeeeCCchhhhhhHHHHhccCCc
Q 018216           14 HKRTIEVDEEQKQNQFLQWARLILPWLTP   42 (359)
Q Consensus        14 ~~~~~~~~~~~~~~~~~~~~~~v~~~l~~   42 (359)
                      ....-+|.|.||||+.||+..+|.+||..
T Consensus        21 ~~~~~dV~Y~aPcGKklRs~~ev~~yL~~   49 (73)
T cd01397          21 GRIQGEVAYYAPCGKKLRQYPEVIKYLSK   49 (73)
T ss_pred             CCccceEEEECCCCcccccHHHHHHHHHh
Confidence            45667999999999999999999999985


No 23 
>PF01429 MBD:  Methyl-CpG binding domain;  InterPro: IPR001739 Methylation at CpG dinucleotide, the most common DNA modification in eukaryotes, has been correlated with gene silencing associated with various phenomena such as genomic imprinting, transposon and chromosome X inactivation, differentiation, and cancer. Effects of DNA methylation are mediated through proteins which bind to symmetrically methylated CpGs. Such proteins contain a specific domain of ~70 residues, the methyl-CpG-binding domain (MBD), which is linked to additional domains associated with chromatin, such as the bromodomain, the AT hook motif,the SET domain, or the PHD finger. MBD-containing proteins appear to act as structural proteins, which recruit a variety of histone deacetylase (HDAC) complexes and chromatin remodelling factors, leading to chromatin compaction and, consequently, to transcriptional repression. The MBD of MeCP2, MBD1, MBD2, MBD4 and BAZ2 mediates binding to DNA, in case of MeCP2, MBD1 and MBD2 preferentially to methylated CpG. In case of human MBD3 and SETDB1 the MBD has been shown to mediate protein-protein interactions [, ]. The MBD folds into an alpha/beta sandwich structure comprising a layer of twisted beta sheet, backed by another layer formed by the alpha1 helix and a hairpin loop at the C terminus. These layers are both amphipathic, with the alpha1 helix and the beta sheet lying parallel and the hydrophobic faces tightly packed against each other. The beta sheet is composed of two long inner strands (beta2 and beta3) sandwiched by two shorter outer strands (beta1 and beta4) [].; GO: 0003677 DNA binding, 0005634 nucleus; PDB: 2KY8_A 1UB1_A 1D9N_A 1IG4_A 1QK9_A 3C2I_A.
Probab=94.72  E-value=0.011  Score=45.64  Aligned_cols=31  Identities=6%  Similarity=-0.002  Sum_probs=27.5

Q ss_pred             ccceeeeeeCCchhhhhhHHHHhccCCcchh
Q 018216           15 KRTIEVDEEQKQNQFLQWARLILPWLTPGEL   45 (359)
Q Consensus        15 ~~~~~~~~~~~~~~~~~~~~~v~~~l~~~~~   45 (359)
                      ....+|.|.+|||+.+|++.||++||..+..
T Consensus        28 ~~~~dv~Y~sP~Gk~~RS~~eV~~yL~~~~~   58 (77)
T PF01429_consen   28 AGKKDVYYYSPCGKRFRSKKEVVRYLKENPS   58 (77)
T ss_dssp             TTSEEEEEEETTSEEESSHHHHHHHHTTSS-
T ss_pred             CCceEEEEECCCCCEEeCHHHHHHHHHhCCC
Confidence            3579999999999999999999999987763


No 24 
>KOG2084 consensus Predicted histone tail methylase containing SET domain [Chromatin structure and dynamics]
Probab=91.28  E-value=0.36  Score=48.68  Aligned_cols=40  Identities=38%  Similarity=0.639  Sum_probs=30.0

Q ss_pred             ccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCC-eEEEeCCCC
Q 018216          288 FINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGE-ELAFSYGEI  334 (359)
Q Consensus       288 FiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GE-ELT~dYg~~  334 (359)
                      ++||||.| |+. +.  ..+.   ...+.+...+.+++ ||+..|.+.
T Consensus       208 ~~~hsC~p-n~~-~~--~~~~---~~~~~~~~~~~~~~~~l~~~y~~~  248 (482)
T KOG2084|consen  208 LFNHSCFP-NIS-VI--FDGR---GLALLVPAGIDAGEEELTISYTDP  248 (482)
T ss_pred             hcccCCCC-CeE-EE--ECCc---eeEEEeecccCCCCCEEEEeeccc
Confidence            78999999 554 22  2333   57777888888887 999999866


No 25 
>cd01396 MeCP2_MBD MeCP2, MBD1, MBD2, MBD3, and MBD4 are members of a protein family that share the methyl-CpG-binding domain (MBD). The MBD, consists of about 70 residues and is defined as the minimal region required for binding to methylated DNA by a methyl-CpG-binding protein which binds specifically to methylated DNA. The MBD can recognize a single symmetrically methylated CpG either as naked DNA or within chromatin.  MeCP2, MBD1 and MBD2 (and likely MBD3) form complexes with histone deacetylase and are involved in histone deacetylase-dependent repression of transcription. MBD4 is an endonuclease that forms a complex with the DNA mismatch-repair protein MLH1.
Probab=90.07  E-value=0.11  Score=39.98  Aligned_cols=29  Identities=7%  Similarity=-0.006  Sum_probs=26.4

Q ss_pred             ccceeeeeeCCchhhhhhHHHHhccCCcc
Q 018216           15 KRTIEVDEEQKQNQFLQWARLILPWLTPG   43 (359)
Q Consensus        15 ~~~~~~~~~~~~~~~~~~~~~v~~~l~~~   43 (359)
                      ....+|.|.+|+|..+|+..+|.+||...
T Consensus        23 ~~k~DvyY~sP~Gkk~RS~~ev~~yL~~~   51 (77)
T cd01396          23 AGKFDVYYISPTGKKFRSKVELARYLEKN   51 (77)
T ss_pred             CCcceEEEECCCCCEEECHHHHHHHHHhC
Confidence            56789999999999999999999999763


No 26 
>KOG1337 consensus N-methyltransferase [General function prediction only]
Probab=82.14  E-value=1.1  Score=46.22  Aligned_cols=41  Identities=24%  Similarity=0.386  Sum_probs=31.4

Q ss_pred             ccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          288 FINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       288 FiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                      +.||++.+   ....+...+   +.+.+++.++|.+|||++++||..
T Consensus       239 ~~NH~~~~---~~~~~~~~d---~~~~l~~~~~v~~geevfi~YG~~  279 (472)
T KOG1337|consen  239 LLNHSPEV---IKAGYNQED---EAVELVAERDVSAGEEVFINYGPK  279 (472)
T ss_pred             hhccCchh---ccccccCCC---CcEEEEEeeeecCCCeEEEecCCC
Confidence            67999887   122333333   389999999999999999999964


No 27 
>PF03638 TCR:  Tesmin/TSO1-like CXC domain, cysteine-rich domain;  InterPro: IPR005172 This entry includes proteins that have two copies of a cysteine rich motif as follows: C-X-C-X4-C-X3-YC-X-C-X6-C-X3-C-X-C-X2-C. The family includes Tesmin Q9Y4I5 from SWISSPROT [] and TSO1 Q9LE32 from SWISSPROT []. This group of proteins is called a CXC domain in [].
Probab=79.03  E-value=0.96  Score=30.72  Aligned_cols=37  Identities=32%  Similarity=0.913  Sum_probs=27.1

Q ss_pred             CCCCcCCCccCCCCCCCcCCCCCCCcccccceeecCCCCCCCCCCCCcc
Q 018216          149 SGCDCEECFEVGLGDGVFGCPCFSGLEDVGIVSECGPSCGCGSECGNRL  197 (359)
Q Consensus       149 ~gC~C~~~~~~~~~~~~~~c~C~~~~~~~~~i~EC~~~C~C~~~C~Nrv  197 (359)
                      .||.|...-...     .||.|++      ....|++.|.| ..|.|+.
T Consensus         4 ~gC~Ckks~Clk-----~YC~Cf~------~g~~C~~~C~C-~~C~N~~   40 (42)
T PF03638_consen    4 KGCNCKKSKCLK-----LYCECFQ------AGRFCTPNCKC-QNCKNTE   40 (42)
T ss_pred             CCCcccCcChhh-----hhCHHHH------CcCcCCCCccc-CCCCCcC
Confidence            578886433233     6899987      36789999999 6788864


No 28 
>PF12937 F-box-like:  F-box-like; PDB: 1P22_A 2OVP_B 2OVR_B 2OVQ_B 1FS1_A 1FS2_C 1FQV_I 1LDK_E 2AST_B 2ASS_B.
Probab=61.39  E-value=5.4  Score=26.98  Aligned_cols=31  Identities=29%  Similarity=0.555  Sum_probs=25.6

Q ss_pred             hhHHHHhccCCcchhhhHhhhhhhhhhhccc
Q 018216           31 QWARLILPWLTPGELANVSLTCRTLSQISKS   61 (359)
Q Consensus        31 ~~~~~v~~~l~~~~~~~~~~~~~~~~~~~~~   61 (359)
                      +-..+|+.||.+.|+..++.+|+.+.++..+
T Consensus         7 Eil~~If~~L~~~dl~~~~~vcr~w~~~~~~   37 (47)
T PF12937_consen    7 EILLEIFSYLDPRDLLRLSLVCRRWRRIAND   37 (47)
T ss_dssp             HHHHHHHTTS-HHHHHHHTTSSHHHHHHHTC
T ss_pred             HHHHHHHhcCCHHHHHHHHHHHHHHHHHHCC
Confidence            4567899999999999999999999877653


No 29 
>PF08666 SAF:  SAF domain;  InterPro: IPR013974  This entry includes a range of different proteins, such as antifreeze proteins, flagellar FlgA proteins, and CpaB pilus proteins. ; PDB: 1C89_A 3NLA_A 3RDN_A 1C8A_A 3FRN_A 1WVO_A 3K3S_H 3G8R_B 1XUU_A 1XUZ_A ....
Probab=57.50  E-value=6.4  Score=28.11  Aligned_cols=15  Identities=27%  Similarity=0.501  Sum_probs=11.4

Q ss_pred             EEEEccCCCCCCeEE
Q 018216          314 CFFASKDIKEGEELA  328 (359)
Q Consensus       314 ~~fA~rdI~~GEELT  328 (359)
                      .+.|.|||++|+.|+
T Consensus         3 vvVA~~di~~G~~i~   17 (63)
T PF08666_consen    3 VVVAARDIPAGTVIT   17 (63)
T ss_dssp             EEEESSTB-TT-BEC
T ss_pred             EEEEeCccCCCCEEc
Confidence            478999999999995


No 30 
>PF00646 F-box:  F-box domain;  InterPro: IPR001810 The F-box domain was first described as a sequence motif found in cyclin-F that interacts with the protein SKP1 [, ]. This relatively conserved structural motif is present in numerous proteins and serves as a link between a target protein and a ubiquitin-conjugating enzyme. The SCF complex (e.g., Skp1-Cullin-F-box) plays a similar role as an E3 ligase in the ubiquitin protein degradation pathway [, ]. Different F-box proteins as a part of SCF complex recruit particular substrates for ubiquitination through specific protein-protein interaction domains.  Many mammalian F-box domains contain leucine-rich or WD-40 repeats (IPR001680 from INTERPRO). However, several F-box proteins either have other previously described domains such as Sec7 domain found in FBS protein or do not contain defined protein-protein interaction domains or motifs.; GO: 0005515 protein binding; PDB: 2E32_A 2E31_A 3V7D_B 1NEX_B 3MKS_D 3L2O_B.
Probab=49.47  E-value=9.1  Score=25.67  Aligned_cols=33  Identities=24%  Similarity=0.345  Sum_probs=26.5

Q ss_pred             hhhHHHHhccCCcchhhhHhhhhhhhhhhcccc
Q 018216           30 LQWARLILPWLTPGELANVSLTCRTLSQISKSI   62 (359)
Q Consensus        30 ~~~~~~v~~~l~~~~~~~~~~~~~~~~~~~~~~   62 (359)
                      ...+.+|+.||.+.|+..++.+++.+.++....
T Consensus         8 ~~il~~Il~~l~~~~~~~l~~vsk~~~~~~~~~   40 (48)
T PF00646_consen    8 DEILQEILSYLDPKDLLRLSLVSKRWRSLVDSP   40 (48)
T ss_dssp             HHHHHHHHHTS-HHHHHHHCTT-HHHHHHHTTH
T ss_pred             HHHHHHHHHHCcHHHHHHHHHHhhHHHHHHcCC
Confidence            356889999999999999999999998876654


No 31 
>PF11403 Yeast_MT:  Yeast metallothionein;  InterPro: IPR022710  Metallothioneins are characterised by an abundance of cysteine residues and a lack of generic secondary structure motifs. This protein functions in primary metal storage, transport and detoxification []. For the first 40 residues in the protein the polypeptide wraps around the metal by forming two large parallel loops separated by a deep cleft containing the metal cluster []. ; PDB: 1AQS_A 1AQR_A 1RJU_V 1FMY_A 1AOO_A 1AQQ_A.
Probab=45.04  E-value=22  Score=22.83  Aligned_cols=18  Identities=28%  Similarity=0.748  Sum_probs=8.2

Q ss_pred             CCCCcCCCccCCCCCCCcCCCCC
Q 018216          149 SGCDCEECFEVGLGDGVFGCPCF  171 (359)
Q Consensus       149 ~gC~C~~~~~~~~~~~~~~c~C~  171 (359)
                      ..|+|..+|...     .+|+|.
T Consensus        22 kscscptgcnsd-----dkcpcg   39 (40)
T PF11403_consen   22 KSCSCPTGCNSD-----DKCPCG   39 (40)
T ss_dssp             TS-SS-TTTTSS-----TT--TT
T ss_pred             hcCCCCCCCCCC-----CcCCCC
Confidence            567787776222     567774


No 32 
>smart00256 FBOX A Receptor for Ubiquitination Targets.
Probab=42.78  E-value=18  Score=22.87  Aligned_cols=31  Identities=23%  Similarity=0.267  Sum_probs=25.9

Q ss_pred             hhHHHHhccCCcchhhhHhhhhhhhhhhccc
Q 018216           31 QWARLILPWLTPGELANVSLTCRTLSQISKS   61 (359)
Q Consensus        31 ~~~~~v~~~l~~~~~~~~~~~~~~~~~~~~~   61 (359)
                      ..+.+|+.||.+.|+..+...++.+..+...
T Consensus         4 ~ll~~I~~~l~~~d~~~~~~vc~~~~~~~~~   34 (41)
T smart00256        4 EILEEILSKLPPKDLLRLRKVSRRWRSLIDS   34 (41)
T ss_pred             HHHHHHHHcCCHHHHHHHHHHHHHHHHHhcC
Confidence            3467889999999999999999998776654


No 33 
>KOG2155 consensus Tubulin-tyrosine ligase-related protein [Posttranslational modification, protein turnover, chaperones]
Probab=39.91  E-value=17  Score=37.09  Aligned_cols=49  Identities=14%  Similarity=0.186  Sum_probs=33.3

Q ss_pred             cccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEeCCCC
Q 018216          285 IARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFSYGEI  334 (359)
Q Consensus       285 ~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~dYg~~  334 (359)
                      ++..+.||-+| |.........-.....-.++-+|++..|||+|-|+-..
T Consensus       204 fGsrvrHsdeP-nf~~aPf~fmPq~vaYsimwp~k~~~tgeE~trDfasg  252 (631)
T KOG2155|consen  204 FGSRVRHSDEP-NFRIAPFMFMPQNVAYSIMWPTKPVNTGEEITRDFASG  252 (631)
T ss_pred             hhhhhccCCCC-cceeeeheecchhcceeEEeeccCCCCchHHHHHHhhc
Confidence            45568999999 66443322222333456688999999999999887543


No 34 
>smart00858 SAF This domain family includes a range of different proteins. Such as antifreeze proteins and flagellar FlgA proteins, and CpaB pilus proteins.
Probab=30.46  E-value=34  Score=24.17  Aligned_cols=16  Identities=25%  Similarity=0.474  Sum_probs=13.6

Q ss_pred             EEEEccCCCCCCeEEE
Q 018216          314 CFFASKDIKEGEELAF  329 (359)
Q Consensus       314 ~~fA~rdI~~GEELT~  329 (359)
                      .+.|.++|++|+.|+-
T Consensus         3 v~va~~~i~~G~~i~~   18 (64)
T smart00858        3 VVVAARDLPAGEVITA   18 (64)
T ss_pred             EEEEeCccCCCCCcch
Confidence            4678999999999984


No 35 
>PF01847 VHL:  von Hippel-Lindau disease tumour suppressor protein;  InterPro: IPR022772 Von Hippel-Lindau disease tumor suppressor (VHL) has two domains: a roughly 100-residue N-terminal domain rich in beta sheet (beta domain) and a smaller alpha-helical domain (alpha domain), held together by two linkers and a polar interface. A large portion of the alpha domain surface, and a small portion of the beta domain, interact with ElonginC. About half of the tumorigenic mutations map to the alpha domain and its residues that contact ElonginC. The remaining mutations map to the beta domain, and significantly, to a beta domain surface patch uninvolved in ElonginC binding. This suggests that two intact macromolecular binding sites may be required for the tumor suppressor effects of VHL []. This entry represents both beta and alpha domains of VHL.; PDB: 3ZRF_L 3ZRC_C 1LM8_V 1LQB_C 1VCB_F.
Probab=24.41  E-value=84  Score=27.56  Aligned_cols=38  Identities=18%  Similarity=0.273  Sum_probs=20.5

Q ss_pred             cccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEe
Q 018216          285 IARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFS  330 (359)
Q Consensus       285 ~aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~d  330 (359)
                      ..+|+|+|..+  +...+++..|.    -..|  ..++||++..++
T Consensus        16 ~V~F~N~s~r~--V~v~Wldy~G~----~~~Y--~~L~Pg~~~~~~   53 (156)
T PF01847_consen   16 FVRFVNRSPRT--VDVYWLDYDGK----PVPY--GTLKPGQGRRQN   53 (156)
T ss_dssp             EEEEEE-SSS---EEEEEE-TTS-----EEE-----B-TTEEEEEE
T ss_pred             EEEEEECCCCE--EEEEEEcCCCc----Eeec--cccCCCCeEEcc
Confidence            35799999775  56677777774    2333  347888877543


No 36 
>smart00317 SET SET (Su(var)3-9, Enhancer-of-zeste, Trithorax) domain. Putative methyl transferase, based on outlier plant homologues
Probab=24.03  E-value=86  Score=24.30  Aligned_cols=16  Identities=25%  Similarity=0.250  Sum_probs=13.8

Q ss_pred             ccceeccccCCCcEEE
Q 018216          216 WGLYADQFIKQGQFIC  231 (359)
Q Consensus       216 wGlrA~~~I~~G~fI~  231 (359)
                      ..++|.++|++|+=|.
T Consensus        98 ~~~~a~r~I~~GeEi~  113 (116)
T smart00317       98 IVIFALRDIKPGEELT  113 (116)
T ss_pred             EEEEECCCcCCCCEEe
Confidence            6788999999999764


No 37 
>cd05468 pVHL von Hippel-Landau (pVHL) tumor suppressor protein. von Hippel-Landau (pVHL) protein, the gene product of VHL, is a critical regulator of the ubiquitous oxygen-sensing pathway. It is conserved throughout evolution, as its homologs are found in organisms ranging from mammals to the Drosophila melanogaster, Anopheles gambiae insects and the Caenorhabditis elegans nematode. pVHL acts as the substrate recognition component of an E3 ubiquitin ligase complex.  Several proteins have been identified as pVHL-binding proteins that are subject to ubiquitin-mediated proteolysis; the best characterized putative substrates are the alpha subunits of the hypoxia-inducible factor (HIF1alpha, HIF2alpha, and HIF3alpha). In addition to HIF degradation, pVHL has been implicated to be involved in HIF independent cellular processes. Germline VHL mutations cause renal cell carcinomas, hemangioblastomas and pheochromocytomas in humans. pVHL can bind to and direct the proper deposition of fibronecti
Probab=23.22  E-value=1.1e+02  Score=26.15  Aligned_cols=37  Identities=16%  Similarity=0.324  Sum_probs=25.0

Q ss_pred             ccccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCeEEEe
Q 018216          286 ARFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEELAFS  330 (359)
Q Consensus       286 aRFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEELT~d  330 (359)
                      .+|+|++-.|  ++..+++..|.    -..|+  .|+||++.+.+
T Consensus        11 v~F~N~t~~~--v~~~Wid~~G~----~~~Y~--~l~pg~~~~~~   47 (141)
T cd05468          11 VRFVNRTDRP--VELYWIDYDGK----PVSYG--TLQPGETVRQN   47 (141)
T ss_pred             EEEEeCCCCe--EEEEEECCCCC----EEEee--eeCCCCEEeec
Confidence            4689999654  66667777774    33444  48888887654


No 38 
>KOG1081 consensus Transcription factor NSD1 and related SET domain proteins [Transcription]
Probab=22.37  E-value=26  Score=36.14  Aligned_cols=138  Identities=10%  Similarity=-0.048  Sum_probs=72.0

Q ss_pred             eCCCCcc---ceeccccCCCcEEEEecCHHHHHHHHHHhhcccCCCCCCc-hhhhhhhcCCCCCcceeEEEeccccCCcc
Q 018216          211 SVNKGWG---LYADQFIKQGQFICELLTTKEARRRQQIYDGLASSPRNSS-ALLVIREHLPSGKACLRMNIDATRIGNIA  286 (359)
Q Consensus       211 t~~kGwG---lrA~~~I~~G~fI~E~i~~~e~~~r~~~~~~~~~~~~~~~-~~~~~~~~~~~~~~~~~~~IDa~~~GN~a  286 (359)
                      .+..+|+   .+|.+.+..|++|...++.........-.+..-....... -.+..      .........++...|+..
T Consensus       120 ~~~c~~~~~d~~~~~~~~~~~~vw~~vg~~~~~~c~vc~~~~~~~~~~~~~~~f~~------~~~~~~~~~~~~~~g~~~  193 (463)
T KOG1081|consen  120 LEKCSKRCTDCRAFKKREVGDLVWSKVGEYPWWPCMVCHDPLLPKGMKHDHVNFFG------CYAWTHEKRVFPYEGQSS  193 (463)
T ss_pred             CcccccCCcceeeeccccceeEEeEEcCcccccccceecCcccchhhccccceecc------chhhHHHhhhhhccchHH
Confidence            4455666   8888899999999976665322111000000000000000 00000      000111234445599999


Q ss_pred             cccccCCCCCceEEEEEEeCCCcccEEEEEEccCCCCCCe------EEEeCCCCCCCCCCceeecCCCCCccccCCCC
Q 018216          287 RFINHSCDGGNLSTTLVRSSGSILPRLCFFASKDIKEGEE------LAFSYGEIRARPRGLPCYCGSTSCFGILPSEN  358 (359)
Q Consensus       287 RFiNHSC~PNN~~~v~~~~~g~~~p~i~~fA~rdI~~GEE------LT~dYg~~~~~~~~~~C~CGs~~Crg~l~~~~  358 (359)
                      ++++|++.|+-........ +  .+++..++.+.++-++-      ++.+|....+ .....+.+.+..|...+....
T Consensus       194 ~~l~~~~~~~s~~~~~~~~-~--~~r~~~~~~q~~~~~~~~e~k~~~~~~~~~~~~-~~~~~~~~~~~~~~~k~~~~~  267 (463)
T KOG1081|consen  194 KLIPHSKKPASTMSEKIKE-A--KARFGKLKAQWEAGIKQKELKPEEYKRIKVVCP-IGDQQIYSAAVSCIKKLLAKP  267 (463)
T ss_pred             Hhhhhccccchhhhhhhhc-c--cchhhhcccchhhccchhhcccccccccccccC-cCcccccchhhhhhhhccccc
Confidence            9999999995333222222 1  24677777777666655      5555554322 233447888888877666544


No 39 
>KOG4289 consensus Cadherin EGF LAG seven-pass G-type receptor [Signal transduction mechanisms]
Probab=21.75  E-value=1.5e+02  Score=35.08  Aligned_cols=78  Identities=18%  Similarity=0.357  Sum_probs=40.3

Q ss_pred             CCCCCcCCCccCCCCCCCcCCCCCCCcccccceeecCCC--CCCCC--------CCCCccc----------cCCCcccEE
Q 018216          148 ESGCDCEECFEVGLGDGVFGCPCFSGLEDVGIVSECGPS--CGCGS--------ECGNRLT----------QRGISVRLK  207 (359)
Q Consensus       148 ~~gC~C~~~~~~~~~~~~~~c~C~~~~~~~~~i~EC~~~--C~C~~--------~C~Nrv~----------Q~g~~~~le  207 (359)
                      .--|.|+..-.... .+.-.|.|.-    .....||...  |+|.+        .|.|-..          ..+-.+.+|
T Consensus      1789 ~G~CqCKe~hy~~~-~~Cl~CdC~~----Gs~Sr~C~adGqC~C~pgaiGRqCdrCd~pfaevttlgCrv~y~~CPk~ie 1863 (2531)
T KOG4289|consen 1789 NGQCQCKENHYRPI-GSCLPCDCYF----GSDSRECDADGQCPCKPGAIGRQCDRCDNPFAEVTTLGCRVVYDKCPKAIE 1863 (2531)
T ss_pred             CcceeeccccccCC-Ccceeecccc----CCCcccccCCCcCCCCCccccccccccCChhhhccccCcEEccccCcHHHh
Confidence            34577764332221 1135677762    1356789865  88874        4665422          111111111


Q ss_pred             EE-EeCCCCccceeccccCCCcEE
Q 018216          208 IV-RSVNKGWGLYADQFIKQGQFI  230 (359)
Q Consensus       208 V~-~t~~kGwGlrA~~~I~~G~fI  230 (359)
                      -. .=+..-+|+-|..+-|+|++=
T Consensus      1864 ~giwWP~tkfG~~a~~pCP~G~~G 1887 (2531)
T KOG4289|consen 1864 AGIWWPQTKFGSPAAVPCPKGSSG 1887 (2531)
T ss_pred             hcccccccccCCcccccCCCCcch
Confidence            11 012334688888888888873


No 40 
>TIGR03569 NeuB_NnaB N-acetylneuraminate synthase. This family is a subset of the Pfam model pfam03102 and is believed to include only authentic NeuB N-acetylneuraminate (sialic acid) synthase enzymes. The majority of the genes identified by this model are observed adjacent to both the NeuA and NeuC genes which together effect the biosynthesis of CMP-N-acetylneuraminate from UDP-N-acetylglucosamine.
Probab=20.32  E-value=55  Score=32.19  Aligned_cols=20  Identities=40%  Similarity=0.499  Sum_probs=17.7

Q ss_pred             EEEEEEccCCCCCCeEEEeC
Q 018216          312 RLCFFASKDIKEGEELAFSY  331 (359)
Q Consensus       312 ~i~~fA~rdI~~GEELT~dY  331 (359)
                      |-.++|.+||++||.||.+-
T Consensus       277 rrsl~a~~di~~G~~lt~~~  296 (329)
T TIGR03569       277 RKSLVAAKDIKKGEIFTEDN  296 (329)
T ss_pred             ceEEEEccCcCCCCEecHHh
Confidence            67899999999999999764


Done!