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:            172063333.fa
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query sequence: gi|172063333|ref|YP_001810984.1|
Accession:      [none]
Description:    FHA domain-containing protein [Burkholderia ambifaria MC40-6]

Scores for sequence family classification (score includes all domains):
Model           Description                             Score    E-value  N
--------        -----------                             -----    ------- ---
LysM            LysM domain                              23.4    8.5e-06   1
Peptidase_M16_C Peptidase M16 inactive domain             8.9       0.05   2
DUF2390         Protein of unknown function (DUF2390)     8.1      0.052   1
DUF1013         Protein of unknown function (DUF1013)     6.1      0.061   1
PT              PT repeat                                 8.5       0.13   2
SSV1_ORF_D-335  ORF D-335-like protein                    6.5       0.17   1
DUF1974         Domain of unknown function (DUF1974)      3.9       0.22   1
Type_III_YscX   Type III secretion system YscX (type_     5.5       0.32   1
LR8             LR8 protein                               5.3       0.52   1
DUF2125         Uncharacterized protein conserved in      1.9        0.7   1
PSCyt1          Planctomycete cytochrome C                3.7        1.4   1
PNPase          Polyribonucleotide nucleotidyltransfe     2.0        1.4   1
Phage_tail_T    Minor tail protein T                      2.5        1.5   1
DUF1094         Protein of unknown function (DUF1094)     3.9        1.7   1
XPC-binding     XPC-binding domain                        4.7        1.8   1
GumN            GumN protein                              2.8          2   1
Chlam_PMP       Chlamydia polymorphic membrane protei     6.2        2.2   1
Apidaecin       Apidaecin                                 2.9        2.3   1
HyaE            Hydrogenase-1 expression protein HyaE     2.3        2.4   1
ABM             Antibiotic biosynthesis monooxygenase     4.4        2.6   1
EST1_DNA_bind   Est1 DNA/RNA binding domain               3.1        2.8   1
Peptidase_C26   Peptidase C26                            -1.9        2.8   1
DUF1112         Protein of unknown function (DUF1112)     1.9        2.9   1
Aminotran_4     Aminotransferase class IV                 2.1        2.9   1
Opi1            Transcription factor Opi1                 1.4          3   1
HTS             Homoserine O-succinyltransferase         -0.1        3.2   1
DUF438          Family of unknown function (DUF438)       2.4        3.2   1
MTP18           Mitochondrial 18 KDa protein (MTP18)      2.3        3.4   1
Hydrolase       haloacid dehalogenase-like hydrolase      2.1        3.7   1
DUF2550         Protein of unknown function (DUF2550)     2.8        3.8   1
Ribonuc_2-5A    Ribonuclease 2-5A                         2.2        3.9   1
DUF1970         Domain of unknown function (DUF1970)      0.4          4   1
Methyltransf_11 Methyltransferase domain                  3.2          4   2
CARD            Caspase recruitment domain                3.2        4.3   1
Mu-transpos_C   Mu transposase, C-terminal                3.2        4.3   1
Peptidase_M15   D-ala-D-ala dipeptidase                   1.4        4.5   1
Glyco_hydro_79n Glycosyl hydrolase family 79, N-termi    -0.1        4.6   1
Drf_FH3         Diaphanous FH3 Domain                     0.8        4.6   1
Ribosomal_S3_N  Ribosomal protein S3, N-terminal doma     2.2        4.7   1
AIRS_C          AIR synthase related protein, C-termi     2.1        4.8   1
IPPT            IPP transferase                           1.8        4.9   1
EspF            EspF protein repeat                       1.4        4.9   1
DUF2164         Uncharacterized conserved protein (DU     3.0        4.9   1
DUF2190         Uncharacterized conserved protein (DU     2.1          5   1
Toxin_13        Neurotoxin B-IV-like protein              1.9        5.1   1
Cbl_N2          CBL proto-oncogene N-terminus, EF han     2.6        5.1   1
Not1            CCR4-Not complex component, Not1          0.7        5.1   1
MutS_IV         MutS family domain IV                     2.8        5.1   1
FTHFS           Formate--tetrahydrofolate ligase         -0.4        5.3   1
DUF2134         Predicted membrane protein (DUF2134)      1.9        5.4   1
Rgp1            Rgp1                                      0.4        5.5   1
EIAV_GP90       EIAV coat protein, gp90                  -0.6        5.7   1
GerE            Bacterial regulatory proteins, luxR f     3.0        5.8   1
Oxidored_nitro  Nitrogenase component 1 type Oxidored    -0.7        5.8   1
YjeF_N          YjeF-related protein N-terminus           1.0        5.8   1
Transposase_Tc5 Tc5 transposase                           0.5        5.8   1
Leader_Thr      Threonine leader peptide                  3.6        5.9   1
SRR1            SRR1                                      3.5        5.9   1
PhyH            Phytanoyl-CoA dioxygenase (PhyH)          0.4        5.9   1
DUF952          Protein of unknown function (DUF952)      2.5        6.2   1
DUF526          Protein of unknown function (DUF526)      2.1        6.2   1
Type_III_YscP   Type III secretion system  (type_III_     1.3        6.2   1
Peptidase_M4_C  Thermolysin metallopeptidase, alpha-h     0.7        6.3   1
Plasmid_RAQPRD  Plasmid protein of unknown function (     1.2        6.5   1
DUF1990         Domain of unknown function (DUF1990)      1.4        6.7   1
CRISPR_Cse1     CRISPR-associated protein Cse1 (CRISP     0.1        6.7   1
HisG            ATP phosphoribosyltransferase             1.5        6.8   1
ELFV_dehydrog   Glutamate/Leucine/Phenylalanine/Valin     0.6        6.8   1
DUF702          Domain of unknown function (DUF702)       0.5        6.8   1
Pas_Saposin     Pas factor saposin fold                   2.0          7   1
SUR7            SUR7 protein                              1.0          7   1
DUF494          Protein of unknown function (DUF494)      0.8        7.1   1
DnaJ_C          DnaJ C terminal region                    1.8        7.1   1
CysG_dimeriser  Sirohaem synthase dimerisation region     2.4        7.2   1
DUF2296         Predicted integral membrane metal-bin     2.2        7.2   1
PRCH            Photosynthetic reaction centre, H-cha     0.8        7.2   1
tRNA-synt_2c    tRNA synthetases class II (A)            -0.3        7.3   1
Glyco_hydro_42M Beta-galactosidase trimerisation doma     0.5        7.3   1
DUF2036         Uncharacterized conserved protein (DU     0.6        7.6   1
APS_kinase      Adenylylsulphate kinase                   0.1        7.6   1
BRK             BRK domain                                2.6        7.6   1
Nha1_C          Alkali metal cation/H+ antiporter Nha     1.0        7.7   1
DUF1175         Protein of unknown function (DUF1175)    -0.1        7.9   1
Flagellin_N     Bacterial flagellin N-terminus            1.6          8   1
DUF632          Protein of unknown function (DUF632)      0.1          8   1
RAG2            Recombination activating protein 2       -1.7        8.2   1
LEA_3           Late embryogenesis abundant protein       0.9        8.3   1
PEN-2           Presenilin enhancer-2 subunit of gamm     0.3        8.4   1
OpcA            Outer membrane protein OpcA              -0.7        8.4   1
Sipho_tail      Siphovirus tail component protein        -1.5        8.6   1
CMD             Carboxymuconolactone decarboxylase fa     2.3        8.7   1
DUF151          Uncharacterised ACR, COG1259              1.2        8.8   1
Metallophos     Calcineurin-like phosphoesterase          1.5          9   1
CPT             Chloramphenicol phosphotransferase-li     0.2          9   1
NDUF_B8         NADH-ubiquinone oxidoreductase ASHI s    -0.5        9.1   1
DUF2064         Uncharacterized protein conserved in      2.0        9.1   1
gpW             gpW                                       1.6        9.2   1
Cytomega_TRL10  Cytomegalovirus TRL10 protein             0.3        9.2   1
Phage_attach    Phage Head-Tail Attachment               -1.1        9.3   1
Beta-lactamase  Beta-lactamase                            0.2        9.5   1
DUF1349         Protein of unknown function (DUF1349)     0.2        9.6   1
DUF1035         Protein of unknown function (DUF1035)     1.1        9.7   1
DUF1493         Protein of unknown function (DUF1493)     0.1        9.8   1
HTH_psq         helix-turn-helix, Psq domain              2.4        9.9   1
DUF2029         Protein of unknown function (DUF2029)     0.6         10   1

Parsed for domains:
Model           Domain  seq-f seq-t    hmm-f hmm-t      score  E-value
--------        ------- ----- -----    ----- -----      -----  -------
Plasmid_RAQPRD    1/1      31    47 ..     1    17 [.     1.2      6.5
Peptidase_M16_C   1/2      46    84 ..   158   198 .]     2.4      3.4
DUF1974           1/1      61    77 ..   114   130 ..     3.9     0.22
Phage_attach      1/1      61    76 ..     1    16 [.    -1.1      9.3
DUF702            1/1      67    75 ..   161   169 .]     0.5      6.8
DUF1013           1/1      73    88 ..   127   143 .]     6.1    0.061
Toxin_13          1/1      84   100 ..     1    17 [.     1.9      5.1
DnaJ_C            1/1      97   105 ..     1     9 [.     1.8      7.1
CARD              1/1     205   229 ..    64    89 .]     3.2      4.3
APS_kinase        1/1     433   446 ..   149   162 .]     0.1      7.6
Phage_tail_T      1/1     452   458 ..    97   103 .]     2.5      1.5
DUF1493           1/1     527   534 ..   112   119 .]     0.1      9.8
Type_III_YscX     1/1     631   659 ..    80   114 ..     5.5     0.32
DUF1112           1/1     650   658 ..   230   238 .]     1.9      2.9
DUF2036           1/1     691   700 ..     1    10 [.     0.6      7.6
OpcA              1/1     719   739 ..    95   116 ..    -0.7      8.4
DUF1035           1/1     801   816 ..    41    57 ..     1.1      9.7
DUF1990           1/1     877   908 ..     6    39 ..     1.4      6.7
SUR7              1/1     922   944 ..   188   210 .]     1.0        7
DUF2190           1/1     928   946 ..    87   107 ..     2.1        5
AIRS_C            1/1     971  1019 ..     1    62 [.     2.1      4.8
Cbl_N2            1/1    1072  1081 ..    77    86 .]     2.6      5.1
YjeF_N            1/1    1099  1109 ..     1    11 [.     1.0      5.8
Oxidored_nitro    1/1    1127  1137 ..   455   465 .]    -0.7      5.8
ELFV_dehydrog     1/1    1131  1139 ..   270   278 .]     0.6      6.8
DUF1349           1/1    1175  1188 ..   168   181 ..     0.2      9.6
DUF2164           1/1    1275  1289 ..    62    81 .]     3.0      4.9
SSV1_ORF_D-335    1/1    1281  1294 ..    17    31 ..     6.5     0.17
CMD               1/1    1508  1538 ..    58    88 .]     2.3      8.7
DUF2064           1/1    1513  1555 ..     1    48 [.     2.0      9.1
Methyltransf_11   1/2    1530  1553 ..     1    28 [.     1.3       13
NDUF_B8           1/1    1612  1632 ..     1    21 [.    -0.5      9.1
EIAV_GP90         1/1    1622  1633 ..   381   392 .]    -0.6      5.7
gpW               1/1    1646  1653 ..    61    68 .]     1.6      9.2
LR8               1/1    1647  1657 ..   220   230 .]     5.3     0.52
DUF1175           1/1    1649  1658 ..   233   242 .]    -0.1      7.9
XPC-binding       1/1    1705  1722 ..    35    59 .]     4.7      1.8
RAG2              1/1    1818  1835 ..   519   536 .]    -1.7      8.2
Pas_Saposin       1/1    1830  1839 ..    67    76 .]     2.0        7
EspF              1/1    1902  1916 ..    33    47 .]     1.4      4.9
DUF151            1/1    1929  1939 ..    54    64 ..     1.2      8.8
HTS               1/1    1997  2003 ..   327   333 .]    -0.1      3.2
DUF952            1/1    2011  2019 ..     1     9 [.     2.5      6.2
CysG_dimeriser    1/1    2014  2038 ..    40    64 .]     2.4      7.2
Ribosomal_S3_N    1/1    2038  2054 ..    50    71 .]     2.2      4.7
PSCyt1            1/1    2092  2100 ..    54    62 .]     3.7      1.4
CPT               1/1    2139  2151 ..   210   222 .]     0.2        9
PRCH              1/1    2164  2174 ..   135   145 .]     0.8      7.2
tRNA-synt_2c      1/1    2222  2233 ..   598   609 .]    -0.3      7.3
Apidaecin         1/1    2251  2284 ..     1    28 []     2.9      2.3
PT                1/2    2265  2287 ..     1    23 [.     2.5      7.3
Beta-lactamase    1/1    2312  2340 ..     1    35 [.     0.2      9.5
Methyltransf_11   2/2    2316  2336 ..    76    97 ..     1.9      9.2
DUF1094           1/1    2347  2363 ..    81    97 .]     3.9      1.7
Peptidase_M4_C    1/1    2356  2366 ..   157   167 .]     0.7      6.3
Rgp1              1/1    2369  2377 ..   564   572 .]     0.4      5.5
Peptidase_M15     1/1    2424  2446 ..    31    53 ..     1.4      4.5
Type_III_YscP     1/1    2467  2484 ..     1    25 [.     1.3      6.2
Glyco_hydro_79n   1/1    2630  2635 ..   351   356 .]    -0.1      4.6
PNPase            1/1    2640  2646 ..    78    84 .]     2.0      1.4
Aminotran_4       1/1    2697  2718 ..   236   266 .]     2.1      2.9
BRK               1/1    2793  2816 ..     1    23 [.     2.6      7.6
HTH_psq           1/1    2842  2857 ..    32    48 .]     2.4      9.9
DUF2125           1/1    2851  2860 ..   356   365 .]     1.9      0.7
Not1              1/1    2901  2925 ..   391   415 .]     0.7      5.1
SRR1              1/1    2932  2948 ..    40    56 .]     3.5      5.9
DUF632            1/1    2939  2957 ..   329   347 .]     0.1        8
Chlam_PMP         1/1    2946  2969 ..     1    19 []     6.2      2.2
Ribonuc_2-5A      1/1    2949  2963 ..    64    78 ..     2.2      3.9
Mu-transpos_C     1/1    2973  2985 ..    65    81 .]     3.2      4.3
IPPT              1/1    3000  3015 ..   248   264 .]     1.8      4.9
Peptidase_C26     1/1    3001  3014 ..   176   195 ..    -1.9      2.8
DUF526            1/1    3053  3059 ..    64    70 .]     2.1      6.2
ABM               1/1    3064  3083 ..    44    65 .]     4.4      2.6
DUF2390           1/1    3067  3104 ..    38    75 ..     8.1    0.052
HyaE              1/1    3084  3095 ..     1    14 [.     2.3      2.4
MutS_IV           1/1    3091  3106 ..     1    16 [.     2.8      5.1
PhyH              1/1    3128  3158 ..     1    32 [.     0.4      5.9
LEA_3             1/1    3241  3246 ..     1     6 [.     0.9      8.3
FTHFS             1/1    3246  3256 ..   547   557 ..    -0.4      5.3
Hydrolase         1/1    3247  3311 ..    95   169 ..     2.1      3.7
Opi1              1/1    3286  3296 ..     1    11 [.     1.4        3
DUF494            1/1    3301  3312 ..   146   157 .]     0.8      7.1
DUF2550           1/1    3334  3353 ..   117   136 .]     2.8      3.8
Cytomega_TRL10    1/1    3373  3380 ..     1     9 [.     0.3      9.2
Sipho_tail        1/1    3382  3389 ..   684   691 .]    -1.5      8.6
Drf_FH3           1/1    3412  3425 ..     1    14 [.     0.8      4.6
Glyco_hydro_42M   1/1    3413  3428 ..   215   231 .]     0.5      7.3
DUF2029           1/1    3440  3460 ..   179   199 ..     0.6       10
EST1_DNA_bind     1/1    3464  3472 ..    67    75 .]     3.1      2.8
Peptidase_M16_C   2/2    3493  3520 ..   172   198 .]     6.5     0.24
DUF438            1/1    3504  3517 ..    42    56 ..     2.4      3.2
CRISPR_Cse1       1/1    3531  3546 ..   542   557 .]     0.1      6.7
PEN-2             1/1    3558  3581 ..    80   101 .]     0.3      8.4
Leader_Thr        1/1    3655  3673 ..     1    23 []     3.6      5.9
MTP18             1/1    3716  3727 ..   176   187 .]     2.3      3.4
DUF2296           1/1    3718  3729 ..     1    12 [.     2.2      7.2
GumN              1/1    3729  3739 ..   276   286 .]     2.8        2
HisG              1/1    3779  3795 ..   171   189 .]     1.5      6.8
GerE              1/1    3869  3880 ..    47    58 .]     3.0      5.8
DUF1970           1/1    3874  3893 ..    98   117 .]     0.4        4
Transposase_Tc5   1/1    3921  3932 ..     1    13 [.     0.5      5.8
Nha1_C            1/1    3933  3938 ..     1     6 [.     1.0      7.7
DUF2134           1/1    3946  3963 ..    94   115 .]     1.9      5.4
Metallophos       1/1    3968  4082 ..    98   124 .]     1.5        9
PT                2/2    4029  4050 ..     1    23 [.     6.0     0.67
LysM              1/1    4054  4123 .]     1    45 []    23.4  8.5e-06
Flagellin_N       1/1    4084  4094 ..     1    11 [.     1.6        8

Alignments of top-scoring domains:
Plasmid_RAQPRD: domain 1 of 1, from 31 to 47: score 1.2, E = 6.5
                   *->aaaAdaaaERerLaavl<-*
                      +a  + ++E+++Laa++
  gi|1720633    31    GAPPAQTPEQQALAAAQ    47

Peptidase_M16_C: domain 1 of 2, from 46 to 84: score 2.4, E = 3.4
                CS    E ECCCCHHHHHHHHHHHHHHHHHHHHTT-SHHHHHHHHHHH
                   *->a.tdpegenldeRvlelileeleklaeeGiteeelerakaql<-*
                      a+++++  +l + + ++ +++l++l ++G ++++l++  +
  gi|1720633    46    AqAETQ--RLLQ-LAQAEMATLKRLHDQGASQADLDKQVKKA    84

DUF1974: domain 1 of 1, from 61 to 77: score 3.9, E = 0.22
                   *->SAvLKRyeDeGrpeaDl<-*
                      +A+LKR +D+G  +aDl
  gi|1720633    61    MATLKRLHDQGASQADL    77

Phage_attach: domain 1 of 1, from 61 to 76: score -1.1, E = 9.3
                   *->MadfdNLFDeAisrAD<-*
                      Ma    L D+  s+AD
  gi|1720633    61    MATLKRLHDQGASQAD    76

DUF702: domain 1 of 1, from 67 to 75: score 0.5, E = 6.8
                   *->LYDQGpehq<-*
                      L+DQG  ++
  gi|1720633    67    LHDQGASQA    75

DUF1013: domain 1 of 1, from 73 to 88: score 6.1, E = 0.061
                   *->sQidLDaaVekkAakrr<-*
                      sQ+dLD+ V+k A+
  gi|1720633    73    SQADLDKQVKK-AETAW    88

Toxin_13: domain 1 of 1, from 84 to 100: score 1.9, E = 5.1
                CS    --S--SSS-THHHHHHH
                   *->AsatWGaaYPACEnnCR<-*
                      A + WGa   A En  R
  gi|1720633    84    AETAWGAVQAAVENQMR    100

DnaJ_C: domain 1 of 1, from 97 to 105: score 1.8, E = 7.1
                   *->qrlRlaGeG<-*
                      +++R+a+eG
  gi|1720633    97    NQMRVAAEG    105

CARD: domain 1 of 1, from 205 to 229: score 3.2, E = 4.3
                CS    CCHHHHHHHHHHHTT-HHHHHHHHCC
                   *->peafqiFleaLretgqphLadllels<-*
                      ++a+ ++  aLr+ ++  +a  l+l+
  gi|1720633   205    NNAYTALTTALRH-EYSTSAATLRLD    229

APS_kinase: domain 1 of 1, from 433 to 446: score 0.1, E = 7.6
                CS    SS-SEEEETTTS-H
                   *->enPElvldttkqsv<-*
                      ++PE++l  tkq v
  gi|1720633   433    AAPEVRLQYTKQQV    446

Phage_tail_T: domain 1 of 1, from 452 to 458: score 2.5, E = 1.5
                   *->GPDSepg<-*
                      GPDS+pg
  gi|1720633   452    GPDSNPG    458

DUF1493: domain 1 of 1, from 527 to 534: score 0.1, E = 9.8
                   *->KAGRWLYD<-*
                      KAG+W  D
  gi|1720633   527    KAGKWMQD    534

Type_III_YscX: domain 1 of 1, from 631 to 659: score 5.5, E = 0.32
                   *->ldtTlaRLrsLLqesgsseLqaAaeLLestqeDeR<-*
                      +d+    +r++L  + +++L +A  LL ++q D+R
  gi|1720633   631    FDG----VRQALSSGVDPALSKA--LLSKMQNDDR    659

DUF1112: domain 1 of 1, from 650 to 658: score 1.9, E = 2.9
                   *->LLskLqsrD<-*
                      LLsk+q++D
  gi|1720633   650    LLSKMQNDD    658

DUF2036: domain 1 of 1, from 691 to 700: score 0.6, E = 7.6
                   *->dYkLLELded<-*
                      +YkL+++d+d
  gi|1720633   691    SYKLFQADPD    700

OpcA: domain 1 of 1, from 719 to 739: score -0.7, E = 8.4
                CS    EEEEEEEETTEEEEEEEEEEEE
                   *->YtYtqtLsDsLnFRvGAGLGle<-*
                      Yt+t t sD +   vGAGLGl+
  gi|1720633   719    YTFT-TGSDQMTNFVGAGLGLQ    739

DUF1035: domain 1 of 1, from 801 to 816: score 1.1, E = 9.7
                   *->vGssnAtivsLVplFYl<-*
                      vG  +  +vsLVp++Y+
  gi|1720633   801    VGG-DHAVVSLVPIYYV    816

DUF1990: domain 1 of 1, from 877 to 908: score 1.4, E = 6.7
                   *->GatraegatpleedsdlvpsGyrildvrvrlGsG<-*
                      Gat   g++ +e+  + ++++y++l++ v+ G+G
  gi|1720633   877    GAT--GGYSQYESIDAHKTEWYQHLERYVEMGTG    908

SUR7: domain 1 of 1, from 922 to 944: score 1.0, E = 7
                   *->GikaslGtkmfAftWaAVaflLl<-*
                      Gi a +G  +++++W+AV++ ++
  gi|1720633   922    GIGAPVGGLLLGAAWTAVGAGMA    944

DUF2190: domain 1 of 1, from 928 to 946: score 2.1, E = 5
                   *->gNtliGvAwaaaAagsGaatv<-*
                      g+ l G+Aw+a++ag+  a++
  gi|1720633   928    GGLLLGAAWTAVGAGM--AVG    946

AIRS_C: domain 1 of 1, from 971 to 1019: score 2.1, E = 4.8
                CS    STTSEEEEEECTT..BTT----TTSHHHHHTT---S--G.......G
                   *->kpGDvlillgsttkktGrsGlhGnglslvrklleesglpdsdepdrp
                      + GD + l+gs+       G+ G+gl  + k l+e+  ++     r
  gi|1720633   971    ARGDWISLIGSG------AGIAGGGLGVAAKSLAEASSLL-----RT 1006

                CS G....---HHHHHHH
                   avqgtlgdplleptl<-*
                    +q  ++d+l +  +
  gi|1720633  1007 TGQ--VSDLLTVASR    1019

Cbl_N2: domain 1 of 1, from 1072 to 1081: score 2.6, E = 5.1
                CS    GHHHHHHHHC
                   *->tLlrNWqtLa<-*
                      +L++NW+ L
  gi|1720633  1072    SLVKNWDALS    1081

YjeF_N: domain 1 of 1, from 1099 to 1109: score 1.0, E = 5.8
                   *->lGipglvLMEr<-*
                      +G+ +++LMEr
  gi|1720633  1099    MGLGSHTLMER    1109

Oxidored_nitro: domain 1 of 1, from 1127 to 1137: score -0.7, E = 5.8
                CS    HHHHHHHHHHH
                   *->lareianains<-*
                      +a+++ana+++
  gi|1720633  1127    VAFDVANAMFA    1137

ELFV_dehydrog: domain 1 of 1, from 1131 to 1139: score 0.6, E = 6.8
                CS    HHHHHHHCS
                   *->VAdAMkarG<-*
                      VA+AM a+G
  gi|1720633  1131    VANAMFANG    1139

DUF1349: domain 1 of 1, from 1175 to 1188: score 0.2, E = 9.6
                   *->SPereGfkVtFsef<-*
                      +P ++G+ ++Fs++
  gi|1720633  1175    APADKGLRARFSDW    1188

DUF2164: domain 1 of 1, from 1275 to 1289: score 3.0, E = 4.9
                   *->kledLDieieedlyeLEkKp<-*
                      +led    +e+ +y+LE
  gi|1720633  1275    RLED----VEYYVYALE-DA    1289

SSV1_ORF_D-335: domain 1 of 1, from 1281 to 1294: score 6.5, E = 0.17
                   *->YYVYkIEkDneGkvr<-*
                      YYVY +E   eG+vr
  gi|1720633  1281    YYVYALED-AEGNVR    1294

CMD: domain 1 of 1, from 1508 to 1538: score 2.3, E = 8.7
                CS    --HHHHHHHHHHHHHHHHHHHHHHHHHHHHH
                   *->vteeeiaealawaaayaggpaalaAlaaaee<-*
                      +++e++ae +++  +++ +p +l+A  a
  gi|1720633  1508    MSAEDLAEMVRQSGWDGKSPVILYACGAGSR    1538

DUF2064: domain 1 of 1, from 1513 to 1555: score 2.0, E = 9.1
                   *->Tlavaaalqaag..dvvlaysgdldgaaraadawalearkaaalggv
                      +++++++++++g+++v+l++ g++ +  ++a+ +       a++ gv
  gi|1720633  1513    LAEMVRQSGWDGksPVILYACGAGSRVSPLAQLL-------ANVLGV 1552

                   lvy<-*
                    vy
  gi|1720633  1553 DVY    1555

Methyltransf_11: domain 1 of 2, from 1530 to 1553: score 1.3, E = 13
                   *->LDvGcGtGllaraLarrvgpgarvtGvD<-*
                      L + cG G  + +La++    a+v GvD
  gi|1720633  1530    LYA-CGAGSRVSPLAQL---LANVLGVD    1553

NDUF_B8: domain 1 of 1, from 1612 to 1632: score -0.5, E = 9.1
                   *->daviraevlgVsflalkkriv<-*
                       avi ++  +V ++  + r v
  gi|1720633  1612    LAVIDQGAHPVNWSGVAPRPV    1632

EIAV_GP90: domain 1 of 1, from 1622 to 1633: score -0.6, E = 5.7
                   *->vktsGitPLPis<-*
                      v+ sG+ P P+s
  gi|1720633  1622    VNWSGVAPRPVS    1633

gpW: domain 1 of 1, from 1646 to 1653: score 1.6, E = 9.2
                CS    --------
                   *->RrPlGvrl<-*
                      R PlGv++
  gi|1720633  1646    RVPLGVYV    1653

LR8: domain 1 of 1, from 1647 to 1657: score 5.3, E = 0.52
                   *->tsLglylrkmf<-*
                      ++Lg+y+++m
  gi|1720633  1647    VPLGVYVWRML    1657

DUF1175: domain 1 of 1, from 1649 to 1658: score -0.1, E = 7.9
                   *->LGVYRlnfLk<-*
                      LGVY ++ L+
  gi|1720633  1649    LGVYVWRMLT    1658

XPC-binding: domain 1 of 1, from 1705 to 1722: score 4.7, E = 1.8
                   *->NPqLlqlIqqNqeeFlrLLnepaeg<-*
                      NP+Llq       eF +L+++p+ g
  gi|1720633  1705    NPELLQ-------EFADLIQAPVFG    1722

RAG2: domain 1 of 1, from 1818 to 1835: score -1.7, E = 8.2
                   *->tmkkltpaKKsFLRRLFe<-*
                      +m+  t a K F RRL e
  gi|1720633  1818    IMNGKTIAPKEFARRLLE    1835

Pas_Saposin: domain 1 of 1, from 1830 to 1839: score 2.0, E = 7
                   *->avalLEdpey<-*
                      a++lLEdp y
  gi|1720633  1830    ARRLLEDPAY    1839

EspF: domain 1 of 1, from 1902 to 1916: score 1.4, E = 4.9
                   *->PDvAqRLvqHLAEHG<-*
                      PDvA RL  H+ E
  gi|1720633  1902    PDVATRLYLHVPEAP    1916

DUF151: domain 1 of 1, from 1929 to 1939: score 1.2, E = 8.8
                   *->gvkfPRPLTHD<-*
                      g ++PRPL HD
  gi|1720633  1929    GGPPPRPLAHD    1939

HTS: domain 1 of 1, from 1997 to 2003: score -0.1, E = 3.2
                   *->VYQeTPY<-*
                      VYQ TP
  gi|1720633  1997    VYQPTPL    2003

DUF952: domain 1 of 1, from 2011 to 2019: score 2.5, E = 6.2
                   *->IataseWea<-*
                      I+ta+eW+a
  gi|1720633  2011    IFTADEWAA    2019

CysG_dimeriser: domain 1 of 1, from 2014 to 2038: score 2.4, E = 7.2
                   *->aelvlaGdeeeAealleqaLagaaa<-*
                      a+ ++a +e+ Aea++++ L+g ++
  gi|1720633  2014    ADEWAARTEQIAEAEADDSLSGKTE    2038

Ribosomal_S3_N: domain 1 of 1, from 2038 to 2054: score 2.2, E = 4.7
                   *->nknlknaSnegisrieIeRkid<-*
                      +k ++++     sr++  Rki+
  gi|1720633  2038    EKIYPPT-----SRVYTDRKIA    2054

PSCyt1: domain 1 of 1, from 2092 to 2100: score 3.7, E = 1.4
                   *->ddgeMPPpd<-*
                      d+++MPP+d
  gi|1720633  2092    DTLLMPPDD    2100

CPT: domain 1 of 1, from 2139 to 2151: score 0.2, E = 9
                CS    HHHHHHHHHTT--
                   *->peECArrIrerle<-*
                      p+E A+r++ +le
  gi|1720633  2139    PDEVAERVAPHLE    2151

PRCH: domain 1 of 1, from 2164 to 2174: score 0.8, E = 7.2
                CS    SSBSEEETTT-
                   *->GhnKIkPlrva<-*
                      G nK  Plrva
  gi|1720633  2164    GSNKHDPLRVA    2174

tRNA-synt_2c: domain 1 of 1, from 2222 to 2233: score -0.3, E = 7.3
                CS    XXXXXXXXXXXX
                   *->GdtVtavvDedr<-*
                      G+  t++vD+dr
  gi|1720633  2222    GTRLTGQVDRDR    2233

Apidaecin: domain 1 of 1, from 2251 to 2284: score 2.9, E = 2.3
                   *->lrakPEAEPG......nnRPvYiPqPRPPHPRlr<-*
                      l a P AEP ++++++  RP+  P PRP  PR+
  gi|1720633  2251    LPAGPHAEPAptttrtASRPIADPPPRPTVPRVA    2284

PT: domain 1 of 2, from 2265 to 2287: score 2.5, E = 7.3
                CS    --S--------SS-----SS-SS
                   *->qPTaqPTaqPTvqPTgqPtaqPT<-*
                       ++++P a+P++ PT++ +a PT
  gi|1720633  2265    RTASRPIADPPPRPTVPRVADPT    2287

Beta-lactamase: domain 1 of 1, from 2312 to 2340: score 0.2, E = 9.5
                CS    HHHHHHHHHH.HHHTTESEEEEEEEETTEEEEEEE
                   *->ldrlieelmaddveggipGvavavvrdGkvvvdag<-*
                      ++ +++++++      + G a+ v++ G +v  ++
  gi|1720633  2312    FRTALADRPG------VIGDAARVTQPGGTVSQHW    2340

Methyltransf_11: domain 2 of 2, from 2316 to 2336: score 1.9, E = 9.2
                   *->laedperalrEiaRVLKPGGkl<-*
                      l +d++ ++  +aRV  PGG++
  gi|1720633  2316    L-ADRPGVIGDAARVTQPGGTV    2336

DUF1094: domain 1 of 1, from 2347 to 2363: score 3.9, E = 1.7
                   *->GreaeeIaanLteAFde<-*
                      G e+  ++++ ++AFd+
  gi|1720633  2347    GLEPTRVIDDVMQAFDA    2363

Peptidase_M4_C: domain 1 of 1, from 2356 to 2366: score 0.7, E = 6.3
                CS    HHHHHHHHTT-
                   *->dVkqAfsaVGV<-*
                      dV qAf+a+G
  gi|1720633  2356    DVMQAFDAAGL    2366

Rgp1: domain 1 of 1, from 2369 to 2377: score 0.4, E = 5.5
                   *->LrFeFVtss<-*
                      Lr+eFVt++
  gi|1720633  2369    LRVEFVTHD    2377

Peptidase_M15: domain 1 of 1, from 2424 to 2446: score 1.4, E = 4.5
                CS    --TTB-SSS-EEEHHHHHHHHHH
                   *->PvdlyqvaRcLvhkeaaeaLaka<-*
                      P +l+++aR+++h e+  aL +a
  gi|1720633  2424    PLALPEGARTVAHDELGQALDAA    2446

Type_III_YscP: domain 1 of 1, from 2467 to 2484: score 1.3, E = 6.2
                   *->tqGDllsPLHEPRYLagmeAastaq<-*
                      +qGDl+        + +++A++ta+
  gi|1720633  2467    AQGDLT-------PQERVAASATAN    2484

Glyco_hydro_79n: domain 1 of 1, from 2630 to 2635: score -0.1, E = 4.6
                   *->YLDqLG<-*
                      +LD+LG
  gi|1720633  2630    WLDTLG    2635

PNPase: domain 1 of 1, from 2640 to 2646: score 2.0, E = 1.4
                CS    B-TTS-B
                   *->pRIDGRe<-*
                      +R+DGR+
  gi|1720633  2640    TRVDGRH    2646

Aminotran_4: domain 1 of 1, from 2697 to 2718: score 2.1, E = 2.9
                CS    EEEEC.TCE.CT..TTSS-HHH.HHHHHHHH
                   *->PVveidDgrlpifeGdGkpGpvdtrrLqely<-*
                      PV+ + Dg+ ++      p ++ t++L  +
  gi|1720633  2697    PVRML-DGK-RV------PPAL-TKALENAG    2718

BRK: domain 1 of 1, from 2793 to 2816: score 2.6, E = 7.6
                   *->sldgeeR.VpVinrktGkrLtGde<-*
                       ++ge R ++V+  +tGk+L+++e
  gi|1720633  2793    AVSGEQRwIHVVDDETGKVLGYAE    2816

HTH_psq: domain 1 of 1, from 2842 to 2857: score 2.4, E = 9.9
                   *->rSTLwlrrlrggyslkr<-*
                      r+ ++ rr+rgg+++ r
  gi|1720633  2842    RTVFR-RRRRGGIPVGR    2857

DUF2125: domain 1 of 1, from 2851 to 2860: score 1.9, E = 0.7
                   *->GpiPLGqiPp<-*
                      G iP+G+ P+
  gi|1720633  2851    GGIPVGRSPR    2860

Not1: domain 1 of 1, from 2901 to 2925: score 0.7, E = 5.1
                   *->LFqsvakycmglkssqqvmqsegae<-*
                      LFqs++    +l+ +++++   +a+
  gi|1720633  2901    LFQSIVPVRRPLSRAEEGDDANAAN    2925

SRR1: domain 1 of 1, from 2932 to 2948: score 3.5, E = 5.9
                   *->YDPvFnalekefleelG<-*
                      +DP + +++k++l+ ++
  gi|1720633  2932    VDPHLLPSDKSALNAVD    2948

DUF632: domain 1 of 1, from 2939 to 2957: score 0.1, E = 8
                   *->pdkevleAIksFasvvhai<-*
                      +dk +l+A++   ++v ++
  gi|1720633  2939    SDKSALNAVDFSGNSVSHL    2957

Chlam_PMP: domain 1 of 1, from 2946 to 2969: score 6.2, E = 2.2
                   *->ditFsgNsa.....aggGGAiyas<-*
                      +++FsgNs+++  +a ++GA+ +
  gi|1720633  2946    AVDFSGNSVshllsAIRNGALDGG    2969

Ribonuc_2-5A: domain 1 of 1, from 2949 to 2963: score 2.2, E = 3.9
                   *->YkgsSvrDLLRaiRN<-*
                      ++g+Sv +LL aiRN
  gi|1720633  2949    FSGNSVSHLLSAIRN    2963

Mu-transpos_C: domain 1 of 1, from 2973 to 2985: score 3.2, E = 4.3
                   *->VlVyrqdGNeerflceA<-*
                      ++V+r dG     ++eA
  gi|1720633  2973    IYVFRTDG----AVAEA    2985

IPPT: domain 1 of 1, from 3000 to 3015: score 1.8, E = 4.9
                   *->lTWFRkdsqPvhWldig<-*
                      ++WF k++q v++l+
  gi|1720633  3000    MRWFDKIPQ-VNNLHHR    3015

Peptidase_C26: domain 1 of 1, from 3001 to 3014: score -1.9, E = 2.8
                CS    TS-HHHHHHHS.--EEEEEE
                   *->rllelLglaeeeieVNSlHh<-*
                      r++      +  ++VN+lHh
  gi|1720633  3001    RWF---D--KI-PQVNNLHH    3014

DUF526: domain 1 of 1, from 3053 to 3059: score 2.1, E = 6.2
                   *->RlaaLEa<-*
                      R+aaLEa
  gi|1720633  3053    RVAALEA    3059

ABM: domain 1 of 1, from 3064 to 3083: score 4.4, E = 2.6
                CS    EEEEESSHHHHHHHHTSHHHHH
                   *->ilevWedeaAfeahlqsphfka<-*
                      ++ +W+d++A  + l+s+ f+a
  gi|1720633  3064    LHDQWKDAEA--RSLNSAPFRA    3083

DUF2390: domain 1 of 1, from 3067 to 3104: score 8.1, E = 0.052
                   *->lWLaqrGlaldaerlqaLqaaaepWreevvqPLRalRr<-*
                      +W +++ ++l+ + ++aL a++++  +ev++ LR+lRr
  gi|1720633  3067    QWKDAEARSLNSAPFRALAARTRHGYNEVLGDLRRLRR    3104

HyaE: domain 1 of 1, from 3084 to 3095: score 2.3, E = 2.4
                   *->LiqRLltrfGwpvV<-*
                      L++R  tr+G+++V
  gi|1720633  3084    LAAR--TRHGYNEV    3095

MutS_IV: domain 1 of 1, from 3091 to 3106: score 2.8, E = 5.1
                   *->gfdeeLDelrdlleel<-*
                      g++e L  lr+l++++
  gi|1720633  3091    GYNEVLGDLRRLRRNG    3106

PhyH: domain 1 of 1, from 3128 to 3158: score 0.4, E = 5.9
                CS    HHHHHHSEEEETT..-S-H HHHHHHHHHHHHH
                   *->affredGylviegErflsp.eevaalraeaerl<-*
                      +f  ++G++  e    ++ +e+++a+++ a+r+
  gi|1720633  3128    QFGLRNGAAKPEQ--YFAThEQIEAYNKQAQRP    3158

LEA_3: domain 1 of 1, from 3241 to 3246: score 0.9, E = 8.3
                   *->MARSls<-*
                      MARS+s
  gi|1720633  3241    MARSFS    3246

FTHFS: domain 1 of 1, from 3246 to 3256: score -0.4, E = 5.3
                   *->SlSdDPkllGa<-*
                      SlSdDP+l Ga
  gi|1720633  3246    SLSDDPALMGA    3256

Hydrolase: domain 1 of 1, from 3247 to 3311: score 2.1, E = 3.7
                CS    EESSEH HHHHHHHHTTHG GTSEEEECHHHTEEGSTTSTTTSHHHH
                   *->lTggdr.naeallealgla.lfdvivdsdevggvgpivvgKPkpeif
                      l++++  + + + ++ g+ + +  ++d + ++   ++  gK ++  f
  gi|1720633  3247    LSDDPAlMGALVERYGGVMpVVTLVIDPHSLA---AVL-GKAHDPDF 3289

                CS HHHHHHHTSSGGGE.....EEEESSCHHHH
                   llalerlgvkpeevakpevlmVGDgvnDap<-*
                   l+a+++lg +         +++ Dg+n+++
  gi|1720633  3290 LRAVRALGDQKN-------AVF-DGANGMI    3311

Opi1: domain 1 of 1, from 3286 to 3296: score 1.4, E = 3
                   *->DvdvlraaSaL<-*
                      D+d+lra+ aL
  gi|1720633  3286    DPDFLRAVRAL    3296

DUF494: domain 1 of 1, from 3301 to 3312: score 0.8, E = 7.1
                   *->eLlFdeeegllH<-*
                      + +Fd ++g++H
  gi|1720633  3301    NAVFDGANGMIH    3312

DUF2550: domain 1 of 1, from 3334 to 3353: score 2.8, E = 3.8
                   *->mdpdAlTGFlSWLESaPpgr<-*
                      +d++A TGF   LE aP++
  gi|1720633  3334    LDEAARTGFAIRLETAPSRA    3353

Cytomega_TRL10: domain 1 of 1, from 3373 to 3380: score 0.3, E = 9.2
                   *->mYprvmhav<-*
                       Y+rvm+a+
  gi|1720633  3373    -YVRVMSAL    3380

Sipho_tail: domain 1 of 1, from 3382 to 3389: score -1.5, E = 8.6
                   *->sWntkpPd<-*
                      sW++++Pd
  gi|1720633  3382    SWSNAHPD    3389

Drf_FH3: domain 1 of 1, from 3412 to 3425: score 0.8, E = 4.6
                CS    THHHHHHHHHHHHH
                   *->gGhelvleaiinfk<-*
                      gGhe++l+a+++ k
  gi|1720633  3412    GGHEKLLRAVLERK    3425

Glyco_hydro_42M: domain 1 of 1, from 3413 to 3428: score 0.5, E = 7.3
                CS    HHHHHHHHHHHH.TT--
                   *->fldallekllaevagle<-*
                      ++++ll+++l++ +gle
  gi|1720633  3413    GHEKLLRAVLER-KGLE    3428

DUF2029: domain 1 of 1, from 3440 to 3460: score 0.6, E = 10
                   *->plafagNQSlrGvLaRLgged<-*
                      +++f+gNQ l+G+La L +e+
  gi|1720633  3440    GADFVGNQDLTGALAKLIVEN    3460

EST1_DNA_bind: domain 1 of 1, from 3464 to 3472: score 3.1, E = 2.8
                   *->edipDrLvG<-*
                      +di++rL+G
  gi|1720633  3464    PDIINRLHG    3472

Peptidase_M16_C: domain 2 of 2, from 3493 to 3520: score 6.5, E = 0.24
                CS    HHHHHHHHHHHHT T-SHHHHHHHHHHH
                   *->elileeleklaee.Giteeelerakaql<-*
                      +l+++e+e++   +G+t++e++++++ +
  gi|1720633  3493    QLLFAEVERIGIAdGMTPAEVAALRQKI    3520

DUF438: domain 1 of 1, from 3504 to 3517: score 2.4, E = 3.2
                   *->iaeGgmppeEiqkLC<-*
                      ia+G m+p+E++ L+
  gi|1720633  3504    IADG-MTPAEVAALR    3517

CRISPR_Cse1: domain 1 of 1, from 3531 to 3546: score 0.1, E = 6.7
                   *->tArkalrrlAlrlFda<-*
                      +Arka  ++A +l++a
  gi|1720633  3531    QARKATIEHATQLWEA    3546

PEN-2: domain 1 of 1, from 3558 to 3581: score 0.3, E = 8.4
                   *->GeRleWgatgDllsFiiPDk..LG<-*
                      G R+++g + D+   ++PD+ ++G
  gi|1720633  3558    GARTFGGRWLDAYGGNVPDRvtTG    3581

Leader_Thr: domain 1 of 1, from 3655 to 3673: score 3.6, E = 5.9
                   *->mkrislnttiittieittGyGAG<-*
                       +ri+    i+  +  t G GAG
  gi|1720633  3655    -RRIA--GDITAATF-TGGAGAG    3673

MTP18: domain 1 of 1, from 3716 to 3727: score 2.3, E = 3.4
                   *->pypFDepVEeav<-*
                      p +FD+ V+++v
  gi|1720633  3716    PRVFDRVVDRLV    3727

DUF2296: domain 1 of 1, from 3718 to 3729: score 2.2, E = 7.2
                   *->wyDRilDiLvGe<-*
                      ++DR++D Lv++
  gi|1720633  3718    VFDRVVDRLVRG    3729

GumN: domain 1 of 1, from 3729 to 3739: score 2.8, E = 2
                   *->GvldlLRarGy<-*
                      G++++L+a G+
  gi|1720633  3729    GLQAQLKANGF    3739

HisG: domain 1 of 1, from 3779 to 3795: score 1.5, E = 6.8
                CS    HCHH.HHHHHHHHHHHHCH
                   *->ekedlelidelvtrLegVi<-*
                      ++++  + ++++++++gV+
  gi|1720633  3779    HTKA--IAADMLAQMQGVV    3795

GerE: domain 1 of 1, from 3869 to 3880: score 3.0, E = 5.8
                CS    T-SSHHHHHHHH
                   *->nvhsrvelirla<-*
                      n+ +r +l+r a
  gi|1720633  3869    NAENRSQLARIA    3880

DUF1970: domain 1 of 1, from 3874 to 3893: score 0.4, E = 4
                   *->sqLnnYAesLnDYLAsqAGv<-*
                      sqL   A++  DYL s   +
  gi|1720633  3874    SQLARIASATGDYLTSAVSI    3893

Transposase_Tc5: domain 1 of 1, from 3921 to 3932: score 0.5, E = 5.8
                   *->ifeLrkkFPdplp<-*
                       + L+ +FP+ +p
  gi|1720633  3921    -GRLHSHFPNATP    3932

Nha1_C: domain 1 of 1, from 3933 to 3938: score 1.0, E = 7.7
                   *->WMqRLP<-*
                      WM+RLP
  gi|1720633  3933    WMERLP    3938

DUF2134: domain 1 of 1, from 3946 to 3963: score 1.9, E = 5.4
                   *->agilGilfGlppsrtlsAtAtA<-*
                      ++i+Gi+    ++ +l ++AtA
  gi|1720633  3946    VAIFGIM----QGINLYTSATA    3963

Metallophos: domain 1 of 1, from 3968 to 4082: score 1.5, E = 9
                CS    SSSSSCC
                   *->Hgplsps........................................
                      H p+ +s+ +++++ + + ++ +  ++++  + ++ +  +  ++++
  gi|1720633  3968    HKPPGSSpqpsgsnagtpapdasfapspsatsapsqptasaapspsa 4014

                CS                                                HHH
                   ...............................................lee
                   + +++ +  +  +  +++ ++ +   + +++++ + +   +   +++l+
  gi|1720633  4015 saspsasatstagaqpssapsqptqaptnspsptptpltvnvvagdsLWT 4064

                CS HHH TTSSEEEEESSSSE
                   llk.ngvdlvirGHtHvp<-*
                   ++++++  l++++H+  +
  gi|1720633  4065 IAErHEGSLLAAAHVPAD    4082

PT: domain 2 of 2, from 4029 to 4050: score 6.0, E = 0.67
                CS    --S--------SS-----SS-SS
                   *->qPTaqPTaqPTvqPTgqPtaqPT<-*
                      qP+++P +qPT++PT+ P++ PT
  gi|1720633  4029    QPSSAP-SQPTQAPTNSPSPTPT    4050

LysM: domain 1 of 1, from 4054 to 4123: score 23.4, E = 8.5e-06
                CS    EE--TT--HHHHHHHT                        TS-IIII
                   *->YtVkkGDTLwsIArry........................gvsvseL
                      + V  GD+Lw IA r++++     + + +++++ +++++ +v ++e+
  gi|1720633  4054    VNVVAGDSLWTIAERHegsllaaahvpadeqqrmsdgqriNVALNEI 4100

                CS IITS-TT-SS. .--TT-EEEET
                   aelNpgLsspn.tlyvGQvLkiP<-*
                   ++lN++L+++ +++ +GQ+ ++
  gi|1720633  4101 LQLNRNLAGNPgEIAPGQQVIVG    4123

Flagellin_N: domain 1 of 1, from 4084 to 4094: score 1.6, E = 8
                CS    HHHHTS--TTT
                   *->qeRLSSGlRIN<-*
                      q+R+S G RIN
  gi|1720633  4084    QQRMSDGQRIN    4094

//