BLASTP 2.2.25+


Reference:
Stephen F. Altschul, Thomas L. Madden, Alejandro A. Schäffer,
Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997),
"Gapped BLAST and PSI-BLAST: a new generation of protein database
search programs", Nucleic Acids Res. 25:3389-3402.



Reference for composition-based statistics:
Alejandro A. Schäffer, L. Aravind, Thomas L. Madden, Sergei
Shavirin, John L. Spouge, Yuri I. Wolf, Eugene V. Koonin, and
Stephen F. Altschul (2001), "Improving the accuracy of PSI-BLAST
protein database searches with composition-based statistics and
other refinements", Nucleic Acids Res. 29:2994-3005.



Database: All non-redundant GenBank CDS translations+PDB+SwissProt+PIR+PRF
excluding environmental samples from WGS projects
           15,229,318 sequences; 5,219,829,388 total letters



Query= Rv2761c

Length=364
                                                                      Score     E
Sequences producing significant alignments:                          (Bits)  Value

gi|15609898|ref|NP_217277.1|  type I restriction/modification sys...   737    0.0   
gi|306807435|ref|ZP_07444103.1|  type I restriction/modification ...   735    0.0   
gi|306798669|ref|ZP_07436971.1|  type I restriction/modification ...   734    0.0   
gi|254365412|ref|ZP_04981457.1|  type I restriction/modification ...   733    0.0   
gi|340627762|ref|YP_004746214.1|  putative type I restriction/mod...   733    0.0   
gi|339295611|gb|AEJ47722.1|  type I restriction/modification syst...   575    5e-162
gi|240171166|ref|ZP_04749825.1|  putative type I restriction/modi...   488    9e-136
gi|183981974|ref|YP_001850265.1|  type I restriction/modification...   423    2e-116
gi|342858379|ref|ZP_08715034.1|  type I restriction/modification ...   423    3e-116
gi|296171678|ref|ZP_06852892.1|  type I restriction/modification ...   355    7e-96 
gi|343518295|ref|ZP_08755288.1|  type I restriction modification ...   154    1e-35 
gi|294624820|ref|ZP_06703480.1|  type I restriction-modification ...   151    2e-34 
gi|331681144|ref|ZP_08381781.1|  type I restriction-modification ...   130    3e-28 
gi|94986115|ref|YP_605479.1|  restriction modification system DNA...   126    5e-27 
gi|258546307|ref|ZP_05706541.1|  restriction modification system ...   119    8e-25 
gi|319946656|ref|ZP_08020890.1|  putative restriction modificatio...   114    2e-23 
gi|322386250|ref|ZP_08059882.1|  type I restriction system specif...   113    5e-23 
gi|217968470|ref|YP_002353704.1|  restriction modification system...   112    8e-23 
gi|315612675|ref|ZP_07887587.1|  HsdS family protein [Streptococc...   110    5e-22 
gi|319777746|ref|YP_004137397.1|  restriction modification system...   108    1e-21 
gi|238810194|dbj|BAH69984.1|  hypothetical protein [Mycoplasma fe...   105    9e-21 
gi|308190349|ref|YP_003923280.1|  hypothetical protein MFE_08350 ...   105    1e-20 
gi|332704541|ref|ZP_08424629.1|  restriction modification system ...   104    3e-20 
gi|78484678|ref|YP_390603.1|  restriction modification system DNA...   103    4e-20 
gi|53805024|ref|YP_113333.1|  type I restriction-modification sys...   103    5e-20 
gi|329114036|ref|ZP_08242800.1|  Putative type-1 restriction enzy...   102    8e-20 
gi|91773198|ref|YP_565890.1|  restriction modification system DNA...  99.0    1e-18 
gi|206896558|ref|YP_002247704.1|  type I restriction/modification...  97.1    5e-18 
gi|71065438|ref|YP_264165.1|  type I restriction modification sys...  96.3    8e-18 
gi|239827072|ref|YP_002949696.1|  restriction modification system...  95.9    9e-18 
gi|291288456|ref|YP_003505272.1|  restriction modification system...  93.6    5e-17 
gi|189499314|ref|YP_001958784.1|  restriction modification system...  90.5    4e-16 
gi|284048511|ref|YP_003398850.1|  restriction modification system...  90.1    5e-16 
gi|261366731|ref|ZP_05979614.1|  putative type I restriction modi...  88.6    1e-15 
gi|242280198|ref|YP_002992327.1|  restriction modification system...  88.6    2e-15 
gi|281421790|ref|ZP_06252789.1|  type I restriction-modification ...  87.8    3e-15 
gi|150006174|ref|YP_001300918.1|  type I restriction endonuclease...  87.0    4e-15 
gi|312875033|ref|ZP_07735051.1|  type I restriction modification ...  86.3    7e-15 
gi|336115613|ref|YP_004570380.1|  putative type I restriction mod...  85.9    1e-14 
gi|268609384|ref|ZP_06143111.1|  putative type I restriction-modi...  85.5    1e-14 
gi|310826741|ref|YP_003959098.1|  hypothetical protein ELI_1147 [...  84.7    2e-14 
gi|338767478|gb|EGP22297.1|  Type I restriction-modification syst...  84.3    3e-14 
gi|229002234|dbj|BAH57700.1|  hypothetical protein [Staphylococcu...  82.8    8e-14 
gi|300821374|ref|ZP_07101522.1|  type I restriction modification ...  82.0    2e-13 
gi|89093019|ref|ZP_01165970.1|  putative type I restriction enzym...  81.6    2e-13 
gi|158315561|ref|YP_001508069.1|  restriction modification system...  81.3    2e-13 
gi|85711748|ref|ZP_01042804.1|  type I restriction-modification s...  80.9    3e-13 
gi|229195092|ref|ZP_04321867.1|  type I restriction-modification ...  80.9    3e-13 
gi|38505922|ref|NP_942540.1|  type I site-specific deoxyribonucle...  80.5    4e-13 
gi|345342870|gb|EGW75262.1|  type I restriction modification DNA ...  80.1    5e-13 


>gi|15609898|ref|NP_217277.1| type I restriction/modification system specificity determinant 
HsdS [Mycobacterium tuberculosis H37Rv]
 gi|15842299|ref|NP_337336.1| type I restriction system specificity protein [Mycobacterium 
tuberculosis CDC1551]
 gi|31793935|ref|NP_856428.1| type I restriction/modification system specificity determinant 
HsdS [Mycobacterium bovis AF2122/97]
 69 more sequence titles
 Length=364

 Score =  737 bits (1902),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 364/364 (100%), Positives = 364/364 (100%), Gaps = 0/364 (0%)

Query  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60
            MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV
Sbjct  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60

Query  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS
Sbjct  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180
            QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS
Sbjct  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180

Query  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240
            LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR
Sbjct  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240

Query  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300
            IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR
Sbjct  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300

Query  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360
            IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL
Sbjct  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360

Query  361  PAGR  364
            PAGR
Sbjct  361  PAGR  364


>gi|306807435|ref|ZP_07444103.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu007]
 gi|306968824|ref|ZP_07481485.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu009]
 gi|308346165|gb|EFP35016.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu007]
 gi|308353575|gb|EFP42426.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu009]
Length=364

 Score =  735 bits (1898),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 363/364 (99%), Positives = 363/364 (99%), Gaps = 0/364 (0%)

Query  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60
            MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV
Sbjct  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60

Query  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS
Sbjct  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180
            QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS
Sbjct  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180

Query  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240
            LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR
Sbjct  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240

Query  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300
            IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR
Sbjct  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300

Query  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360
            IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL
Sbjct  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360

Query  361  PAGR  364
            P GR
Sbjct  361  PTGR  364


>gi|306798669|ref|ZP_07436971.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu006]
 gi|308341049|gb|EFP29900.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis SUMu006]
Length=364

 Score =  734 bits (1894),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 363/364 (99%), Positives = 363/364 (99%), Gaps = 0/364 (0%)

Query  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60
            MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV
Sbjct  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60

Query  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS
Sbjct  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180
            QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS
Sbjct  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180

Query  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240
            LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR
Sbjct  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240

Query  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300
            IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR
Sbjct  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300

Query  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360
            IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLP LMSGIDGL
Sbjct  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPPLMSGIDGL  360

Query  361  PAGR  364
            PAGR
Sbjct  361  PAGR  364


>gi|254365412|ref|ZP_04981457.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis str. Haarlem]
 gi|134150925|gb|EBA42970.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis str. Haarlem]
Length=364

 Score =  733 bits (1893),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 363/364 (99%), Positives = 363/364 (99%), Gaps = 0/364 (0%)

Query  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60
            MSRVEKVEKVRLGDHLDFSNGHTS HTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV
Sbjct  1    MSRVEKVEKVRLGDHLDFSNGHTSVHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60

Query  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS
Sbjct  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180
            QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS
Sbjct  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180

Query  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240
            LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR
Sbjct  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240

Query  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300
            IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR
Sbjct  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300

Query  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360
            IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL
Sbjct  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360

Query  361  PAGR  364
            PAGR
Sbjct  361  PAGR  364


>gi|340627762|ref|YP_004746214.1| putative type I restriction/modification system specificity determinant 
HSDS [Mycobacterium canettii CIPT 140010059]
 gi|340005952|emb|CCC45119.1| putative type I restriction/modification system specificity determinant 
HSDS (S protein) [Mycobacterium canettii CIPT 140010059]
Length=364

 Score =  733 bits (1893),  Expect = 0.0, Method: Compositional matrix adjust.
 Identities = 363/364 (99%), Positives = 363/364 (99%), Gaps = 0/364 (0%)

Query  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60
            MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV
Sbjct  1    MSRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIV  60

Query  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS
Sbjct  61   VGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180
            QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS
Sbjct  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSS  180

Query  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240
            LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR
Sbjct  181  LASRSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPR  240

Query  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300
            IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR
Sbjct  241  IPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQR  300

Query  301  IQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360
            IQPTQLLRVWV DVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL
Sbjct  301  IQPTQLLRVWVCDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGL  360

Query  361  PAGR  364
            PAGR
Sbjct  361  PAGR  364


>gi|339295611|gb|AEJ47722.1| type I restriction/modification system specificity determinant 
hsdS [Mycobacterium tuberculosis CCDC5079]
Length=287

 Score =  575 bits (1481),  Expect = 5e-162, Method: Compositional matrix adjust.
 Identities = 286/287 (99%), Positives = 287/287 (100%), Gaps = 0/287 (0%)

Query  78   VWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDR  137
            +WVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDR
Sbjct  1    MWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDR  60

Query  138  PRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACLDAQHFDST  197
            PRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACLDAQHFDST
Sbjct  61   PRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACLDAQHFDST  120

Query  198  VAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNITSLPSEMALA  257
            VAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNITSLPSEMALA
Sbjct  121  VAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNITSLPSEMALA  180

Query  258  STEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWVRDVRRL  317
            STEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWVRDVRRL
Sbjct  181  STEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWVRDVRRL  240

Query  318  TPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGLPAGR  364
            TPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGLPAGR
Sbjct  241  TPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSGIDGLPAGR  287


>gi|240171166|ref|ZP_04749825.1| putative type I restriction/modification system specificity determinant 
[Mycobacterium kansasii ATCC 12478]
Length=366

 Score =  488 bits (1255),  Expect = 9e-136, Method: Compositional matrix adjust.
 Identities = 247/349 (71%), Positives = 281/349 (81%), Gaps = 4/349 (1%)

Query  8    EKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSY  67
            E V+LGDHLDF++G       P     G +PV GANG+IGY+ Q NARGPLIV+GRVGSY
Sbjct  4    EHVKLGDHLDFAHGRAL----PTRLSDGGFPVCGANGIIGYADQPNARGPLIVIGRVGSY  59

Query  68   CGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNV  127
            CGS+ Y D DVWVTDNALACRAK+PEETRYWYYAL G GLNRYRAGSGQPLLSQ +LR+V
Sbjct  60   CGSVHYHDGDVWVTDNALACRAKRPEETRYWYYALRGCGLNRYRAGSGQPLLSQTILRDV  119

Query  128  SVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTA  187
            S  A AAPDR RIGE+LGA DDKIAAN RVIEAAEALMLAIV  +S  VPLS+LA +STA
Sbjct  120  STRAAAAPDRSRIGEVLGALDDKIAANKRVIEAAEALMLAIVHNISDRVPLSTLARKSTA  179

Query  188  CLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNI  247
              + Q FD +VA+YSF AFD GAQP+ V  RTI+S K V+++PCVLF KLNPRIPRIWN+
Sbjct  180  SRNPQEFDGSVAYYSFPAFDDGAQPTVVNSRTIKSVKFVLTRPCVLFSKLNPRIPRIWNV  239

Query  248  TSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLL  307
             SLP EMAL STEFVVL PVGVDTSALW+ALRQ DV A LR+ V G+TGSRQRIQP +LL
Sbjct  240  ISLPVEMALTSTEFVVLHPVGVDTSALWSALRQSDVGARLRRRVAGVTGSRQRIQPAELL  299

Query  308  RVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
               VRDVRRLT   A AI+ LGALC+ RR+ESA+LA+ RDALLPLL+ G
Sbjct  300  HEQVRDVRRLTAAQATAISGLGALCHVRRVESAQLAAGRDALLPLLVCG  348


>gi|183981974|ref|YP_001850265.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium marinum M]
 gi|183175300|gb|ACC40410.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium marinum M]
Length=361

 Score =  423 bits (1088),  Expect = 2e-116, Method: Compositional matrix adjust.
 Identities = 215/347 (62%), Positives = 254/347 (74%), Gaps = 4/347 (1%)

Query  10   VRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCG  69
            +RLG+HL+F     SG  SP     GR PVYGANG IG++A+HN  GPLIVVGRVGSYCG
Sbjct  1    MRLGEHLEFR----SGSCSPPRRARGRVPVYGANGPIGFAAEHNVSGPLIVVGRVGSYCG  56

Query  70   SLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSV  129
            S+RYCDSDVWVTDNA  CRA  P ETRYWYYAL    LN +RAGSGQPLL+Q  LR VSV
Sbjct  57   SVRYCDSDVWVTDNAYVCRANDPAETRYWYYALQTCRLNEHRAGSGQPLLNQRTLREVSV  116

Query  130  SAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACL  189
                AP+R RI E+LGA DDKIA N+RVIEAAEALM+A+VG + A V LS LA RST  +
Sbjct  117  HVAQAPERRRIAEVLGALDDKIANNERVIEAAEALMVAMVGSVDARVALSGLARRSTKLV  176

Query  190  DAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNITS  249
            +   FD  VAH+S  AFD GA    V G +++S K  +S+PCVLF KLNPR+PRIWN+  
Sbjct  177  NPADFDDVVAHFSLPAFDAGAHARPVAGASVKSGKFHLSEPCVLFAKLNPRVPRIWNVVR  236

Query  250  LPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRV  309
            LP +MALAS EFVVL P+GVDTS LW+ALRQP+V   L QLV G +GS QRI P  LL +
Sbjct  237  LPPQMALASCEFVVLSPLGVDTSVLWSALRQPEVSTSLAQLVAGTSGSHQRIGPKDLLDL  296

Query  310  WVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
             V DVRRL    +A I +LGALC+ RR + A+LA+ RDALLP L++G
Sbjct  297  QVPDVRRLGAAQSATITDLGALCHARRGQCAQLAALRDALLPGLVTG  343


>gi|342858379|ref|ZP_08715034.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium colombiense CECT 3035]
 gi|342134083|gb|EGT87263.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium colombiense CECT 3035]
Length=375

 Score =  423 bits (1087),  Expect = 3e-116, Method: Compositional matrix adjust.
 Identities = 221/348 (64%), Positives = 261/348 (75%), Gaps = 6/348 (1%)

Query  9    KVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYC  68
            K  LG+HLDFS G  SG  +P  +  G  PVYGANG IGY+A HNA GPLIV+GRVGSYC
Sbjct  2    KSVLGEHLDFSTG--SG--APPRDSDGSVPVYGANGAIGYAAAHNATGPLIVLGRVGSYC  57

Query  69   GSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVS  128
            GSLRYCDSDVWVT+NAL CRAK P ETRYWYYA+    L  +R+GSGQPLL+Q +L +V 
Sbjct  58   GSLRYCDSDVWVTENALVCRAKDPRETRYWYYAMQTCRLADHRSGSGQPLLNQRILHDVP  117

Query  129  VSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTAC  188
            V  VAA +R RI E+LGAFDDKIAANDRVI+ AE LM+A+       V LSSLA+RSTA 
Sbjct  118  VWTVAARERRRIAELLGAFDDKIAANDRVIDTAERLMVALAESAPDRVRLSSLATRSTAS  177

Query  189  LDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNIT  248
            +    FD+ VAH+S  AFD GA+P  V    +RSAKL+++Q CVLF KLNPRIPRIWN+ 
Sbjct  178  V--SEFDARVAHFSLPAFDDGAKPCLVPAAEVRSAKLLLTQSCVLFAKLNPRIPRIWNVV  235

Query  249  SLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLR  308
             LPS+MA+ASTEFVVL PVG+D+SALW+ALRQPDV   L++ V G +GS QRIQP  LL 
Sbjct  236  GLPSDMAVASTEFVVLVPVGIDSSALWSALRQPDVSRRLQRRVAGTSGSHQRIQPRDLLD  295

Query  309  VWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
            V V DVRRL+P     +  LGALC+ RR ESARL+  RDALLPLL+ G
Sbjct  296  VLVSDVRRLSPEALRTLTGLGALCHARRTESARLSGLRDALLPLLIRG  343


>gi|296171678|ref|ZP_06852892.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium parascrofulaceum ATCC BAA-614]
 gi|295894039|gb|EFG73802.1| type I restriction/modification system specificity determinant 
HsdS (S protein) [Mycobacterium parascrofulaceum ATCC BAA-614]
Length=301

 Score =  355 bits (911),  Expect = 7e-96, Method: Compositional matrix adjust.
 Identities = 181/293 (62%), Positives = 214/293 (74%), Gaps = 0/293 (0%)

Query  64   VGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGV  123
            +G+ CGSLRYCDSDVWVT+NAL CRAK P ETRYWYYAL    L  +R+GSGQPLL+Q +
Sbjct  1    MGAQCGSLRYCDSDVWVTENALVCRAKDPRETRYWYYALQTCRLGDHRSGSGQPLLNQRI  60

Query  124  LRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLAS  183
            LR+V+  +V AP R  I E+LGA DDK++AN+RVI AAE LM+A+     A+ PLSSLA 
Sbjct  61   LRDVAAPSVGAPRRAAIAELLGALDDKVSANERVIGAAENLMVALAESAPAHAPLSSLAR  120

Query  184  RSTACLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPR  243
            RSTA  D   FD  VAH+S  AFD GA P  V G +I+S K VV  PCVLF KLNPRIPR
Sbjct  121  RSTAFRDVTEFDDVVAHFSLPAFDAGATPRVVRGESIKSGKFVVPGPCVLFAKLNPRIPR  180

Query  244  IWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQP  303
            IWN+  LP EMALAS+EFVVL P+ VD SALWAA+RQP V A L++   G +GS QRI+P
Sbjct  181  IWNVVGLPPEMALASSEFVVLTPMAVDASALWAAVRQPGVSARLQRRAAGTSGSHQRIRP  240

Query  304  TQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
              L+ V V DVRRL       I  LGALC+ RR E+ARLA+ RDALLPLL+SG
Sbjct  241  RDLMEVAVPDVRRLPAAALRTITGLGALCHARRTENARLAAGRDALLPLLISG  293


>gi|343518295|ref|ZP_08755288.1| type I restriction modification DNA specificity domain protein 
[Haemophilus pittmaniae HK 85]
 gi|343393786|gb|EGV06338.1| type I restriction modification DNA specificity domain protein 
[Haemophilus pittmaniae HK 85]
Length=419

 Score =  154 bits (390),  Expect = 1e-35, Method: Compositional matrix adjust.
 Identities = 86/247 (35%), Positives = 133/247 (54%), Gaps = 16/247 (6%)

Query  11   RLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGS  70
             LG HL FSNG TS    P     G YPVYGANG+IGYS+++N+    I++GRVGSYCGS
Sbjct  3    NLGQHLKFSNGKTS----PNRTDTGFYPVYGANGIIGYSSEYNSDENTIIIGRVGSYCGS  58

Query  71   LRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVS  130
            + + D   WVTDNA++ R+   +E+ +W+Y L    LN  R GSGQPLL+Q +L+++S+ 
Sbjct  59   VHFSDKKCWVTDNAISARSMNLKESIFWFYYLKSLNLNSMRVGSGQPLLNQTILKSISLP  118

Query  131  AVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACLD  190
             +   ++ ++GEIL +FD KI  N ++ +  E +  A+    S +V    + ++  A  D
Sbjct  119  DLNI-NKIKVGEILLSFDQKIDLNTQINQTLEQIAQALFK--SWFVDFDPVRAKVQALSD  175

Query  191  AQHFDSTVAHYSFAAFDGGAQPSRVGG-------RTIRSAKLVVSQPCVLFPKLNPRIPR  243
                +   A  +      G  P  +         R    A+   + PC +       +P+
Sbjct  176  GLSLEQ--AELAAMQAISGKTPEELTALSQTQPDRYAELAETAKAFPCEIVEIDGVEVPK  233

Query  244  IWNITSL  250
             WN+ +L
Sbjct  234  GWNLKAL  240


>gi|294624820|ref|ZP_06703480.1| type I restriction-modification system specificity determinant 
[Xanthomonas fuscans subsp. aurantifolii str. ICPB 11122]
 gi|292600884|gb|EFF44961.1| type I restriction-modification system specificity determinant 
[Xanthomonas fuscans subsp. aurantifolii str. ICPB 11122]
Length=389

 Score =  151 bits (381),  Expect = 2e-34, Method: Compositional matrix adjust.
 Identities = 114/356 (33%), Positives = 161/356 (46%), Gaps = 44/356 (12%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G+Y V+G+NG IG++    A GP +++GR G+Y G +R+     WV D A     KK  +
Sbjct  30   GKYRVFGSNGAIGWTENALAEGPGVILGRKGAYRG-VRFWREPFWVIDTAYYVVPKKKLD  88

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
             R+ YYA+    L     GS  P  ++  +    ++  +  ++  I  +LG  DDKI  N
Sbjct  89   MRWLYYAIKHHKLGEIDDGSPIPSTTRAAVYVRELTVPSLKEQGEISYVLGVLDDKIELN  148

Query  155  DRVIEAAEALMLAIVGRLSAYVPLSSLASRSTACLD----------------AQHFDSTV  198
             R+ +  EA++ A+    S +V    +A       +                AQ    TV
Sbjct  149  RRMNQTLEAMVHALFK--SWFVDFDGVAPEDMQESELGFIPKGWQVIAFGDVAQQVKGTV  206

Query  199  ----------AHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNIT  248
                       HYS  AFD    P R  G  I+S K  V   CVL  KLNP IPRIW I 
Sbjct  207  NPMTSPEETFTHYSLPAFDVAQLPVRELGEAIKSNKTPVPNECVLVSKLNPHIPRIWLIG  266

Query  249  SLPSEMALASTEFVVLRPVGVDTSALWAAL-RQPDVLAELRQLVGGMTGSRQRIQPTQLL  307
                  A+ STEF+V  P     SA    L    +  + LRQLV G + S QR++P QL 
Sbjct  267  G-AGHNAVCSTEFIVWMPKKPANSAFVYVLASSSEFNSALRQLVTGTSNSHQRVKPEQLA  325

Query  308  RVWVRDVRRLTPGHAAAIANLGALCN-------ERRIESARLASCRDALLPLLMSG  356
             +      R+   +  AI+   A            R++S +LA  RD LLP L+SG
Sbjct  326  NI------RVIAVNDEAISKFSAQSKPLMEKLLHHRLQSQQLAQLRDTLLPKLISG  375


>gi|331681144|ref|ZP_08381781.1| type I restriction-modification system specificity determinant 
[Escherichia coli H299]
 gi|331081365|gb|EGI52526.1| type I restriction-modification system specificity determinant 
[Escherichia coli H299]
Length=434

 Score =  130 bits (328),  Expect = 3e-28, Method: Compositional matrix adjust.
 Identities = 64/154 (42%), Positives = 91/154 (60%), Gaps = 5/154 (3%)

Query  12   LGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSL  71
            +G+HL   NG +S    P+    G YPVYG+NG+IGYS ++NA    I++GRVGSYCGS+
Sbjct  10   IGEHLLIRNGKSS----PSRAITGEYPVYGSNGIIGYSDEYNANENTIIIGRVGSYCGSV  65

Query  72   RYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSA  131
                   WVTDNA+   AK   E+ +W+Y L    LN Y  GSGQPL++Q ++  +SV+ 
Sbjct  66   YISGKKCWVTDNAIIGTAKNENESHFWFYLLKKIDLNNYSTGSGQPLINQTIINTISVTI  125

Query  132  VA-APDRPRIGEILGAFDDKIAANDRVIEAAEAL  164
               +  R  IG  L  FD KI  +  + ++ E +
Sbjct  126  PKLSEKRVSIGHFLRHFDQKINLSLNINQSLEQM  159


>gi|94986115|ref|YP_605479.1| restriction modification system DNA specificity subunit [Deinococcus 
geothermalis DSM 11300]
 gi|94556396|gb|ABF46310.1| restriction modification system DNA specificity domain [Deinococcus 
geothermalis DSM 11300]
Length=417

 Score =  126 bits (317),  Expect = 5e-27, Method: Compositional matrix adjust.
 Identities = 115/376 (31%), Positives = 154/376 (41%), Gaps = 59/376 (15%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G  PVYG+NG++G+       GP IV+GR G+  G++ Y     W  D            
Sbjct  31   GSVPVYGSNGIVGFHDSALTGGPTIVIGRKGT-VGAVHYSPIPCWPIDTTFFVSDSDRSL  89

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVL--RNVSVSAVAAPDRPRIGEILGAFDDKIA  152
             RY YY L   GL    A S  P L++     R V +    A  R  I  ILG  DDKI 
Sbjct  90   VRYSYYLLKSLGLENMNADSAVPGLNRDAAHARIVLIPRDKAEQR-AIAHILGTLDDKIE  148

Query  153  ANDRVIEAAEALMLAIV---------------GRLSAYVPLSSLASRSTACLDAQHFDST  197
             N +  E  EA+  A+                GR      L  L +        +  DS 
Sbjct  149  LNRKQSETLEAMARALFKAWFVDFEPVRAKMEGRWQRGQSLPGLPAHLYDLFPDRLVDSE  208

Query  198  V---------------------------------AHYSFAAFDGGAQPSRVGGRTIRSAK  224
            +                                  HYS  AFD    PS   G  I+S K
Sbjct  209  LGEIPEGWRVFAFGDVAQQGKGVVNPGNSPQDLFTHYSLPAFDSAHCPSIEPGHAIKSNK  268

Query  225  LVVSQPCVLFPKLNPRIPRIWNI-TSLPSEMALASTEFVVLRPVGVDTSA-LWAALRQPD  282
              V    VL  KLNP IPR+W++ T+ P+  A+ STEF+V  P     SA L+     P+
Sbjct  269  TPVPDGAVLVSKLNPHIPRVWHVGTAGPN--AVCSTEFIVWAPKAPANSAFLYCLASSPE  326

Query  283  VLAELRQLVGGMTGSRQRIQPTQL--LRVWVRDVRRLTPGHAAAIANLGALCNERRIESA  340
                + QLV G + S QR++P QL  +RV+      +        + L  +   R+ +S 
Sbjct  327  FSGAMHQLVTGTSNSHQRVKPDQLREIRVFAATENAIEAFSEWVRSPLEKILQNRQ-QSR  385

Query  341  RLASCRDALLPLLMSG  356
             LA  RDALLP L+SG
Sbjct  386  TLAQLRDALLPRLISG  401


>gi|258546307|ref|ZP_05706541.1| restriction modification system DNA specificity domain protein 
[Cardiobacterium hominis ATCC 15826]
 gi|258518451|gb|EEV87310.1| restriction modification system DNA specificity domain protein 
[Cardiobacterium hominis ATCC 15826]
Length=391

 Score =  119 bits (298),  Expect = 8e-25, Method: Compositional matrix adjust.
 Identities = 63/154 (41%), Positives = 89/154 (58%), Gaps = 8/154 (5%)

Query  5    EKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRV  64
            E+++   L D +DF NG +     P+S   G  P+YG NG++GYS ++N    LI+ GRV
Sbjct  3    EELKIKFLSDLIDFKNGKS---IKPSS---GDIPIYGGNGILGYSEKYNYNNILII-GRV  55

Query  65   GSYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVL  124
            G+YCGS+ Y     WV+DNA+A   K      Y YY +    LN  + GS QPLL+QGVL
Sbjct  56   GAYCGSIHYHKEKCWVSDNAIAGEVKSDYSIDYLYYLMKSLNLNDRQVGSSQPLLTQGVL  115

Query  125  RNVSVSAV-AAPDRPRIGEILGAFDDKIAANDRV  157
             N+SV    ++  +  I  +L A D KIA N ++
Sbjct  116  NNISVKIYESSQTQQSIAAVLSALDKKIALNKQI  149


>gi|319946656|ref|ZP_08020890.1| putative restriction modification system DNA specificity subunit 
[Streptococcus australis ATCC 700641]
 gi|319746704|gb|EFV98963.1| putative restriction modification system DNA specificity subunit 
[Streptococcus australis ATCC 700641]
Length=386

 Score =  114 bits (286),  Expect = 2e-23, Method: Compositional matrix adjust.
 Identities = 58/157 (37%), Positives = 91/157 (58%), Gaps = 8/157 (5%)

Query  9    KVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYC  68
            K++LGD +DF NG +        +  G  P+YG NG++GY+ + N     IVVGRVG+YC
Sbjct  3    KIKLGDIIDFKNGKS------VKKSDGVIPIYGGNGILGYTDKSNFSHT-IVVGRVGAYC  55

Query  69   GSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVS  128
            GS+   ++  WV+DNA+A   K+ ++  Y YY L    LN  + GS QPL++Q +LR++ 
Sbjct  56   GSIYVEENSCWVSDNAIAGVPKEGQDLTYLYYVLKSLNLNSKQIGSSQPLITQSMLRDMV  115

Query  129  VS-AVAAPDRPRIGEILGAFDDKIAANDRVIEAAEAL  164
            V   +    + RI   +   D KI  N+++ +  EA+
Sbjct  116  VDIEINIEKQKRIANSISIIDQKIQINNQINQELEAM  152


>gi|322386250|ref|ZP_08059882.1| type I restriction system specificity protein [Streptococcus 
cristatus ATCC 51100]
 gi|321269712|gb|EFX52640.1| type I restriction system specificity protein [Streptococcus 
cristatus ATCC 51100]
Length=394

 Score =  113 bits (282),  Expect = 5e-23, Method: Compositional matrix adjust.
 Identities = 73/252 (29%), Positives = 121/252 (49%), Gaps = 50/252 (19%)

Query  9    KVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYC  68
            K++LGD +DF NG +        +  G  P+YG NG++GY+ + N     IVVGRVG+YC
Sbjct  3    KIKLGDVIDFKNGKS------IKKSDGNIPIYGGNGILGYTDKSNFSHT-IVVGRVGAYC  55

Query  69   GSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVS  128
            GS+   ++  WV+DNA+A   K+ ++  Y YY L    LN  + GS QPL++Q +L+++ 
Sbjct  56   GSIHVEENLCWVSDNAIAGIPKEGQDLTYLYYVLKSLNLNSKQIGSSQPLITQSMLKDMV  115

Query  129  VSA-VAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTA  187
            V   +    + RI + +   D KI  N+++ +  EA+                       
Sbjct  116  VDVEIDNEKQKRIAKSILIIDQKIQINNQINQELEAMA----------------------  153

Query  188  CLDAQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSA--KLVVSQPCVLFPKLNPRIPRIW  245
                     T+  Y F  FD   Q     G+  +S+  K+V +      P+L   IP  W
Sbjct  154  --------KTLYDYWFVQFDFPDQ----NGKPYKSSGGKMVYN------PELKREIPEGW  195

Query  246  NITSLPSEMALA  257
             +  L  +++++
Sbjct  196  GVEKLKDKLSVS  207


>gi|217968470|ref|YP_002353704.1| restriction modification system DNA specificity domain-containing 
protein [Thauera sp. MZ1T]
 gi|217505797|gb|ACK52808.1| restriction modification system DNA specificity domain protein 
[Thauera sp. MZ1T]
Length=378

 Score =  112 bits (281),  Expect = 8e-23, Method: Compositional matrix adjust.
 Identities = 65/157 (42%), Positives = 86/157 (55%), Gaps = 11/157 (7%)

Query  10   VRLGDHLDFSNGHTSGHTSPASEPG--GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSY  67
            V LG+ +DF NG        A +PG  G YP YG+NG+IG +         I++GRVG+Y
Sbjct  7    VTLGEVVDFFNGK-------AIKPGQDGEYPAYGSNGLIGGAPDWKYENS-IIIGRVGAY  58

Query  68   CGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNV  127
            CGS+ YC S  W +DN +  R K   +  Y+YY L    LNRY  G+ QPL++Q VL+ V
Sbjct  59   CGSVAYCKSRFWASDNTIVARPKS-GDVGYFYYLLKALELNRYAGGAAQPLVTQTVLKGV  117

Query  128  SVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEAL  164
                   P + RI  IL A+DD I  N R I   E +
Sbjct  118  PARVPDIPTQRRIASILSAYDDLIENNTRRIAILEEM  154


>gi|315612675|ref|ZP_07887587.1| HsdS family protein [Streptococcus sanguinis ATCC 49296]
 gi|315315262|gb|EFU63302.1| HsdS family protein [Streptococcus sanguinis ATCC 49296]
Length=388

 Score =  110 bits (274),  Expect = 5e-22, Method: Compositional matrix adjust.
 Identities = 92/396 (24%), Positives = 162/396 (41%), Gaps = 60/396 (15%)

Query  6    KVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVG  65
            K+E+  LGD ++F  G+    +       G+YPV  +NG++GY  ++   GP I +GR G
Sbjct  3    KLEECILGDLVEFQRGYDLPKSKFVE---GKYPVQSSNGILGYHNEYKVEGPGITIGRSG  59

Query  66   SYCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLR  125
            +  G+      + +  + +L  +  K  +  Y YY L    L   ++GSG P +++  L 
Sbjct  60   T-VGNPHLIRENFFPHNTSLFVKEFKGNDIEYIYYLLQYLDLGNQKSGSGVPTMNRNHLH  118

Query  126  NVSVSAVAAPD-RPRIGEILGAFDDKIAANDRVIEAAEALMLAIV---------------  169
             + + A      + R  +IL   D KI  N+++ +  E +   +                
Sbjct  119  PIKIRAYRDKTCQQRTIKILSLIDKKIQINNQINQELEVMAKTLYDYWFVQFDFPDQNGK  178

Query  170  -----GRLSAYVP--------------LSSLASRSTACLDAQHF-DSTVAHYSFAAFDGG  209
                 G    Y P              LSSL       ++          +YS   +D  
Sbjct  179  PYKSSGGKMVYNPELKREIPVRWGVEKLSSLLEIGRETINPMKTPKEEFKYYSIPEYDVS  238

Query  210  AQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPR-IWNITSLPSEMALASTEFVVLRPVG  268
               S   G TIRS K +V +  +L  KLNP   R ++N+     E A++STEF+V +   
Sbjct  239  GSFSYELGETIRSNKFIVEKSDLLVSKLNPWFNRVVYNL----EENAISSTEFIVWKTFN  294

Query  269  -VDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWVRDVRRLTPGHAAAIAN  327
              + + L+      + +    +   G + S +R+ P  ++   +       P     I  
Sbjct  295  RFEKNFLYQVATSKEFIEYCTRFTTGTSNSHKRVSPDIMVGFQI-------PFEKTYIQK  347

Query  328  LGALCN-------ERRIESARLASCRDALLPLLMSG  356
             G + +       +  +++  L   RD LLP+LM+G
Sbjct  348  FGEITDSIRTQVLQNNVQNQELTQLRDWLLPMLMNG  383


>gi|319777746|ref|YP_004137397.1| restriction modification system DNA specificity domain [Mycoplasma 
fermentans M64]
 gi|318038821|gb|ADV35020.1| Restriction modification system DNA specificity domain [Mycoplasma 
fermentans M64]
Length=372

 Score =  108 bits (270),  Expect = 1e-21, Method: Compositional matrix adjust.
 Identities = 47/137 (35%), Positives = 80/137 (59%), Gaps = 0/137 (0%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G +P +G NG+I Y  + N     I++GRVG+ CG +   +   +VTDNAL  + K+   
Sbjct  22   GNFPAFGGNGIITYVNKWNVDEDAIIIGRVGANCGCVNITNKKSFVTDNALIFKPKEKNM  81

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
             R+++Y LL   LN++  GS QPLL+QG+L N+ ++  +     +I +IL   D++I  N
Sbjct  82   ARFYFYFLLHLNLNKFHIGSSQPLLTQGILGNIKINIPSLNKCQKISKILDNIDNQIERN  141

Query  155  DRVIEAAEALMLAIVGR  171
            + +++  + +  AI  R
Sbjct  142  NSMVQKLQVMGQAIFNR  158


>gi|238810194|dbj|BAH69984.1| hypothetical protein [Mycoplasma fermentans PG18]
Length=172

 Score =  105 bits (263),  Expect = 9e-21, Method: Compositional matrix adjust.
 Identities = 45/135 (34%), Positives = 80/135 (60%), Gaps = 0/135 (0%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G +P +G NG+I Y  + N     I++GRVG+ CG +   +   +VTDNAL  + K+   
Sbjct  27   GNFPAFGGNGIITYVNKWNVDEDAIIIGRVGANCGCVNITNKKSFVTDNALIFKPKEKNM  86

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
             R+++Y LL   LN++  GS QPLL+QG+L N+ ++  +     +I +IL   D++I  N
Sbjct  87   ARFYFYFLLHLNLNKFHIGSSQPLLTQGILGNIKINIPSLNKCQKISKILDNIDNQIERN  146

Query  155  DRVIEAAEALMLAIV  169
            + +++  ++   A++
Sbjct  147  NSMVQKLQSFEQALI  161


>gi|308190349|ref|YP_003923280.1| hypothetical protein MFE_08350 [Mycoplasma fermentans JER]
 gi|307625091|gb|ADN69396.1| hypothetical protein MFE_08350 [Mycoplasma fermentans JER]
Length=167

 Score =  105 bits (261),  Expect = 1e-20, Method: Compositional matrix adjust.
 Identities = 45/135 (34%), Positives = 79/135 (59%), Gaps = 0/135 (0%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G +P +G NG+I Y  + N     I++GRVG+ CG +   +   +VTDNAL  + K+   
Sbjct  22   GNFPAFGGNGIITYVNKWNVDEDAIIIGRVGANCGCVNITNKKSFVTDNALIFKPKEKNM  81

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
             R+++Y LL   LN++  GS QPLL+QG+L N+ ++  +     +I +IL   D++I  N
Sbjct  82   ARFYFYFLLHLNLNKFHIGSSQPLLTQGILGNIKINIPSLNKCQKISKILDNIDNQIERN  141

Query  155  DRVIEAAEALMLAIV  169
            + +++  +    A++
Sbjct  142  NSMVQKLQCFEQALI  156


>gi|332704541|ref|ZP_08424629.1| restriction modification system protein with DNA specificity 
domain [Desulfovibrio africanus str. Walvis Bay]
 gi|332554690|gb|EGJ51734.1| restriction modification system protein with DNA specificity 
domain [Desulfovibrio africanus str. Walvis Bay]
Length=442

 Score =  104 bits (259),  Expect = 3e-20, Method: Compositional matrix adjust.
 Identities = 87/290 (30%), Positives = 126/290 (44%), Gaps = 46/290 (15%)

Query  113  GSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAI----  168
            GS  P  + G LR++        ++ RI  +LGA DDKI  N + + A E L  AI    
Sbjct  145  GSSVPGFNLGQLRSIRFPLPPLSEQSRISRVLGALDDKIEQNQQAVRALERLAQAIFCAW  204

Query  169  ----------VGRLSAY----------------------VP-------LSSLASRSTACL  189
                      V   +++                      VP       +S LA+ S   +
Sbjct  205  FVDFEPIKAKVAGATSFPSMPQPVFDALSIRLIDSKIGPVPEGWKVGTVSDLATLSKTQI  264

Query  190  DAQHF-DSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNIT  248
              Q + D    ++S  AFD G +     G+ I+S K VV + CVL  KLNPRIPRIW   
Sbjct  265  KPQDYPDELFDYFSIPAFDTGKRAFLELGKAIKSNKFVVVEGCVLLSKLNPRIPRIWLPP  324

Query  249  SLPSEMALASTEFVVLRPV-GVDTSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLL  307
                +  + STEF+V  P   +D   L+   +Q      L Q   G + S QR++P  LL
Sbjct  325  PPNGKRQITSTEFLVFVPCSSIDRHYLYCQFQQSSFRENLAQGASGTSSSHQRVRPNDLL  384

Query  308  -RVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
             +  +   + +    A  I  L +  +   +ES +LA  RD LLP L+SG
Sbjct  385  GKAVIVPPKPIRMEFAHLIDPLFSFASACLLESTKLAEMRDYLLPKLLSG  434


>gi|78484678|ref|YP_390603.1| restriction modification system DNA specificity subunit [Thiomicrospira 
crunogena XCL-2]
 gi|78362964|gb|ABB40929.1| Type I restriction enzyme, S subunit [Thiomicrospira crunogena 
XCL-2]
Length=371

 Score =  103 bits (257),  Expect = 4e-20, Method: Compositional matrix adjust.
 Identities = 54/154 (36%), Positives = 80/154 (52%), Gaps = 7/154 (4%)

Query  10   VRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCG  69
            + L   L+  NG    +TS      G Y ++G+NGVIG +  +N    LI+ GRVG+YCG
Sbjct  6    IELSKALNLKNGKALKNTS-----NGIYQIFGSNGVIGTTELNNNENALII-GRVGAYCG  59

Query  70   SLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSV  129
            S+       W +DN +    K      YWYY L  F L ++  G+ QPLL+Q  L+ + +
Sbjct  60   SIELSQEKFWASDNTIVAEPKPDNVLHYWYYRLKSFPLRKFAGGAAQPLLTQNTLKPLKI  119

Query  130  SA-VAAPDRPRIGEILGAFDDKIAANDRVIEAAE  162
            +A     ++ +I  IL  +DD I  N+R I   E
Sbjct  120  AAHTDYLEQDKIANILKVYDDLIENNNRRIALLE  153


>gi|53805024|ref|YP_113333.1| type I restriction-modification system S subunit [Methylococcus 
capsulatus str. Bath]
 gi|53758785|gb|AAU93076.1| type I restriction-modification system, S subunit, EcoA family 
[Methylococcus capsulatus str. Bath]
Length=416

 Score =  103 bits (257),  Expect = 5e-20, Method: Compositional matrix adjust.
 Identities = 106/376 (29%), Positives = 150/376 (40%), Gaps = 59/376 (15%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G  P+  ++GV    A+   +GP +V GR G+  G + Y + D W  +  L  R  K  +
Sbjct  30   GDVPLVSSSGVTDTHAKAMVKGPGVVTGRYGT-LGQVFYVEQDFWPLNTTLYVRDFKGND  88

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPR-IGEILGAFDDKIAA  153
             R+  Y L     + Y   +  P L++  L    V   + P+  R I  ILG  DDKI  
Sbjct  89   PRFISYFLRDVDFHAYSDKAAVPGLNRNHLHQAKVRIPSDPNEQRAIAHILGTLDDKIEL  148

Query  154  NDRVIEAAEALMLAIVGRLSAYVPLSSLASR--------STACLDAQHF-----DSTVAH  200
            N R  E  EA+  A+    + +V    + ++         T     QH      D  V  
Sbjct  149  NRRQNETLEAMARALFK--AWFVDFEPVRAKCRGDRPVAPTGWQWPQHILDLFPDRLVES  206

Query  201  -----------YSFAAF----DGGAQPSRVGG----------------------RTIRSA  223
                       +SF        G   PSR  G                       +I+S 
Sbjct  207  ELGEIPEGWRVFSFGDVAEQGKGFVNPSREPGERFTHYSLPAFDAGKMPVIEPGESIKSN  266

Query  224  KLVVSQPCVLFPKLNPRIPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAAL-RQPD  282
            K  V    VL  KLNP IPRIW +       A+ STEF+V  P     SA    L   P+
Sbjct  267  KTPVPDGAVLVSKLNPHIPRIW-LVGEAGNRAVCSTEFIVWTPKSPAQSAFVYCLASSPE  325

Query  283  VLAELRQLVGGMTGSRQRIQPTQL--LRVWVRDVRRLTPGHAAAIANLGALCNERRIESA  340
             +  + QLV G + S QR++P QL  +RV+  +   +      A   +       R +S 
Sbjct  326  FVGAMCQLVTGTSNSHQRVKPDQLREIRVFAGNENVVETFSKTAEPLMDQFLQNTR-QSR  384

Query  341  RLASCRDALLPLLMSG  356
             LA  RD LLP L+SG
Sbjct  385  ILAQLRDTLLPKLISG  400


>gi|329114036|ref|ZP_08242800.1| Putative type-1 restriction enzyme MjaXP specificity protein 
[Acetobacter pomorum DM001]
 gi|326696575|gb|EGE48252.1| Putative type-1 restriction enzyme MjaXP specificity protein 
[Acetobacter pomorum DM001]
Length=439

 Score =  102 bits (255),  Expect = 8e-20, Method: Compositional matrix adjust.
 Identities = 87/260 (34%), Positives = 112/260 (44%), Gaps = 44/260 (16%)

Query  140  IGEILGAFDDKIAANDRVIEAAEALMLAIV------------------------------  169
            I  ILG+ DDKI  N R  E  EA+  A+                               
Sbjct  166  IASILGSLDDKIDLNRRTNETLEAMARALFRDWFVDFGPTRAKMAGEAPYLAPELWELFP  225

Query  170  GRLS--------AYVPLSSLASRSTACLD-AQHFDSTVAHYSFAAFDGGAQPSRVGGRTI  220
             RL            PL+ LA  S + ++ AQ  D    H+S  AFD G  P  V G  I
Sbjct  226  DRLDDEGNPEGWQSWPLADLAILSKSSINPAQFSDEYFLHFSLPAFDKGMMPDLVKGEEI  285

Query  221  RSAKLVVSQPCVLFPKLNPRIPRIWNITSLPSEMA-LASTEFVVLRPVGVDTSALW--AA  277
            +S K  VS   +L  KLNP  PR+W +T+       + STEF+VL P+  D  AL   A 
Sbjct  286  KSGKFSVSSNSILLSKLNPETPRVWMVTAHEEPYQRICSTEFMVLNPLQKDWLALIYCAC  345

Query  278  LRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWVRDVRRLTPGHAAAIAN-LGALCNERR  336
            L QP     L+ +V G + S QR+QP  +++  +     +        A  L A  N  R
Sbjct  346  LSQP-FRETLQGMVTGTSKSHQRVQPLAVMQTHLLHATDILMRQFDLTAQPLLAKMNFNR  404

Query  337  IESARLASCRDALLPLLMSG  356
             ES  LA  RD LLP LMSG
Sbjct  405  NESNTLAQLRDLLLPKLMSG  424


>gi|91773198|ref|YP_565890.1| restriction modification system DNA specificity subunit [Methanococcoides 
burtonii DSM 6242]
 gi|91712213|gb|ABE52140.1| Restriction modification system DNA specificity subunit [Methanococcoides 
burtonii DSM 6242]
Length=351

 Score = 99.0 bits (245),  Expect = 1e-18, Method: Compositional matrix adjust.
 Identities = 55/164 (34%), Positives = 84/164 (52%), Gaps = 10/164 (6%)

Query  12   LGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSL  71
            LGD   F NG       P     G  P+YG NG++GYS + N  G  +++GRVG+YCGSL
Sbjct  181  LGDVASFHNGKKR----PDDIIEGNIPIYGGNGILGYSDKSNNEGVTVIIGRVGAYCGSL  236

Query  72   RYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSA  131
                + VW++DNAL  +    E + + ++ L    LN    GS  PLL+Q +L+++    
Sbjct  237  YIERNPVWISDNALVAKPINKEHSSFLFFLLKSLQLNEIAEGSSHPLLTQNLLKSIQ---  293

Query  132  VAAPDRPRIGEILGAFD---DKIAANDRVIEAAEALMLAIVGRL  172
            +  P   RI   +   D   +KI  N++ I   E L   ++ +L
Sbjct  294  IILPPEHRIEPFVYQADTWFNKIDKNNKQIRTLEKLRDTLLSKL  337


 Score = 57.0 bits (136),  Expect = 5e-06, Method: Compositional matrix adjust.
 Identities = 44/163 (27%), Positives = 71/163 (44%), Gaps = 11/163 (6%)

Query  6    KVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVG  65
            K +K +LGD L  + G +           G+ PVY + G+ GY  +       +++GR G
Sbjct  4    KWKKCKLGDVLVLNYGKSLPERKRVE---GKIPVYSSAGLTGYHNETLVNSEGLIIGRKG  60

Query  66   SYCGSLRYCDSDVWVTDNALACRAKKPEETRYW----YYALLGFGLNRYRAGSGQPLLSQ  121
            +  G + Y  +  +  D A       PEET+Y+    YY L   GL      S  P L++
Sbjct  61   T-VGKIYYSKTPFFCIDTAYYIL---PEETKYYLNFIYYLLKTIGLEELNEDSAVPRLNR  116

Query  122  GVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEAL  164
                +  +     P++  I  +L + DDKI    R  +  EA+
Sbjct  117  NTAYSQDILLPPLPEQRAIASVLSSLDDKIDLLHRQNKTLEAM  159


>gi|206896558|ref|YP_002247704.1| type I restriction/modification enzyme [Coprothermobacter proteolyticus 
DSM 5265]
 gi|206739175|gb|ACI18253.1| type I restriction/modification enzyme [Coprothermobacter proteolyticus 
DSM 5265]
Length=678

 Score = 97.1 bits (240),  Expect = 5e-18, Method: Compositional matrix adjust.
 Identities = 57/168 (34%), Positives = 92/168 (55%), Gaps = 6/168 (3%)

Query  8    EKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNA--RGPLIVVGRVG  65
            E V LG+ +   NG +  +T  ++   G+YPV G+NG+ GY+  ++    G  IVVGRVG
Sbjct  508  EWVELGELIALRNGISISNTLVSNR--GKYPVCGSNGIYGYTDNNDKLLFGETIVVGRVG  565

Query  66   SYCGSLRYCDSDVWVTDNALACRA--KKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGV  123
            +YCG++ Y D  +WVTDNA+      K   +T+Y YY LL   L +Y   +GQP +SQ +
Sbjct  566  AYCGNVHYYDVPIWVTDNAIVVTVTNKDKLKTKYLYYFLLSKDLGKYANVTGQPYISQSI  625

Query  124  LRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGR  171
            + ++ V       + +I + L    + +    R+ E A+  +  I+ R
Sbjct  626  ISSLKVPLPPIEKQQKIVDFLNVQFETLTNIRRLKENAKQTIKMILDR  673


 Score = 43.1 bits (100),  Expect = 0.077, Method: Compositional matrix adjust.
 Identities = 45/173 (27%), Positives = 72/173 (42%), Gaps = 22/173 (12%)

Query  5    EKVEKVRLGDHLDFSNGH--TSGHTSPASEP---GGRYPVYGANGVIGYSAQHNARGPLI  59
            +K E V L +  D   G   T   T   + P   GG+ P Y       Y  Q N  G +I
Sbjct  336  QKWEVVSLREICDIQKGTSITKADTVEGNVPVIAGGQEPAY-------YHNQSNRDGNII  388

Query  60   VVGRVGSYCGSLRYCDSDVWVTD-NALACRAKKPEETRYWYYALLGFGLNRYRA--GSGQ  116
             V   G+Y G + Y D  ++ +D   +    ++   T+Y +Y L     + Y+   G+GQ
Sbjct  389  TVSASGAYAGFVNYFDIPIFASDCTTIKSNDEEKALTKYIFYILKSRQEDLYKLQRGAGQ  448

Query  117  PLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIV  169
            P +    L N+    +  P  P   E++   D + A    +IE   A+   I+
Sbjct  449  PHVYPNDLANIQ---IPLPPLPVQQELVARLDKQQA----IIEQCNAMEKTIL  494


>gi|71065438|ref|YP_264165.1| type I restriction modification system methylase [Psychrobacter 
arcticus 273-4]
 gi|71038423|gb|AAZ18731.1| probable type I restriction modification system methylase [Psychrobacter 
arcticus 273-4]
Length=424

 Score = 96.3 bits (238),  Expect = 8e-18, Method: Compositional matrix adjust.
 Identities = 56/145 (39%), Positives = 75/145 (52%), Gaps = 6/145 (4%)

Query  12   LGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSL  71
            LGD     NG TS   +      G  PV+G+NG IG   + N R   I++GRVGSYCGSL
Sbjct  253  LGDFAILGNGKTSPDRAV-----GDIPVFGSNGKIGDCDESN-RDNTIIIGRVGSYCGSL  306

Query  72   RYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSA  131
            +Y     W+TDNA++   K  +   Y +  L    LN  R GSGQPLL+Q +LR++    
Sbjct  307  QYYPFKCWITDNAMSAEMKNKDHNIYLFQLLSRDNLNDRRTGSGQPLLNQSILRSIKTIT  366

Query  132  VAAPDRPRIGEILGAFDDKIAANDR  156
             + P       I  +F  KI   +R
Sbjct  367  PSVPLIDEYSRIANSFYKKINKANR  391


 Score = 43.1 bits (100),  Expect = 0.071, Method: Compositional matrix adjust.
 Identities = 36/134 (27%), Positives = 58/134 (44%), Gaps = 20/134 (14%)

Query  48   YSAQHNARGPLI---VVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPE------ETRYW  98
            Y      R P I   V+ RVGSY G + Y +S      N   C  +          +R+ 
Sbjct  57   YEDFSRGRDPQIGDLVISRVGSY-GVVSYVNS------NEKFCLGQNTAFIIPKINSRFL  109

Query  99   YYALLG----FGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
            YY L+     + + ++  G+ Q  +S   +R   +      ++  I  ILG+ DDKI  N
Sbjct  110  YYQLISPFVKWQIEQFVVGAVQKTISLKSIRQFQIKLPPVTEQKAIAHILGSLDDKIELN  169

Query  155  DRVIEAAEALMLAI  168
             ++ E  EA+  A+
Sbjct  170  RQMNETLEAMAQAL  183


>gi|239827072|ref|YP_002949696.1| restriction modification system DNA specificity domain-containing 
protein [Geobacillus sp. WCH70]
 gi|239807365|gb|ACS24430.1| restriction modification system DNA specificity domain protein 
[Geobacillus sp. WCH70]
Length=428

 Score = 95.9 bits (237),  Expect = 9e-18, Method: Compositional matrix adjust.
 Identities = 93/345 (27%), Positives = 142/345 (42%), Gaps = 53/345 (15%)

Query  59   IVVGRVGSYCGSLRYCDSDVWVTDNALACRAK------KPEETRYWYYALLGFG-LNRYR  111
            + + R G+  G     + DV     +   R K       P   +Y + + + +G +N   
Sbjct  71   LCIARTGATTGINTVIEEDVDAVFASYLVRFKLNKEIVDPTFIKYIFKSNMWYGYVNSII  130

Query  112  AGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGR  171
            +GS QP  +   + N  +S     ++ +I  +L   D KI  N+++ +  E +  AI  R
Sbjct  131  SGSAQPGANAQQMSNFKMSIPDLDEQKKIASVLSVLDKKIVLNNKINKTLEEMAQAIFKR  190

Query  172  L-------------------------SAYVP-------LSSLASRSTACLD-AQHFDSTV  198
                                      S  +P       L +L   +TA +D  ++ +   
Sbjct  191  WFVDFEFPNENGKPYKSSGGKFVESESGMIPEGWKEGTLDNLVVINTASVDPKENPEILY  250

Query  199  AHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWNITSLPSEMALAS  258
             HYS  AFD    P    GR I+S K +V     L  KLNP   R+W+   + +E A++S
Sbjct  251  EHYSIPAFDEQKYPKFEYGREIKSNKYLVRPNSFLVSKLNPTTKRVWDPLCI-TENAISS  309

Query  259  TEFVVLRPVGVD-TSALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLL--RVWVRDVR  315
            TEF+   P  +   S L+  L        L +   G TGSRQR++P + L   V + D  
Sbjct  310  TEFINYLPKDISYQSYLYCMLNSERFSEHLIKHATGSTGSRQRVKPAETLTFNVILPDTE  369

Query  316  RLTPGHAAAIANLGALCNER----RIESARLASCRDALLPLLMSG  356
             L         NL     E+    +I SA L   RD LLP LMSG
Sbjct  370  TL-----KKFDNLIRPIREKLKINQINSAVLKDVRDILLPKLMSG  409


>gi|291288456|ref|YP_003505272.1| restriction modification system DNA specificity domain-containing 
protein [Denitrovibrio acetiphilus DSM 12809]
 gi|290885616|gb|ADD69316.1| restriction modification system DNA specificity domain protein 
[Denitrovibrio acetiphilus DSM 12809]
Length=405

 Score = 93.6 bits (231),  Expect = 5e-17, Method: Compositional matrix adjust.
 Identities = 52/161 (33%), Positives = 83/161 (52%), Gaps = 2/161 (1%)

Query  8    EKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSY  67
            ++ +L D  +  NG   G  +      G+  V+GANG I  + +     P+IV+GRVG+Y
Sbjct  5    KEYKLADLANLRNGK--GLNNKFYTDFGKSGVWGANGQIASTNEVLNSDPVIVIGRVGAY  62

Query  68   CGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNV  127
            CGS+   + + WVTDNA+    K   +  + YY L    ++R   GS QPL++Q  +  +
Sbjct  63   CGSIHMAEGNNWVTDNAIQATPKNDTDLNFLYYLLKSLNVSRAATGSAQPLITQSGIGVL  122

Query  128  SVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAI  168
               A +   +  +  IL + DDKI  N ++ E  E +  AI
Sbjct  123  ECKAPSPKIQKEVASILSSLDDKIELNRKMNETLEEMARAI  163


>gi|189499314|ref|YP_001958784.1| restriction modification system DNA specificity domain [Chlorobium 
phaeobacteroides BS1]
 gi|189494755|gb|ACE03303.1| restriction modification system DNA specificity domain [Chlorobium 
phaeobacteroides BS1]
Length=430

 Score = 90.5 bits (223),  Expect = 4e-16, Method: Compositional matrix adjust.
 Identities = 55/148 (38%), Positives = 76/148 (52%), Gaps = 9/148 (6%)

Query  5    EKVEKVRLGDHLDFSNGHTS-GHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGR  63
            E  E V +GD   F++G T    T+P        PVYG NGV+GYSAQ      ++++GR
Sbjct  246  ESWEVVEIGDVFKFTSGKTKPKDTAPEPSVERTVPVYGGNGVLGYSAQSLLNEDVLILGR  305

Query  64   VGSYCGSLRYCDSDVWVTDNALACRAKKPEETRYW---YYALLGFGLNRYRAGSGQPLLS  120
            VG YCG         WVTDNAL  + +K    R +   ++A L   LN+Y    GQPL++
Sbjct  306  VGEYCGCAHLTKPVSWVTDNALYAKEEKRSVNRSYARTHFAHL--NLNQYSNKMGQPLIT  363

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFD  148
            QG++  V       P R    E+  AF+
Sbjct  364  QGIINRVK---FGLPSREEQDELANAFE  388


>gi|284048511|ref|YP_003398850.1| restriction modification system DNA specificity domain protein 
[Acidaminococcus fermentans DSM 20731]
 gi|283952732|gb|ADB47535.1| restriction modification system DNA specificity domain protein 
[Acidaminococcus fermentans DSM 20731]
Length=168

 Score = 90.1 bits (222),  Expect = 5e-16, Method: Compositional matrix adjust.
 Identities = 49/155 (32%), Positives = 85/155 (55%), Gaps = 3/155 (1%)

Query  8    EKVRLGD-HLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGS  66
            E+ +LG+  ++F +G   G  +   E  G YPVYG NG+ GY+ Q+N  G   ++GR G+
Sbjct  15   EQRKLGELTIEFKSG--KGIKANEIEARGTYPVYGGNGLRGYAEQYNHDGEYALIGRQGA  72

Query  67   YCGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRN  126
             CG++ +     + T++A+A +A +  +T++ YY L    L +Y   S QP L+   L  
Sbjct  73   LCGNMTFSMGKAYFTEHAVAVKANEINDTKFLYYILCNMNLGQYSGQSAQPGLAVNKLIA  132

Query  127  VSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAA  161
            +        ++ +I   LGAFD+ I  + R+ ++A
Sbjct  133  LKAFVPGKQEQLKISSYLGAFDNLITLHQRMQKSA  167


>gi|261366731|ref|ZP_05979614.1| putative type I restriction modification system methylase [Subdoligranulum 
variabile DSM 15176]
 gi|282571558|gb|EFB77093.1| putative type I restriction modification system methylase [Subdoligranulum 
variabile DSM 15176]
Length=350

 Score = 88.6 bits (218),  Expect = 1e-15, Method: Compositional matrix adjust.
 Identities = 51/166 (31%), Positives = 81/166 (49%), Gaps = 7/166 (4%)

Query  8    EKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSY  67
            E  +LG+ + F NG +   T       G  PVYG NGV+ Y+  HN     +++GRVG+Y
Sbjct  191  EHSKLGELMTFGNGKSRPKTD------GPIPVYGGNGVLSYTDHHNIENA-VLIGRVGAY  243

Query  68   CGSLRYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNV  127
            CGS+       WV+DNA+  ++K  ++  + Y+ L    L  +  G+GQ LL+Q +L N+
Sbjct  244  CGSVYLEQGICWVSDNAIFAKSKITKDEFFDYFLLKRLNLFNHHVGTGQQLLTQEILNNI  303

Query  128  SVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLS  173
             V                +  + I  N R I   + L   ++ RL+
Sbjct  304  EVPKPVTEQIELFNRKATSIFETIFTNSREIIRLQELSDLLLSRLA  349


>gi|242280198|ref|YP_002992327.1| restriction modification system DNA specificity domain protein 
[Desulfovibrio salexigens DSM 2638]
 gi|242123092|gb|ACS80788.1| restriction modification system DNA specificity domain protein 
[Desulfovibrio salexigens DSM 2638]
Length=415

 Score = 88.6 bits (218),  Expect = 2e-15, Method: Compositional matrix adjust.
 Identities = 47/134 (36%), Positives = 66/134 (50%), Gaps = 0/134 (0%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G  P+ G++GV GY  +    GP I VGR G+  G + Y D D W  + AL  +  +   
Sbjct  32   GDIPILGSSGVTGYHNESKVAGPGITVGRSGASIGVVTYSDIDFWPLNTALYVKDFRGNH  91

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAAN  154
             R+ YY L  F   RY +GS QP L++  +    +       +  I  ILG  D+KI  N
Sbjct  92   PRFAYYFLKQFDFKRYNSGSAQPSLNRNFVHPTKIRIPPLKTQQAIAHILGTIDEKIDLN  151

Query  155  DRVIEAAEALMLAI  168
             R+ E  EA+  AI
Sbjct  152  RRMNETLEAMAQAI  165


>gi|281421790|ref|ZP_06252789.1| type I restriction-modification enzyme S subunit [Prevotella 
copri DSM 18205]
 gi|281404148|gb|EFB34828.1| type I restriction-modification enzyme S subunit [Prevotella 
copri DSM 18205]
Length=450

 Score = 87.8 bits (216),  Expect = 3e-15, Method: Compositional matrix adjust.
 Identities = 39/107 (37%), Positives = 63/107 (59%), Gaps = 1/107 (0%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G+YP+YG NG+ GY+  +N  G  I++GRVG YCGS+ Y ++ +WVTDNA   +      
Sbjct  324  GKYPIYGGNGITGYNESYNVDGINIIIGRVGFYCGSVHYVNNKIWVTDNAFVTKIMGNVY  383

Query  95   T-RYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRI  140
            T ++ YY L  + L +Y   + QP++S   +  ++V      ++ RI
Sbjct  384  TPKFLYYLLQQYDLQQYSNSTAQPVISGKTVYPINVMLPPLSEQYRI  430


>gi|150006174|ref|YP_001300918.1| type I restriction endonuclease S subunit [Bacteroides vulgatus 
ATCC 8482]
 gi|149934598|gb|ABR41296.1| type I restriction endonuclease S subunit [Bacteroides vulgatus 
ATCC 8482]
Length=358

 Score = 87.0 bits (214),  Expect = 4e-15, Method: Compositional matrix adjust.
 Identities = 61/197 (31%), Positives = 97/197 (50%), Gaps = 15/197 (7%)

Query  24   SGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDN  83
            SG   P S   G   +YG+ G+IG +   +    +++V RVG+  G L+       VTDN
Sbjct  39   SGKDKPTS--NGTVALYGSTGIIGMTRLASYNEEIVLVARVGANAGQLQITTIPCGVTDN  96

Query  84   ALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEI  143
             L   AK  E  RY YY L  + LNR   GSGQPL++  +L+ + +      +R +I  +
Sbjct  97   TLIINAK--EWNRYIYYYLQHYNLNRLVFGSGQPLITGSMLKKLKIIYGEESERNKIVNL  154

Query  144  LGAFDDKIAANDRVIEAAEALMLAIVGRLSAYVPLSSLASRSTAC-------LDAQHFDS  196
            L   D++IA  +++IE  + L  AI  RL   V  S++   S  C       ++ ++   
Sbjct  155  LCLLDERIATQNKIIEDLKKLKSAISERLFKSVKGSTVL-LSDLCDIVKGKQINGENLSD  213

Query  197  TVAHYSFAAFDGGAQPS  213
            +  +Y     +GG +PS
Sbjct  214  SGNYY---VMNGGTEPS  227


>gi|312875033|ref|ZP_07735051.1| type I restriction modification DNA specificity domain protein 
[Lactobacillus iners LEAF 2053A-b]
 gi|311089428|gb|EFQ47854.1| type I restriction modification DNA specificity domain protein 
[Lactobacillus iners LEAF 2053A-b]
Length=146

 Score = 86.3 bits (212),  Expect = 7e-15, Method: Compositional matrix adjust.
 Identities = 52/146 (36%), Positives = 77/146 (53%), Gaps = 7/146 (4%)

Query  12   LGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSL  71
            L  +++F NG       P  +  G  PVYG NG++ Y+   N +   +V+GRVG YCGS+
Sbjct  4    LSKYVEFKNGKKR----PTLK--GTIPVYGGNGILDYTNTANMQSG-VVIGRVGVYCGSV  56

Query  72   RYCDSDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSA  131
                 + WV+DNA+    K+  +  Y YY L    LN  R G+ QPLL+Q +L N+ V  
Sbjct  57   FLVREECWVSDNAIKAMCKENIDLGYLYYLLSSLHLNERRIGTSQPLLTQNILNNIEVEI  116

Query  132  VAAPDRPRIGEILGAFDDKIAANDRV  157
                 + +I  IL   D+KI  N+ +
Sbjct  117  PELAIQKKISSILELLDEKIKLNNEI  142


>gi|336115613|ref|YP_004570380.1| putative type I restriction modification system, specificity 
protein [Bacillus coagulans 2-6]
 gi|335369043|gb|AEH54994.1| Putative Type I restriction modification system, specificity 
protein [Bacillus coagulans 2-6]
Length=436

 Score = 85.9 bits (211),  Expect = 1e-14, Method: Compositional matrix adjust.
 Identities = 75/277 (28%), Positives = 111/277 (41%), Gaps = 54/277 (19%)

Query  124  LRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRLSA---------  174
            L+ + V      ++  I  IL   D+KI  N+++ E  E +  A+               
Sbjct  155  LKQIKVRLPELNEQKAIANILSTLDEKIETNNQINEKLEEMAQALFKHWFVDFEFPNENG  214

Query  175  ----------------YVP-------LSSLASRSTACLDAQHF-DSTVAHYSFAAFDGGA  210
                             +P       L  L   S   +  Q + +    HYS  AFD G 
Sbjct  215  KPYKSSGGEMVESELGMIPKGWEVGKLGELVEISNKSIKPQEYPEKFFEHYSIPAFDNGR  274

Query  211  QPSRVGGRTIRSAKLVVSQPCVLFPKLNPRIPRIWN-ITSLPSEMALASTEFVVLRPVGV  269
             P    G  I+S+K V++   VL  KLNP   RIW  IT   +  A+ STEF+V      
Sbjct  275  TPEIQKGEEIKSSKYVINNHMVLVSKLNPSTKRIWKPITQ--TNNAICSTEFIVYMAKKE  332

Query  270  DT-SALWAALRQPDVLAELRQLVGGMTGSRQRIQPTQLLRVWV-----RDVRRLT----P  319
            +  S ++  L        L     G TGSRQR++P + L   V     + +R+ +    P
Sbjct  333  NILSYVFEFLNNEKFQGILTANATGSTGSRQRVKPKETLNYKVVLPPLKLMRQFSKIIEP  392

Query  320  GHAAAIANLGALCNERRIESARLASCRDALLPLLMSG  356
             H    +N+        IE+ +L   RD LLP L+SG
Sbjct  393  MHLKIGSNM--------IENQKLKGLRDILLPKLISG  421


>gi|268609384|ref|ZP_06143111.1| putative type I restriction-modification system specificity subunit 
[Ruminococcus flavefaciens FD-1]
Length=367

 Score = 85.5 bits (210),  Expect = 1e-14, Method: Compositional matrix adjust.
 Identities = 46/153 (31%), Positives = 81/153 (53%), Gaps = 10/153 (6%)

Query  36   RYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEET  95
             Y VYG NG+ GY++ +N  G   ++GR G+YCG++RY     ++T++A+   A     T
Sbjct  33   EYAVYGGNGIRGYTSDYNFEGDCAIIGRQGAYCGNVRYFSGKAYMTEHAVIACANSEHNT  92

Query  96   RYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPD---RPRIGEILGAFDDKI-  151
            RY  Y L    L R    S QP +S   ++ +S+  V  P    + +I  ++ + ++KI 
Sbjct  93   RYLSYVLTAMDLGRLSGQSAQPGIS---VKTLSIQKVKMPSLNLQRKIVAVISSLEEKIE  149

Query  152  ---AANDRVIEAAEALMLAIVGRLSAYVPLSSL  181
               A N+ + + A+AL   ++  +   VP +S+
Sbjct  150  LNAAINENLEQQAQALFKDMISDVQEQVPFTSV  182


>gi|310826741|ref|YP_003959098.1| hypothetical protein ELI_1147 [Eubacterium limosum KIST612]
 gi|308738475|gb|ADO36135.1| hypothetical protein ELI_1147 [Eubacterium limosum KIST612]
Length=415

 Score = 84.7 bits (208),  Expect = 2e-14, Method: Compositional matrix adjust.
 Identities = 47/142 (34%), Positives = 75/142 (53%), Gaps = 4/142 (2%)

Query  28   SPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALAC  87
            S   +  G YPVYG NG+ GY +++N  G  ++VGR G+ CG++       W T++A+  
Sbjct  258  SDQIDEAGEYPVYGGNGIRGYYSKYNYEGEYLLVGRQGALCGNVHKIKGCFWATEHAVVT  317

Query  88   RAKKPEETRYWYYALLGFGLNRYRAGS-GQPLLSQGVLRNVSVSAVAAPDRPRIGEILGA  146
            +  +  E  + YY L G  LNRY + S  QP LS   ++N  +  V  P   +I EI   
Sbjct  318  KNVEGVELSFLYYLLNGMNLNRYASNSAAQPGLSVNTIQN--IKTVFPPIEEQI-EISTY  374

Query  147  FDDKIAANDRVIEAAEALMLAI  168
             +D   + D +IE  ++L+  +
Sbjct  375  LNDICHSIDSIIENKQSLIFEL  396


>gi|338767478|gb|EGP22297.1| Type I restriction-modification system methyltransferase subunit 
[Escherichia coli PCN033]
Length=635

 Score = 84.3 bits (207),  Expect = 3e-14, Method: Compositional matrix adjust.
 Identities = 45/129 (35%), Positives = 69/129 (54%), Gaps = 2/129 (1%)

Query  33   PGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKP  92
            P G  PVYG NGVIGY+ +    G  I++G+VG+ CG++RY   D+WVT+N+L  +   P
Sbjct  487  PTGPIPVYGGNGVIGYNNEAMENGENIIIGKVGALCGNVRYVQGDIWVTNNSLIIKNYSP  546

Query  93   EETRYWYYALLGFGLN--RYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDK  150
             +    Y A L   LN  R   G+ Q  L+   ++NV VS      + ++   L   D+ 
Sbjct  547  NKVLTPYLAKLLSTLNLRRLAVGTAQQYLTTTQIKNVEVSIPPLTTQYKLNIWLDELDNT  606

Query  151  IAANDRVIE  159
            +   + +IE
Sbjct  607  LEQYNNLIE  615


>gi|229002234|dbj|BAH57700.1| hypothetical protein [Staphylococcus aureus]
 gi|238768520|dbj|BAH66832.1| type I restriction-modification system endonuclease [Staphylococcus 
aureus]
 gi|343098331|dbj|BAK57490.1| type I restriction-modification system endonuclease [Staphylococcus 
aureus]
Length=433

 Score = 82.8 bits (203),  Expect = 8e-14, Method: Compositional matrix adjust.
 Identities = 55/185 (30%), Positives = 83/185 (45%), Gaps = 10/185 (5%)

Query  178  LSSLASRSTACLDAQHFDS-TVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFPK  236
            L  +AS      + +  +  TV H+S  A+D   Q        I+S K +++  CVLF K
Sbjct  236  LKDIASHKKETFNPKKSEEVTVKHFSLPAYDNEEQAIEEEVNKIKSNKWIINNNCVLFSK  295

Query  237  LNPRIPRIWNITSLP----SEMALASTEFVVLR-PVGVDTSALWAALRQPDVLAELRQLV  291
            +NP   RIW    LP     ++ +AS+EFVV+  P     S ++        +  L+   
Sbjct  296  MNPDTKRIW----LPVIDNKKLNVASSEFVVMESPNNKINSFIYNICLNSQFIDYLKANT  351

Query  292  GGMTGSRQRIQPTQLLRVWVRDVRRLTPGHAAAIANLGALCNERRIESARLASCRDALLP  351
             G T SRQR++PT  +   +     +   ++  I          R E  +L   RD LLP
Sbjct  352  TGSTNSRQRVKPTIAVNYKLAIEDSIVKKYSEIITPYMEEMKILRSEIGKLTQLRDTLLP  411

Query  352  LLMSG  356
             LMSG
Sbjct  412  KLMSG  416


>gi|300821374|ref|ZP_07101522.1| type I restriction modification DNA specificity domain protein 
[Escherichia coli MS 119-7]
 gi|331680404|ref|ZP_08381063.1| type I restriction enzyme EcoR124II specificity protein (S protein)(S.EcoR124II) 
[Escherichia coli H591]
 gi|300526263|gb|EFK47332.1| type I restriction modification DNA specificity domain protein 
[Escherichia coli MS 119-7]
 gi|331071867|gb|EGI43203.1| type I restriction enzyme EcoR124II specificity protein (S protein)(S.EcoR124II) 
[Escherichia coli H591]
Length=388

 Score = 82.0 bits (201),  Expect = 2e-13, Method: Compositional matrix adjust.
 Identities = 51/145 (36%), Positives = 76/145 (53%), Gaps = 6/145 (4%)

Query  5    EKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRV  64
            E   ++R+GD  DF  G  +G+T P S  GG+YPVYG+NG++G   ++N+     V+G +
Sbjct  210  ENTREMRIGDIYDFQYG--TGNTIPKS--GGQYPVYGSNGIVGSHDKYNSEDS-PVIGHI  264

Query  65   GSYCGSLRYCDSDVWVTDNALACRAKKPEE-TRYWYYALLGFGLNRYRAGSGQPLLSQGV  123
            G+Y G + +     +VT N + CR K  E   +Y YY LL          + QP +S  +
Sbjct  265  GAYAGIVNWGQGKHFVTYNGVICRHKSKEVLQKYAYYLLLLQDFGSKSNSASQPFVSYNI  324

Query  124  LRNVSVSAVAAPDRPRIGEILGAFD  148
            L    V      ++ RI EIL  FD
Sbjct  325  LNAPIVLVPPLQEQARIVEILDKFD  349


>gi|89093019|ref|ZP_01165970.1| putative type I restriction enzyme, S subunit [Oceanospirillum 
sp. MED92]
 gi|89082669|gb|EAR61890.1| putative type I restriction enzyme, S subunit [Oceanospirillum 
sp. MED92]
Length=394

 Score = 81.6 bits (200),  Expect = 2e-13, Method: Compositional matrix adjust.
 Identities = 41/147 (28%), Positives = 73/147 (50%), Gaps = 4/147 (2%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G YPVYG NG+ GY++++N  G   ++GR G+ CG++ + +   + T++A+A +A +  +
Sbjct  46   GAYPVYGGNGLRGYTSEYNHEGSYALIGRQGALCGNMNFSNGKAFFTEHAIAVQANEKND  105

Query  95   TRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKI---  151
            T + YY L    L +Y   S QP LS   L  +        ++  IG      D  I   
Sbjct  106  TLFLYYKLGSMNLGQYSGQSAQPGLSVNKLSELETFTAGKVEQTAIGNYFHKLDTLINQH  165

Query  152  -AANDRVIEAAEALMLAIVGRLSAYVP  177
               +D++    +A++  +  +    VP
Sbjct  166  QQKHDKLSNLKKAMLEKMFPKAGETVP  192


>gi|158315561|ref|YP_001508069.1| restriction modification system DNA specificity subunit [Frankia 
sp. EAN1pec]
 gi|158110966|gb|ABW13163.1| restriction modification system DNA specificity domain [Frankia 
sp. EAN1pec]
Length=374

 Score = 81.3 bits (199),  Expect = 2e-13, Method: Compositional matrix adjust.
 Identities = 50/144 (35%), Positives = 72/144 (50%), Gaps = 1/144 (0%)

Query  25   GHTSPASEP-GGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDN  83
            G   PA E   G +PV G+ GV G+  +    GP I +GR GS  G++ Y  SD W  + 
Sbjct  18   GFDLPAPERRAGCFPVVGSAGVSGWHDRGPIAGPGITLGRSGSSIGTVTYVPSDYWPLNT  77

Query  84   ALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEI  143
             L     +  + R+ Y+ L      R+ +GS QP L++  +  V + A   P++  I  +
Sbjct  78   VLFVEDFQGNDPRFLYFLLRTIDFARFNSGSAQPSLNRNYIAAVELRAPEYPEQRAIAAV  137

Query  144  LGAFDDKIAANDRVIEAAEALMLA  167
            LGA DDKIA N R+   A  L  A
Sbjct  138  LGALDDKIALNHRLASTARELAEA  161


>gi|85711748|ref|ZP_01042804.1| type I restriction-modification system, S subunit, EcoA family 
protein [Idiomarina baltica OS145]
 gi|85694363|gb|EAQ32305.1| type I restriction-modification system, S subunit, EcoA family 
protein [Idiomarina baltica OS145]
Length=663

 Score = 80.9 bits (198),  Expect = 3e-13, Method: Compositional matrix adjust.
 Identities = 63/187 (34%), Positives = 90/187 (49%), Gaps = 11/187 (5%)

Query  177  PLSSLASRSTACLD-AQHFDSTVAHYSFAAFDGGAQPSRVGGRTIRSAKLVVSQPCVLFP  235
            P S++A   T  +  A+  +    HYS  AFD G  P+   G  I+S K  V    +L  
Sbjct  277  PFSNIARLDTTSVKPAKEPEKIWEHYSIPAFDDGMSPAFDLGVDIKSNKYRVFPASILVS  336

Query  236  KLNPRIPRIWNITSLPSEMALASTEFVVLRPVGVDTSALWAALRQPDVLAE-LRQLVGGM  294
            KLNP  PR W      S  A+ STEF+   P+  +  A  A + + +     +   V G 
Sbjct  337  KLNPHFPRTWLPDVFDSNAAICSTEFMQFVPIKPNQRAFVAGVVKSESFQNGIMMRVTGS  396

Query  295  TGSRQRIQPTQLLRVWV----RDVRRLTPGHAAAIA-NLGALCNERRIESARLASCRDAL  349
            TGSRQR QP Q+  + V     ++R +   ++  IA  L +  +  R E+  LA  RD L
Sbjct  397  TGSRQRAQPKQVAEMEVLLPSEELRNI---YSVLIAPQLESQASNIR-EALNLADVRDTL  452

Query  350  LPLLMSG  356
            LP L+SG
Sbjct  453  LPKLLSG  459


>gi|229195092|ref|ZP_04321867.1| type I restriction-modification enzyme, S subunit [Bacillus cereus 
m1293]
 gi|228588321|gb|EEK46364.1| type I restriction-modification enzyme, S subunit [Bacillus cereus 
m1293]
Length=475

 Score = 80.9 bits (198),  Expect = 3e-13, Method: Compositional matrix adjust.
 Identities = 51/162 (32%), Positives = 75/162 (47%), Gaps = 11/162 (6%)

Query  2    SRVEKVEKVRLGDHLDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVV  61
            +R + V K++ GD L          T       G  PVYG NG+ G   + N     IV+
Sbjct  33   TRFDSVLKIKSGDSL----------TKAKMNEQGMIPVYGGNGITGTHDKSNVETETIVI  82

Query  62   GRVGSYCGSLRYCDSDVWVTDNALACR-AKKPEETRYWYYALLGFGLNRYRAGSGQPLLS  120
            GRVG YCGS+     + WVTDNA      +K  + ++ Y+ L    L +Y   S QP++S
Sbjct  83   GRVGYYCGSVHLTSEEAWVTDNAFVLSFPEKIIDKKFIYWNLKHCNLGQYSKSSAQPVIS  142

Query  121  QGVLRNVSVSAVAAPDRPRIGEILGAFDDKIAANDRVIEAAE  162
               +  V ++     ++ RI E +     KI     +IE AE
Sbjct  143  GKTIGPVGINVPPYLEQKRIVEKVERLLGKIEEAKALIEEAE  184


>gi|38505922|ref|NP_942540.1| type I site-specific deoxyribonuclease [Synechocystis sp. PCC 
6803]
 gi|38423946|dbj|BAD02154.1| type I site-specific deoxyribonuclease [Synechocystis sp. PCC 
6803]
Length=394

 Score = 80.5 bits (197),  Expect = 4e-13, Method: Compositional matrix adjust.
 Identities = 45/157 (29%), Positives = 77/157 (50%), Gaps = 6/157 (3%)

Query  16   LDFSNGHTSGHTSPASEPGGRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCD  75
            L++ N  T  + S     GG YPVYG+NG+IG+   +    P I+VGR GS  G + + +
Sbjct  16   LEYGNTLTEENRS-----GGDYPVYGSNGIIGFHKAYLLDSPNIIVGRKGS-VGEVVWAN  69

Query  76   SDVWVTDNALACRAKKPEETRYWYYALLGFGLNRYRAGSGQPLLSQGVLRNVSVSAVAAP  135
             + W  D       K+    R+ Y+ L  F L +  + +G P L++  +  +  +  + P
Sbjct  70   KNCWAIDTTYYVTLKQDNSLRFIYWLLKSFDLRKLDSSTGVPGLNRNDVYRIKCNLPSLP  129

Query  136  DRPRIGEILGAFDDKIAANDRVIEAAEALMLAIVGRL  172
            ++ +I EIL   D+ IA  +  I   + +   +V  L
Sbjct  130  EQEKIAEILDTMDEAIAKTEECIAKLKKIKAGLVHDL  166


>gi|345342870|gb|EGW75262.1| type I restriction modification DNA specificity domain protein 
[Escherichia coli STEC_C165-02]
Length=633

 Score = 80.1 bits (196),  Expect = 5e-13, Method: Compositional matrix adjust.
 Identities = 41/127 (33%), Positives = 69/127 (55%), Gaps = 2/127 (1%)

Query  35   GRYPVYGANGVIGYSAQHNARGPLIVVGRVGSYCGSLRYCDSDVWVTDNALACRAKKPEE  94
            G  PVYG NGVIGY+ +    G  I++G+VG+ CG++RY   D+WVT+N+L  +   P +
Sbjct  487  GPIPVYGGNGVIGYNYEAMENGENIIIGKVGALCGNVRYVQGDIWVTNNSLIIKNYLPNK  546

Query  95   TRYWYYALLGFGLN--RYRAGSGQPLLSQGVLRNVSVSAVAAPDRPRIGEILGAFDDKIA  152
                Y A L   LN  +   G+ Q  L+   ++NV ++      + ++   L   D+ + 
Sbjct  547  VLTPYLAKLLSTLNLRKLAVGTAQQYLTTTQIKNVEINIPPLTTQHKLNMWLDELDNTLE  606

Query  153  ANDRVIE  159
              +++IE
Sbjct  607  QYNKLIE  613



Lambda     K      H
   0.321    0.136    0.409 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 

Effective search space used: 685166507650


  Database: All non-redundant GenBank CDS translations+PDB+SwissProt+PIR+PRF
excluding environmental samples from WGS projects
    Posted date:  Sep 5, 2011  4:36 AM
  Number of letters in database: 5,219,829,388
  Number of sequences in database:  15,229,318



Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Neighboring words threshold: 11
Window for multiple hits: 40