Query         031239
Match_columns 163
No_of_seqs    44 out of 46
Neff          4.0 
Searched_HMMs 46136
Date          Fri Mar 29 11:16:26 2013
Command       hhsearch -i /work/01045/syshi/csienesis_hhblits_a3m/031239.a3m -d /work/01045/syshi/HHdatabase/Cdd.hhm -o /work/01045/syshi/hhsearch_cdd/031239hhsearch_cdd -cpu 12 -v 0 

 No Hit                             Prob E-value P-value  Score    SS Cols Query HMM  Template HMM
  1 PF01726 LexA_DNA_bind:  LexA D  42.3      39 0.00085   23.2   3.4   48  112-159     6-55  (65)
  2 PRK15345 type III secretion sy  40.0      35 0.00075   31.1   3.5   32  105-147   218-254 (326)
  3 COG0278 Glutaredoxin-related p  36.8      24 0.00052   27.3   1.8   32  128-159    31-63  (105)
  4 PF14003 YlbE:  YlbE-like prote  36.5      45 0.00098   23.7   3.0   43  100-142    11-53  (65)
  5 PF04539 Sigma70_r3:  Sigma-70   31.7      54  0.0012   22.0   2.7   35  115-149     4-39  (78)
  6 PF04031 Las1:  Las1-like ;  In  31.4      30 0.00066   27.5   1.6   47  104-151    75-122 (154)
  7 PF14549 P22_Cro:  DNA-binding   25.1      38 0.00082   23.3   1.0   28  122-150     2-29  (60)
  8 KOG2244 Highly conserved prote  25.1      60  0.0013   32.2   2.7   28   14-42    575-602 (786)
  9 PF12528 DUF3728:  Prepilin pep  24.0      31 0.00067   25.0   0.4   18    4-21     15-32  (84)
 10 PF03826 OAR:  OAR domain;  Int  23.1      55  0.0012   18.6   1.2    9  102-110     8-16  (21)
 11 PRK01964 4-oxalocrotonate taut  22.3      61  0.0013   21.3   1.6   26  128-153    20-45  (64)
 12 PF09339 HTH_IclR:  IclR helix-  22.1      21 0.00045   22.8  -0.7   48  110-160     1-48  (52)
 13 PRK02289 4-oxalocrotonate taut  22.1      66  0.0014   21.1   1.7   22  133-154    25-46  (60)
 14 PF13961 DUF4219:  Domain of un  21.8 1.1E+02  0.0024   17.7   2.4   15   99-113     5-19  (27)
 15 PRK06683 hypothetical protein;  21.5      90  0.0019   22.4   2.4   34  112-146    39-75  (82)
 16 PF00570 HRDC:  HRDC domain Blo  21.2 1.5E+02  0.0033   19.2   3.3   63   78-156     4-66  (68)
 17 TIGR02049 gshA_ferroox glutama  20.5      70  0.0015   30.0   2.0   37  103-153   183-219 (403)

No 1  
>PF01726 LexA_DNA_bind:  LexA DNA binding domain;  InterPro: IPR006199 This is the DNA binding domain of the LexA SOS regulon repressor which prevents expression of DNA repair proteins in bacteria. The aligned region contains a variant form of the helix-turn-helix DNA binding motif []. This domain usually at the N terminus is found associated with IPR006198 from INTERPRO the auto-proteolytic domain of LexA 3.4.21.88 from EC.; GO: 0004252 serine-type endopeptidase activity, 0006508 proteolysis; PDB: 1LEA_A 1JHH_A 3JSP_A 1JHF_A 3JSO_B 1LEB_A 3K2Z_A.
Probab=42.34  E-value=39  Score=23.21  Aligned_cols=48  Identities=27%  Similarity=0.431  Sum_probs=31.0

Q ss_pred             HHHHHhHHHHHHHHH-hhchhhHHHHHHHhCCC-CCcccchhhhhhhhcc
Q 031239          112 DRARRLTEEYSEALL-QIFPQKVQLLGKALGIP-ENSVRTYTEAEIRAGY  159 (163)
Q Consensus       112 DRaRRl~e~~s~~~~-~~f~~~a~~LG~aLGi~-~~av~~FtEa~IRAgv  159 (163)
                      +|=++.++--.+.+. ..++|+...+++.||+. ..+|.-.-++-.+-|.
T Consensus         6 ~rQ~~vL~~I~~~~~~~G~~Pt~rEIa~~~g~~S~~tv~~~L~~Le~kG~   55 (65)
T PF01726_consen    6 ERQKEVLEFIREYIEENGYPPTVREIAEALGLKSTSTVQRHLKALERKGY   55 (65)
T ss_dssp             HHHHHHHHHHHHHHHHHSS---HHHHHHHHTSSSHHHHHHHHHHHHHTTS
T ss_pred             HHHHHHHHHHHHHHHHcCCCCCHHHHHHHhCCCChHHHHHHHHHHHHCcC
Confidence            444455555455444 68999999999999998 7777766666666554


No 2  
>PRK15345 type III secretion system protein SsaL; Provisional
Probab=40.03  E-value=35  Score=31.08  Aligned_cols=32  Identities=34%  Similarity=0.374  Sum_probs=25.2

Q ss_pred             HHHHhHHHHHHHhH-----HHHHHHHHhhchhhHHHHHHHhCCCCCcc
Q 031239          105 LRLKATLDRARRLT-----EEYSEALLQIFPQKVQLLGKALGIPENSV  147 (163)
Q Consensus       105 LRlKAtLDRaRRl~-----e~~s~~~~~~f~~~a~~LG~aLGi~~~av  147 (163)
                      -|+.||+|+-||++     |+||+.           +|+++|+|.+++
T Consensus       218 ~rLaa~i~dLrRLLlFLgLEd~C~r-----------~a~a~~l~gd~i  254 (326)
T PRK15345        218 SRLAAALVRLRRLLLFLGLEKECQR-----------EEWICQLPPNTL  254 (326)
T ss_pred             chHHHHHHHHHHHHHHhchHHHHHH-----------HHHHhCCChhHH
Confidence            79999999999985     777765           567777776654


No 3  
>COG0278 Glutaredoxin-related protein [Posttranslational modification, protein turnover, chaperones]
Probab=36.77  E-value=24  Score=27.33  Aligned_cols=32  Identities=25%  Similarity=0.296  Sum_probs=28.6

Q ss_pred             hchhhHHHHHHHhC-CCCCcccchhhhhhhhcc
Q 031239          128 IFPQKVQLLGKALG-IPENSVRTYTEAEIRAGY  159 (163)
Q Consensus       128 ~f~~~a~~LG~aLG-i~~~av~~FtEa~IRAgv  159 (163)
                      .|+.++-.+=..+| ++-+.|+|-.+++||+++
T Consensus        31 GFS~~~vqiL~~~g~v~~~~vnVL~d~eiR~~l   63 (105)
T COG0278          31 GFSAQAVQILSACGVVDFAYVDVLQDPEIRQGL   63 (105)
T ss_pred             CccHHHHHHHHHcCCcceeEEeeccCHHHHhcc
Confidence            57778888889999 899999999999999986


No 4  
>PF14003 YlbE:  YlbE-like protein
Probab=36.45  E-value=45  Score=23.75  Aligned_cols=43  Identities=16%  Similarity=0.305  Sum_probs=35.5

Q ss_pred             hhHHHHHHHhHHHHHHHhHHHHHHHHHhhchhhHHHHHHHhCC
Q 031239          100 KIIWALRLKATLDRARRLTEEYSEALLQIFPQKVQLLGKALGI  142 (163)
Q Consensus       100 ~~~wALRlKAtLDRaRRl~e~~s~~~~~~f~~~a~~LG~aLGi  142 (163)
                      -|.|..+|---=+.--..-.++-+.+-+.||++++.+.+.++.
T Consensus        11 ~P~WYR~LsR~P~~l~~fe~~a~~~y~kT~p~rVek~~n~lqM   53 (65)
T PF14003_consen   11 QPIWYRILSRNPEELEAFEKEAKHFYKKTIPHRVEKFSNQLQM   53 (65)
T ss_pred             CcHHHHHHccCHHHHHHHHHHHHHHHhccccHHHHHHHhHHHH
Confidence            3578888877777777788888999999999999999887754


No 5  
>PF04539 Sigma70_r3:  Sigma-70 region 3;  InterPro: IPR007624 The bacterial core RNA polymerase complex, which consists of five subunits, is sufficient for transcription elongation and termination but is unable to initiate transcription. Transcription initiation from promoter elements requires a sixth, dissociable subunit called a sigma factor, which reversibly associates with the core RNA polymerase complex to form a holoenzyme []. RNA polymerase recruits alternative sigma factors as a means of switching on specific regulons. Most bacteria express a multiplicity of sigma factors. Two of these factors, sigma-70 (gene rpoD), generally known as the major or primary sigma factor, and sigma-54 (gene rpoN or ntrA) direct the transcription of a wide variety of genes. The other sigma factors, known as alternative sigma factors, are required for the transcription of specific subsets of genes.  With regard to sequence similarity, sigma factors can be grouped into two classes, the sigma-54 and sigma-70 families. Sequence alignments of the sigma70 family members reveal four conserved regions that can be further divided into subregions eg. sub-region 2.2, which may be involved in the binding of the sigma factor to the core RNA polymerase; and sub-region 4.2, which seems to harbor a DNA-binding 'helix-turn-helix' motif involved in binding the conserved -35 region of promoters recognised by the major sigma factors [, ].  Region 3 forms a discrete compact three helical domain within the sigma-factor. Region is not normally involved in the recognition of promoter DNA, but in some specific bacterial promoters containing an extended -10 promoter element, residues within region 3 play an important role. Region 3 primarily is involved in binding the core RNA polymerase in the holoenzyme [].; GO: 0003677 DNA binding, 0003700 sequence-specific DNA binding transcription factor activity, 0016987 sigma factor activity, 0006352 transcription initiation, DNA-dependent, 0006355 regulation of transcription, DNA-dependent; PDB: 1L0O_C 1KU2_B 1RP3_C 1SC5_A 1TTY_A 2BE5_F 2A6E_F 2CW0_F 2A69_P 2A6H_P ....
Probab=31.73  E-value=54  Score=22.01  Aligned_cols=35  Identities=20%  Similarity=0.286  Sum_probs=20.8

Q ss_pred             HHhHHHHHHHHHhh-chhhHHHHHHHhCCCCCcccc
Q 031239          115 RRLTEEYSEALLQI-FPQKVQLLGKALGIPENSVRT  149 (163)
Q Consensus       115 RRl~e~~s~~~~~~-f~~~a~~LG~aLGi~~~av~~  149 (163)
                      +++-..+.++..++ -.|+.+.++..|||+..-|.-
T Consensus         4 ~~i~~a~~~L~~~lgr~Pt~eEiA~~lgis~~~v~~   39 (78)
T PF04539_consen    4 RKIERARRELEQELGREPTDEEIAEELGISVEEVRE   39 (78)
T ss_dssp             HHHHHHHHHHHHHHSS--BHHHHHHHHTS-HHHHHH
T ss_pred             HHHHHHHHHHHHHhCCCCCHHHHHHHHcccHHHHHH
Confidence            33333444444444 689999999999999776653


No 6  
>PF04031 Las1:  Las1-like ;  InterPro: IPR007174 Las1 is an essential nuclear protein involved in cell morphogenesis and cell surface growth [].
Probab=31.36  E-value=30  Score=27.54  Aligned_cols=47  Identities=19%  Similarity=0.312  Sum_probs=33.7

Q ss_pred             HHHHHhHHHHHHHhHHHHHHHHH-hhchhhHHHHHHHhCCCCCcccchh
Q 031239          104 ALRLKATLDRARRLTEEYSEALL-QIFPQKVQLLGKALGIPENSVRTYT  151 (163)
Q Consensus       104 ALRlKAtLDRaRRl~e~~s~~~~-~~f~~~a~~LG~aLGi~~~av~~Ft  151 (163)
                      .+|+.-+.-=. |.+-..+|..+ ..+..+...+++.+|||.|-|++=-
T Consensus        75 ~lr~~ysmaiv-RfVNgl~D~~Q~~~~a~si~~~A~~iglP~~lVdlRH  122 (154)
T PF04031_consen   75 ELRLSYSMAIV-RFVNGLVDPSQQGKYARSIASLAKEIGLPSWLVDLRH  122 (154)
T ss_pred             HHHHHHHHhHH-HHHHHhhhHhhccchhhhHHHHHHHcCCCHHHHHHHH
Confidence            45655554443 45666667666 6777888999999999999987633


No 7  
>PF14549 P22_Cro:  DNA-binding transcriptional regulator Cro; PDB: 1RZS_A 3BD1_A 3QWS_A 2HIN_B.
Probab=25.11  E-value=38  Score=23.25  Aligned_cols=28  Identities=25%  Similarity=0.444  Sum_probs=19.9

Q ss_pred             HHHHHhhchhhHHHHHHHhCCCCCcccch
Q 031239          122 SEALLQIFPQKVQLLGKALGIPENSVRTY  150 (163)
Q Consensus       122 s~~~~~~f~~~a~~LG~aLGi~~~av~~F  150 (163)
                      .+.+.+.|+. ...|+++|||.+.+|.-.
T Consensus         2 k~~aI~~~G~-~~~lAkalGVs~~aVs~W   29 (60)
T PF14549_consen    2 KKDAIKYFGG-QSKLAKALGVSPQAVSQW   29 (60)
T ss_dssp             HHHHHHHHSS-HHHHHHHHTS-HHHHHHH
T ss_pred             HHHHHHHHCC-HHHHHHHHCCCHHHHHHh
Confidence            3456677874 558999999998887544


No 8  
>KOG2244 consensus Highly conserved protein containing a thioredoxin domain [General function prediction only]
Probab=25.07  E-value=60  Score=32.23  Aligned_cols=28  Identities=39%  Similarity=0.735  Sum_probs=23.3

Q ss_pred             cchhHHHHHHHHHHHHHhhccCcchhhhh
Q 031239           14 GLEDYLFVLLSRFLNALETKGGAHWLAEN   42 (163)
Q Consensus        14 ~lEdyaFVLLSR~ln~le~~gg~~~~~~~   42 (163)
                      =||||||+ .+-+|..-++.|+.+||.-.
T Consensus       575 FldDYAFl-I~gLLDlYea~~~~e~LkwA  602 (786)
T KOG2244|consen  575 FLDDYAFL-ISGLLDLYEAGGGIEWLKWA  602 (786)
T ss_pred             hhhhHHHH-HHHHHHHHHccCchHHHHHH
Confidence            58999996 57788999999999999433


No 9  
>PF12528 DUF3728:  Prepilin peptidase dependent protein C (DUF3728);  InterPro: IPR022204  This family of proteins is found in bacteria. Proteins in this family are typically between 106 and 121 amino acids in length. The family is found in association with PF07963 from PFAM. There are two completely conserved C residues that may be functionally important. This family is frequently annotated as prepilin peptidase dependent protein C however there is little accompanying literature to confirm this. 
Probab=23.96  E-value=31  Score=25.03  Aligned_cols=18  Identities=33%  Similarity=0.447  Sum_probs=15.5

Q ss_pred             chhhHHHhhhcchhHHHH
Q 031239            4 LFNQWRLCEIGLEDYLFV   21 (163)
Q Consensus         4 ~rQ~WRlaEi~lEdyaFV   21 (163)
                      +||.||+|.=-||-|..-
T Consensus        15 qrQAWr~a~q~le~~~~~   32 (84)
T PF12528_consen   15 QRQAWRLAHQQLESYEAG   32 (84)
T ss_pred             HHHHHHHHHHHHhhcccc
Confidence            799999999998888765


No 10 
>PF03826 OAR:  OAR domain;  InterPro: IPR003654 This 14 amino acid motif has been identified within the C-terminal region of several Paired-like homeodomain (HD) containing proteins [, ]. It was named OAR domain after the initials of otp, aristaless, and rax []. Although it has been proposed that this domain could be important for transactivation and be involved in protein-protein interactions or DNA binding [, ], is function is not yet known. Some proteins known to contain a OAR domain include human RIEG, defects in which are the cause of Rieger syndrome []; human OG12X and Mus musculus (Mouse) Og12x, whose function is not yet known []; vertebrate Rax, which plays a role in the proliferation and/or differentiation of retinal cells []; Drosophila DRX, which appears to be important in brain development []; and human SHOX, encoded by the short stature homeobox-containing gene. Defects or lack of this protein are the cause of short stature associated with the Turner syndrome [].; GO: 0003677 DNA binding, 0007275 multicellular organismal development, 0005634 nucleus
Probab=23.05  E-value=55  Score=18.62  Aligned_cols=9  Identities=78%  Similarity=0.785  Sum_probs=7.2

Q ss_pred             HHHHHHHhH
Q 031239          102 IWALRLKAT  110 (163)
Q Consensus       102 ~wALRlKAt  110 (163)
                      +-+||+||-
T Consensus         8 Ia~LRlKAk   16 (21)
T PF03826_consen    8 IAALRLKAK   16 (21)
T ss_pred             HHHHHHHHH
Confidence            678999973


No 11 
>PRK01964 4-oxalocrotonate tautomerase; Provisional
Probab=22.35  E-value=61  Score=21.28  Aligned_cols=26  Identities=15%  Similarity=0.409  Sum_probs=18.0

Q ss_pred             hchhhHHHHHHHhCCCCCcccchhhh
Q 031239          128 IFPQKVQLLGKALGIPENSVRTYTEA  153 (163)
Q Consensus       128 ~f~~~a~~LG~aLGi~~~av~~FtEa  153 (163)
                      ++..-.+.+.+.||+|++.+.++-++
T Consensus        20 l~~~it~~l~~~lg~p~~~v~V~i~e   45 (64)
T PRK01964         20 LIREVTEAISATLDVPKERVRVIVNE   45 (64)
T ss_pred             HHHHHHHHHHHHhCcChhhEEEEEEE
Confidence            33334456777899999998887544


No 12 
>PF09339 HTH_IclR:  IclR helix-turn-helix domain;  InterPro: IPR005471 The many bacterial transcription regulation proteins which bind DNA through a 'helix-turn-helix' motif can be classified into subfamilies on the basis of sequence similarities. One of these subfamilies, called 'iclR', groups several proteins including:  gylR, a possible activator protein for the gylABX glycerol operon in Streptomyces.   iclR, the repressor of the acetate operon (also known as glyoxylate bypass operon) in Escherichia coli and Salmonella typhimurium.    These proteins have a Helix-Turn-Helix motif at the N terminus that is similar to that of other DNA-binding proteins [].; GO: 0003677 DNA binding, 0006355 regulation of transcription, DNA-dependent; PDB: 1MKM_A 3MQ0_A 3R4K_A 2G7U_C 2O0Y_C 2XRO_F 2XRN_B 2IA2_D.
Probab=22.15  E-value=21  Score=22.77  Aligned_cols=48  Identities=31%  Similarity=0.567  Sum_probs=31.5

Q ss_pred             HHHHHHHhHHHHHHHHHhhchhhHHHHHHHhCCCCCcccchhhhhhhhccc
Q 031239          110 TLDRARRLTEEYSEALLQIFPQKVQLLGKALGIPENSVRTYTEAEIRAGYV  160 (163)
Q Consensus       110 tLDRaRRl~e~~s~~~~~~f~~~a~~LG~aLGi~~~av~~FtEa~IRAgvv  160 (163)
                      +++|+-+.++-.++.=.   +-+...|.+.+|+|..++.-+-..-.+.|.+
T Consensus         1 sl~ral~iL~~l~~~~~---~~t~~eia~~~gl~~stv~r~L~tL~~~g~v   48 (52)
T PF09339_consen    1 SLERALRILEALAESGG---PLTLSEIARALGLPKSTVHRLLQTLVEEGYV   48 (52)
T ss_dssp             HHHHHHHHHHCHHCTBS---CEEHHHHHHHHTS-HHHHHHHHHHHHHTTSE
T ss_pred             CHHHHHHHHHHHHcCCC---CCCHHHHHHHHCcCHHHHHHHHHHHHHCcCe
Confidence            46677666655444311   1268899999999998887776666666654


No 13 
>PRK02289 4-oxalocrotonate tautomerase; Provisional
Probab=22.07  E-value=66  Score=21.14  Aligned_cols=22  Identities=5%  Similarity=0.336  Sum_probs=16.8

Q ss_pred             HHHHHHHhCCCCCcccchhhhh
Q 031239          133 VQLLGKALGIPENSVRTYTEAE  154 (163)
Q Consensus       133 a~~LG~aLGi~~~av~~FtEa~  154 (163)
                      .+.+.+.||+|+.+|.++-++.
T Consensus        25 t~a~~~~~~~p~~~v~V~i~ev   46 (60)
T PRK02289         25 TEVVSRIAKAPKEAIHVFINDM   46 (60)
T ss_pred             HHHHHHHhCcCcceEEEEEEEe
Confidence            3456678899999999887664


No 14 
>PF13961 DUF4219:  Domain of unknown function (DUF4219)
Probab=21.82  E-value=1.1e+02  Score=17.66  Aligned_cols=15  Identities=33%  Similarity=0.671  Sum_probs=12.7

Q ss_pred             chhHHHHHHHhHHHH
Q 031239           99 GKIIWALRLKATLDR  113 (163)
Q Consensus        99 ~~~~wALRlKAtLDR  113 (163)
                      ..+.|+.|+|+.|+.
T Consensus         5 NY~~W~~~M~~~L~~   19 (27)
T PF13961_consen    5 NYSTWKIRMKAYLES   19 (27)
T ss_pred             CHHHHHHHHHHHHHH
Confidence            467999999999875


No 15 
>PRK06683 hypothetical protein; Provisional
Probab=21.52  E-value=90  Score=22.38  Aligned_cols=34  Identities=21%  Similarity=0.141  Sum_probs=22.1

Q ss_pred             HHHHHhHHHHHHHHH---hhchhhHHHHHHHhCCCCCc
Q 031239          112 DRARRLTEEYSEALL---QIFPQKVQLLGKALGIPENS  146 (163)
Q Consensus       112 DRaRRl~e~~s~~~~---~~f~~~a~~LG~aLGi~~~a  146 (163)
                      ++++.-+.++++.+.   -.|+ +-..||+++|||..+
T Consensus        39 ~~~~~~i~~~~~~~~Vpv~~~~-t~~eLG~A~G~~~~~   75 (82)
T PRK06683         39 MRLTHVIIRTALQHNIPITKVE-SVRKLGKVAGIQVGA   75 (82)
T ss_pred             HHHHHHHHHHHHhcCCCEEEEC-CHHHHHHHhCCcccE
Confidence            445555555555432   2333 779999999999764


No 16 
>PF00570 HRDC:  HRDC domain Bloom syndrome. Werner syndrome.;  InterPro: IPR002121 The HRDC (Helicase and RNase D C-terminal) domain has a putative role in nucleic acid binding. Mutations in the HRDC domain associated with the human BLM gene result in Bloom Syndrome (BS), an autosomal recessive disorder characterised by proportionate pre- and postnatal growth deficiency; sun-sensitive, telangiectatic, hypo- and hyperpigmented skin; predisposition to malignancy; and chromosomal instability [].; GO: 0003676 nucleic acid binding, 0005622 intracellular; PDB: 3SAG_B 3SAH_B 2CPR_A 3SAF_B 3CYM_A 1WUD_A 2HBK_A 2HBJ_A 2HBM_A 2HBL_A ....
Probab=21.20  E-value=1.5e+02  Score=19.21  Aligned_cols=63  Identities=24%  Similarity=0.367  Sum_probs=41.7

Q ss_pred             HHHHHHHHHHhccccccCCCcchhHHHHHHHhHHHHHHHhHHHHHHHHHhhchhhHHHHHHHhCCCCCcccchhhhhhh
Q 031239           78 IGNELFAWQEKGLSEKEGSEDGKIIWALRLKATLDRARRLTEEYSEALLQIFPQKVQLLGKALGIPENSVRTYTEAEIR  156 (163)
Q Consensus        78 i~~eL~~W~~~~~~~~~~~~d~~~~wALRlKAtLDRaRRl~e~~s~~~~~~f~~~a~~LG~aLGi~~~av~~FtEa~IR  156 (163)
                      +..+|..|.......    .|            ..+..-+..+--..+....|.+.+.|-+.-|+.++-++-|-++.|.
T Consensus         4 ~~~~L~~~R~~~A~~----~~------------~~~~~Il~~~~L~~ia~~~P~s~~~L~~i~g~~~~~~~~~g~~il~   66 (68)
T PF00570_consen    4 LLKALKEWREELARE----ED------------VPPYRILSDEALLEIAKRLPTSIEELLQIPGMGKRKVRKYGDEILE   66 (68)
T ss_dssp             HHHHHHHHHHHHHHH----HT------------S-HHHHS-HHHHHHHHHH--SSHHHHHTSTTCGHHHHHHCHHHHHH
T ss_pred             HHHHHHHHHHHHHHH----cC------------cCcccccCHHHHHHHHHhCCCCHHHHHHccCCCHHHHHHHHHHHHh
Confidence            456889998665211    11            1223334444555677778999999999999999999999998875


No 17 
>TIGR02049 gshA_ferroox glutamate--cysteine ligase, T. ferrooxidans family. This family consists of a rare family of glutamate--cysteine ligases, demonstrated first in Thiobacillus ferrooxidans and present in a few other Proteobacteria. It is the first of two enzymes for glutathione biosynthesis. It is also called gamma-glutamylcysteine synthetase.
Probab=20.47  E-value=70  Score=29.95  Aligned_cols=37  Identities=22%  Similarity=0.409  Sum_probs=30.7

Q ss_pred             HHHHHHhHHHHHHHhHHHHHHHHHhhchhhHHHHHHHhCCCCCcccchhhh
Q 031239          103 WALRLKATLDRARRLTEEYSEALLQIFPQKVQLLGKALGIPENSVRTYTEA  153 (163)
Q Consensus       103 wALRlKAtLDRaRRl~e~~s~~~~~~f~~~a~~LG~aLGi~~~av~~FtEa  153 (163)
                      |+-|-|+.              .-+.|+.-+..+++.+|||+|-++.|.+.
T Consensus       183 W~~RrKS~--------------hF~~Y~~va~eF~~~~~IDPWlInp~f~~  219 (403)
T TIGR02049       183 WHSRRKSN--------------HFEAYQEVASEFAKLIGIDPWLINPYFEK  219 (403)
T ss_pred             cccccccc--------------hhHHHHHHHHHHHHHhCCCcccccHhhhc
Confidence            88888864              44567778999999999999999988763


Done!