Query         031868
Match_columns 151
No_of_seqs    119 out of 702
Neff          4.7 
Searched_HMMs 46136
Date          Fri Mar 29 06:28:51 2013
Command       hhsearch -i /work/01045/syshi/csienesis_hhblits_a3m/031868.a3m -d /work/01045/syshi/HHdatabase/Cdd.hhm -o /work/01045/syshi/hhsearch_cdd/031868hhsearch_cdd -cpu 12 -v 0 

 No Hit                             Prob E-value P-value  Score    SS Cols Query HMM  Template HMM
  1 PRK00802 3-methyladenine DNA g 100.0 7.4E-72 1.6E-76  450.1  14.7  150    1-150     2-153 (188)
  2 TIGR00567 3mg DNA-3-methyladen 100.0 1.7E-71 3.7E-76  449.2  15.0  147    4-150     1-154 (192)
  3 PF02245 Pur_DNA_glyco:  Methyl 100.0 3.3E-71 7.2E-76  444.8  11.0  144    5-149     1-149 (184)
  4 cd00540 AAG Alkyladenine DNA g 100.0 1.8E-69 3.8E-74  433.3  14.7  142    8-150     1-147 (179)
  5 COG2094 Mpg 3-methyladenine DN 100.0 8.5E-67 1.8E-71  422.8  13.0  148    2-150     4-157 (200)
  6 KOG4486 3-methyladenine DNA gl 100.0 6.9E-59 1.5E-63  380.6   9.4  148    1-150    56-209 (261)
  7 PF12006 DUF3500:  Protein of u  74.3    0.43 9.4E-06   41.4  -2.3   94    6-99     56-154 (313)
  8 PF03451 HELP:  HELP motif;  In  40.6      23 0.00051   25.0   2.1   25   47-72     43-67  (77)
  9 PF05268 GP38:  Phage tail fibr  33.0      21 0.00046   30.6   1.0   52   30-90     41-92  (260)
 10 PHA02097 hypothetical protein   28.3      33 0.00071   23.2   1.1   17   81-97     11-27  (59)
 11 KOG3297 DNA-directed RNA polym  25.6      50  0.0011   27.5   2.0   47    4-50     13-69  (202)
 12 PF14549 P22_Cro:  DNA-binding   23.9      42 0.00092   22.5   1.0   16  119-134     8-23  (60)
 13 PRK04057 30S ribosomal protein  23.5      15 0.00033   30.3  -1.4   35   16-51     20-57  (203)
 14 PF09857 DUF2084:  Uncharacteri  21.7      57  0.0012   23.8   1.4   30   22-51     18-52  (85)
 15 PF01015 Ribosomal_S3Ae:  Ribos  21.4      85  0.0019   25.6   2.6   36   16-52     26-64  (194)
 16 PF02011 Glyco_hydro_48:  Glyco  20.2      38 0.00083   32.5   0.3   10   72-81    140-149 (619)

No 1  
>PRK00802 3-methyladenine DNA glycosylase; Reviewed
Probab=100.00  E-value=7.4e-72  Score=450.06  Aligned_cols=150  Identities=43%  Similarity=0.659  Sum_probs=141.2

Q ss_pred             CCCCCccccCCCHHhHHhhhcCCeEEEc-CeeEEEEeeeccC-CCCCCccCCCCCCCccccccCCCCeEEEEEeccceee
Q 031868            1 MTILPHHFFQIDALDLAPRLLGKFLRRD-DVLLQITEVEAYR-PNDSACHGRFGITARTAPVFGPGGLAYVYLCYGLHTM   78 (151)
Q Consensus         1 m~~L~~~ff~~~~~~vA~~LLG~~Lv~~-~~~grIVETEAY~-~~D~AsHa~~grT~R~~~mf~~~G~~YVY~~YGmh~c   78 (151)
                      +++|+++||++|+++||++||||+|+++ .++|||||||||. ++|||||||+|||+||++||++|||+|||++||||||
T Consensus         2 ~~~l~~~Ff~~~~~~vA~~LLGk~Lv~~~~~~grIvETEAY~G~~D~A~Ha~~grT~Rn~~mf~~~G~~YVY~~YG~h~~   81 (188)
T PRK00802          2 GMPLPREFFARDALEVARDLLGKVLVHEGGVSGRIVETEAYIGADDPASHSYRGRTPRTEVMFGPPGHAYVYFIYGMHHC   81 (188)
T ss_pred             CCccCHHHHCCCHHHHHHHhCCCEEEECCEEEEEEEEEecccCCCCccccccCCCChhhHHHcCCCceEEEEEecCceeE
Confidence            3568999999999999999999999995 6789999999999 8999999999999999999999999999999999999


Q ss_pred             eeEEeCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCCceee
Q 031868           79 LNVVADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGELNCH  150 (151)
Q Consensus        79 lNivt~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~l~~~  150 (151)
                      |||||+++|+|+||||||+||++|++.|++||+.+...++|||||||||||||||+++||+||+.+++|||+
T Consensus        82 lNvVt~~~g~~~aVLIRA~ep~~g~~~m~~~R~~~~~~~~L~nGPGkL~~AlgI~~~~~g~~l~~~~~l~i~  153 (188)
T PRK00802         82 LNVVCGPEGTGAAVLIRALEPLEGIALMRRRRGGKRPEKNLCNGPGKLCKALGITLADNGADLFDASPLYIE  153 (188)
T ss_pred             EEEEECCCCCccEEEEEeccccccHHHHHHhcccCCcccccccCHHHHHHHhCCCHHHCCCcccCCCCEEEe
Confidence            999999999999999999999999999999995444457899999999999999999999999877889985


No 2  
>TIGR00567 3mg DNA-3-methyladenine glycosylase (3mg). This families are based on the phylogenomic analysis of JA Eisen (1999, Ph.D. Thesis, Stanford University). All proteins in this family for which the function is known are involved in the base excision repair of alkylation damage to DNA. The exact specificty of the type of alkylation damage repaired by each of these varies somewhat between species. Substrates include 3-methyl adenine, 7-methyl-guanaine, and 3-methyl-guanine.
Probab=100.00  E-value=1.7e-71  Score=449.16  Aligned_cols=147  Identities=50%  Similarity=0.761  Sum_probs=138.0

Q ss_pred             CCccccCCCHHhHHhhhcCCeEEEcC----e-eEEEEeeeccC-CCCCCccCCCCC-CCccccccCCCCeEEEEEeccce
Q 031868            4 LPHHFFQIDALDLAPRLLGKFLRRDD----V-LLQITEVEAYR-PNDSACHGRFGI-TARTAPVFGPGGLAYVYLCYGLH   76 (151)
Q Consensus         4 L~~~ff~~~~~~vA~~LLG~~Lv~~~----~-~grIVETEAY~-~~D~AsHa~~gr-T~R~~~mf~~~G~~YVY~~YGmh   76 (151)
                      |+++||++|+++||++||||+|+++.    . +|+|||||||. ++|||||||+|| |+||++||++||++|||++||||
T Consensus         1 l~~~Ff~~~~~~vA~~LLGk~Lv~~~~~g~~~~grIvEtEAY~G~~D~A~Ha~~grrT~Rn~~mfg~~G~~YVY~~YGmh   80 (192)
T TIGR00567         1 MPPEFFQQDAVTLAPRLLGQLLVRRLDKGTEVRGRIVETEAYMGPEDSAAHSRGGRQTPRTRGMFGPPGRLYVYIIYGMH   80 (192)
T ss_pred             CCHHHhcCCHHHHHHHhCCCEEEEECCCCcEEEEEEEEEecccCCCCccccccCCCCCcchHHHccCCceEEEEEecCCE
Confidence            68899999999999999999999862    3 69999999999 999999999996 99999999999999999999999


Q ss_pred             eeeeEEeCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCCceee
Q 031868           77 TMLNVVADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGELNCH  150 (151)
Q Consensus        77 ~clNivt~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~l~~~  150 (151)
                      |||||||+++|.|+||||||+||++|++.|++||+.....++|||||||||||||||+++||+||+++++||++
T Consensus        81 ~~lNvVt~~~g~p~aVLIRA~ep~~G~~~m~~~R~~~~~~~~L~nGPGkL~~ALgI~~~~ng~~l~~~~~l~i~  154 (192)
T TIGR00567        81 YMLNVVAAPEGVPAAVLIRALEPLEGLELMRERRGRSLKDRELTNGPGKLCQALGITMSFNGRDLIQPSALWLE  154 (192)
T ss_pred             EEEEEEECCCCCccEEEEEeccccCchHHHHHhcCCCccccccccCHHHHHHHhCCCHHHCCCcccCCCceEEe
Confidence            99999999999999999999999999999999997654346899999999999999999999999887789985


No 3  
>PF02245 Pur_DNA_glyco:  Methylpurine-DNA glycosylase (MPG);  InterPro: IPR003180 Methylpurine-DNA glycosylase is a base excision-repair protein. It is responsible for the hydrolysis of the deoxyribose N-glycosidic bond, excising 3-methyladenine and 3-methylguanine from damaged DNA []. Its action is induced by alkylating chemotherapeutics, as well as deaminated and lipid peroxidation-induced purine adducts []. MPG without an N-terminal extension excises hypoxanthine with one-third of the efficiency of full-length MPG under similar conditions, suggesting that is function may largely be attributable to the N-terminal extension [].; GO: 0003677 DNA binding, 0003905 alkylbase DNA N-glycosylase activity, 0006284 base-excision repair; PDB: 1BNK_A 1F6O_A 3QI5_A 3UBY_B 1F4R_A 1EWN_A.
Probab=100.00  E-value=3.3e-71  Score=444.81  Aligned_cols=144  Identities=44%  Similarity=0.683  Sum_probs=122.7

Q ss_pred             CccccCCCHHhHHhhhcCCeEEEc----CeeEEEEeeeccC-CCCCCccCCCCCCCccccccCCCCeEEEEEeccceeee
Q 031868            5 PHHFFQIDALDLAPRLLGKFLRRD----DVLLQITEVEAYR-PNDSACHGRFGITARTAPVFGPGGLAYVYLCYGLHTML   79 (151)
Q Consensus         5 ~~~ff~~~~~~vA~~LLG~~Lv~~----~~~grIVETEAY~-~~D~AsHa~~grT~R~~~mf~~~G~~YVY~~YGmh~cl   79 (151)
                      |++||++|+++||++||||.||++    .++|||||||||. ++|||||||+|+|+||++||++|||+|||++|||||||
T Consensus         1 p~~Ff~r~~~~vA~~LLG~~Lv~~~~~~~~~grIvEtEAY~g~~D~AsHa~~g~T~Rn~~mf~~~G~~YVY~~YGmh~cl   80 (184)
T PF02245_consen    1 PRDFFDRDTVEVARDLLGKVLVRRIPGGELSGRIVETEAYLGPEDPASHAYRGRTPRNEVMFGPPGHAYVYLIYGMHHCL   80 (184)
T ss_dssp             --HHHSSBHHHHHHHCTT-EEEEE-TTS-EEEEEEEEEEE-STT-TTSTTGGG--STTGGGGSSTTBEEEEEETTTEEEE
T ss_pred             CcchhccCHHHHHHHhCCCEEEEEeCCCeEEEEEEEEeeccCCCCCcccccCCCCcccHHHcCCCCEEEEEEecCCeeEE
Confidence            689999999999999999999995    4899999999999 89999999999999999999999999999999999999


Q ss_pred             eEEeCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCCcee
Q 031868           80 NVVADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGELNC  149 (151)
Q Consensus        80 Nivt~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~l~~  149 (151)
                      ||||+++|.|+||||||+||++|++.|++||+.++ .++|||||||||||||||+++||+||++++.|||
T Consensus        81 Nivt~~~g~~~aVLIRAieP~~G~~~m~~~R~~~~-~~~L~~GPGkL~~ALgI~~~~~g~~l~~~~~l~i  149 (184)
T PF02245_consen   81 NIVTGPEGVPEAVLIRAIEPLEGIELMRARRGGKR-DKNLTNGPGKLCQALGIDRELNGLDLCGDPSLWI  149 (184)
T ss_dssp             EEE-BSTT---EEEEEEEEEEE-HHHHHHHHHHHH-GGGSSSSHHHHHHHTT--GGGTT-BTTT-SSEEE
T ss_pred             eEEecCCCceEEEEEEeccccCCHHHHHHhcCCCC-ccccccCHHHHHHHhCCCHHHCCcccCCCCceEE
Confidence            99999999999999999999999999999998754 5789999999999999999999999999988998


No 4  
>cd00540 AAG Alkyladenine DNA glycosylase (AAG), also known as 3-methyladenine DNA glycosylase, catalyzes the first step in base excision repair (BER) by cleaving damaged DNA bases within double-stranded DNA to produce an abasic site. AAG bends DNA by intercalating between the base pairs, causing the damaged base to flip out of the double helix and into the enzyme active site for cleavage. Although AAG represents one of six DNA glycosylase classes, it lacks the helix-hairpin-helix active site motif associated with the other BER glycosylases and is structurally quite distinct from them.
Probab=100.00  E-value=1.8e-69  Score=433.33  Aligned_cols=142  Identities=44%  Similarity=0.648  Sum_probs=134.9

Q ss_pred             ccCCCHHhHHhhhcCCeEEEc----CeeEEEEeeeccC-CCCCCccCCCCCCCccccccCCCCeEEEEEeccceeeeeEE
Q 031868            8 FFQIDALDLAPRLLGKFLRRD----DVLLQITEVEAYR-PNDSACHGRFGITARTAPVFGPGGLAYVYLCYGLHTMLNVV   82 (151)
Q Consensus         8 ff~~~~~~vA~~LLG~~Lv~~----~~~grIVETEAY~-~~D~AsHa~~grT~R~~~mf~~~G~~YVY~~YGmh~clNiv   82 (151)
                      ||++|+++||++||||+|+++    .++|||||||||. ++|||||||+|||+||.+||++|||+|||++||||||||||
T Consensus         1 f~~~~~~~vA~~LLGk~Lv~~~~~~~~~grIvEtEAY~G~~D~A~Ha~~gr~~~~~~mfg~~G~~YVY~~YG~h~~lNvV   80 (179)
T cd00540           1 FFNRDTVEVARDLLGKVLVRRLPGGILSGRIVETEAYLGPEDPASHAYRGRTTRRTVMFGPPGHAYVYLIYGMHWCLNVV   80 (179)
T ss_pred             CCCCCHHHHHHHhCCCEEEEECCCCEEEEEEEEEeccCCCCCccccccCCCCccchhhccCCceEEEEeecCCEEEEEEE
Confidence            899999999999999999997    3789999999999 89999999999988888999999999999999999999999


Q ss_pred             eCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCCceee
Q 031868           83 ADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGELNCH  150 (151)
Q Consensus        83 t~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~l~~~  150 (151)
                      |+++|+|+||||||+||++|++.|++||+.++. ++|||||||||||||||+++||.||++++.|||+
T Consensus        81 t~~~g~~~aVLIRAiEp~~G~~~m~~~R~~~~~-~~L~nGPGkL~~AlgI~~~~ng~~l~~~~~l~i~  147 (179)
T cd00540          81 TGPEGEPAAVLIRALEPLEGLELMRERRGGKRK-RDLTNGPGKLCQALGIDKSLNGADLTDSSPLWLE  147 (179)
T ss_pred             ECCCCCccEEEEEeeccccchhhHHhccCCCcc-CeeccChHHHHHHhCCcHHHCCCccCCCCcEEEe
Confidence            999999999999999999999999999976643 7899999999999999999999999887789985


No 5  
>COG2094 Mpg 3-methyladenine DNA glycosylase [DNA replication, recombination, and repair]
Probab=100.00  E-value=8.5e-67  Score=422.83  Aligned_cols=148  Identities=46%  Similarity=0.687  Sum_probs=140.2

Q ss_pred             CCCCccccCCCHHhHHhhhcCCeEEEc--C--eeEEEEeeeccC-CCCCCccCCCCCCCccccccCCCCeEEEEEeccce
Q 031868            2 TILPHHFFQIDALDLAPRLLGKFLRRD--D--VLLQITEVEAYR-PNDSACHGRFGITARTAPVFGPGGLAYVYLCYGLH   76 (151)
Q Consensus         2 ~~L~~~ff~~~~~~vA~~LLG~~Lv~~--~--~~grIVETEAY~-~~D~AsHa~~grT~R~~~mf~~~G~~YVY~~YGmh   76 (151)
                      +.+|++||.+|+++||++||||.|+|+  +  ++|+|||||||. ++|||||||+|||+||++||++|||+|||++||||
T Consensus         4 ~~~p~~F~~rd~~~vAr~LLG~~lv~~~~g~~~~g~IVEtEAY~G~~D~A~Ha~~GrT~R~~aMf~~~G~~YvY~~~G~h   83 (200)
T COG2094           4 MPLPRSFFARDTLVVARELLGKTLVRRIGGLTTSGRIVETEAYLGPDDPACHAYRGRTKRNRAMFGPPGHAYVYRIYGMH   83 (200)
T ss_pred             ccCCHhhhccCHHHHHHHhcCcEEEEecCCcEEEEEEEEEeEecCCCchhhhhccCcCcccHhhccCCceEEEEEEeccE
Confidence            358999999999999999999999986  3  789999999999 99999999999999999999999999999999999


Q ss_pred             eeeeEEeCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCC-ceee
Q 031868           77 TMLNVVADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGE-LNCH  150 (151)
Q Consensus        77 ~clNivt~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~-l~~~  150 (151)
                      |||||||.+||.|+||||||+||++|++.|++||+.+. .++|||||||||||||||+++||.|+.+++. ||+.
T Consensus        84 ~~lNvv~~~eG~p~aVLIRAieP~eg~~~m~~rR~~~~-~~~l~nGPgkLc~ALgI~~~~~g~~l~~~~~~l~l~  157 (200)
T COG2094          84 YCLNVVTGPEGIPAAVLIRAIEPEEGIEPMRQRRGVRV-LRLLTNGPGKLCKALGITREDNGADLINPSEPLWLE  157 (200)
T ss_pred             EEEEEEecCCCCcceEEEEeeccccccchhhhccCccc-cchhccCchHHHHHhCCCHHHcCCcccCCCCceEEe
Confidence            99999999999999999999999999999999998774 4689999999999999999999999888766 8875


No 6  
>KOG4486 consensus 3-methyladenine DNA glycosylase [Replication, recombination and repair]
Probab=100.00  E-value=6.9e-59  Score=380.63  Aligned_cols=148  Identities=40%  Similarity=0.581  Sum_probs=141.8

Q ss_pred             CCCCCccccCCCHHhHHhhhcCCeEEEcC----eeEEEEeeeccC-CCCCCccCCCC-CCCccccccCCCCeEEEEEecc
Q 031868            1 MTILPHHFFQIDALDLAPRLLGKFLRRDD----VLLQITEVEAYR-PNDSACHGRFG-ITARTAPVFGPGGLAYVYLCYG   74 (151)
Q Consensus         1 m~~L~~~ff~~~~~~vA~~LLG~~Lv~~~----~~grIVETEAY~-~~D~AsHa~~g-rT~R~~~mf~~~G~~YVY~~YG   74 (151)
                      |+.++.+||+.|++++|+.|||+.|+|+.    +.|||||||||+ |+|+|||+.+| +||||+.||++|||+|||+|||
T Consensus        56 m~~mgpEFFqipA~~LArrlLG~~l~RR~~g~~~~gRIvEtEAYlgPeD~AcHsRgGr~TPRn~~mFmk~Gt~YVY~iYG  135 (261)
T KOG4486|consen   56 MMDMGPEFFQIPASQLARRLLGKMLCRRIEGRTTKGRIVETEAYLGPEDKACHSRGGRRTPRNSAMFMKAGTCYVYRIYG  135 (261)
T ss_pred             HhhcChhhhcCchHHHHHHHHHHHHhhhcCCceeeeeEEeeeeecCcchhhhhhcCCccCCCcccccccCceEEEEEeee
Confidence            56789999999999999999999999973    679999999999 99999999998 6999999999999999999999


Q ss_pred             ceeeeeEEeCCCCCccEEEecccccCCCHHHHHHHhcCCCCCCccCCChhHHHHHhCCCcccCCCCCCCCCCceee
Q 031868           75 LHTMLNVVADKEGVGAAVLIRSCAPVSGLKTIQQRRAQLTEKPVLLTGPGKVGQALGISTEWSNHPLYMPGELNCH  150 (151)
Q Consensus        75 mh~clNivt~~~g~~~aVLIRA~eP~~G~~~m~~rR~~~~~~~~L~~GPgkL~~ALgI~~~~~g~~l~~~~~l~~~  150 (151)
                      ||+|||||+  +|.+++|||||+||++|+|+|+.||+.+.+++.|+|||+||||||||+++|+.++|..++++|++
T Consensus       136 ~y~c~NIss--~G~GA~VLiRalePl~G~ETmreRrG~~~kdr~L~nGPsKlcQALgi~~e~d~~~La~~e~~wLe  209 (261)
T KOG4486|consen  136 RYECFNISS--VGAGAGVLVRALEPLCGVETMRERRGGRVKDRDLANGPSKLCQALGITREIDKEWLAGSEKIWLE  209 (261)
T ss_pred             ehheeeeec--cccceEEeeeccCcccchHHHHHhcCCCcCChhhhcCcHHHHHHhccccccccccccCCCcEeee
Confidence            999999999  78999999999999999999999999988889999999999999999999999999999999985


No 7  
>PF12006 DUF3500:  Protein of unknown function (DUF3500);  InterPro: IPR021889  This family of proteins is functionally uncharacterised. This protein is found in bacteria and eukaryotes. Proteins in this family are typically between 335 to 438 amino acids in length. This protein has a conserved GHH sequence motif. This protein has two completely conserved G residues that may be functionally important. 
Probab=74.27  E-value=0.43  Score=41.44  Aligned_cols=94  Identities=14%  Similarity=0.052  Sum_probs=64.6

Q ss_pred             ccccCCCHHhHHhhhcCCeEEEcC--eeEEEEeeeccC-CCCCCccCCCCCCCccccccCCCCe--EEEEEeccceeeee
Q 031868            6 HHFFQIDALDLAPRLLGKFLRRDD--VLLQITEVEAYR-PNDSACHGRFGITARTAPVFGPGGL--AYVYLCYGLHTMLN   80 (151)
Q Consensus         6 ~~ff~~~~~~vA~~LLG~~Lv~~~--~~grIVETEAY~-~~D~AsHa~~grT~R~~~mf~~~G~--~YVY~~YGmh~clN   80 (151)
                      .+.+..+..+.|..||=..|--++  ..-.|.+.|.|+ .-....=..++.-.-+-++||.|+.  -+-+.++|-|-|+|
T Consensus        56 l~~lt~~Qr~~a~~lL~~~LS~~Gy~k~~~im~~d~~L~~~~~~~~~~~~~~~Y~~~ifG~Ps~~~~Wg~~f~GHHlsln  135 (313)
T PF12006_consen   56 LDDLTDEQRKAALALLKAALSPEGYEKARGIMRLDDVLGELEGGSPTIRDPESYYFAIFGTPSTTGPWGWQFEGHHLSLN  135 (313)
T ss_pred             hHhCCHHHHHHHHHHHHHhcCHHHHHHHHHHHHHHHHHHHhhcCCcccCCCCceEEEEecCCCCCCCeeEEecccEEEEE
Confidence            345556666788888877776443  234577888887 3332222224456778899999998  78899999999999


Q ss_pred             EEeCCCCCccEEEeccccc
Q 031868           81 VVADKEGVGAAVLIRSCAP   99 (151)
Q Consensus        81 ivt~~~g~~~aVLIRA~eP   99 (151)
                      ++........+=..-|.||
T Consensus       136 ~~~~~~~v~~tP~F~Ga~P  154 (313)
T PF12006_consen  136 CTFVGGQVSITPTFFGAEP  154 (313)
T ss_pred             EEEeCCCEEECceeeCCCC
Confidence            9996655544456677777


No 8  
>PF03451 HELP:  HELP motif;  InterPro: IPR005108  The HELP (Hydrophobic ELP) domain is found in EMAP and EMAP-like proteins (ELPs) [, ]. Although called a domain it contains a predicted transmembrane helix and may not form a globular domain. It is also not clear if these proteins localize to membranes.
Probab=40.59  E-value=23  Score=25.05  Aligned_cols=25  Identities=24%  Similarity=0.367  Sum_probs=21.9

Q ss_pred             ccCCCCCCCccccccCCCCeEEEEEe
Q 031868           47 CHGRFGITARTAPVFGPGGLAYVYLC   72 (151)
Q Consensus        47 sHa~~grT~R~~~mf~~~G~~YVY~~   72 (151)
                      .|+|+|+.-||-..|-+.|-+ ||++
T Consensus        43 VyGYrg~d~R~Nl~y~~~gei-vY~~   67 (77)
T PF03451_consen   43 VYGYRGHDCRNNLFYNATGEI-VYFT   67 (77)
T ss_pred             EeccCCccccccEEECCCCCE-EEEe
Confidence            499999989999999999987 7776


No 9  
>PF05268 GP38:  Phage tail fibre adhesin Gp38;  InterPro: IPR007932 This entry is represented by Bacteriophage T2, Gp38. The characteristics of the protein distribution suggest prophage matches in addition to the phage matches. This family contains several Gp38 proteins from T-even-like phages. Gp38, together with a second phage protein, Gp57, catalyses the organisation of Gp37 but is absent from the phage particle. Gp37 is responsible for receptor recognition [].
Probab=32.99  E-value=21  Score=30.55  Aligned_cols=52  Identities=23%  Similarity=0.219  Sum_probs=39.2

Q ss_pred             eeEEEEeeeccCCCCCCccCCCCCCCccccccCCCCeEEEEEeccceeeeeEEeCCCCCcc
Q 031868           30 VLLQITEVEAYRPNDSACHGRFGITARTAPVFGPGGLAYVYLCYGLHTMLNVVADKEGVGA   90 (151)
Q Consensus        30 ~~grIVETEAY~~~D~AsHa~~grT~R~~~mf~~~G~~YVY~~YGmh~clNivt~~~g~~~   90 (151)
                      +.||=+|.  |. .=-|.|.|++.|-|+. ||..++.--|+-|-|     |+|+...++|.
T Consensus        41 m~GrS~ei--~~-tiGA~hn~n~~~~~~~-~~~~GsaPiV~~ITG-----~lVS~S~~vp~   92 (260)
T PF05268_consen   41 MAGRSKEI--IH-TIGADHNFNGQWFRDR-CFEAGSAPIVFNITG-----NLVSYSKDVPC   92 (260)
T ss_pred             hhccceeE--EE-ecccCcccchhHHHHH-HHHcCCCCEEEEecc-----ceeeccCCcee
Confidence            55666663  22 2347899999999886 899999999988876     78888777764


No 10 
>PHA02097 hypothetical protein
Probab=28.34  E-value=33  Score=23.16  Aligned_cols=17  Identities=29%  Similarity=0.372  Sum_probs=13.6

Q ss_pred             EEeCCCCCccEEEeccc
Q 031868           81 VVADKEGVGAAVLIRSC   97 (151)
Q Consensus        81 ivt~~~g~~~aVLIRA~   97 (151)
                      +||..-..|.+||||-.
T Consensus        11 ~vt~amntp~gv~iri~   27 (59)
T PHA02097         11 VVTSAMNTPGGVIIRIA   27 (59)
T ss_pred             eEEEEeeCCCcEEEEEE
Confidence            56666778999999964


No 11 
>KOG3297 consensus DNA-directed RNA polymerase subunit E' [Transcription]
Probab=25.65  E-value=50  Score=27.51  Aligned_cols=47  Identities=30%  Similarity=0.382  Sum_probs=32.5

Q ss_pred             CCccccCCCHH-----hHHhhhcCCeEEEcCe---eEEEEeeeccC--CCCCCccCC
Q 031868            4 LPHHFFQIDAL-----DLAPRLLGKFLRRDDV---LLQITEVEAYR--PNDSACHGR   50 (151)
Q Consensus         4 L~~~ff~~~~~-----~vA~~LLG~~Lv~~~~---~grIVETEAY~--~~D~AsHa~   50 (151)
                      +|++-|.++..     ++-+.|.+|+|---++   .-.|.|.|-+.  +.|-|||+-
T Consensus        13 I~P~qf~~~~~~ai~~eL~~k~anKvl~nvGLCI~vyDi~~v~e~~v~pGDGas~~~   69 (202)
T KOG3297|consen   13 IPPSQFEKPLEDAIKEELNRKLANKVLPNVGLCICVYDILEVEEGIVLPGDGASYAR   69 (202)
T ss_pred             cChHHhCchHHHHHHHHHHHHHHhhhcccccEEEEEeEeeeecceEEecCCCceEEE
Confidence            56667777643     4556677777754432   34799998875  899999974


No 12 
>PF14549 P22_Cro:  DNA-binding transcriptional regulator Cro; PDB: 1RZS_A 3BD1_A 3QWS_A 2HIN_B.
Probab=23.89  E-value=42  Score=22.47  Aligned_cols=16  Identities=38%  Similarity=0.690  Sum_probs=12.8

Q ss_pred             cCCChhHHHHHhCCCc
Q 031868          119 LLTGPGKVGQALGIST  134 (151)
Q Consensus       119 L~~GPgkL~~ALgI~~  134 (151)
                      ...|..+|+++|||+.
T Consensus         8 ~~G~~~~lAkalGVs~   23 (60)
T PF14549_consen    8 YFGGQSKLAKALGVSP   23 (60)
T ss_dssp             HHSSHHHHHHHHTS-H
T ss_pred             HHCCHHHHHHHHCCCH
Confidence            3568899999999986


No 13 
>PRK04057 30S ribosomal protein S3Ae; Validated
Probab=23.48  E-value=15  Score=30.35  Aligned_cols=35  Identities=11%  Similarity=0.115  Sum_probs=26.4

Q ss_pred             HHhhhcCCeEEEc--CeeEEEEeeeccC-CCCCCccCCC
Q 031868           16 LAPRLLGKFLRRD--DVLLQITEVEAYR-PNDSACHGRF   51 (151)
Q Consensus        16 vA~~LLG~~Lv~~--~~~grIVETEAY~-~~D~AsHa~~   51 (151)
                      .-..-+|++++..  .+.||++|+=.+. .+|. +|++.
T Consensus        20 F~~~~iG~T~a~~~~~l~GRv~EvsL~DL~~d~-~~~~~   57 (203)
T PRK04057         20 FGGVEIGETPADDPEKLIGRVVETTLGDLTGDF-SKQNV   57 (203)
T ss_pred             cCCceEEEEEccChhhcCCcEEEEEHHHhcCCh-hhceE
Confidence            3344578888865  4889999999998 7786 66664


No 14 
>PF09857 DUF2084:  Uncharacterized protein conserved in bacteria (DUF2084);  InterPro: IPR018654  This domain is found in various hypothetical bacterial proteins that have no known function. 
Probab=21.74  E-value=57  Score=23.79  Aligned_cols=30  Identities=27%  Similarity=0.235  Sum_probs=20.8

Q ss_pred             CCeEEEcCeeEEEEeeeccC-C----CCCCccCCC
Q 031868           22 GKFLRRDDVLLQITEVEAYR-P----NDSACHGRF   51 (151)
Q Consensus        22 G~~Lv~~~~~grIVETEAY~-~----~D~AsHa~~   51 (151)
                      |.++..++-.|+|+++|.|. .    .|=.-.-|+
T Consensus        18 G~I~~~rd~~gri~~v~C~TReG~~l~dctl~vF~   52 (85)
T PF09857_consen   18 GRIRHERDDSGRITAVECYTREGWLLSDCTLAVFR   52 (85)
T ss_pred             CeEEEEECCCCCEEEEEEEccCCeeeCCCCHHHHH
Confidence            56666667899999999997 3    354444454


No 15 
>PF01015 Ribosomal_S3Ae:  Ribosomal S3Ae family;  InterPro: IPR001593 Ribosomes are the particles that catalyse mRNA-directed protein synthesis in all organisms. The codons of the mRNA are exposed on the ribosome to allow tRNA binding. This leads to the incorporation of amino acids into the growing polypeptide chain in accordance with the genetic information. Incoming amino acid monomers enter the ribosomal A site in the form of aminoacyl-tRNAs complexed with elongation factor Tu (EF-Tu) and GTP. The growing polypeptide chain, situated in the P site as peptidyl-tRNA, is then transferred to aminoacyl-tRNA and the new peptidyl-tRNA, extended by one residue, is translocated to the P site with the aid the elongation factor G (EF-G) and GTP as the deacylated tRNA is released from the ribosome through one or more exit sites [, ]. About 2/3 of the mass of the ribosome consists of RNA and 1/3 of protein. The proteins are named in accordance with the subunit of the ribosome which they belong to - the small (S1 to S31) and the large (L1 to L44). Usually they decorate the rRNA cores of the subunits.  Many ribosomal proteins, particularly those of the large subunit, are composed of a globular, surfaced-exposed domain with long finger-like projections that extend into the rRNA core to stabilise its structure. Most of the proteins interact with multiple RNA elements, often from different domains. In the large subunit, about 1/3 of the 23S rRNA nucleotides are at least in van der Waal's contact with protein, and L22 interacts with all six domains of the 23S rRNA. Proteins S4 and S7, which initiate assembly of the 16S rRNA, are located at junctions of five and four RNA helices, respectively. In this way proteins serve to organise and stabilise the rRNA tertiary structure. While the crucial activities of decoding and peptide transfer are RNA based, proteins play an active role in functions that may have evolved to streamline the process of protein synthesis. In addition to their function in the ribosome, many ribosomal proteins have some function 'outside' the ribosome [, ]. A number of eukaryotic and archaebacterial ribosomal proteins can be grouped on the basis of sequence similarities. One of these families consists of proteins that have from 220 to 250 amino acids and represents Rps1 (eukaryotic) and Rps3Ae (archaeal and eukaryotic).; GO: 0003735 structural constituent of ribosome, 0006412 translation, 0005622 intracellular, 0005840 ribosome; PDB: 2XZN_4 2XZM_4 3U5C_B 3U5G_B.
Probab=21.44  E-value=85  Score=25.61  Aligned_cols=36  Identities=19%  Similarity=0.143  Sum_probs=24.8

Q ss_pred             HHhhhcCCeEEEc--CeeEEEEeeeccC-CCCCCccCCCC
Q 031868           16 LAPRLLGKFLRRD--DVLLQITEVEAYR-PNDSACHGRFG   52 (151)
Q Consensus        16 vA~~LLG~~Lv~~--~~~grIVETEAY~-~~D~AsHa~~g   52 (151)
                      .-..-+|++++-+  .+.||++|+=.+. .+|.. |+|+.
T Consensus        26 F~~~~iG~T~~~~~~~l~gRv~Evsl~DL~~d~~-~~~~K   64 (194)
T PF01015_consen   26 FGNRNIGKTPANKPEKLKGRVFEVSLADLTNDFS-KAYRK   64 (194)
T ss_dssp             SSSSEECEEEEE-CCCCCC-EEEEECHCCCSTTT-TSS-E
T ss_pred             hCcceeeEEEcCCcccccCeEEEEEHHHhcCchh-hhcEE
Confidence            3444589999876  5889999999998 66644 77764


No 16 
>PF02011 Glyco_hydro_48:  Glycosyl hydrolase family 48;  InterPro: IPR000556 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 48 GH48 from CAZY comprises enzymes with several known activities; endoglucanase (3.2.1.4 from EC); cellobiohydrolase (3.2.1.91 from EC). The largest cellulase gene sequenced to date is one of the cellulases (celA) from the genome of the thermophilic anaerobic bacterium Caldocellum saccharolyticum. The celA gene product is a polypeptide of 1751 amino acids; this has a multidomain structure comprising two catalytic domains and two cellulose-binding domains, linked by Pro-Thr-rich regions. The N-terminal domain encodes an endoglucanase activity on carboxymethylcellulose, consistent with its similarity to several endo-1, 4-beta-D-glucanase sequences. The C-terminal domain shows similarity to a cellulase from Clostridium thermocellum (CelS), which acts synergistically with a second component to hydrolyse crystalline cellulose [].; GO: 0004553 hydrolase activity, hydrolyzing O-glycosyl compounds, 0005975 carbohydrate metabolic process; PDB: 1FAE_A 1FBO_A 1FBW_A 1FCE_A 1F9D_A 1F9O_A 1G9G_A 1G9J_A 2QNO_A 1L1Y_E ....
Probab=20.17  E-value=38  Score=32.48  Aligned_cols=10  Identities=50%  Similarity=1.032  Sum_probs=6.4

Q ss_pred             eccceeeeeE
Q 031868           72 CYGLHTMLNV   81 (151)
Q Consensus        72 ~YGmh~clNi   81 (151)
                      +|+|||++.|
T Consensus       140 ~Y~MHWL~DV  149 (619)
T PF02011_consen  140 IYGMHWLLDV  149 (619)
T ss_dssp             --S-BSEEET
T ss_pred             ceeeeeeeec
Confidence            7999999987


Done!