Query         016586
Match_columns 387
No_of_seqs    288 out of 766
Neff          4.0 
Searched_HMMs 46136
Date          Fri Mar 29 08:11:16 2013
Command       hhsearch -i /work/01045/syshi/csienesis_hhblits_a3m/016586.a3m -d /work/01045/syshi/HHdatabase/Cdd.hhm -o /work/01045/syshi/hhsearch_cdd/016586hhsearch_cdd -cpu 12 -v 0 

 No Hit                             Prob E-value P-value  Score    SS Cols Query HMM  Template HMM
  1 KOG2480 3-hydroxy-3-methylglut 100.0  3E-102  6E-107  790.2  25.8  355   24-378    42-423 (602)
  2 cd00643 HMG-CoA_reductase_clas 100.0 8.9E-74 1.9E-78  574.4  23.3  239  137-377     1-247 (403)
  3 TIGR00920 2A060605 3-hydroxy-3 100.0 1.2E-73 2.6E-78  609.8  23.9  246  133-378   458-709 (886)
  4 TIGR00533 HMG_CoA_R_NADP 3-hyd 100.0 2.6E-71 5.6E-76  556.7  24.0  231  134-367     1-238 (402)
  5 COG1257 HMG1 Hydroxymethylglut 100.0 3.9E-68 8.4E-73  535.5  17.6  237  139-376     2-246 (436)
  6 PF00368 HMG-CoA_red:  Hydroxym 100.0 2.5E-65 5.4E-70  510.2  19.0  215  161-377     1-221 (373)
  7 cd00644 HMG-CoA_reductase_clas 100.0 1.7E-62 3.7E-67  495.2  18.9  207  159-369     3-221 (417)
  8 TIGR00532 HMG_CoA_R_NAD hydrox 100.0 5.5E-61 1.2E-65  481.3  19.7  223  144-371     3-238 (393)
  9 cd00365 HMG-CoA_reductase Hydr 100.0 1.3E-58 2.8E-63  462.2  20.9  213  151-369     1-220 (376)
 10 PF14150 YesK:  YesK-like prote  74.7     3.4 7.3E-05   34.5   3.2   57   35-95      1-58  (81)
 11 PF11947 DUF3464:  Protein of u  72.1     8.7 0.00019   35.5   5.5   70   26-115    46-129 (153)
 12 PRK09509 fieF ferrous iron eff  71.4      16 0.00035   35.8   7.6   68  283-350   213-280 (299)
 13 TIGR01297 CDF cation diffusion  60.4      63  0.0014   30.4   9.0   66  283-348   192-258 (268)
 14 PF14715 FixP_N:  N-terminal do  58.4      12 0.00025   28.7   3.0   25   31-55     18-42  (51)
 15 PF14012 DUF4229:  Protein of u  53.7      20 0.00043   28.7   3.8   30   33-62     30-59  (69)
 16 PF01545 Cation_efflux:  Cation  47.3 1.1E+02  0.0024   28.9   8.4   66  284-349   206-272 (284)
 17 COG0053 MMT1 Predicted Co/Zn/C  45.1      98  0.0021   31.0   8.0   78  282-359   214-291 (304)
 18 COG1862 YajC Preprotein transl  41.6      21 0.00046   30.6   2.4   27   34-60      6-32  (97)
 19 PRK13454 F0F1 ATP synthase sub  38.9      29 0.00064   32.0   3.0   21   36-56     30-50  (181)
 20 COG5548 Small integral membran  38.0      46 0.00099   29.0   3.8   48   36-99     30-77  (105)
 21 PF09551 Spore_II_R:  Stage II   37.7      69  0.0015   28.9   5.0   58  267-328    38-112 (130)
 22 PF09580 Spore_YhcN_YlaJ:  Spor  37.6 2.5E+02  0.0054   25.1   8.7   70  284-357    73-142 (177)
 23 PF01578 Cytochrom_C_asm:  Cyto  37.5      67  0.0015   29.6   5.2   50   39-99    162-211 (214)
 24 TIGR03063 srtB_target sortase   36.1      37 0.00081   23.4   2.4   26   25-56      2-27  (29)
 25 PF07045 DUF1330:  Protein of u  35.5      33 0.00071   26.4   2.3   43  237-282    17-59  (65)
 26 PF11360 DUF3110:  Protein of u  34.0      57  0.0012   27.3   3.7   34  245-281    11-44  (86)
 27 PF13858 DUF4199:  Protein of u  33.7 1.4E+02  0.0031   26.1   6.4   49   39-91     34-85  (163)
 28 TIGR00782 ccoP cytochrome c ox  33.0      39 0.00084   33.2   3.0   25   31-55     24-48  (285)
 29 TIGR02837 spore_II_R stage II   29.1      91   0.002   29.4   4.5   56  269-328    75-147 (168)
 30 PF07074 TRAP-gamma:  Transloco  29.1      36 0.00079   32.1   1.9   49   27-82    117-167 (170)
 31 PF11834 DUF3354:  Domain of un  27.9 1.4E+02  0.0031   24.1   4.9   34  282-315    25-68  (69)
 32 PRK13108 prolipoprotein diacyl  27.1 1.2E+02  0.0027   32.4   5.7   43   46-88    228-272 (460)
 33 PF06129 Chordopox_G3:  Chordop  25.7      57  0.0012   28.8   2.4   21   36-56      2-22  (109)
 34 TIGR02206 intg_mem_TP0381 cons  24.9 1.4E+02  0.0031   28.5   5.2   41   60-102    63-109 (222)
 35 PRK14645 hypothetical protein;  23.3 5.8E+02   0.013   23.4   8.8   60  283-342     6-69  (154)
 36 PF09531 Ndc1_Nup:  Nucleoporin  21.5      78  0.0017   34.1   3.0   26   34-62    140-165 (602)
 37 TIGR00739 yajC preprotein tran  20.5      75  0.0016   26.3   2.1   22   37-58      3-24  (84)

No 1  
>KOG2480 consensus 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) reductase [Lipid transport and metabolism]
Probab=100.00  E-value=3e-102  Score=790.25  Aligned_cols=355  Identities=55%  Similarity=0.824  Sum_probs=328.1

Q ss_pred             CCCCCCCccchhhhhhhHHHHHHHHHHHHHHHHhH--hhhhcCCCccccchHHHHHHHHH-HHHHHHHHhhcCccccccc
Q 016586           24 TPKASDALPLPLYLTNAIFFTLFFSVAYYLLHRWR--EKIRNSTPLHVVTLSEIAAIVSL-IASFIYLLGFFGIDFVQSF  100 (387)
Q Consensus        24 ~~~~~~~lplp~~~~n~~~~~~f~~~~~~l~~~wr--~kir~s~plh~~~~~e~~~~~~~-~as~~yl~~ffgi~~~~~~  100 (387)
                      ..++++.+|||+|+++..|+.|||+++|+++++|+  +|+|.|+|||++++++..+.... |++++|+++|++|+++|++
T Consensus        42 ~~~~~~~~~l~~~~~i~~~~~l~~~~~~~~~~~~~~~~k~~~s~~~~~~~~~~~~~~~~~~i~~~~~~~~~~~~~~~~~~  121 (602)
T KOG2480|consen   42 VCKASKYLPLSLYLTIESPVLLLLSVVYIGFENWRLADKIRISTPLHVVDLSELAAHIGESIAENIDLLNFAKIDLIKSV  121 (602)
T ss_pred             cchhhhhcccchHHhhhhHHHHHHHHHHHHHhhhhhhhhhhccCceeEEEccchhcceehhhhhhhhhhccccccchhhc
Confidence            45789999999999999999999999999999999  99999999999999999998888 9999999999999999999


Q ss_pred             cCCCCCccCCcccCCCc----cc--cCCCCC-c---------cccc-CCCCCHHHHHHHHHcCCCCccccccccCChHHH
Q 016586          101 ISRATPEAWDLEEDDSD----II--SRPPAP-I---------SITT-LSSAQDEDVVQSVIDGSIPSYALESKLGDCRRA  163 (387)
Q Consensus       101 ~~~~~~~~~~~~~~~~~----~~--~~~~~~-~---------~~~~-~~~~~deei~~~v~~g~i~~y~LE~~l~d~~rA  163 (387)
                      +++...+.|..+..-.+    ..  +.++.| .         +... .-+.+|+||++.|..|++|+|+||++|+|++||
T Consensus       122 ~s~~~~~~~~~~~~~~~~~~~~~~~~~~~~pi~~l~~~~~~~~~g~~~~~~~d~Eiv~lv~~G~~p~y~lEs~l~d~~Ra  201 (602)
T KOG2480|consen  122 ISDDDVESWVTDRTILCVSGRLVTFPSIPKPIRPLLEILNLFPSGQEVKSLSDEEIVQLVIAGKIPLYALESKLGDAERA  201 (602)
T ss_pred             cccccccccccccccccCCCcccccCCCCCCccchhhhcccCccccccccCChHHHHHHhhcCcccceeccccccchhhh
Confidence            98887778876544211    11  222211 1         0111 245689999999999999999999999999999


Q ss_pred             HHHHHHHHHHHhCCccC---CCCCCCCChhhhhhcccccceEEeeeceeeecceEEcCeEEeeeccccchhHHHHhhhhh
Q 016586          164 AAIRREALQKMTGRSLQ---GLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLLDGFEYSVPMATTEGCLVASTNRGC  240 (387)
Q Consensus       164 v~iRR~~l~~~tg~sL~---~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlInG~ey~VPMATtEgsLVASanRGa  240 (387)
                      +.+||++|.+.++.+..   .||+++|||+.+.|+||||+|||+||||||||||+|||++|||||||||||||||+||||
T Consensus       202 v~iRR~~l~~~~~~~~~~~~~lP~~~yDY~~Vl~aCCENvIGY~piPVGVaGPLlldG~~y~VPMATTEGaLVAStnRGc  281 (602)
T KOG2480|consen  202 VSIRRTALSRNAREPTGTSDKLPYEGYDYSRVLGACCENVIGYMPIPVGVAGPLLLDGTEYYVPMATTEGALVASTNRGC  281 (602)
T ss_pred             HHHHHHHHHHhhcCcccccccCCccCcCHHHHHHHhhhccccceecccccccceEECCeEEEeeeccccchhhHhhccch
Confidence            99999999999998776   999999999999999999999999999999999999999999999999999999999999


Q ss_pred             hhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhhcCCcceeeEEEEecCCeEEEEEEee
Q 016586          241 KAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVFNRSSRFARLQHIQCSIAGKNLYIRFCCT  320 (387)
Q Consensus       241 KaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~van~tSR~GgL~~I~~~i~G~~l~lrf~~d  320 (387)
                      |+|+.+||+++++.+|||||+|++||+++.||.++|.||+.||||+.|+++||++|||+|||+|++.++|+++|+||...
T Consensus       282 KaI~a~GGa~sv~~~dGMTRgPvVRFps~~rA~~~k~WLe~~E~~~~lk~~FnstSRFaRLq~i~~~iaG~~lyiRF~~~  361 (602)
T KOG2480|consen  282 KAINAGGGATSVLEKDGMTRGPVVRFPSARRAAEAKIWLESPENQEVLKKAFNSTSRFARLQSIHTTIAGRNLYIRFCTS  361 (602)
T ss_pred             hheecCCCeeEEEeecCccccceeecchHHHHHHHHHHhcChhhHHHHHHHhcccchhhhhheeeEEeecceeEEEEEec
Confidence            99999999999999999999999999999999999999999999999999999999999999999999999999999999


Q ss_pred             ccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCccccccce----eEEEeeeeeeeeh
Q 016586          321 TGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSC----SYILLSVTNEIWL  378 (387)
Q Consensus       321 TGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A----~~~~~~v~~e~~~  378 (387)
                      ||||||||||+|++|+++.+|.++||||+|++||||||||+++.|    +=|+.||-.|+-|
T Consensus       362 TGDAMGMNMISKgve~~l~~l~~~fpdM~vi~iSGNyCtDKKpAAiNWieGRGKsVV~Ea~I  423 (602)
T KOG2480|consen  362 TGDAMGMNMISKGVENALRFLSEEFPDMQVISISGNYCTDKKPAAINWIEGRGKSVVAEATI  423 (602)
T ss_pred             cccchhhhhHHHHHHHHHHHHHHhCCCceEEEeccCcccCCchHhhhhhccCCceEEEEEec
Confidence            999999999999999999999999999999999999999999987    4578899888866


No 2  
>cd00643 HMG-CoA_reductase_classI Class I hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR). Hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR), class I enzyme, homotetramer. Catalyzes the synthesis of coenzyme A and mevalonate in isoprenoid synthesis. In mammals this is the rate limiting committed step in cholesterol biosynthesis. Class I enzymes are found predominantly in eukaryotes and contain N-terminal membrane regions. With the exception of Archaeoglobus fulgidus, most archeae are assigned to class I, based on sequence similarity of the active site, even though they lack membrane regions. Yeast and human HMGR are divergent in their N-terminal regions, but are conserved in their active site. In contrast, human and bacterial HMGR differ in their active site architecture.
Probab=100.00  E-value=8.9e-74  Score=574.42  Aligned_cols=239  Identities=51%  Similarity=0.780  Sum_probs=229.3

Q ss_pred             HHHHHHHHcCCCCccccccccCChHHHHHHHHHHHHHHhCCccCCCCCCCCChhhhhhcccccceEEeeeceeeecceEE
Q 016586          137 EDVVQSVIDGSIPSYALESKLGDCRRAAAIRREALQKMTGRSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLL  216 (387)
Q Consensus       137 eei~~~v~~g~i~~y~LE~~l~d~~rAv~iRR~~l~~~tg~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlI  216 (387)
                      +||++++.+|++++|+||+.+.|++||+.+||+++++.+|.+++++|++++||+++.++||||+||++|||+||||||+|
T Consensus         1 ~ei~~~~~~~~~~~~~le~~~~~~~~a~~~Rr~~~~~~~~~~l~~~~~~~~d~~~~~~~~iEN~IG~~~vPlGvAgpl~I   80 (403)
T cd00643           1 EEIIDLLSAGHIKLYKLEKSLEDAERAVRIRRLYLEKSTGKSLEHLPYTTYDYSEVLGRNIENVIGYVQVPVGVAGPLLI   80 (403)
T ss_pred             ChHHHhhhcCCCCchhccccccChHHHHHHHHHHHHhccccccccCCCCCcCHHHHhccccccceeeEeeceEEecceEE
Confidence            58999999999999999999999999999999999999999999999999999999999999999999999999999999


Q ss_pred             c----CeEEeeeccccchhHHHHhhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHH
Q 016586          217 D----GFEYSVPMATTEGCLVASTNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVF  292 (387)
Q Consensus       217 n----G~ey~VPMATtEgsLVASanRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~va  292 (387)
                      |    |++|+|||||||||||||+|||||+++.+||+++.+++|+|||+|+|+|+++.+|.++++|+++  |+++|+++|
T Consensus        81 nG~~a~~~~~VPmATtEgslVAS~srGak~i~~~GG~~~~v~~~~Mtrap~~~~~~~~~a~~~~~wi~~--~~~~i~~~a  158 (403)
T cd00643          81 NGEYAGGEFYVPMATTEGALVASTNRGCKAINLSGGATTRVLGDGMTRAPVFRFPSAREAAEFKAWIEE--NFEAIKEVA  158 (403)
T ss_pred             ecccCCcEEEEEeeechhHHHHHHHHHHHHHHhcCCceEEEcCCcceeccEEEECchhhHHHHHHHHHH--HHHHHHHHH
Confidence            9    7799999999999999999999999999999999999999999999999999999999999986  999999999


Q ss_pred             hhhcCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCccccccceeE----E
Q 016586          293 NRSSRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSCSY----I  368 (387)
Q Consensus       293 n~tSR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A~~----~  368 (387)
                      |++||||||++|++++.|+++|+||.|+||||||||||||++|++|+||++++|++.+++|||||||||+++|.-    +
T Consensus       159 ~stsr~g~l~~i~~~~~g~~v~lrf~~~TgDAMG~NMv~~~~e~v~~~i~~~~~~~~~~~i~gN~~tdKK~sa~n~~~gR  238 (403)
T cd00643         159 ESTSRHARLQSIKPYIAGRSVYLRFEYTTGDAMGMNMVTKATEAACDWIEENFPDMEVISLSGNFCTDKKPSAINWIEGR  238 (403)
T ss_pred             HhcCCCcccceEEEEecCCEEEEEEEEEcCcchhhhHHHHHHHHHHHHHHHhCCCceEEEEeccccccccchhhccccCC
Confidence            999999999999999999999999999999999999999999999999999999998999999999999877765    3


Q ss_pred             Eeeeeeeee
Q 016586          369 LLSVTNEIW  377 (387)
Q Consensus       369 ~~~v~~e~~  377 (387)
                      +..|.+|+=
T Consensus       239 Gk~V~Ae~~  247 (403)
T cd00643         239 GKSVVAEAT  247 (403)
T ss_pred             ceEEEEEEE
Confidence            346665553


No 3  
>TIGR00920 2A060605 3-hydroxy-3-methylglutaryl-coenzyme A reductase.
Probab=100.00  E-value=1.2e-73  Score=609.78  Aligned_cols=246  Identities=50%  Similarity=0.813  Sum_probs=236.2

Q ss_pred             CCCHHHHHHHHHcCCCCccccccccCChHHHHHHHHHHHHHHhC--CccCCCCCCCCChhhhhhcccccceEEeeeceee
Q 016586          133 SAQDEDVVQSVIDGSIPSYALESKLGDCRRAAAIRREALQKMTG--RSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGI  210 (387)
Q Consensus       133 ~~~deei~~~v~~g~i~~y~LE~~l~d~~rAv~iRR~~l~~~tg--~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGV  210 (387)
                      ..+|+||++++.+|+++.|+||+.++|++||+++||++|++.+|  .+++++|++++||+++.++||||+|||+|||+||
T Consensus       458 ~l~~~ei~~~~~~g~~~~~~le~~~~~~~~av~~Rr~~l~~~~~~~~~l~~l~~~~~d~~~v~~~~iEN~IG~~~vPvGV  537 (886)
T TIGR00920       458 ALSDAEVISLVNAKHIPAYKLETVLDNPERGVAIRRQILSKKLPMPDALDVLPYKNYDYSKVMGACCENVIGYMPIPVGV  537 (886)
T ss_pred             cCCHHHHHHHHhcCCCCccchhhhccCHHHHHHHHHHHHHhhcCCchhhhcCCcCCcCHHHHHhcccceeeeeeecceee
Confidence            35899999999999999999999999999999999999999987  6899999999999999999999999999999999


Q ss_pred             ecceEEcCeEEeeeccccchhHHHHhhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHH
Q 016586          211 AGPLLLDGFEYSVPMATTEGCLVASTNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAV  290 (387)
Q Consensus       211 AGPLlInG~ey~VPMATtEgsLVASanRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~  290 (387)
                      ||||+|||++|+|||||||||||||+|||||+++.+||+++++++|+|||+|+|+|+++.+|.++++|++++++++.+++
T Consensus       538 AGpL~InG~~y~VPMATtEgsLVAS~nRGaK~i~~~GG~~a~v~~d~MtRapvv~f~~~~~A~~~~~wi~~~e~~~~ik~  617 (886)
T TIGR00920       538 AGPLLLDGKEYQVPMATTEGCLVASTNRGCRALMLGGGVRSRVLADGMTRGPVVRLPSACRAAEAKAWLEVPENFAVIKD  617 (886)
T ss_pred             eccEEECCeEEEEeeeeccchHHHHHhHHHHHHHhcCCeEEEEccCceeeeeeEEeCcHHhHHHHHHHHhChHhHHHHHH
Confidence            99999999999999999999999999999999999999999999999999999999999999999999998789999999


Q ss_pred             HHhhhcCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCccccccce----e
Q 016586          291 VFNRSSRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSC----S  366 (387)
Q Consensus       291 van~tSR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A----~  366 (387)
                      ++++++|||||++|++++.|+++|+||.++||||||||||||++|+++++|+++||++++++|||||||||+++|    +
T Consensus       618 ~~~StsR~g~L~~I~~~i~g~~v~lrf~~~TGDAMG~NMVn~atEa~~~~i~~~~~~~~~~sIsgN~ctDKK~sAiN~i~  697 (886)
T TIGR00920       618 AFDSTSRFARLKKIHIAMAGRNLYIRFQAKTGDAMGMNMISKGTEQALAELQEHFPDMQILSLSGNYCTDKKPAAINWIE  697 (886)
T ss_pred             hhhccCcccceeeEEEEeeCCEEEEEEEEEcCchhcchHHHHHHHHHHHHHHHhCCCCeEEEEecccccchhhhhhcccc
Confidence            999999999999999999999999999999999999999999999999999988999999999999999999999    4


Q ss_pred             EEEeeeeeeeeh
Q 016586          367 YILLSVTNEIWL  378 (387)
Q Consensus       367 ~~~~~v~~e~~~  378 (387)
                      -|+.+|.+|+=|
T Consensus       698 gRGk~Vvae~~I  709 (886)
T TIGR00920       698 GRGKSVVCEATI  709 (886)
T ss_pred             cCCcEEEEEEEE
Confidence            466677777533


No 4  
>TIGR00533 HMG_CoA_R_NADP 3-hydroxy-3-methylglutaryl Coenzyme A reductase, hydroxymethylglutaryl-CoA reductase (NADP). This model represents archaeal examples of the enzyme hydroxymethylglutaryl-CoA reductase (NADP) (EC 1.1.1.34) and the catalytic domain of eukaryotic examples, which also contain a hydrophobic N-terminal domain. This enzyme synthesizes mevalonate, a precursor of isopentenyl pyrophosphate (IPP), a building block for the synthesis of cholesterol, isoprenoids, and other molecules. A related hydroxymethylglutaryl-CoA reductase, typified by an example from Pseudomonas mevalonii, is NAD-dependent and catabolic.
Probab=100.00  E-value=2.6e-71  Score=556.74  Aligned_cols=231  Identities=41%  Similarity=0.684  Sum_probs=221.6

Q ss_pred             CCHHHHHHHHHcCCCCccccccccCChHHHHHHHHHHHHHHhCCccCCCCCCCCChhhhhhcccccceEEeeeceeeecc
Q 016586          134 AQDEDVVQSVIDGSIPSYALESKLGDCRRAAAIRREALQKMTGRSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGP  213 (387)
Q Consensus       134 ~~deei~~~v~~g~i~~y~LE~~l~d~~rAv~iRR~~l~~~tg~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGP  213 (387)
                      ++|+||++++.+|+++.|+||+++ |++||+++||++|++.+|.++++|+...+|++++.++||||+||++|||+|||||
T Consensus         1 ~~~~~~~~~~~~~~~~~~~l~~~~-~~~~a~~~R~~~l~~~~~~~~~~l~~~~~d~~~~~~~~iEN~IG~~~vPlGvagp   79 (402)
T TIGR00533         1 MENNEILELVLNGKIKLYQLEKKL-GTTRAVEIRRKFIEKLAGLESEHLPNYSIDYERAFGANIENVIGYMQIPLGVAGP   79 (402)
T ss_pred             CChHHHHHHHhcCCCCchhhhhhc-CcHHHHHHHHHHHHHhcCCChhhccccCCCHHHHhhcCceeccceeeeceeEecc
Confidence            479999999999999999999987 6999999999999999999999999999999999999999999999999999999


Q ss_pred             eEEcC----eEEeeeccccchhHHHHhhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHH
Q 016586          214 LLLDG----FEYSVPMATTEGCLVASTNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLA  289 (387)
Q Consensus       214 LlInG----~ey~VPMATtEgsLVASanRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~  289 (387)
                      |+|||    ++|+|||||||||||||+|||||+++.+||+++++++++|+|||+|+|+++.|+.++.+|+++  |+++|+
T Consensus        80 l~InG~~a~~~~~VPmAT~E~slVAS~srGak~i~~~GG~~~~v~~~~M~r~~~~~~~~~~da~~~~~~i~~--~~~~i~  157 (402)
T TIGR00533        80 LKIDGEYAKGEYYIPLATTEGALVASVNRGCSAITAGGGATVRVTKDGMTRAPVVRTPSVVRAGACRIWIDE--NQNAIK  157 (402)
T ss_pred             EEEcCcccCceEEEeeeechhhHHHHHhHHHHHHHhcCCeEEEECCCcEEccEEEEeCChhhHHHHHHHHHH--HHHHHH
Confidence            99999    799999999999999999999999999999999999999999999999999999999999985  999999


Q ss_pred             HHHhhhcCCcceeeEEE-EecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCC--CCeEEEEeCCCcccccccee
Q 016586          290 VVFNRSSRFARLQHIQC-SIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFP--DMDVIGISGMCLVIYLKSCS  366 (387)
Q Consensus       290 ~van~tSR~GgL~~I~~-~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp--~~~visIsgN~ctDrl~~A~  366 (387)
                      ++||+++|||||++|++ .+.|+++|+||.|||||||||||||+++|++|++|++++|  ++++++||||||+||++++.
T Consensus       158 ~~a~~tsr~g~l~~i~~~~~~g~~v~l~f~~~tgDAMG~Nmvn~a~e~v~~~i~~~~~~~~~~~~~I~sN~~tdkk~sa~  237 (402)
T TIGR00533       158 EAAESTTRHGKLQKIQPICLAGDLLYPRFVTTTGDAMGMNMVTIATEYALKQMVEEYGWEGMEVVAVSGNYCTDKKPAAI  237 (402)
T ss_pred             HHHHhhCCCCceeeeEEEeccCCEEEEEEEEEccccccchHHHHHHHHHHHHHHHHcCCCCceEEEEeccCccchhhhhh
Confidence            99999999999999998 5789999999999999999999999999999999999886  56789999999999995554


Q ss_pred             E
Q 016586          367 Y  367 (387)
Q Consensus       367 ~  367 (387)
                      .
T Consensus       238 n  238 (402)
T TIGR00533       238 N  238 (402)
T ss_pred             c
Confidence            3


No 5  
>COG1257 HMG1 Hydroxymethylglutaryl-CoA reductase [Lipid metabolism]
Probab=100.00  E-value=3.9e-68  Score=535.49  Aligned_cols=237  Identities=38%  Similarity=0.519  Sum_probs=223.2

Q ss_pred             HHHHHHcCCCCccccccccCChHHHHHHHHHHHHHHhCCccCCCCCCCCChh---hhhhcccccceEEeeeceeeecceE
Q 016586          139 VVQSVIDGSIPSYALESKLGDCRRAAAIRREALQKMTGRSLQGLPLDGFDYD---SILGQCCEMPVGYVQIPVGIAGPLL  215 (387)
Q Consensus       139 i~~~v~~g~i~~y~LE~~l~d~~rAv~iRR~~l~~~tg~sL~~L~~~~~Dy~---~i~g~~iEN~IG~vqIPvGVAGPLl  215 (387)
                      +++++.+|++++|++|++.+ .++++++||++|++.++.+++++.....|++   +++++||||+||++|+|+||||||+
T Consensus         2 ~~~~~~~~~~~~~~~~~f~~-~~~a~~~Rr~~l~~~~~l~~~~~~~~~~d~~~~~~~~~~~iENvIG~~~lPlGia~~l~   80 (436)
T COG1257           2 VVQKVLAGEIKLSELEKFYK-ANKAVEERRQALERFTGLSLEEIGSLSIDGSLPIDVANRNIENVIGTVQLPLGIAGPLL   80 (436)
T ss_pred             hhhhhhcccchhhhhhhhcc-ccccHHHHHHHHHHHhCCChhhhhhhccCcchhHHHHHHhhhhhhceeeeecccccceE
Confidence            67899999999999999984 7899999999999999999999998888876   9999999999999999999999999


Q ss_pred             EcCeEEeeeccccchhHHHHhhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhh
Q 016586          216 LDGFEYSVPMATTEGCLVASTNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVFNRS  295 (387)
Q Consensus       216 InG~ey~VPMATtEgsLVASanRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~van~t  295 (387)
                      |||++|+|||||||||||||+|||||+|+++||+++++++|+|||||+|+|.+..++.+++.|++.+++++.++++++++
T Consensus        81 InG~~Y~iPmATtEgalVAs~~rgaK~i~~~GG~~~~~~~~~Mtr~~v~~~~~~~~a~~~~~w~~~~e~~e~~k~~~~s~  160 (436)
T COG1257          81 INGKEYYIPMATTEGALVASASRGAKLITASGGATARVTEDGMTRQPVFRFVSLPEAAKFAIWVKKKEIIELAKEAAPST  160 (436)
T ss_pred             EcCceeEEeeeccchHHHHHHHhHHHHHHhcCCcEEEEeccccccccEEecCchHHHHHHHHHhhhHHHHHHHHHhhhcc
Confidence            99999999999999999999999999999999999999999999999999999999999999999888999999999999


Q ss_pred             cCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhC-CCCeEEEEeCCCcccc--ccce--eEEEe
Q 016586          296 SRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDF-PDMDVIGISGMCLVIY--LKSC--SYILL  370 (387)
Q Consensus       296 SR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~f-p~~~visIsgN~ctDr--l~~A--~~~~~  370 (387)
                      +||||+++|++.+.|+++|+||.+||||||||||||+++|+++++|+++| ++..+++||||||||+  .+.+  ..+..
T Consensus       161 ~~~~kl~~i~~~~~g~~l~l~~~~dT~DAMGmNMvn~~~E~v~~~i~~~~~~~~~~~~issN~ctdkk~~a~~~~~grgk  240 (436)
T COG1257         161 GRGGKLIHIEPFVEGNLLYLRFYVDTGDAMGMNMVNIATEAVAPFIEEETFGGAVLLAISSNLCTDKKPAAINSIEGRGK  240 (436)
T ss_pred             CCcceeeeeeeeccCcEEEEEEEEecchhhhhHHHHHHHHHHHHHHHHhccCCceEEEEecccccccchhheEEEEeccc
Confidence            99999999999999999999999999999999999999999999999996 6888899999999999  4443  44555


Q ss_pred             eeeeee
Q 016586          371 SVTNEI  376 (387)
Q Consensus       371 ~v~~e~  376 (387)
                      +|..|+
T Consensus       241 ~v~ae~  246 (436)
T COG1257         241 TVVAEA  246 (436)
T ss_pred             eeeccc
Confidence            666554


No 6  
>PF00368 HMG-CoA_red:  Hydroxymethylglutaryl-coenzyme A reductase;  InterPro: IPR002202 Synonym(s): 3-hydroxy-3-methylglutaryl-coenzyme A reductase, HMG-CoA reductase.  There are two distinct classes of hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase enzymes: class I consists of eukaryotic and most archaeal enzymes (1.1.1.34 from EC), while class II consists of prokaryotic enzymes (1.1.1.88 from EC) [, ]. Class I HMG-CoA reductases catalyse the NADP-dependent synthesis of mevalonate from 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA). In vertebrates, membrane-bound HMG-CoA reductase is the rate-limiting enzyme in the biosynthesis of cholesterol and other isoprenoids. In plants, mevalonate is the precursor of all isoprenoid compounds []. The reduction of HMG-CoA to mevalonate is regulated by feedback inhibition by sterols and non-sterol metabolites derived from mevalonate, including cholesterol. In archaea, HMG-CoA reductase is a cytoplasmic enzyme involved in the biosynthesis of the isoprenoids side chains of lipids []. Class I HMG-CoA reductases consist of an N-terminal membrane domain (lacking in archaeal enzymes), and a C-terminal catalytic region. The catalytic region can be subdivided into three domains: an N-domain (N-terminal), a large L-domain, and a small S-domain (inserted within the L-domain). The L-domain binds the substrate, while the S-domain binds NADP.  Class II HMG-CoA reductases catalyse the reverse reaction of class I enzymes, namely the NAD-dependent synthesis of HMG-CoA from mevalonate and CoA []. Some bacteria, such as Pseudomonas mevalonii, can use mevalonate as the sole carbon source. Class II enzymes lack a membrane domain. Their catalytic region is structurally related to that of class I enzymes, but it consists of only two domains: a large L-domain and a small S-domain (inserted within the L-domain). As with class I enzymes, the L-domain binds substrate, but the S-domain binds NAD (instead of NADP in class I).; GO: 0004420 hydroxymethylglutaryl-CoA reductase (NADPH) activity, 0050662 coenzyme binding, 0015936 coenzyme A metabolic process, 0055114 oxidation-reduction process; PDB: 2Q6B_C 2Q6C_D 1HWI_C 1HWJ_C 3CD5_B 2R4F_C 1HWL_B 1HWK_C 1DQA_B 3CCT_D ....
Probab=100.00  E-value=2.5e-65  Score=510.18  Aligned_cols=215  Identities=40%  Similarity=0.575  Sum_probs=193.8

Q ss_pred             HHHHHHHHHHHHHHhCCc--cCCCCCCCCChhhhhhcccccceEEeeeceeeecceEEcCeEEeeeccccchhHHHHhhh
Q 016586          161 RRAAAIRREALQKMTGRS--LQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLLDGFEYSVPMATTEGCLVASTNR  238 (387)
Q Consensus       161 ~rAv~iRR~~l~~~tg~s--L~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlInG~ey~VPMATtEgsLVASanR  238 (387)
                      +|++++||++|++.++.+  +++++..++|+++++++||||+||+++||+||||||+|||++|+|||||||||||||+||
T Consensus         1 ~ka~~~Rr~~l~~~~~~~~~~~~l~~~~~~~~~~~~~~iEN~IG~~~lPlGvagpl~InG~~y~VPmATeE~svVAsasr   80 (373)
T PF00368_consen    1 KKAIEERREFLEELTGTSEQLEHLPNYSLDPEEVADQMIENVIGYVQLPLGVAGPLLINGKDYYVPMATEEPSVVASASR   80 (373)
T ss_dssp             -CHHHHHHHHHHTTSSCCCHHHHCTTTTS-HHTHHHHHSSSEEEEEEEEEEEEEEEEETTEEEEEEEE-SSTTHHHHHHH
T ss_pred             CchHHHHHHHHHhhhCcchhHHhCCccCCCHHHHhccCCceeeeEEEEeeEEEccEEEcCEEEEeeeceechhHHHhhhh
Confidence            579999999999999988  999999999999999999999999999999999999999999999999999999999999


Q ss_pred             hhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhhcCCcceeeEEEEecCCeEEEEEE
Q 016586          239 GCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVFNRSSRFARLQHIQCSIAGKNLYIRFC  318 (387)
Q Consensus       239 GaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~van~tSR~GgL~~I~~~i~G~~l~lrf~  318 (387)
                      |||+++.+||+++.+.++.|+|||+|.|+++.|+.++++|+++  |+++|+++||.+||||||++|++++.|+++++||.
T Consensus        81 Gak~i~~sGG~~t~v~~~~m~rq~~~~f~~~~~a~~~~~~i~~--~~~~l~~~a~~~sr~ggl~~i~~~~~g~~v~l~~~  158 (373)
T PF00368_consen   81 GAKAINESGGFTTSVLEDGMTRQPVFVFDDVEDAAAFKQWIEE--NFEELKEIANSTSRGGGLRDIEVRIVGRFVHLRFH  158 (373)
T ss_dssp             HHHHHHHTTSBEEEEEEEEEEEEEEEE-SSHSHHHHHHHHHHH--THHHHHHHHHHHCCTEEEEEEEEEECSTEEEEEEE
T ss_pred             HHHHHhhcCCeEEEEcCCceeeeeEEEeCCHhhHHHHHHHHHH--HHHHHHHhhhhcCCceeeeEEEEEecCCEEEEEEE
Confidence            9999999999999999999999999999999999999999995  99999999999999999999999999999999999


Q ss_pred             eeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCcccccccee----EEEeeeeeeee
Q 016586          319 CTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSCS----YILLSVTNEIW  377 (387)
Q Consensus       319 ~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A~----~~~~~v~~e~~  377 (387)
                      |||||||||||||+++|++|+||++++++..+++|+||||+||++++.    -+++.|++|+.
T Consensus       159 ~dtgDAMGaNmvn~~~e~v~~~i~~~~~~~~~~~I~sN~~~dKK~s~~n~~~grG~~v~A~~~  221 (373)
T PF00368_consen  159 YDTGDAMGANMVNTATEAVCPWIEEQFGGMVLMSILSNLCTDKKPSAINWIEGRGKLVTAEAT  221 (373)
T ss_dssp             EEETTC--HHHHHHHHHHHHHHHHHHHTSEEEEEEE-SSSCCSS--HHHHHH--ECEEEEEEE
T ss_pred             EEccccchhhHHHHHHHHHHHHHHHhcccceEEeeecccchhhhhhhhheeccCCEEEEEEEE
Confidence            999999999999999999999999999888889999999999554443    33344444443


No 7  
>cd00644 HMG-CoA_reductase_classII Class II hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR). Hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR), class II, prokaryotic enzyme is a homodimer. Class II enzymes are found primarily in prokaryotes and Archaeoglobus fulgidus and are soluble as they lack the membrane region. Enzymes catalyze the synthesis of coenzyme A and mevalonate in isoprenoid synthesis. Bacteria, such as Pseudomonas mevalonii, which rely solely on mevalonate for their carbon source, catalyze the reverse reaction, using an NAD-dependent HMGR to deacetylate mevalonate into 3-hydroxy-3-methylglutaryl-CoA. Human and bacterial HMGR differ in their active site architecture.
Probab=100.00  E-value=1.7e-62  Score=495.22  Aligned_cols=207  Identities=21%  Similarity=0.232  Sum_probs=197.3

Q ss_pred             ChHHHHHHHHHHHHHH---hCCccCCCCCCCCChhhhhhcccccceEEeeeceeeecceEEcCeEEeeeccccchhHHHH
Q 016586          159 DCRRAAAIRREALQKM---TGRSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLLDGFEYSVPMATTEGCLVAS  235 (387)
Q Consensus       159 d~~rAv~iRR~~l~~~---tg~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlInG~ey~VPMATtEgsLVAS  235 (387)
                      -+++++++||++|++.   ++..++.++.++.|+++++++||||+||+++||+||||||+|||++|+|||||||||||||
T Consensus         3 fy~~s~~eR~~~l~~~~~lt~~~~~~l~~~~~~~~~~~~~~iEN~IG~~~vP~Gvag~l~InG~~y~VPMATeE~svVAs   82 (417)
T cd00644           3 FYKLSPEERLQILAEFAGLSEEDVQLLKSGGALPLELADQMIENVIGTFSLPLGVATNFLVNGKDYLVPMATEEPSVVAA   82 (417)
T ss_pred             ccccCHHHHHHHHHHhcCCCHHHHHHHhccCCChHHHHhcCcccccceeecceeEeccEEECCeEEEEeeeeccchHHHH
Confidence            4678899999999996   6667888999999999999999999999999999999999999999999999999999999


Q ss_pred             hhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhhc-----CCcceeeEEEEecC
Q 016586          236 TNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVFNRSS-----RFARLQHIQCSIAG  310 (387)
Q Consensus       236 anRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~van~tS-----R~GgL~~I~~~i~G  310 (387)
                      +|||||+++.+||+++.+.+++|+|||+|  .++.|+.++.+|+++  |+++|++++|+++     |||||++|++++.|
T Consensus        83 ~srGak~i~~~GG~~~~~~~~~m~~q~~~--~~~~~a~~~~~~i~~--~~~~l~~~a~~~~ps~~~rgggl~~i~~~~~~  158 (417)
T cd00644          83 ASNAAKIARKSGGFKTSSSDRLMIGQIQL--VDVSDPAKARAFILA--HKDEILEIANEAHPSLVKRGGGARDIEVRVLD  158 (417)
T ss_pred             HhHHHHHHHhcCCEEEEECCCceEeccEE--CcccCHHHHHHHHHH--hHHHHHHHHHhhChhHHhcCCceeEEEEEecC
Confidence            99999999999999999999999999986  777899999999986  9999999999984     99999999999998


Q ss_pred             ----CeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCccccccceeEEE
Q 016586          311 ----KNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSCSYIL  369 (387)
Q Consensus       311 ----~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A~~~~  369 (387)
                          +++|+||.+||||||||||||+++|+++++|++++|+..+++||||||+||+++|++++
T Consensus       159 ~~~~~~v~l~~~~dtgDAMGaNmvn~~~eav~~~l~~~~~~~~~~~IlsN~~tdk~v~A~~~i  221 (417)
T cd00644         159 ADLGDFLSVHLLVDTKDAMGANIVNTMLEAVAPLLEEITGGEVLLRILSNYATERLVRAKVSI  221 (417)
T ss_pred             CCCCCeEEEEEEEEccchhcchhHHHHHHHHHHHHHHhCCcceeEEEeccCCCCceEEEEEEE
Confidence                99999999999999999999999999999999999988889999999999999999876


No 8  
>TIGR00532 HMG_CoA_R_NAD hydroxymethylglutaryl-CoA reductase, degradative. Most known examples of hydroxymethylglutaryl-CoA reductase are NADP-dependent (EC 1.1.1.34) from eukaryotes and archaea, involved in the biosynthesis of mevalonate from 3-hydroxy-3-methylglutaryl-CoA. This model, in contrast, is built from the two examples in completed genomes of sequences closely related to the degradative, NAD-dependent hydroxymethylglutaryl-CoA reductase of Pseudomonas mevalonii, a bacterium that can use mevalonate as its sole carbon source.
Probab=100.00  E-value=5.5e-61  Score=481.34  Aligned_cols=223  Identities=19%  Similarity=0.225  Sum_probs=201.6

Q ss_pred             HcCCCCccccccccCChHHHHHHHHHHHHHHhC---CccCCCCCCCCChhhhhhcccccceEEeeeceeeecceEEcCeE
Q 016586          144 IDGSIPSYALESKLGDCRRAAAIRREALQKMTG---RSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLLDGFE  220 (387)
Q Consensus       144 ~~g~i~~y~LE~~l~d~~rAv~iRR~~l~~~tg---~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlInG~e  220 (387)
                      .+|++++|++++..+-+++++++||++|++.++   .+++.|+...+|++...++||||+||+++||+||||||+|||++
T Consensus         3 ~~~~~~~~~~~~~~~f~~~s~~eR~~~l~~~~~l~~~~~~~l~~~~l~~~~a~~~~iEN~Ig~~~vPlGvAg~l~InG~~   82 (393)
T TIGR00532         3 KEGKIRFAELSKIFGFYHKSVEEKLKEIAEFAELSDEEVKAFFSNGANEDFAFDRMIENVIGTFEFPIGIAKNFKIDGKD   82 (393)
T ss_pred             cccccchhhhccccChhhcCHHHHHHHHHHhcCCCHHHHHHHhhCCCCHHHHhccCcccccceeeeceeEeccEEECCeE
Confidence            479999999999775599999999999999999   56778888778888777999999999999999999999999999


Q ss_pred             EeeeccccchhHHHHhhhhhhhhhccCCeEEEeccCeeecce-EEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhhc---
Q 016586          221 YSVPMATTEGCLVASTNRGCKAIYASGGAASMLLRDGMTRAP-IVRFASAMRASELKFFLEDPNNFETLAVVFNRSS---  296 (387)
Q Consensus       221 y~VPMATtEgsLVASanRGaKaI~~sGG~~t~v~~d~MTRap-v~~f~s~~dA~~l~~wie~p~n~~~L~~van~tS---  296 (387)
                      |+|||||||||||||+|||||+++.+||+++.+.+++|+||+ ++.++|..+   +.+|+.+  ++++|++++|+++   
T Consensus        83 y~VPMATeE~svVAs~srGak~~~~~GG~~~~v~~~~M~gqi~~~~~~~~~~---a~~~~~~--~~~~i~~~a~~~~~~~  157 (393)
T TIGR00532        83 YLIPIAIEEPSVVAAANFAAKIAEEADGFTSDGEGLGIIGQIQQIKIKNEKA---AKFEFLD--LGDEIIERAEECDPML  157 (393)
T ss_pred             EEEEeeeccccHHHHHHHHHHHHHhcCCeEEEEcCCceEEEEEEEecCCHHH---HHHHHHH--HHHHHHHHHHHhCHHH
Confidence            999999999999999999999999999999999999999887 566676544   3445543  8999999999986   


Q ss_pred             --CCcceeeEEEE----ecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCccccccceeEEEe
Q 016586          297 --RFARLQHIQCS----IAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLVIYLKSCSYILL  370 (387)
Q Consensus       297 --R~GgL~~I~~~----i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ctDrl~~A~~~~~  370 (387)
                        ||||+++|+++    +.|+++|+||.|||||||||||||+++|+++++|++++|+..+++||||||+||+++++++..
T Consensus       158 ~~rggg~~~i~~r~~~~~~g~~v~l~~~~dtgDAMGaNmvn~~~Eav~~~i~~~~~~~~~~~IlsN~~~dk~~~arg~~~  237 (393)
T TIGR00532       158 NNLGGGCKDIEARVIDIIEGGILILHIIVDTCDAMGANALNSIAEKVAEFIELEFGGECVLKIISNDAAEFTAKARAKAD  237 (393)
T ss_pred             HhhcCCeEEEEEEeeecccCCEEEEEEEEecccccccHHHHHHHHHHHHHHHHhCCCceEEEEecCccccceEEEEEEEE
Confidence              99999999986    568999999999999999999999999999999999998877799999999999998877764


Q ss_pred             e
Q 016586          371 S  371 (387)
Q Consensus       371 ~  371 (387)
                      -
T Consensus       238 ~  238 (393)
T TIGR00532       238 F  238 (393)
T ss_pred             h
Confidence            3


No 9  
>cd00365 HMG-CoA_reductase Hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR). Hydroxymethylglutaryl-coenzyme A (HMG-CoA) reductase (HMGR) is a tightly regulated enzyme, which catalyzes the synthesis of coenzyme A and mevalonate in isoprenoid synthesis. In mammals, this is the rate limiting committed step in cholesterol biosynthesis. Bacteria, such as Pseudomonas mevalonii, which rely solely on mevalonate for their carbon source, catalyze the reverse reaction, using an NAD-dependent HMGR to deacetylate mevalonate into 3-hydroxy-3-methylglutaryl-CoA. There are two classes of HMGR: class I enzymes which are found predominantly in eukaryotes and contain N-terminal membrane regions and class II enzymes which are found primarily in prokaryotes and are soluble as they lack the membrane region. With the exception of Archaeoglobus fulgidus, most archeae are assigned to class I, based on sequence similarity of the active site, even though they lack membrane regions. Yeast and human HMGR
Probab=100.00  E-value=1.3e-58  Score=462.24  Aligned_cols=213  Identities=19%  Similarity=0.194  Sum_probs=198.8

Q ss_pred             cccccccCChHHHHHHHHHHHHHHhCCccCCCCCCCCChhhhhhcccccceEEeeeceeeecceEEcCeEEeeeccccch
Q 016586          151 YALESKLGDCRRAAAIRREALQKMTGRSLQGLPLDGFDYDSILGQCCEMPVGYVQIPVGIAGPLLLDGFEYSVPMATTEG  230 (387)
Q Consensus       151 y~LE~~l~d~~rAv~iRR~~l~~~tg~sL~~L~~~~~Dy~~i~g~~iEN~IG~vqIPvGVAGPLlInG~ey~VPMATtEg  230 (387)
                      |.+|+.. |++|++++|++++. .++.+++.|+...+|++.. ++||||+||+++||+||||||+|||++|+||||||||
T Consensus         1 ~~~~~~~-~~~~~~~~~~~~~~-lt~~~~~~l~~~~l~~~~~-~~~iEN~IG~~~vPlGva~~l~InG~~y~VPmATeE~   77 (376)
T cd00365           1 PAFRTLS-PHAARLDHIGQLLG-LSHDDVQLLANAALPMDIA-NGMIENVIGTFELPYAVASNFQIDGRDVLVPLVTEEP   77 (376)
T ss_pred             Cchhhhc-CHHHHHHHHHHHhc-cChHHHhhccccCCCHHHH-hcccccccceeeeeeEEecCEEEcCcEEEEEeEecch
Confidence            4567776 79999999999985 5888899999988898876 9999999999999999999999999999999999999


Q ss_pred             hHHHHhhhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCCCCHHHHHHHHhhhcCC-----cceeeEE
Q 016586          231 CLVASTNRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDPNNFETLAVVFNRSSRF-----ARLQHIQ  305 (387)
Q Consensus       231 sLVASanRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p~n~~~L~~van~tSR~-----GgL~~I~  305 (387)
                      |||||+|||||+++.+||+++.+.+++|+|||+|.|.+  |+.++++|+++ +|+++|++++|+++++     ||+++|+
T Consensus        78 svVAs~srgak~i~~~GG~~~~~~~~~m~~q~~~~~~~--d~~~~~~~i~~-~~~~~i~~~a~~~~~~~~~rgggl~~i~  154 (376)
T cd00365          78 SIVAAASYMAKLARAGGGFTTSSSAPLMHAQVQIVLIQ--DPLNAKLSLLR-SGKDEIIELANRKDQLLNSLGGGCRDIE  154 (376)
T ss_pred             hHHHHHHHHHHHHHhcCCEEEEECCCcEEeeeEEEeCC--CHHHHHHHHHh-hhHHHHHHHHhhhCcchhccCccceEEE
Confidence            99999999999999999999999999999999999987  99999999983 3999999999999666     9999999


Q ss_pred             EEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCe--EEEEeCCCccccccceeEEE
Q 016586          306 CSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMD--VIGISGMCLVIYLKSCSYIL  369 (387)
Q Consensus       306 ~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~--visIsgN~ctDrl~~A~~~~  369 (387)
                      +++.|+++++||.+||||||||||||+++|+++++|++++++.+  .++|++|+++.|+++|++++
T Consensus       155 ~~~~~~~v~l~~~~dtgDAMGaNmvn~~~eav~~~i~~~~~~~~~~~~si~~n~~~g~~v~A~~~I  220 (376)
T cd00365         155 VHTFGPMLVAHLIVDVGDAMGANMINTMAEAVAPLMEAYTGGMQVRLRSLSNLTGDGRLARAQARI  220 (376)
T ss_pred             EEeeCCEEEEEEEEEccchhhchhHHHHHHHHHHHHHHHcCCCceEEEEEecccCCCcEEEEEEEE
Confidence            99999999999999999999999999999999999999888764  48899999999999999986


No 10 
>PF14150 YesK:  YesK-like protein
Probab=74.68  E-value=3.4  Score=34.50  Aligned_cols=57  Identities=21%  Similarity=0.475  Sum_probs=35.9

Q ss_pred             hhhhhhHHHHHHHHHHHHHHHHhHhhhhcCCCccccc-hHHHHHHHHHHHHHHHHHhhcCcc
Q 016586           35 LYLTNAIFFTLFFSVAYYLLHRWREKIRNSTPLHVVT-LSEIAAIVSLIASFIYLLGFFGID   95 (387)
Q Consensus        35 ~~~~n~~~~~~f~~~~~~l~~~wr~kir~s~plh~~~-~~e~~~~~~~~as~~yl~~ffgi~   95 (387)
                      +|+.-++++.+.|+++|++=+||++|-    |--.++ ..=+..+..++-|.++.=|+.|++
T Consensus         1 f~llg~~~~ii~f~~S~~lr~r~p~k~----~~~il~~ililis~~~v~~S~f~vGGweGmg   58 (81)
T PF14150_consen    1 FYLLGIVTFIIVFGVSVLLRKRFPKKQ----PEIILPLILILISLLTVLISIFLVGGWEGMG   58 (81)
T ss_pred             CcHHHHHHHHHHHHHHHHHHHhCCCcc----hhHHHHHHHHHHHHHHHHHHHheEcchhhhh
Confidence            367778999999999999999999984    443333 122233333444444444555544


No 11 
>PF11947 DUF3464:  Protein of unknown function (DUF3464);  InterPro: IPR021855  This family of proteins are functionally uncharacterised. This protein is found in bacteria and eukaryotes. Proteins in this family are typically between 137 to 196 amino acids in length. 
Probab=72.13  E-value=8.7  Score=35.45  Aligned_cols=70  Identities=29%  Similarity=0.639  Sum_probs=42.0

Q ss_pred             CCCCCccchhhhhhh------HHH-------HHHHHHHHHHHHHhHhhhhcCCCccccchHHHHHHHH-HHHHHHHHHhh
Q 016586           26 KASDALPLPLYLTNA------IFF-------TLFFSVAYYLLHRWREKIRNSTPLHVVTLSEIAAIVS-LIASFIYLLGF   91 (387)
Q Consensus        26 ~~~~~lplp~~~~n~------~~~-------~~f~~~~~~l~~~wr~kir~s~plh~~~~~e~~~~~~-~~as~~yl~~f   91 (387)
                      ...+...+|.+++|=      +|.       ..+|.+.|||..+  .         .+   |+-..+. ++.-+.|.+|+
T Consensus        46 ~~~~~~~IP~~Vs~RM~rRm~~~~GiP~~lG~~~f~~~y~l~~~--~---------~~---dvP~~~~~~~S~~~Fg~gl  111 (153)
T PF11947_consen   46 RDEDDSAIPEVVSNRMLRRMAVFVGIPTALGVAVFVVFYYLKSR--Q---------IV---DVPPWAVLLVSLVFFGLGL  111 (153)
T ss_pred             ccccccccCHHHHHHHHHHHHHHhchHHHHHHHHHHHHHHHHhc--c---------cc---ccCchHHHHHHHHHHHHHH
Confidence            356777888899983      222       2367777888764  1         12   2323333 34445688999


Q ss_pred             cCccccccccCCCCCccCCcccCC
Q 016586           92 FGIDFVQSFISRATPEAWDLEEDD  115 (387)
Q Consensus        92 fgi~~~~~~~~~~~~~~~~~~~~~  115 (387)
                      .||.|-  +++    .+||.+.+.
T Consensus       112 lGisYG--ilS----aSWD~~r~G  129 (153)
T PF11947_consen  112 LGISYG--ILS----ASWDPEREG  129 (153)
T ss_pred             Hhhhhh--hcc----cccCCCCCC
Confidence            999984  222    248876544


No 12 
>PRK09509 fieF ferrous iron efflux protein F; Reviewed
Probab=71.43  E-value=16  Score=35.76  Aligned_cols=68  Identities=13%  Similarity=0.192  Sum_probs=51.3

Q ss_pred             CCHHHHHHHHhhhcCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeE
Q 016586          283 NNFETLAVVFNRSSRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDV  350 (387)
Q Consensus       283 ~n~~~L~~van~tSR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~v  350 (387)
                      +..+++++..++....-++.+++.+..|+..++.+.+...+.|-..=.....+++-+.|++++|+.++
T Consensus       213 ~~~~~I~~~i~~~~~v~~v~~l~~~~~G~~~~v~v~i~v~~~~~~~e~h~i~~~ie~~l~~~~~~~~v  280 (299)
T PRK09509        213 EERQEIIDIVTSWPGVSGAHDLRTRQSGPTRFIQLHLEMEDNLPLVQAHMIADQVEQALLRRFPGSDV  280 (299)
T ss_pred             HHHHHHHHHHHhCCCCcCceeeeeEeeCCeEEEEEEEEECCCCCHHHHHHHHHHHHHHHHHhCCCCEE
Confidence            34567777777765556778899999999988888888877775444556678888888888887765


No 13 
>TIGR01297 CDF cation diffusion facilitator family transporter. This model describes a broadly distributed family of transporters, a number of which have been shown to transport divalent cations of cobalt, cadmium and/or zinc. The family has six predicted transmembrane domains. Members of the family are variable in length because of variably sized inserts, often containing low-complexity sequence.
Probab=60.37  E-value=63  Score=30.43  Aligned_cols=66  Identities=6%  Similarity=0.145  Sum_probs=49.6

Q ss_pred             CCHHHHHHHHhhhcCCcceeeEEEEecCC-eEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCC
Q 016586          283 NNFETLAVVFNRSSRFARLQHIQCSIAGK-NLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDM  348 (387)
Q Consensus       283 ~n~~~L~~van~tSR~GgL~~I~~~i~G~-~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~  348 (387)
                      +..+++++..++....-++.+++.+..|+ .+++.+.+...+.|-..-..+..+.+-+.|++++|+.
T Consensus       192 ~~~~~i~~~i~~~~~v~~v~~~~~~~~G~~~~~v~~~v~v~~~~~~~~ah~i~~~i~~~i~~~~~~v  258 (268)
T TIGR01297       192 EDLEEIKKAILSIPGVKGVHDLHIWRIGPGKLFLDVHVVVDPDLDLKQAHDIALEIEREILKRHPGI  258 (268)
T ss_pred             ccHHHHHHHHhcCCCcccceEeEEEEcCCCCEEEEEEEEECCCCChhHHHHHHHHHHHHHHHhcCCC
Confidence            45677777777554455778899999998 8888888888776655555677778888787777776


No 14 
>PF14715 FixP_N:  N-terminal domain of cytochrome oxidase-cbb3, FixP 
Probab=58.38  E-value=12  Score=28.66  Aligned_cols=25  Identities=32%  Similarity=0.547  Sum_probs=22.3

Q ss_pred             ccchhhhhhhHHHHHHHHHHHHHHH
Q 016586           31 LPLPLYLTNAIFFTLFFSVAYYLLH   55 (387)
Q Consensus        31 lplp~~~~n~~~~~~f~~~~~~l~~   55 (387)
                      =|||-...-.++.|..|++.|+++.
T Consensus        18 nplP~ww~~~f~~tivfa~~Y~~~y   42 (51)
T PF14715_consen   18 NPLPRWWLWLFYGTIVFAVGYLVLY   42 (51)
T ss_pred             CCCCHHHHHHHHHHHHHHHHHHHHH
Confidence            4899999999999999999999874


No 15 
>PF14012 DUF4229:  Protein of unknown function (DUF4229)
Probab=53.71  E-value=20  Score=28.71  Aligned_cols=30  Identities=23%  Similarity=0.557  Sum_probs=22.3

Q ss_pred             chhhhhhhHHHHHHHHHHHHHHHHhHhhhh
Q 016586           33 LPLYLTNAIFFTLFFSVAYYLLHRWREKIR   62 (387)
Q Consensus        33 lp~~~~n~~~~~~f~~~~~~l~~~wr~kir   62 (387)
                      .|+.+.=.+-...=+..+|++++|||++.-
T Consensus        30 ~p~~~~~l~A~vis~~lS~~ll~~~R~~~~   59 (69)
T PF14012_consen   30 VPLLVAALLALVISMPLSYVLLRRLRDRAS   59 (69)
T ss_pred             chHHHHHHHHHHHHHHHHHHHHHHHHHHHH
Confidence            566666555556667789999999999863


No 16 
>PF01545 Cation_efflux:  Cation efflux family;  InterPro: IPR002524 Members of this family are integral membrane proteins, that are found to increase tolerance to divalent metal ions such as cadmium, zinc, and cobalt. These proteins are considered to be efflux pumps that remove these ions from cells [, ], however others are implicated in ion uptake []. The family has six predicted transmembrane domains. Members of the family are variable in length because of variably sized inserts, often containing low-complexity sequence.; GO: 0008324 cation transmembrane transporter activity, 0006812 cation transport, 0055085 transmembrane transport, 0016021 integral to membrane; PDB: 3BYP_B 3BYR_A 3H90_B 2QFI_B 2ZZT_A.
Probab=47.30  E-value=1.1e+02  Score=28.90  Aligned_cols=66  Identities=9%  Similarity=0.196  Sum_probs=52.9

Q ss_pred             CHHHHHHHHhhhcCCcceeeEEEEecCC-eEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCe
Q 016586          284 NFETLAVVFNRSSRFARLQHIQCSIAGK-NLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMD  349 (387)
Q Consensus       284 n~~~L~~van~tSR~GgL~~I~~~i~G~-~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~  349 (387)
                      ..+++++..++....-++.+++.+..|+ .+++.+.+..-..|..-=.....+++.+.+++++|+..
T Consensus       206 ~~~~i~~~i~~~~~v~~v~~~~~~~~g~~~~~v~i~v~v~~~~~v~~~~~i~~~i~~~l~~~~~~i~  272 (284)
T PF01545_consen  206 LVEKIRRIIESVPGVIEVHDLRVWQVGRNKYVVEIHVQVDPDMSVEEAHEIRERIEKRLREKFPGIY  272 (284)
T ss_dssp             HHHHHHHHHHHTSS-SEEEEEEEEEETT-EEEEEEEEEETTTSBHHHHHHHHHHHHHHHHHHSTTCE
T ss_pred             chhHHHHhhccCCceEeccceEEEEecCCcEEEEEEEEeCCCCCHHHHHHHHHHHHHHHHHhCCCcE
Confidence            3467777777767777999999999998 89999999988877655566788899999999998876


No 17 
>COG0053 MMT1 Predicted Co/Zn/Cd cation transporters [Inorganic ion transport and metabolism]
Probab=45.13  E-value=98  Score=30.98  Aligned_cols=78  Identities=12%  Similarity=0.253  Sum_probs=60.1

Q ss_pred             CCCHHHHHHHHhhhcCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCCcc
Q 016586          282 PNNFETLAVVFNRSSRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMCLV  359 (387)
Q Consensus       282 p~n~~~L~~van~tSR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~ct  359 (387)
                      |+..++++++..+...--++.+++.+..|+..+++++....+.|.-==.-..++.+-+.|++.+|+..-+-|.-+=+.
T Consensus       214 ~~~~~~i~~~i~~~~~V~~v~~lr~R~~G~~~~id~~i~v~~~ls~~eah~I~~~ie~~i~~~~~~~~~v~IhveP~~  291 (304)
T COG0053         214 PEDLEKIRAIILSVPGVKGVHDLRTRKSGSRIFIDVHIEVDPDLSLEEAHEIADEVEKRIKKEFPKVADVTIHVEPLG  291 (304)
T ss_pred             HHHHHHHHHHHhcCCcceeeecceeeeeCCeEEEEEEEEECCCCChHHHHHHHHHHHHHHHHhcCCCceEEEEecCCc
Confidence            445667777777777667889999999999999999999999888777777888888888888885444444444333


No 18 
>COG1862 YajC Preprotein translocase subunit YajC [Intracellular trafficking and secretion]
Probab=41.63  E-value=21  Score=30.59  Aligned_cols=27  Identities=26%  Similarity=0.363  Sum_probs=21.7

Q ss_pred             hhhhhhhHHHHHHHHHHHHHHHHhHhh
Q 016586           34 PLYLTNAIFFTLFFSVAYYLLHRWREK   60 (387)
Q Consensus        34 p~~~~n~~~~~~f~~~~~~l~~~wr~k   60 (387)
                      +-.++..+.|.++|.+.|||+.|=..|
T Consensus         6 ~~~~~~ll~~vl~~~ifyFli~RPQrK   32 (97)
T COG1862           6 GSGLVLLLPLVLIFAIFYFLIIRPQRK   32 (97)
T ss_pred             cccHHHHHHHHHHHHHHHHhhcCHHHH
Confidence            457788899999999999998765433


No 19 
>PRK13454 F0F1 ATP synthase subunit B'; Provisional
Probab=38.93  E-value=29  Score=32.03  Aligned_cols=21  Identities=14%  Similarity=0.501  Sum_probs=17.0

Q ss_pred             hhhhhHHHHHHHHHHHHHHHH
Q 016586           36 YLTNAIFFTLFFSVAYYLLHR   56 (387)
Q Consensus        36 ~~~n~~~~~~f~~~~~~l~~~   56 (387)
                      |..+.++|.++|.+.||+|.+
T Consensus        30 ~~~q~~~~lI~F~iL~~ll~k   50 (181)
T PRK13454         30 FPNQIFWLLVTLVAIYFVLTR   50 (181)
T ss_pred             cchHHHHHHHHHHHHHHHHHH
Confidence            445777888889999999987


No 20 
>COG5548 Small integral membrane protein [Function unknown]
Probab=38.05  E-value=46  Score=28.99  Aligned_cols=48  Identities=31%  Similarity=0.410  Sum_probs=34.5

Q ss_pred             hhhhhHHHHHHHHHHHHHHHHhHhhhhcCCCccccchHHHHHHHHHHHHHHHHHhhcCcccccc
Q 016586           36 YLTNAIFFTLFFSVAYYLLHRWREKIRNSTPLHVVTLSEIAAIVSLIASFIYLLGFFGIDFVQS   99 (387)
Q Consensus        36 ~~~n~~~~~~f~~~~~~l~~~wr~kir~s~plh~~~~~e~~~~~~~~as~~yl~~ffgi~~~~~   99 (387)
                      -|-.++||.+-|.++-||++|=|+.               --+.++.+|.. ||+||++-|+-|
T Consensus        30 SL~sG~~~G~~~~~A~yL~~~g~~~---------------Gl~~A~~~s~~-Ll~~~~~R~~~s   77 (105)
T COG5548          30 SLLSGVFSGLLLFVAAYLQLQGQTW---------------GLILATVVSAA-LLVFFALRLVRS   77 (105)
T ss_pred             hhHHHHHHhHHHHHHHHHHHcCccc---------------CeehHHHHHHH-HHHhcchhcccc
Confidence            3567899999999999999875442               11334455543 889999999876


No 21 
>PF09551 Spore_II_R:  Stage II sporulation protein R (spore_II_R);  InterPro: IPR014202  This entry is designated stage II sporulation protein R. A comparative genome analysis of all sequenced genomes of Firmicutes shows that the proteins are strictly conserved among the sub-set of endospore-forming species. SpoIIR is a signalling protein that links the activation of sigma E to the transcriptional activity of sigma F during sporulation [, ].
Probab=37.66  E-value=69  Score=28.94  Aligned_cols=58  Identities=17%  Similarity=0.251  Sum_probs=40.4

Q ss_pred             CChhcHHHHHHHhcCCCCHHHHHHHHhhh-cCCcceeeEEEEe----------------cCCeEEEEEEeeccccccch
Q 016586          267 ASAMRASELKFFLEDPNNFETLAVVFNRS-SRFARLQHIQCSI----------------AGKNLYIRFCCTTGDAMGMN  328 (387)
Q Consensus       267 ~s~~dA~~l~~wie~p~n~~~L~~van~t-SR~GgL~~I~~~i----------------~G~~l~lrf~~dTGDAMGmN  328 (387)
                      ....+..+.++|+++  |.++|.+++++. ..-|.-.++++.+                +|.+=-+|  +.-|++-|.|
T Consensus        38 ~~~~~~~ea~~~i~~--~~~~Ie~~A~~~l~~~G~~y~v~v~~~~~~FPtK~YG~~~~PaG~YeAlr--I~IG~g~G~N  112 (130)
T PF09551_consen   38 SQAKSKEEAREVIRE--NLPEIEQIAEEVLAEEGYDYPVKVELGRFYFPTKTYGDIVLPAGEYEALR--ITIGEGKGHN  112 (130)
T ss_pred             ccCCCHHHHHHHHHH--hHHHHHHHHHHHHHHhCCCCcEEEEEEeeeCCCceECCEeccCCceEEEE--EEecCccCcc
Confidence            445678888999986  999999999977 4444444555543                25554444  5568888988


No 22 
>PF09580 Spore_YhcN_YlaJ:  Sporulation lipoprotein YhcN/YlaJ (Spore_YhcN_YlaJ);  InterPro: IPR019076  This entry contains YhcN and YlaJ, which are predicted lipoproteins that have been detected as spore proteins but not vegetative proteins in Bacillus subtilis. Both appear to be expressed under control of the RNA polymerase sigma-G factor. The YlaJ-like members of this family have a low-complexity, strongly acidic, 40-residue C-terminal domain. 
Probab=37.61  E-value=2.5e+02  Score=25.14  Aligned_cols=70  Identities=11%  Similarity=0.047  Sum_probs=56.5

Q ss_pred             CHHHHHHHHhhhcCCcceeeEEEEecCCeEEEEEEeeccccccchhhhHHHHHHHHHHHhhCCCCeEEEEeCCC
Q 016586          284 NFETLAVVFNRSSRFARLQHIQCSIAGKNLYIRFCCTTGDAMGMNMVSKGVQNVLDFLQNDFPDMDVIGISGMC  357 (387)
Q Consensus       284 n~~~L~~van~tSR~GgL~~I~~~i~G~~l~lrf~~dTGDAMGmNMVnk~~E~v~~~I~~~fp~~~visIsgN~  357 (387)
                      +.+.-..+++...+--++.+..+.+.|+.+||=+..+ .+   .....+.-+.|.+.+++.+|++.-+-+|+|.
T Consensus        73 ~~~~a~~i~~~v~~~~~V~~A~vvv~~~~a~Vav~~~-~~---~~~~~~i~~~V~~~v~~~~p~~~~V~Vs~D~  142 (177)
T PF09580_consen   73 RQQLADRIANRVKKVPGVEDATVVVTDDNAYVAVDLD-FN---RFNTKKIKKKVEKAVKSADPRIYNVYVSTDP  142 (177)
T ss_pred             HHHHHHHHHHHHhcCCCceEEEEEEECCEEEEEEEec-cc---ccchhHHHHHHHHHHHHhCCCccEEEEEcCH
Confidence            3444456667778888999999999999999999888 44   5778888899999999989987777777763


No 23 
>PF01578 Cytochrom_C_asm:  Cytochrome C assembly protein;  InterPro: IPR002541 This entry consists of various proteins involved in cytochrome c assembly from mitochondria and bacteria; CycK from Rhizobium leguminosarum [], CcmC from Escherichia coli and Paracoccus denitrificans [, ] and orf240 from Triticum aestivum (Wheat) mitochondria []. The members of this family are probably integral membrane proteins with six predicted transmembrane helices that may comprise the membrane component of an ABC (ATP binding cassette) transporter complex. This transporter may be necessary for transport of some component needed for cytochrome c assembly. One member, R. leguminosarum CycK, contains a putative haem-binding motif []. Wheat orf240 also contains a putative haem-binding motif and is a proposed ABC transporter with c-type haem as its proposed substrate []. However it seems unlikely that all members of this family transport haem or c-type apocytochromes because P. denitrificans CcmC transports neither [].; GO: 0006461 protein complex assembly, 0008535 respiratory chain complex IV assembly, 0016020 membrane
Probab=37.51  E-value=67  Score=29.60  Aligned_cols=50  Identities=28%  Similarity=0.508  Sum_probs=35.3

Q ss_pred             hhHHHHHHHHHHHHHHHHhHhhhhcCCCccccchHHHHHHHHHHHHHHHHHhhcCcccccc
Q 016586           39 NAIFFTLFFSVAYYLLHRWREKIRNSTPLHVVTLSEIAAIVSLIASFIYLLGFFGIDFVQS   99 (387)
Q Consensus        39 n~~~~~~f~~~~~~l~~~wr~kir~s~plh~~~~~e~~~~~~~~as~~yl~~ffgi~~~~~   99 (387)
                      +.-+++.++-.+|+.+|.++ +-|          ...++..+.++-++-++++||+.++..
T Consensus       162 ~~sli~Wl~y~~~lh~r~~~-~~~----------gr~~a~~~i~gf~~~~~~~~gv~~~~~  211 (214)
T PF01578_consen  162 VWSLITWLVYGAYLHLRSWK-GWR----------GRRAAYLSIIGFLLLLLSYFGVNLLLE  211 (214)
T ss_pred             HHHHHHHHHHHHHHHHHHhh-chh----------hHHHHHHHHHHHHHHHHHHHHHHHhcC
Confidence            34455656666666676555 444          235678888888999999999998764


No 24 
>TIGR03063 srtB_target sortase B cell surface sorting signal. Two different classes of sorting signal, both analogous to the sortase A signal LPXTG, may be recognized by the sortase SrtB. These are given as NXZTN and NPKXZ. Proteins sorted by this class of sortase are less common than the sortase A and LPXTG system. This model describes a number of cell surface protein C-terminal regions from Gram-positive bacteria that appear to be sortase B (SrtB) sorting signals.
Probab=36.12  E-value=37  Score=23.43  Aligned_cols=26  Identities=27%  Similarity=0.671  Sum_probs=17.0

Q ss_pred             CCCCCCccchhhhhhhHHHHHHHHHHHHHHHH
Q 016586           25 PKASDALPLPLYLTNAIFFTLFFSVAYYLLHR   56 (387)
Q Consensus        25 ~~~~~~lplp~~~~n~~~~~~f~~~~~~l~~~   56 (387)
                      ||.||.-|+=+      |.++|-++..||++|
T Consensus         2 PkT~D~a~i~l------y~~l~~~s~~~Li~k   27 (29)
T TIGR03063         2 PKTGDSAQIGL------YAVLFLGSGLFLIRK   27 (29)
T ss_pred             CCCccchhHHH------HHHHHHHHHHHHhhc
Confidence            46788877644      445666666777755


No 25 
>PF07045 DUF1330:  Protein of unknown function (DUF1330);  InterPro: IPR010753 This family consists of several hypothetical bacterial proteins of around 90 residues in length. The function of this family is unknown.; PDB: 2FIU_B 3HHL_A 3DCA_D 3LO3_I.
Probab=35.47  E-value=33  Score=26.40  Aligned_cols=43  Identities=19%  Similarity=0.229  Sum_probs=26.4

Q ss_pred             hhhhhhhhccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcCC
Q 016586          237 NRGCKAIYASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLEDP  282 (387)
Q Consensus       237 nRGaKaI~~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~p  282 (387)
                      ..|.|++..+|...+.--...-.+-.++.|+|.+   .++.|.++|
T Consensus        17 ~~GG~~l~~~~~~~~leG~~~~~~~viieFPs~~---aa~~~~~sp   59 (65)
T PF07045_consen   17 KYGGRVLARGGEPEVLEGDWDPDRVVIIEFPSME---AAKAWYNSP   59 (65)
T ss_dssp             HTT-EEEEECEEEEEEEST-SSSEEEEEEESSHH---HHHHHHCSH
T ss_pred             HcCCEEEEECCceeEEecCCCCCeEEEEECCCHH---HHHHHHCCH
Confidence            4677777664444433332555777899999863   356788875


No 26 
>PF11360 DUF3110:  Protein of unknown function (DUF3110);  InterPro: IPR021503  This family of proteins has no known function. 
Probab=34.00  E-value=57  Score=27.28  Aligned_cols=34  Identities=24%  Similarity=0.281  Sum_probs=30.1

Q ss_pred             ccCCeEEEeccCeeecceEEEeCChhcHHHHHHHhcC
Q 016586          245 ASGGAASMLLRDGMTRAPIVRFASAMRASELKFFLED  281 (387)
Q Consensus       245 ~sGG~~t~v~~d~MTRapv~~f~s~~dA~~l~~wie~  281 (387)
                      ..+|+.|.-.++   |..++.|.+-+||.+...-|+.
T Consensus        11 ~~eGI~si~~~~---~~~Vl~FE~edDA~RYa~lLEA   44 (86)
T PF11360_consen   11 ETEGIYSIQNKD---RNVVLMFEDEDDAERYAGLLEA   44 (86)
T ss_pred             CCCcEEEEEeCC---CCEEEEEccHHHHHHHHHHHHh
Confidence            578999988777   8889999999999999999975


No 27 
>PF13858 DUF4199:  Protein of unknown function (DUF4199)
Probab=33.73  E-value=1.4e+02  Score=26.07  Aligned_cols=49  Identities=18%  Similarity=0.380  Sum_probs=34.8

Q ss_pred             hhHHHHHHHHHHHHHHHHhHhhhhcCCCccccchHHH---HHHHHHHHHHHHHHhh
Q 016586           39 NAIFFTLFFSVAYYLLHRWREKIRNSTPLHVVTLSEI---AAIVSLIASFIYLLGF   91 (387)
Q Consensus        39 n~~~~~~f~~~~~~l~~~wr~kir~s~plh~~~~~e~---~~~~~~~as~~yl~~f   91 (387)
                      |-+.+.+.+...|+.+++-|||..+.    ..++.+-   .-.++++|++++.+.-
T Consensus        34 ~~~~~~~~~~~i~~~i~~~R~~~~~g----~isf~~a~~~g~~~~~ia~li~~v~~   85 (163)
T PF13858_consen   34 GILSMVITIIFIYFAIRRYRKKYNGG----FISFGQAFKVGFLISLIAGLISAVFQ   85 (163)
T ss_pred             HHHHHHHHHHHHHHHHHHHHHHccCC----CeeHHHHHHHHHHHHHHHHHHHHHHH
Confidence            55556677777899999999987665    4777773   4456677777765543


No 28 
>TIGR00782 ccoP cytochrome c oxidase, cbb3-type, subunit III. This model describes a di-heme subunit of approximately 26 kDa of the cbb3 type copper and heme-containing cytochrome oxidase.
Probab=32.99  E-value=39  Score=33.17  Aligned_cols=25  Identities=20%  Similarity=0.442  Sum_probs=22.9

Q ss_pred             ccchhhhhhhHHHHHHHHHHHHHHH
Q 016586           31 LPLPLYLTNAIFFTLFFSVAYYLLH   55 (387)
Q Consensus        31 lplp~~~~n~~~~~~f~~~~~~l~~   55 (387)
                      =|||-+..-.++.|+.|+++|+++-
T Consensus        24 n~~P~ww~~~f~~~i~~~~~y~~~y   48 (285)
T TIGR00782        24 NPLPRWWLWTFYATIVWGFGYLVAY   48 (285)
T ss_pred             CCCCHHHHHHHHHHHHHHHHHHHHh
Confidence            4899999999999999999998875


No 29 
>TIGR02837 spore_II_R stage II sporulation protein R. A comparative genome analysis of all sequenced genomes of shows a number of proteins conserved strictly among the endospore-forming subset of the Firmicutes. This protein, a member of this panel, is designated stage II sporulation protein R.
Probab=29.13  E-value=91  Score=29.40  Aligned_cols=56  Identities=13%  Similarity=0.169  Sum_probs=37.7

Q ss_pred             hhcHHHHHHHhcCCCCHHHHHHHHhhh-cCCcceeeEEEEe----------------cCCeEEEEEEeeccccccch
Q 016586          269 AMRASELKFFLEDPNNFETLAVVFNRS-SRFARLQHIQCSI----------------AGKNLYIRFCCTTGDAMGMN  328 (387)
Q Consensus       269 ~~dA~~l~~wie~p~n~~~L~~van~t-SR~GgL~~I~~~i----------------~G~~l~lrf~~dTGDAMGmN  328 (387)
                      ..+..+.++|+.+  |.++|.+++++. ..-|.=.++++.+                +|++=-  +.+.-|++-|+|
T Consensus        75 ~~s~~ea~~~i~~--~l~~Ie~~a~~~l~~~G~~y~v~v~~~~~~FPtK~YG~~~~PaG~YeA--lrI~IG~g~G~N  147 (168)
T TIGR02837        75 LKSLEEARRVIRE--NLPEIERIAESVIKAEGADYKVRVELGKYSFPTKLYGNIVLPAGEYEA--LRILIGEGAGAN  147 (168)
T ss_pred             CCCHHHHHHHHHH--hhHHHHHHHHHHHHHhCCCCCeEEEEEEEeCCCcccCCEeccCCceEE--EEEEecCcCCcc
Confidence            4466777889986  999999999877 4444444555543                244444  445678888998


No 30 
>PF07074 TRAP-gamma:  Translocon-associated protein, gamma subunit (TRAP-gamma);  InterPro: IPR009779 This family consists of several eukaryotic translocon-associated protein, gamma subunit (TRAP-gamma) sequences. The translocation site (translocon), at which nascent polypeptides pass through the endoplasmic reticulum membrane, contains a component previously called 'signal sequence receptor' that is now renamed as 'translocon-associated protein' (TRAP). The TRAP complex is comprised of four membrane proteins alpha, beta, gamma and delta, which are present in a stoichiometric relation, and are genuine neighbours in intact microsomes. The gamma subunit is predicted to span the membrane four times [].; GO: 0006613 cotranslational protein targeting to membrane, 0005784 Sec61 translocon complex, 0030176 integral to endoplasmic reticulum membrane
Probab=29.06  E-value=36  Score=32.08  Aligned_cols=49  Identities=24%  Similarity=0.644  Sum_probs=34.9

Q ss_pred             CCCCccchhhhhhhHHHHHHHHHHHHHHHHhHhhhhcCCCcc--ccchHHHHHHHHHH
Q 016586           27 ASDALPLPLYLTNAIFFTLFFSVAYYLLHRWREKIRNSTPLH--VVTLSEIAAIVSLI   82 (387)
Q Consensus        27 ~~~~lplp~~~~n~~~~~~f~~~~~~l~~~wr~kir~s~plh--~~~~~e~~~~~~~~   82 (387)
                      .-.|.-+.+..+|++|.++..-.++|+|++|-       |--  +++..=-+.+++++
T Consensus       117 d~Ea~~~SifynNalFl~l~i~~SF~~lk~~~-------p~~Nyi~S~~~asgl~all  167 (170)
T PF07074_consen  117 DYEATTFSIFYNNALFLALVIVFSFYLLKNFS-------PVFNYIFSMSGASGLVALL  167 (170)
T ss_pred             hhhheeeehhhhchHHHHHHHHHHHHHHccCC-------ccceeeehHHHHHHHHHHH
Confidence            34577788899999999999999999998774       433  44444444455544


No 31 
>PF11834 DUF3354:  Domain of unknown function (DUF3354);  InterPro: IPR021789 Potassium channels take part in important processes of higher plants, including opening and closing of stomatal pores and leaf movement. Inward rectifying potassium (K(+)in) channels play an important role in turgor regulation and ion uptake in higher plants. All of them comprise, from their N-terminal to their C-terminal ends: a short hydrophilic region, a hydrophobic region structurally analogous and partially homologous to the transmembrane domain of voltage-gated animal channels from the Shaker superfamily, a putative cyclic nucleotide-binding domain, and a conserved C-terminal KHA domain. Between these last two regions, some of them (AKT1, AKT2 and SKT1) contain an ankyrin-repeat domain with six repeats homologous to those of human erythrocyte ankyrin.  This entry represents the KHA domain which is unique to plant K(+)in channels. The KHA domain contains two high-homology blocks enriched for hydrophobic and acidic residues, respectively. The KHA domain is essential for interaction of plant K(+)in channels. The KHA domain mediates tetramerization and/or stabilisation of the heteromers [, , ]. 
Probab=27.88  E-value=1.4e+02  Score=24.06  Aligned_cols=34  Identities=24%  Similarity=0.436  Sum_probs=27.4

Q ss_pred             CCCHHHHHHHHhhh----------cCCcceeeEEEEecCCeEEE
Q 016586          282 PNNFETLAVVFNRS----------SRFARLQHIQCSIAGKNLYI  315 (387)
Q Consensus       282 p~n~~~L~~van~t----------SR~GgL~~I~~~i~G~~l~l  315 (387)
                      |...++|.++|.+.          .-+|.+.+|.+--.|+.||+
T Consensus        25 P~SleeLl~ia~~kfg~~~~~v~~~dgaeIdDI~~IRDgD~L~~   68 (69)
T PF11834_consen   25 PDSLEELLKIASEKFGFSATKVLNEDGAEIDDIDVIRDGDHLYL   68 (69)
T ss_pred             CccHHHHHHHHHHHhCCCceEEEcCCCCEEeEEEEEEcCCEEEE
Confidence            45778888777653          66889999999999999986


No 32 
>PRK13108 prolipoprotein diacylglyceryl transferase; Reviewed
Probab=27.12  E-value=1.2e+02  Score=32.44  Aligned_cols=43  Identities=14%  Similarity=0.320  Sum_probs=31.9

Q ss_pred             HHHHHHHHHHHhHhhhhcCCCc--cccchHHHHHHHHHHHHHHHH
Q 016586           46 FFSVAYYLLHRWREKIRNSTPL--HVVTLSEIAAIVSLIASFIYL   88 (387)
Q Consensus        46 f~~~~~~l~~~wr~kir~s~pl--h~~~~~e~~~~~~~~as~~yl   88 (387)
                      +|.+.|.+.|=|-|-+|.....  .-++.+.+.+++.+++.++|+
T Consensus       228 lYli~Ygi~RF~iEflR~d~~~~~~gl~~~Q~lSl~~il~gl~~~  272 (460)
T PRK13108        228 FYVAFYCAGRFCVELLRDDPATLIAGIRINSFTSTFVFIGAVVYI  272 (460)
T ss_pred             HHHHHHHHHHHHhhhhccCchhhhcCccHHHHHHHHHHHHHHHHH
Confidence            3456677788888999987632  237889998888888887665


No 33 
>PF06129 Chordopox_G3:  Chordopoxvirus G3 protein;  InterPro: IPR010367 This family consists of several poxvirus specific G3 proteins. The function of this family is unknown.
Probab=25.71  E-value=57  Score=28.78  Aligned_cols=21  Identities=33%  Similarity=0.695  Sum_probs=17.7

Q ss_pred             hhhhhHHHHHHHHHHHHHHHH
Q 016586           36 YLTNAIFFTLFFSVAYYLLHR   56 (387)
Q Consensus        36 ~~~n~~~~~~f~~~~~~l~~~   56 (387)
                      ++.|-+||.+|...+||+..-
T Consensus         2 ~ll~l~ffi~Fl~~~Y~~~y~   22 (109)
T PF06129_consen    2 WLLYLIFFILFLVLCYFFNYY   22 (109)
T ss_pred             cHHHHHHHHHHHHHHHHHhhc
Confidence            567889999999999998754


No 34 
>TIGR02206 intg_mem_TP0381 conserved hypothetical integral membrane protein TIGR02206. This model represents a family of hydrophobic proteins with seven predicted transmembrane alpha helices. Members are found in Bacillus subtilis (ywaF), TP0381 from Treponema pallidum (TP0381), Streptococcus pyogenes, Rhodococcus erythropolis, etc.
Probab=24.91  E-value=1.4e+02  Score=28.55  Aligned_cols=41  Identities=29%  Similarity=0.520  Sum_probs=34.5

Q ss_pred             hhhcCCCccccchHHHHHHHHHH------HHHHHHHhhcCccccccccC
Q 016586           60 KIRNSTPLHVVTLSEIAAIVSLI------ASFIYLLGFFGIDFVQSFIS  102 (387)
Q Consensus        60 kir~s~plh~~~~~e~~~~~~~~------as~~yl~~ffgi~~~~~~~~  102 (387)
                      -+.+|-|||.-+++-+++++.++      -.+.|.+|+.|  -+|.++.
T Consensus        63 ~~~e~LPLhlC~ia~~~~~~~l~~~~~~~~~~~~~~gi~G--a~~Ali~  109 (222)
T TIGR02206        63 TLDESLPLHLCDIAIILAAIMLITKRRWFFQLTYFWGIGG--SFQALLT  109 (222)
T ss_pred             chhhcCChhhccHHHHHHHHHHHcCcHHHHHHHHHHHHHH--HHHHHhc
Confidence            57789999999999999988885      46888899877  7888873


No 35 
>PRK14645 hypothetical protein; Provisional
Probab=23.27  E-value=5.8e+02  Score=23.36  Aligned_cols=60  Identities=8%  Similarity=0.039  Sum_probs=42.7

Q ss_pred             CCHHHHHHHHhhh--cCCcceeeEEEEecCCeEEEEEEeecc--ccccchhhhHHHHHHHHHHH
Q 016586          283 NNFETLAVVFNRS--SRFARLQHIQCSIAGKNLYIRFCCTTG--DAMGMNMVSKGVQNVLDFLQ  342 (387)
Q Consensus       283 ~n~~~L~~van~t--SR~GgL~~I~~~i~G~~l~lrf~~dTG--DAMGmNMVnk~~E~v~~~I~  342 (387)
                      +|.+.|.+..+..  ..+-.|.+|+..-.|+.-++|+.+|..  +.++..-.....+++-+.|.
T Consensus         6 ~~~~~i~~li~~~~~~~G~elvdve~~~~~~~~ilrV~ID~~~~~~v~lddC~~vSr~is~~LD   69 (154)
T PRK14645          6 ENNPDLQQLAEGALEPLGYEVLEVQVQRSGGKRIVLVRIDRKDEQPVTVEDLERASRALEAELD   69 (154)
T ss_pred             ccHHHHHHHHHHHHHHcCCEEEEEEEEeCCCCeEEEEEEECCCCCCcCHHHHHHHHHHHHHHhc
Confidence            4667777776654  567789999988788888899999973  34666666666666655554


No 36 
>PF09531 Ndc1_Nup:  Nucleoporin protein Ndc1-Nup;  InterPro: IPR019049  Ndc1 is a nucleoporin protein that is a component of the Nuclear Pore Complex, and, in fungi, also of the Spindle Pole Body. It consists of six transmembrane segments, three luminal loops, both concentrated at the N terminus and cytoplasmic domains largely at the C terminus, all of which are well conserved. 
Probab=21.55  E-value=78  Score=34.10  Aligned_cols=26  Identities=38%  Similarity=0.710  Sum_probs=13.1

Q ss_pred             hhhhhhhHHHHHHHHHHHHHHHHhHhhhh
Q 016586           34 PLYLTNAIFFTLFFSVAYYLLHRWREKIR   62 (387)
Q Consensus        34 p~~~~n~~~~~~f~~~~~~l~~~wr~kir   62 (387)
                      |+|+.   |+.+|.+++|.+.|-++++-|
T Consensus       140 ~lyl~---~~~~~l~~~ys~~~~~~d~~~  165 (602)
T PF09531_consen  140 PLYLW---FSGLFLGLVYSVQHLYFDRDR  165 (602)
T ss_pred             HHHHH---HHHHHHHHHHHHHHHHcCCCE
Confidence            55553   344555555555555554433


No 37 
>TIGR00739 yajC preprotein translocase, YajC subunit. While this protein is part of the preprotein translocase in Escherichia coli, it is not essential for viability or protein secretion. The N-terminus region contains a predicted membrane-spanning region followed by a region consisting almost entirely of residues with charged (acidic, basic, or zwitterionic) side chains. This small protein is about 100 residues in length, and is restricted to bacteria; however, this protein is absent from some lineages, including spirochetes and Mycoplasmas.
Probab=20.51  E-value=75  Score=26.27  Aligned_cols=22  Identities=32%  Similarity=0.549  Sum_probs=16.6

Q ss_pred             hhhhHHHHHHHHHHHHHHHHhH
Q 016586           37 LTNAIFFTLFFSVAYYLLHRWR   58 (387)
Q Consensus        37 ~~n~~~~~~f~~~~~~l~~~wr   58 (387)
                      ++..+++.++|.+.||++.|=.
T Consensus         3 ~~~l~~~vv~~~i~yf~~~rpq   24 (84)
T TIGR00739         3 LTTLLPLVLIFLIFYFLIIRPQ   24 (84)
T ss_pred             HHHHHHHHHHHHHHHHheechH
Confidence            3456788899999999986533


Done!