hmmpfam - search one or more sequences against HMM database
HMMER 2.3.2 (Oct 2003)
Copyright (C) 1992-2003 HHMI/Washington University School of Medicine
Freely distributed under the GNU General Public License (GPL)
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
HMM file:                 /usr1/pfam-23.0/Pfam_fs
Sequence file:            33597077.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query sequence: gi|33597077|ref|NP_884720.1|
Accession:      [none]
Description:    adhesin [Bordetella parapertussis 12822]

Scores for sequence family classification (score includes all domains):
Model           Description                             Score    E-value  N
--------        -----------                             -----    ------- ---
Fil_haemagg     Haemagluttinin repeat                   162.2    4.9e-46  16
Haemagg_act     haemagglutination activity domain       116.4    2.4e-33   1
Strep_SA_rep    Streptococcal surface antigen repeat      9.2       0.13   2
DUF2190         Uncharacterized conserved protein (DU     7.0       0.17   1
RHS_repeat      RHS Repeat                                8.2       0.18   2
LEA_1           Late embryogenesis abundant (LEA) gro     8.5       0.18   1
P_gingi_FimA    Porphyromonas gingivalis major fimbri     3.1        0.2   1
HutD            HutD                                      6.1       0.22   2
HupH_C          HupH hydrogenase expression protein,      5.1       0.24   1
DUF2604         Protein of Unknown function (DUF2604)     5.2       0.39   1
WWbp            WW-domain ligand protein                  4.1       0.42   1
DUF2600         Protein of unknown function (DUF2600)     4.1       0.45   1
DUF490          Family of unknown function (DUF490)       4.1       0.61   1
Pentapeptide    Pentapeptide repeats (8 copies)           6.4       0.87   1
Phage_Capsid_P3 P3 major capsid protein                   2.0       0.98   1
CDC48_N         Cell division protein 48 (CDC48), N-t     4.4        1.2   1
HMG_box         HMG (high mobility group) box             5.0        1.4   1
AHS2            Allophanate hydrolase subunit 2           2.4        1.4   1
Enolase_N       Enolase, N-terminal domain                3.2        1.5   1
GlcNAc_2-epim   N-acylglucosamine 2-epimerase (GlcNAc     1.5        1.6   1
DUF541          Protein of unknown function (DUF541)      3.6        1.6   1
FlgD            Flagellar hook capping protein            3.1        1.6   1
Androgen_recep  Androgen receptor                         0.9        1.6   1
DPBB_1          Rare lipoprotein A (RlpA)-like double     4.1        1.8   1
COLFI           Fibrillar collagen C-terminal domain      0.1        2.1   1
MmgE_PrpD       MmgE/PrpD family                          2.1        2.3   1
AIRS            AIR synthase related protein, N-termi     2.8        2.3   2
Acylphosphatase Acylphosphatase                           4.3        2.3   1
CbbQ_C          CbbQ/NirQ/NorQ C-terminal                 3.6        2.3   2
Ery_res_leader1 Erythromycin resistance leader peptid     5.8        2.6   1
Nrf1_activ_bdg  Nrf1 activator activation site bindin     3.3        2.9   1
DUF1643         Protein of unknown function (DUF1643)     0.2          3   1
UPF0164         Uncharacterised protein family (UPF01     1.0        3.3   1
IATP            Mitochondrial ATPase inhibitor, IATP      3.7        3.5   1
SPO11_like      SPO11 homologue                           2.7        3.9   1
DUF1840         Domain of unknown function (DUF1840)      2.2        3.9   1
PAD_M           Protein-arginine deiminase (PAD) midd     1.1          4   1
Type_III_YscP   Type III secretion system  (type_III_     1.9        4.1   1
AD              Anticodon-binding domain                  3.1        4.2   1
HNF-1A_C        Hepatocyte nuclear factor 1 (HNF-1),      1.8        4.4   1
PseudoU_synth_1 tRNA pseudouridine synthase               1.7        4.4   1
DUF917          Protein of unknown function (DUF917)      0.9        4.7   1
Orexin          Prepro-orexin                             0.8        4.8   1
LolB            Outer membrane lipoprotein LolB           1.7        5.1   1
DUF2398         Protein of unknown function (DUF2398)    -0.2        5.2   1
DUF1885         Domain of unknown function (DUF1885)      1.3        5.2   1
PACT_coil_coil  Pericentrin-AKAP-450 domain of centro     2.6        5.2   1
BtpA            BtpA family                               0.8        5.3   1
Entericidin     Entericidin EcnA/B family                 3.6        5.4   1
Nairovirus_M    Nairovirus M polyprotein-like            -1.9        5.5   1
TraS            Plasmid conjugative transfer entry ex    -0.0        5.5   1
ATP-synt_ab     ATP synthase alpha/beta family, nucle     0.9        5.5   1
RNR_N           Ribonucleotide reductase N-terminal       2.1        5.7   1
DUF2322         Uncharacterized protein conserved in      1.9        5.8   1
Treacle         Treacher Collins syndrome protein Tre    -1.0        5.9   1
DCD             2'-deoxycytidine 5'-triphosphate deam    -0.8        5.9   1
DUF1098         Protein of unknown function (DUF1098)     2.8        6.1   1
DUF1109         Protein of unknown function (DUF1109)     0.1        6.3   1
DUF166          Uncharacterized ArCR, COG1810            -0.2        6.3   1
FbpA            Fibronectin-binding protein A N-termi    -0.3        6.5   1
Antimicrobial10 Metchnikowin family                       2.1        6.5   1
HSP70           Hsp70 protein                            -0.4        7.4   1
DNA_pol3_theta  DNA polymerase III, theta subunit         1.9        7.4   1
Sigma70_r2      Sigma-70 region 2                         2.7        7.4   1
DUF992          Protein of unknown function (DUF992)      1.7        7.7   1
DGOK            2-keto-3-deoxy-galactonokinase           -1.6        7.9   1
ATP-grasp_2     ATP-grasp domain                         -0.4        8.3   1
YopE            Yersinia virulence determinant (YopE)     2.0        8.4   1
DUF1205         Protein of unknown function (DUF1205)     0.1        8.4   1
Glyco_hydro_70  Glycosyl hydrolase family 70             -1.4        8.5   1
CagZ            CagZ                                     -1.1        8.6   1
Autoind_synth   Autoinducer synthetase                    0.3        8.8   2
Flavi_M         Flavivirus envelope glycoprotein M        1.2        8.8   1
Phage-MuB_C     Mu B transposition protein, C termina     2.0          9   1
Carot_N         Orange carotenoid protein, N-terminal     0.6        9.1   1
Adeno_PIX       Adenovirus hexon-associated protein (     1.0        9.2   1
HTH_DeoR        DeoR-like helix-turn-helix domain         2.1        9.3   1
Birna_VP3       Birnavirus VP3 protein                   -0.2        9.4   1
DUF500          Family of unknown function (DUF500)       1.2        9.4   1
FERM_C          FERM C-terminal PH-like domain            0.9        9.5   1
DUF547          Protein of unknown function, DUF547       1.0        9.6   1
Transposase_20  Transposase IS116/IS110/IS902 family      1.6        9.9   1
Head_binding    Head binding                             -0.0        9.9   1

Parsed for domains:
Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------        ------- ----- -----    ----- -----      -----  -------
DUF992            1/1      58    75 ..     1    18 [.     1.7      7.7
Haemagg_act       1/1      69   202 ..     1   134 [.   116.4  2.4e-33
Orexin            1/1      70    84 ..   117   131 .]     0.8      4.8
COLFI             1/1      76    89 ..     1    15 [.     0.1      2.1
Fil_haemagg       1/16    256   337 ..     1    75 []    28.6  1.6e-07
DUF1098           1/1     266   290 ..     1    25 [.     2.8      6.1
Entericidin       1/1     284   293 ..    34    43 .]     3.6      5.4
DUF2190           1/1     345   367 ..    73    97 ..     7.0     0.17
Fil_haemagg       2/16    346   405 ..     1    75 []     3.8      2.2
Fil_haemagg       3/16    415   475 ..     1    75 []     5.2     0.92
ATP-synt_ab       1/1     461   470 ..   202   211 .]     0.9      5.5
P_gingi_FimA      1/1     501   520 ..     1    21 [.     3.1      0.2
Fil_haemagg       4/16    524   547 ..    25    48 ..     1.7      9.1
Fil_haemagg       5/16    548   629 ..     1    75 []     6.7     0.34
Phage-MuB_C       1/1     591   609 ..    46    63 ..     2.0        9
Fil_haemagg       6/16    633   692 ..     1    75 []    13.6   0.0034
Phage_Capsid_P3   1/1     649   677 ..   366   394 .]     2.0     0.98
Fil_haemagg       7/16    695   766 ..     1    75 []    20.9  2.7e-05
RHS_repeat        1/2     709   723 ..    25    44 .]     1.3       17
DCD               1/1     730   762 ..     1    33 [.    -0.8      5.9
Fil_haemagg       8/16    777   833 ..     1    75 []    13.3   0.0043
YopE              1/1     797   806 ..    61    70 .]     2.0      8.4
DUF1840           1/1     816   831 ..     8    23 ..     2.2      3.9
Fil_haemagg       9/16    838   892 ..     1    75 []     8.1     0.13
Fil_haemagg      10/16    900   950 ..     1    75 []    19.3    8e-05
Fil_haemagg      11/16    953   980 ..     1    33 [.     0.9       16
DUF500            1/1     961   977 ..    26    42 ..     1.2      9.4
DUF490            1/1     962   992 ..     1    39 [.     4.1     0.61
PAD_M             1/1     984   991 ..     1     8 [.     1.1        4
Nrf1_activ_bdg    1/1     995  1010 ..   103   118 .]     3.3      2.9
Fil_haemagg      12/16    996  1053 ..     1    75 []    22.4    1e-05
HutD              1/2    1033  1061 ..   168   200 .]     4.9     0.46
Ery_res_leader1   1/1    1056  1069 ..     1    17 []     5.8      2.6
Fil_haemagg      13/16   1058  1105 ..     1    71 [.     4.4      1.6
CDC48_N           1/1    1082  1100 ..     1    19 [.     4.4      1.2
Androgen_recep    1/1    1126  1159 ..   427   460 .]     0.9      1.6
Fil_haemagg      14/16   1163  1203 ..    34    75 .]     6.2     0.46
RNR_N             1/1    1244  1267 ..    32    55 ..     2.1      5.7
Nairovirus_M      1/1    1308  1318 ..     1    11 [.    -1.9      5.5
HMG_box           1/1    1322  1338 ..    53    69 .]     5.0      1.4
Strep_SA_rep      1/2    1322  1345 ..     1    25 []     5.7      1.3
HSP70             1/1    1323  1392 ..   545   619 .]    -0.4      7.4
PseudoU_synth_1   1/1    1326  1366 ..    72   125 .]     1.7      4.4
MmgE_PrpD         1/1    1328  1357 ..     1    33 [.     2.1      2.3
Birna_VP3         1/1    1350  1382 ..    87   119 ..    -0.2      9.4
Transposase_20    1/1    1351  1365 ..    73    88 .]     1.6      9.9
AD                1/1    1360  1393 ..     9    42 ..     3.1      4.2
Adeno_PIX         1/1    1360  1386 ..    89   115 .]     1.0      9.2
Strep_SA_rep      2/2    1363  1387 ..     1    25 []     3.6      5.3
CbbQ_C            1/2    1377  1391 ..    72    86 .]     0.2       23
Enolase_N         1/1    1381  1402 ..     1    22 [.     3.2      1.5
Head_binding      1/1    1609  1625 ..    96   114 .]    -0.0      9.9
DUF166            1/1    1637  1649 ..   252   264 .]    -0.2      6.3
DUF547            1/1    1687  1697 ..   143   153 .]     1.0      9.6
IATP              1/1    1720  1734 ..    14    28 ..     3.7      3.5
FlgD              1/1    1735  1775 ..    97   155 .]     3.1      1.6
DGOK              1/1    1794  1807 ..   285   298 .]    -1.6      7.9
CbbQ_C            2/2    1795  1810 ..    71    86 .]     3.5      2.6
RHS_repeat        2/2    1823  1836 ..     1    14 [.     6.9     0.41
Pentapeptide      1/1    1912  1926 ..    26    40 .]     6.4     0.87
CagZ              1/1    1963  1976 ..   186   199 .]    -1.1      8.6
DUF2604           1/1    1971  1981 ..    66    76 .]     5.2     0.39
Fil_haemagg      15/16   1975  2011 ..     1    47 [.     1.2       13
DNA_pol3_theta    1/1    2016  2029 ..    49    62 ..     1.9      7.4
Glyco_hydro_70    1/1    2055  2066 ..   457   468 ..    -1.4      8.5
HutD              2/2    2124  2158 ..   164   200 .]     1.1      5.7
DUF541            1/1    2137  2159 ..     1    23 [.     3.6      1.6
Fil_haemagg      16/16   2178  2235 ..     1    75 []     6.0     0.52
Antimicrobial10   1/1    2278  2297 ..     1    20 [.     2.1      6.5
LEA_1             1/1    2291  2315 ..    50    73 .]     8.5     0.18
WWbp              1/1    2321  2332 ..     1    12 [.     4.1     0.42
ATP-grasp_2       1/1    2390  2398 ..   131   139 ..    -0.4      8.3
TraS              1/1    2394  2413 ..   145   164 .]    -0.0      5.5
GlcNAc_2-epim     1/1    2468  2494 ..     1    25 [.     1.5      1.6
Sigma70_r2        1/1    2571  2594 ..    48    71 .]     2.7      7.4
DUF1643           1/1    2616  2621 ..   147   152 .]     0.2        3
FbpA              1/1    2648  2666 ..     1    19 [.    -0.3      6.5
DUF1109           1/1    2695  2706 ..     1    13 [.     0.1      6.3
Flavi_M           1/1    2778  2784 ..    69    75 .]     1.2      8.8
BtpA              1/1    2826  2844 ..   258   276 .]     0.8      5.3
PACT_coil_coil    1/1    2830  2843 ..    75    94 .]     2.6      5.2
DUF2322           1/1    2911  2926 ..    84   101 .]     1.9      5.8
AIRS              1/2    3049  3058 ..   115   124 .]     0.0       16
SPO11_like        1/1    3213  3227 ..    29    43 .]     2.7      3.9
FERM_C            1/1    3322  3332 ..    86    96 .]     0.9      9.5
DUF2600           1/1    3385  3399 ..     1    15 [.     4.1     0.45
UPF0164           1/1    3427  3456 ..   313   345 .]     1.0      3.3
DPBB_1            1/1    3445  3462 ..     1    18 [.     4.1      1.8
Treacle           1/1    3458  3473 ..   495   510 .]    -1.0      5.9
Acylphosphatase   1/1    3494  3521 ..    50    76 ..     4.3      2.3
DUF2398           1/1    3535  3548 ..     1    14 [.    -0.2      5.2
Autoind_synth     1/2    3556  3571 ..   172   187 .]     0.3      9.3
DUF1885           1/1    3641  3646 ..   138   143 .]     1.3      5.2
AHS2              1/1    3745  3764 ..   273   292 .]     2.4      1.4
HupH_C            1/1    3752  3764 ..   111   123 .]     5.1     0.24
HTH_DeoR          1/1    3755  3767 ..     1    13 [.     2.1      9.3
LolB              1/1    3801  3807 ..   163   169 .]     1.7      5.1
AIRS              2/2    3897  3907 ..   114   124 .]     2.7      2.3
HNF-1A_C          1/1    3912  3919 ..    83    90 .]     1.8      4.4
Autoind_synth     2/2    3966  3981 ..   172   187 .]     0.1       10
Type_III_YscP     1/1    3983  4002 ..     1    26 [.     1.9      4.1
DUF1205           1/1    4048  4055 ..   104   111 .]     0.1      8.4
DUF917            1/1    4092  4102 ..     1    11 [.     0.9      4.7
Carot_N           1/1    4202  4212 ..   155   165 .]     0.6      9.1

Alignments of top-scoring domains:
DUF992: domain 1 of 1, from 58 to 75: score 1.7, E = 7.7
                   *->aaaaalaapaaaaaaang<-*
                      a+++a++apa aaaaa+g
  gi|3359707    58    AVTLAVLAPAGAAAAAGG    75

Haemagg_act: domain 1 of 1, from 69 to 202: score 116.4, E = 2.4e-33
                CS    X---EE-TT-EEE-... -TTS--EEE-----TTSEEEEEEEE--B-
                   *->aqslpdgglvtngtvti.qatgvpitggtipqtsgnlfisfssFnVg
                      a+++  g ++++ t   ++ +gvp++ +++p  +g++++++++FnV
  gi|3359707    69    AAAAAGGIAPQGATQVApARNGVPVIRIAAPDATGISHNRYTEFNVR 115

                CS TTEEEEEE-SS EEETTTEEE---TT-.SS--SEEEEEE-SS--EEEEEE
                   qggtvlFnatq.gqSqLgGnltgNpnLapgaasnIlNRVTGgnaSqIqGl
                   q g+vl+n+t +g S L+G + gNp L  g a++IlN+VTG ++  + G+
  gi|3359707   116 QPGVVLNNSTAeGVSALAGRISGNPGL-RGPARAILNEVTGVSPTTLEGA 164

                CS EEEEE.S-EEEEEE-TT-EEEEEEEEESEEEEEEE-SEE
                   levlGganAnVfLaNPNGIifgggAfiNvgglfvtTgai<-*
                   lev+G   A+V  aNPNG + +g ++iN++gl+++Tg +
  gi|3359707   165 LEVFG-PAADVLVANPNGLTANGLSTINIRGLTLSTGRP    202

Orexin: domain 1 of 1, from 70 to 84: score 0.8, E = 4.8
                   *->aAaAsalAPrGrsrv<-*
                      aAaA  +AP+G  +v
  gi|3359707    70    AAAAGGIAPQGATQV    84

COLFI: domain 1 of 1, from 76 to 89: score 0.1, E = 2.1
                   *->lksPeGksrknPARt<-*
                       ++P+G ++ +PAR+
  gi|3359707    76    -IAPQGATQVAPARN    89

Fil_haemagg: domain 1 of 16, from 256 to 337: score 28.6, E = 1.6e-07
                CS    EEESSEEEETT                          TT--EE-SEE
                   *->vaaGaltlnaa..........................alggldlnng
                      v+aGa++++ +++++ ++++++   ++ + ++  +++a+g+++ +++
  gi|3359707   256    VVAGANRYDYGrrtrlqptgrpaaapdasaghaidgsAAGAMYGRHI 302

                CS EEEE-STT-EEEE-S-EEESSE...................EEEEESS-E
                   gtlsagggltltaagllnnggtliaggdlllaagalrnggdltltaaGdL
                    ++s+++gl+++ +gl  ++g                    ++++aaG +
  gi|3359707   303 TLVSSDAGLGVRHRGLVSAAGH-------------------IAIDAAGGV 333

                CS EEEE
                   dnsg<-*
                   ++++
  gi|3359707   334 SVQA    337

DUF1098: domain 1 of 1, from 266 to 290: score 2.8, E = 6.1
                   *->rrnrrrrrptrdvsfdpvdnvnein<-*
                      +r+ r +++ r+ +  ++ ++++i+
  gi|3359707   266    GRRTRLQPTGRPAAAPDASAGHAID    290

Entericidin: domain 1 of 1, from 284 to 293: score 3.6, E = 5.4
                   *->saGeAiedaA<-*
                      saG+Ai+++A
  gi|3359707   284    SAGHAIDGSA    293

DUF2190: domain 1 of 1, from 345 to 367: score 7.0, E = 0.17
                   *->DdTngekvvtttAtgNtliGvAwaa<-*
                      D T+  + +t+ A+g++liG+++a+
  gi|3359707   345    DGTA--GSATIRAGGDALIGSVIAF    367

Fil_haemagg: domain 2 of 16, from 346 to 405: score 3.8, E = 2.2
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EE ESSE...
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaaglln.nggtlia
                       +aG+ t++a+    +d+  g  ++ gg++ ++aa+     +g+l a
  gi|3359707   346    GTAGSATIRAG----GDALIGSVIAFGGDIEVRAAAGRIaTHGALNA 388

                CS ................EEEEESS-EEEEE
                   ggdlllaagalrnggdltltaaGdLdnsg<-*
                     +           + ++l a+++L ++g
  gi|3359707   389 PRG-----------R-VRLQAGDALAVHG    405

Fil_haemagg: domain 3 of 16, from 415 to 475: score 5.2, E = 0.92
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      ++ G+ +++a+           + +++ g++l a+++  +++  ia+
  gi|3359707   415    ATRGDASIAAD-----------LHATRDGVSLAARNATLAHARVIAT 450

                CS ..............    .EEEEESS-EEEEE
                   gdlllaagalrngg....dltltaaGdLdnsg<-*
                            +++++++++++ ++ +G+L+++g
  gi|3359707   451 PP-------AAQARpstpAIDIALSGTLTLRG    475

ATP-synt_ab: domain 1 of 1, from 461 to 470: score 0.9, E = 5.5
                CS    S-B-TTT-B-
                RF    xxxxxxxxxx
                   *->PAIDvllSvS<-*
                      PAID+ lS++
  gi|3359707   461    PAIDIALSGT    470

P_gingi_FimA: domain 1 of 1, from 501 to 520: score 3.1, E = 0.2
                   *->vlktSnpnRafePnGlgdddA<-*
                      ++ t  p+Ra+ Pn l+ +dA
  gi|3359707   501    IVDTRQPQRAV-PNALLASDA    520

Fil_haemagg: domain 4 of 16, from 524 to 547: score 1.7, E = 9.1
                CS    E-STT-EEEE-S-EEESSE.....
                   *->sagggltltaagllnnggtliagg<-*
                      +ag+ ++++aagl n gg+++++g
  gi|3359707   524    AAGQPISIRAAGLRNEGGAIDSTG    547

Fil_haemagg: domain 5 of 16, from 548 to 629: score 6.7, E = 0.34
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      +  G   l+++    + + n g +++ g+l +t  g+l n   l+++
  gi|3359707   548    TGQGHIGLRIG----GPAVNLGYIATHGTLEATVDGTLEN-RMLLSA 589

                CS ..............            .EEEEESS-EEEEE
                   gdlllaagalrngg............dltltaaGdLdnsg<-*
                   ++l++a+  l+n    + +++++++++l l+ + ++dn+g
  gi|3359707   590 AALRVATHDLSNHAamtasgpdtgtpALDLSVQHRVDNAG    629

Phage-MuB_C: domain 1 of 1, from 591 to 609: score 2.0, E = 9
                CS    SHHHHHHHH GGGT-TT-T
                   *->aLRiLtqtL.RlAaMtAkG<-*
                      aLR+ t+ L++ AaMtA+G
  gi|3359707   591    ALRVATHDLsNHAAMTASG    609

Fil_haemagg: domain 6 of 16, from 633 to 692: score 13.6, E = 0.0034
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      +a Gal+l+ + ++ +  n+ +    + g +  aa+ l+n+gtl+++
  gi|3359707   633    AARGALRLRGG-AELSIVNQPAGFMLAHG-QDIAAARLANAGTLSST 677

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                   ++             +tlt aGdL+nsg
  gi|3359707   678 AA-------------GTLTLAGDLENSG    692

Phage_Capsid_P3: domain 1 of 1, from 649 to 677: score 2.0, E = 0.98
                CS    -TTEEEEEEEEEEEESB-----------X
                   *->vnqnARLLmGYEYFtsRtELvnAGtistt<-*
                      vnq A  ++      +   L nAGt+s t
  gi|3359707   649    VNQPAGFMLAHGQDIAAARLANAGTLSST    677

Fil_haemagg: domain 7 of 16, from 695 to 766: score 20.9, E = 2.7e-05
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      ++a+al+++a   +  + ++ +t   g+gltl a   l+ +g ++++
  gi|3359707   695    QSAHALRVQAR--RIASRGRLATATDGAGLTLAADKNLTLAGRTSSA 739

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                   g+ll+ aga+  ++ +++ta+ dL++++
  gi|3359707   740 GALLARAGAALVNE-GIVTARQDLSWRA    766

RHS_repeat: domain 1 of 2, from 709 to 723: score 1.3, E = 17
                   *->aagrltaltlttstdasGtt<-*
                      ++grl ++     td++G t
  gi|3359707   709    SRGRLATA-----TDGAGLT    723

DCD: domain 1 of 1, from 730 to 762: score -0.8, E = 5.9
                   *->lsmTtdadsaGILpdraIeaLiaeekiksaadf<-*
                      l++   + saG L  ra +aL+ e+ ++   d+
  gi|3359707   730    LTLAGRTSSAGALLARAGAALVNEGIVTARQDL    762

Fil_haemagg: domain 8 of 16, from 777 to 833: score 13.3, E = 0.0043
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      v+a +++ +a+       ++ g++++  +l+l aa+++  +g+l a
  gi|3359707   777    VVARSVDMRAG----QGFDHRGAIGSVTDLVLKAARIDS-AGVLRAN 818

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                   +d             +  +a++++++++
  gi|3359707   819 QD-------------IDMHADDAMRLTA    833

YopE: domain 1 of 1, from 797 to 806: score 2.0, E = 8.4
                   *->GaaselVasA<-*
                      G++++lV++A
  gi|3359707   797    GSVTDLVLKA    806

DUF1840: domain 1 of 1, from 816 to 831: score 2.2, E = 3.9
                   *->kaaaDviMlgDlArqL<-*
                      +a +D+ M++D+A++L
  gi|3359707   816    RANQDIDMHADDAMRL    831

Fil_haemagg: domain 9 of 16, from 838 to 892: score 8.1, E = 0.13
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                       a+ +l l+a+     +l ++g  +ag++lt+ta++l n+g   +a
  gi|3359707   838    LAGRDLALAAD-----QLEQSGMAQAGRTLTATAGALENDGLLDAAD 879

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                                  + l ++ ++ n+g
  gi|3359707   880 ---------------AKLRTTRAFVNRG    892

Fil_haemagg: domain 10 of 16, from 900 to 950: score 19.3, E = 8e-05
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      +a G            + +n+g l +g+ l+l+aag l n+g+ +a+
  gi|3359707   900    QAQGP-----------QIRNAGVLRTGALLALQAAGRLENTGGMAAS 935

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                   g              l+++aaG + nsg
  gi|3359707   936 GS-------------LSIAAAGPFANSG    950

Fil_haemagg: domain 11 of 16, from 953 to 980: score 0.9, E = 16
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEE
                   *->vaaGaltlnaaalggldlnnggtlsagggltlt<-*
                       a G+ + + +     + +n+g++s gg+l+l+
  gi|3359707   953    GANGDASFALS-----SFANTGSISVGGDLALR    980

DUF500: domain 1 of 1, from 961 to 977: score 1.2, E = 9.4
                   *->AVdaFaelGniTLGGnl<-*
                      A+ +Fa+ G+i+ GG+l
  gi|3359707   961    ALSSFANTGSISVGGDL    977

DUF490: domain 1 of 1, from 962 to 992: score 4.1, E = 0.61
                   *->ldglsatglgggGtltasGtldlaadlpadLtikgvadn<-*
                      l+++++t     G++++ G+l l + ++++Lt+++  d+
  gi|3359707   962    LSSFANT-----GSISVGGDLAL-RLPDVELTLDA--DH    992

PAD_M: domain 1 of 1, from 984 to 991: score 1.1, E = 4
                CS    EEEEEE--
                   *->VdISLdvD<-*
                      V+++Ld+D
  gi|3359707   984    VELTLDAD    991

Nrf1_activ_bdg: domain 1 of 1, from 995 to 1010: score 3.3, E = 2.9
                   *->avgsegqavqVvtLqd<-*
                      +v++ ++++qV++L++
  gi|3359707   995    PVSQGTTLLQVASLDN    1010

Fil_haemagg: domain 12 of 16, from 996 to 1053: score 22.4, E = 1e-05
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                      v+ G+  l++a     +l+n +   + g l+++a+g++ n++tl+ag
  gi|3359707   996    VSQGTTLLQVA-----SLDNRARSETPGRLSVQARGAIRNQDTLAAG 1037

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                   ++l               +aa+d++n++
  gi|3359707  1038 QGL------------WLESAANDIENGA    1053

HutD: domain 1 of 2, from 1033 to 1061: score 4.9, E = 0.46
                CS    EE-TT-EE..EEESEEEE.EEEEE.EEEEEEEE
                   *->aLeagDtLmvwlegdalllpltgdrgrlllvsl<-*
                      +L+ag++L  wle +a+     g+ g+ll+
  gi|3359707  1033    TLAAGQGL--WLESAAND-IENGA-GALLWSGA    1061

Ery_res_leader1: domain 1 of 1, from 1056 to 1069: score 5.8, E = 2.6
                   *->LrisGaAllhLrLRTdR<-*
                      L+ sGa    LrLR  R
  gi|3359707  1056    LLWSGAD---LRLRGTR    1069

Fil_haemagg: domain 13 of 16, from 1058 to 1105: score 4.4, E = 1.6
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                       ++ +l+l+ +    + ++ ++++ +++g +l+a++ + ng
  gi|3359707  1058    WSGADLRLRGT---RIINREAAIIESAAGMVLDARAEIDNGL----- 1096

                CS ...............EEEEESS-E
                   gdlllaagalrnggdltltaaGdL<-*
                                  +++ a+GdL
  gi|3359707  1097 ---------------GIIRAGGDL    1105

CDC48_N: domain 1 of 1, from 1082 to 1100: score 4.4, E = 1.2
                CS    EEEEEEECCS.-HHCCCCC
                   *->klrvleAdkcpddgrsiar<-*
                      + +vl A+  +d+g++i+r
  gi|3359707  1082    AGMVLDARAEIDNGLGIIR    1100

Androgen_recep: domain 1 of 1, from 1126 to 1159: score 0.9, E = 1.6
                CS    XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
                   *->GggaaGPgsGsPsAaAsssWHtLFtAEEGqLYGP<-*
                      G +a G g+     +A+  WH LFtA    +  P
  gi|3359707  1126    GDAAIGGGTYDHYHSAAVVWHELFTAGAAGIRVP    1159

Fil_haemagg: domain 14 of 16, from 1163 to 1203: score 6.2, E = 0.46
                CS    E-S-EEESSE...................EEEEESS-EEEEE
                   *->aagllnnggtliaggdlllaagalrnggdltltaaGdLdnsg<-*
                      ++ +  ++ + +agg+l l+ g + +++ ++++++G+++ +g
  gi|3359707  1163    GKDVRVAQSVVQAGGNLHLNQGEQKGRQ-ARVSNQGRIEAAG    1203

RNR_N: domain 1 of 1, from 1244 to 1267: score 2.1, E = 5.7
                   *->ehVkPNtvkFhSlkErldYLvend<-*
                      +hV P + kFh+l E ld+L++n+
  gi|3359707  1244    QHVIPAFWKFHTLYEFLDFLLSNN    1267

Nairovirus_M: domain 1 of 1, from 1308 to 1318: score -1.9, E = 5.5
                   *->PTVLDRQSRAS<-*
                      P VLD Q  AS
  gi|3359707  1308    PPVLDPQAKAS    1318

HMG_box: domain 1 of 1, from 1322 to 1338: score 5.0, E = 1.4
                CS    HHHHHHHHHHHHTTTSS
                   *->kAekeKarYekempeYk<-*
                      +A+++ a+Y+k+++eY
  gi|3359707  1322    AAQALVAQYHKDLAEYA    1338

Strep_SA_rep: domain 1 of 2, from 1322 to 1345: score 5.7, E = 1.3
                   *->adYeakLaqYqadLArvqkanaane<-*
                       + +a +aqY +dLA++  a++a++
  gi|3359707  1322    -AAQALVAQYHKDLAEYATALEAAQ    1345

HSP70: domain 1 of 1, from 1323 to 1392: score -0.4, E = 7.4
                CS    HHHHHHHHHHHHCCC..CCCGCHHHHH...HHHHHHHHHHHHHC
                   *->aEslvystekslkeeklgdkvpeadkeklekkveeaiewLkeal...
                      a +lv + +k l+e+  + ++++a+ +   +++++a + +  al+ +
  gi|3359707  1323    AQALVAQYHKDLAEY--ATALEAAQRA---EAIRTARQRVNGALrar 1364

                CS TSSH...HHHHHHHHHHHHHCHHHHHHHHHH
                   egeddAekeeikakteeLqkvvqkigekmyq<-*
                    ge     + ++++t e++++v ++++ + +
  gi|3359707  1365 YGEK---LAQLETRTPEVDAAVAALAQTIFD    1392

PseudoU_synth_1: domain 1 of 1, from 1326 to 1366: score 1.7, E = 4.4
                CS    EEES--...-......-HHHHCHHHHHCS-T...TEEEEEEEE.--T
                   *->fvgpgeltdrdnkndippedikkaLnaklrddladigvltapaegly
                      +v + ++         ++ +++ aL+a  r+     +++ta++ ++
  gi|3359707  1326    LVAQYHK---------DLAEYATALEAAQRA----EAIRTARQRVNG 1359

                CS T--TTTT
                   lfharys<-*
                     +ary
  gi|3359707  1360 ALRARYG    1366

MmgE_PrpD: domain 1 of 1, from 1328 to 1357: score 2.1, E = 2.3
                CS    HHHHHHHHHHHH-----HH....HHHHHHHHHH
                   *->dsvteqlAdyvasvryssippdvealetakkll<-*
                       ++ ++lA+y++++++ +  +  ea++ta+ ++
  gi|3359707  1328    AQYHKDLAEYATALEAAQR-A--EAIRTARQRV    1357

Birna_VP3: domain 1 of 1, from 1350 to 1382: score -0.2, E = 9.4
                   *->iEtqsRevskeqrkklerklrnlEakGpdqdeA<-*
                      i t    v +  r++++ kl  lE + p +d A
  gi|3359707  1350    IRTARQRVNGALRARYGEKLAQLETRTPEVDAA    1382

Transposase_20: domain 1 of 1, from 1351 to 1365: score 1.6, E = 9.9
                   *->lvalrhvedpgsrafy<-*
                      ++a+++v +++ ra+y
  gi|3359707  1351    RTARQRV-NGALRARY    1365

AD: domain 1 of 1, from 1360 to 1393: score 3.1, E = 4.2
                   *->rvkkrlrkaieeakkkaariGkGVSpEGQeiFda<-*
                      ++++r+ ++++++++    + ++V +++Q iFda
  gi|3359707  1360    ALRARYGEKLAQLETRTPEVDAAVAALAQTIFDA    1393

Adeno_PIX: domain 1 of 1, from 1360 to 1386: score 1.0, E = 9.2
                   *->alrslaedkLtvllaeLeeLsqqLaal<-*
                      alr+++ +kL++l+++  e+ +++aal
  gi|3359707  1360    ALRARYGEKLAQLETRTPEVDAAVAAL    1386

Strep_SA_rep: domain 2 of 2, from 1363 to 1387: score 3.6, E = 5.3
                   *->adYeakLaqYqadLArvqkanaane<-*
                      a Y +kLaq ++   +v +a aa
  gi|3359707  1363    ARYGEKLAQLETRTPEVDAAVAALA    1387

CbbQ_C: domain 1 of 2, from 1377 to 1391: score 0.2, E = 23
                   *->pdvreaLmeviaayF<-*
                      p+v +a  ++++++F
  gi|3359707  1377    PEVDAAVAALAQTIF    1391

Enolase_N: domain 1 of 1, from 1381 to 1402: score 3.2, E = 1.5
                CS    EEEEEEEEEE-TTS-EEEEEEE
                   *->aivkihaRqilDsRGnPtvEvd<-*
                      a v+  a  i+D+R  P  Ev+
  gi|3359707  1381    AAVAALAQTIFDARAKPAAEVE    1402

Head_binding: domain 1 of 1, from 1609 to 1625: score -0.0, E = 9.9
                   *->dGaqqfYfanvLKYdPdkf<-*
                      +G+q  Yf ++++Y Pd+
  gi|3359707  1609    YGSQ--YFFDLIRYRPDRP    1625

DUF166: domain 1 of 1, from 1637 to 1649: score -0.2, E = 6.3
                   *->yihkeAvekAlgk<-*
                      ++++e++ +A+g
  gi|3359707  1637    RLIREQIARAMGG    1649

DUF547: domain 1 of 1, from 1687 to 1697: score 1.0, E = 9.6
                   *->eqLdeAareFi<-*
                      eqL +A+r+F+
  gi|3359707  1687    EQLARATRDFV    1697

IATP: domain 1 of 1, from 1720 to 1734: score 3.7, E = 3.5
                CS    XXXXXXXXXXXXHHH
                   *->vggggGSiReAGGaF<-*
                      +++++++ ReAGGa+
  gi|3359707  1720    TRTAASATREAGGAL    1734

FlgD: domain 1 of 1, from 1735 to 1775: score 3.1, E = 1.6
                   *->lvpddgntivyldgggggetGsvsggvevgLpsaAdevtvtikDaaG
                      +++  + t v  d+gg        +++e    +aA  ++ +i+DaaG
  gi|3359707  1735    MAS--AGT-VLADTGG--------ATIE--SGNAAFLGKDVILDAAG 1768

                   qvvnvelVrtid<-*
                          Vr i+
  gi|3359707  1769 GA-----VRLIN    1775

DGOK: domain 1 of 1, from 1794 to 1807: score -1.6, E = 7.9
                   *->GDeAvrAGLseiAr<-*
                      GD A++ GL  +A
  gi|3359707  1794    GDIAIQGGLLDAAQ    1807

CbbQ_C: domain 2 of 2, from 1795 to 1810: score 3.5, E = 2.6
                   *->DpdvreaLmeviaayF<-*
                      D++++ +L+++++ay+
  gi|3359707  1795    DIAIQGGLLDAAQAYL    1810

RHS_repeat: domain 2 of 2, from 1823 to 1836: score 6.9, E = 0.41
                   *->YDanGrLtavtdps<-*
                      YDa+GrL++  d++
  gi|3359707  1823    YDAHGRLVSRHDAR    1836

Pentapeptide: domain 1 of 1, from 1912 to 1926: score 6.4, E = 0.87
                CS    -B-TT-B-TT-B-TC
                   *->AnLsgAnLsganLsg<-*
                      A ++g+ L+ga+Lsg
  gi|3359707  1912    ARFQGTTLEGAILSG    1926

CagZ: domain 1 of 1, from 1963 to 1976: score -1.1, E = 8.6
                   *->ERAKLitqLKLnLE<-*
                      ERAK + qL L L
  gi|3359707  1963    ERAKRVDQLFLRLS    1976

DUF2604: domain 1 of 1, from 1971 to 1981: score 5.2, E = 0.39
                   *->LFLsLKAGvAG<-*
                      LFL L AG++G
  gi|3359707  1971    LFLRLSAGAGG    1981

Fil_haemagg: domain 15 of 16, from 1975 to 2011: score 1.2, E = 13
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                       +aGa+ ++a+            lsa+ggl+++a++  ++g+ ++ag
  gi|3359707  1975    LSAGAGGYEAD----------IDLSAQGGLQAHAGRGQTAGAHTSAG 2011

                CS
                   <-*

  gi|3359707     -     -

DNA_pol3_theta: domain 1 of 1, from 2016 to 2029: score 1.9, E = 7.4
                   *->ReYFReRLiyyRis<-*
                      ReY R R ++yR++
  gi|3359707  2016    REYERSRMTHYRNA    2029

Glyco_hydro_70: domain 1 of 1, from 2055 to 2066: score -1.4, E = 8.5
                   *->RYGKGAmsAsDt<-*
                      RYG+GA +AsD+
  gi|3359707  2055    RYGQGASTASDA    2066

HutD: domain 2 of 2, from 2124 to 2158: score 1.1, E = 5.7
                CS    ..E EEE-TT-EE..EEESEEEE.EEEEE.EEEEEEEE
                   *->yae.iaLeagDtLmvwlegdalllpltgdrgrlllvsl<-*
                      ++e+ +++ag+++ ++++g+a++ ++t+d  ++l +s
  gi|3359707  2124    RQEnSQVHAGQMA-LQMAGEATQ-LATAD-TAALSISA    2158

DUF541: domain 1 of 1, from 2137 to 2159: score 3.6, E = 1.6
                   *->itVsGegevevaPDlAtislgVe<-*
                      +++ Ge+++ +++D+A +s++ +
  gi|3359707  2137    LQMAGEATQLATADTAALSISAT    2159

Fil_haemagg: domain 16 of 16, from 2178 to 2235: score 6.0, E = 0.52
                CS    EEESSEEEETTTT--EE-SEEEEEE-STT-EEEE-S-EEESSE....
                   *->vaaGaltlnaaalggldlnnggtlsagggltltaagllnnggtliag
                        +G+l+++a+  ++  + +g + s +  l+l+a++++n++++ +
  gi|3359707  2178    YFGGNLDIRAT--EQDIDLAGTKFSGADWLALSARRDINARAAASRY 2222

                CS ...............EEEEESS-EEEEE
                   gdlllaagalrnggdltltaaGdLdnsg<-*
                                     + + dL++++
  gi|3359707  2223 ---------------QESSEQHDLQLHQ    2235

Antimicrobial10: domain 1 of 1, from 2278 to 2297: score 2.1, E = 6.5
                   *->mqLnLGAiFLALLGvmAtAt<-*
                      mqL L A  L L   mA At
  gi|3359707  2278    MQLHLQARDLNLDASMAQAT    2297

LEA_1: domain 1 of 1, from 2291 to 2315: score 8.5, E = 0.18
                   *->AKeAqA.eMdlheaKAehaaekesa<-*
                      A +AqA+ Mdlh a+  +a+ ++ +
  gi|3359707  2291    ASMAQAtRMDLHVARDLNAISRQDE    2315

WWbp: domain 1 of 1, from 2321 to 2332: score 4.1, E = 0.42
                   *->pGGGwePslGea<-*
                       GG+we slG+a
  gi|3359707  2321    TGGNWEVSLGAA    2332

ATP-grasp_2: domain 1 of 1, from 2390 to 2398: score -0.4, E = 8.3
                CS    EESS-SS-H
                   *->ASseGGvdI<-*
                      AS +GG+dI
  gi|3359707  2390    ASNQGGMDI    2398

TraS: domain 1 of 1, from 2394 to 2413: score -0.0, E = 5.5
                   *->GGvDsAGntFGsrWqDydDR<-*
                      GG+D AG        Dy D+
  gi|3359707  2394    GGMDIAGRIHAEALHDYRDK    2413

GlcNAc_2-epim: domain 1 of 1, from 2468 to 2494: score 1.5, E = 1.6
                   *->Gf.FerLdad.GqpldadrrrivlqaR<-*
                      G+  +rL+ d+Gq++  +r++++++aR
  gi|3359707  2468    GGvSGRLNTDsGQAVVVQRAERWASAR    2494

Sigma70_r2: domain 1 of 1, from 2571 to 2594: score 2.7, E = 7.4
                CS    HHHHHHHHHHHHHHHHCHHTSTTT
                   *->sfstwLytiarnaiiDalRrrrrr<-*
                      ++ +++ +++ ++++D +  ++++
  gi|3359707  2571    KPDPYANRVIVQLAQDDVATQAAQ    2594

DUF1643: domain 1 of 1, from 2616 to 2621: score 0.2, E = 3
                   *->ePcHPL<-*
                      +P+HPL
  gi|3359707  2616    HPVHPL    2621

FbpA: domain 1 of 1, from 2648 to 2666: score -0.3, E = 6.5
                   *->DGlflhAvveELqetLvgg<-*
                      DG  ++ +v++L+e+L ++
  gi|3359707  2648    DGHAVAGIVQQLRERLPPA    2666

DUF1109: domain 1 of 1, from 2695 to 2706: score 0.1, E = 6.3
                   *->LAadlepvrrrav<-*
                      LAad ++ +++av
  gi|3359707  2695    LAAD-GAQPAPAV    2706

Flavi_M: domain 1 of 1, from 2778 to 2784: score 1.2, E = 8.8
                   *->LvAPAYa<-*
                      +vAP+Y+
  gi|3359707  2778    AVAPSYS    2784

BtpA: domain 1 of 1, from 2826 to 2844: score 0.8, E = 5.3
                   *->inNevDieRvskivelake<-*
                      + + +++eR ++ +++ak
  gi|3359707  2826    FSERISAERIRRVANAAKV    2844

PACT_coil_coil: domain 1 of 1, from 2830 to 2843: score 2.6, E = 5.2
                   *->iAieRqsltdyMrklankWh<-*
                      i++eR      +r+ an++
  gi|3359707  2830    ISAER------IRRVANAAK    2843

DUF2322: domain 1 of 1, from 2911 to 2926: score 1.9, E = 5.8
                   *->DRLLavieggeaLrvkvv<-*
                      D+LLav++   ++++ ++
  gi|3359707  2911    DLLLAVAA--SGYSYPLI    2926

AIRS: domain 1 of 2, from 3049 to 3058: score 0.0, E = 16
                CS    HT-EEEE--E
                   *->aGvplvGGdT<-*
                      +Gvpl+GG+
  gi|3359707  3049    HGVPLLGGNA    3058

SPO11_like: domain 1 of 1, from 3213 to 3227: score 2.7, E = 3.9
                   *->GarEtPaGasRlAss<-*
                      Ga+E PaG  Rl s
  gi|3359707  3213    GAGEIPAGWLRLPSG    3227

FERM_C: domain 1 of 1, from 3322 to 3332: score 0.9, E = 9.5
                CS    HHHHHHCT-XX
                   *->qfylrrrqpdt<-*
                      qf+++rrq d
  gi|3359707  3322    QFHAQRRQGDV    3332

DUF2600: domain 1 of 1, from 3385 to 3399: score 4.1, E = 0.45
                   *->PaVhreLakWrerAg<-*
                      P++ r+L+ Wrer+g
  gi|3359707  3385    PVSARYLQQWRERSG    3399

UPF0164: domain 1 of 1, from 3427 to 3456: score 1.0, E = 3.3
                   *->vaFsssvategnasqlvrasaGaEiqilsvrld<-*
                      v ++s + te  a  l r +aG    +++++ d
  gi|3359707  3427    VQLGSDALTETAARRLFRKHAG---NASWYSQD    3456

DPBB_1: domain 1 of 1, from 3445 to 3462: score 4.1, E = 1.8
                   *->taaGsaslYGndghGrkt<-*
                        aG+as+Y++d  G++t
  gi|3359707  3445    KHAGNASWYSQDASGQLT    3462

Treacle: domain 1 of 1, from 3458 to 3473: score -1.0, E = 5.9
                   *->vGKTsQvRAAlAPAKE<-*
                      +G   QvRA  APA
  gi|3359707  3458    SGQLTQVRAPAAPAIG    3473

Acylphosphatase: domain 1 of 1, from 3494 to 3521: score 4.3, E = 2.3
                CS    E-H HHHHHHHHHHHH--STT-EEEEEE
                   *->Gde.aaleeFlewLkkggPPlarVdkve<-*
                      G+  a+l+++l ++++g P  arV+k++
  gi|3359707  3494    GNNaAQLARMLPAIRQGLPATARVEKIT    3521

DUF2398: domain 1 of 1, from 3535 to 3548: score -0.2, E = 5.2
                   *->lRaLLanPLitrtg<-*
                      lR LL+ +Lit +g
  gi|3359707  3535    LRNLLNHYLITTAG    3548

Autoind_synth: domain 1 of 2, from 3556 to 3571: score 0.3, E = 9.3
                CS    EEEEEE-SHHHHHHHH
                   *->yLlhLPvDddsqqrLl<-*
                      y+++ +vD++ +++++
  gi|3359707  3556    YAGRVDVDAAGHKQIV    3571

DUF1885: domain 1 of 1, from 3641 to 3646: score 1.3, E = 5.2
                   *->MyFykr<-*
                      MyFy+r
  gi|3359707  3641    MYFYSR    3646

AHS2: domain 1 of 1, from 3745 to 3764: score 2.4, E = 1.4
                   *->LAQlrPGdkvRFvpvsleeA<-*
                      LAQ+ PG+ +R  ++  + A
  gi|3359707  3745    LAQAAPGSEIRLSEILQAVA    3764

HupH_C: domain 1 of 1, from 3752 to 3764: score 5.1, E = 0.24
                   *->SaeRLrEileale<-*
                      S+ RL+Eil+a++
  gi|3359707  3752    SEIRLSEILQAVA    3764

HTH_DeoR: domain 1 of 1, from 3755 to 3767: score 2.1, E = 9.3
                   *->RieqIlelLkqqG<-*
                      R+ +Il+ ++q+G
  gi|3359707  3755    RLSEILQAVAQSG    3767

LolB: domain 1 of 1, from 3801 to 3807: score 1.7, E = 5.1
                   *->kLvvdqW<-*
                      kLv++qW
  gi|3359707  3801    KLVMQQW    3807

AIRS: domain 2 of 2, from 3897 to 3907: score 2.7, E = 2.3
                CS    HHT-EEEE--E
                   *->eaGvplvGGdT<-*
                      ++ +pl+GG+T
  gi|3359707  3897    RYSMPLLGGNT    3907

HNF-1A_C: domain 1 of 1, from 3912 to 3919: score 1.8, E = 4.4
                   *->TQMaSSsQ<-*
                      +QMaS+ Q
  gi|3359707  3912    AQMASTIQ    3919

Autoind_synth: domain 2 of 2, from 3966 to 3981: score 0.1, E = 10
                CS    EEEEEE-SHHHHHHHH
                   *->yLlhLPvDddsqqrLl<-*
                      y+++ +vD++ ++r++
  gi|3359707  3966    YAGRVDVDAAGRKRIV    3981

Type_III_YscP: domain 1 of 1, from 3983 to 4002: score 1.9, E = 4.1
                   *->tqGDllsPLHEPRYLagmeAastaqs<-*
                       qG+l     EPR +ag+++a++ ++
  gi|3359707  3983    -QGGLN----EPRAQAGASTAAP-NR    4002

DUF1205: domain 1 of 1, from 4048 to 4055: score 0.1, E = 8.4
                   *->geLPdNVR<-*
                      +eLP+N R
  gi|3359707  4048    TELPANTR    4055

DUF917: domain 1 of 1, from 4092 to 4102: score 0.9, E = 4.7
                   *->LseeDleaiat<-*
                      +s++Dlea+a+
  gi|3359707  4092    ISPTDLEALAR    4102

Carot_N: domain 1 of 1, from 4202 to 4212: score 0.6, E = 9.1
                   *->RdiVvdMGfDp<-*
                      R+++ +MG+ +
  gi|3359707  4202    RNAAAEMGASA    4212

//