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= Rv2755c

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

gi|15842293|ref|NP_337330.1|  HsdS-related protein [Mycobacterium...   184    3e-45
gi|28868302|ref|NP_790921.1|  type I restriction-modification sys...   113    9e-24
gi|158315561|ref|YP_001508069.1|  restriction modification system...   100    8e-20
gi|294339299|emb|CAZ87655.1|  Putative Restriction modification s...  99.8    1e-19
gi|242280198|ref|YP_002992327.1|  restriction modification system...  99.0    2e-19
gi|220908458|ref|YP_002483769.1|  restriction modification system...  97.8    5e-19
gi|198284497|ref|YP_002220818.1|  restriction modification system...  97.1    8e-19
gi|291296825|ref|YP_003508223.1|  restriction modification system...  95.9    2e-18
gi|339834363|gb|EGQ62133.1|  type I restriction-modification syst...  95.9    2e-18
gi|2408222|gb|AAB70708.1|  HsdS [Klebsiella pneumoniae]               95.1    3e-18
gi|53805024|ref|YP_113333.1|  type I restriction-modification sys...  94.7    4e-18
gi|121583503|ref|YP_973929.1|  restriction modification system DN...  94.0    6e-18
gi|331650404|ref|ZP_08351476.1|  type I restriction-modification ...  93.6    1e-17
gi|49658898|emb|CAF28524.1|  putative HsdS-like DNA methylase [Ye...  92.8    2e-17
gi|322691669|ref|YP_004221239.1|  hypothetical protein BLLJ_1480 ...  91.3    4e-17
gi|86146743|ref|ZP_01065063.1|  type I restriction-modification s...  90.9    6e-17
gi|344198531|ref|YP_004782857.1|  restriction modification system...  89.4    2e-16
gi|288937352|ref|YP_003441411.1|  restriction modification system...  89.4    2e-16
gi|253687261|ref|YP_003016451.1|  restriction modification system...  88.6    3e-16
gi|219851731|ref|YP_002466163.1|  restriction modification system...  86.7    9e-16
gi|253563523|ref|ZP_04840980.1|  restriction modification system ...  85.5    2e-15
gi|32455520|ref|NP_862272.1|  HsdS' [Lactobacillus sakei] >gi|244...  84.7    4e-15
gi|301162175|emb|CBW21720.1|  putative type I DNA restriction-mod...  84.7    4e-15
gi|160945577|ref|ZP_02092803.1|  hypothetical protein FAEPRAM212_...  84.0    6e-15
gi|126726006|ref|ZP_01741848.1|  putative type I restriction-modi...  81.6    3e-14
gi|242278888|ref|YP_002991017.1|  restriction modification system...  81.6    4e-14
gi|194364815|ref|YP_002027425.1|  restriction modification system...  81.3    4e-14
gi|260767611|ref|ZP_05876547.1|  type I restriction-modification ...  81.3    5e-14
gi|159904437|ref|YP_001548099.1|  restriction modification system...  81.3    5e-14
gi|344340139|ref|ZP_08771065.1|  restriction modification system ...  80.9    5e-14
gi|119357508|ref|YP_912152.1|  restriction modification system DN...  80.5    9e-14
gi|21226532|ref|NP_632454.1|  type I restriction-modification sys...  80.1    9e-14
gi|343920509|gb|EGV31240.1|  restriction modification system DNA ...  79.7    1e-13
gi|284108343|ref|ZP_06386407.1|  Restriction modification system ...  78.2    4e-13
gi|257051191|ref|YP_003129024.1|  restriction modification system...  78.2    4e-13
gi|329913308|ref|ZP_08275914.1|  Type I restriction-modification ...  78.2    4e-13
gi|313140399|ref|ZP_07802592.1|  type I restriction-modification ...  77.4    6e-13
gi|148262630|ref|YP_001229336.1|  restriction modification system...  75.9    2e-12
gi|315648621|ref|ZP_07901718.1|  restriction modification system ...  75.9    2e-12
gi|315612675|ref|ZP_07887587.1|  HsdS family protein [Streptococc...  74.3    6e-12
gi|323965458|gb|EGB60913.1|  type I restriction modification DNA ...  73.9    8e-12
gi|75909705|ref|YP_324001.1|  restriction modification system DNA...  73.2    1e-11
gi|295136493|ref|YP_003587169.1|  restriction modification system...  71.6    3e-11
gi|339730590|emb|CCC41940.1|  type I site-specific deoxyribonucle...  71.6    4e-11
gi|334123285|ref|ZP_08497313.1|  type I restriction system specif...  71.2    4e-11
gi|262369880|ref|ZP_06063207.1|  sty SBLI [Acinetobacter johnsoni...  71.2    4e-11
gi|227546694|ref|ZP_03976743.1|  type I restriction enzyme HindVI...  70.5    8e-11
gi|301055838|ref|YP_003794049.1|  type I restriction modification...  70.1    9e-11
gi|328947972|ref|YP_004365309.1|  restriction modification system...  69.7    1e-10
gi|261867039|ref|YP_003254961.1|  restriction modification system...  69.7    1e-10


>gi|15842293|ref|NP_337330.1| HsdS-related protein [Mycobacterium tuberculosis CDC1551]
 gi|31793929|ref|NP_856422.1| type I restriction/modification system specificity determinant 
HsdS [Mycobacterium bovis AF2122/97]
 gi|57117020|ref|YP_177904.1| type I restriction/modification system specificity determinant 
HsdS [Mycobacterium tuberculosis H37Rv]
 64 more sequence titles
 Length=91

 Score =  184 bits (468),  Expect = 3e-45, Method: Compositional matrix adjust.
 Identities = 91/91 (100%), Positives = 91/91 (100%), Gaps = 0/91 (0%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG
Sbjct  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
           AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR
Sbjct  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91


>gi|28868302|ref|NP_790921.1| type I restriction-modification system, S subunit [Pseudomonas 
syringae pv. tomato str. DC3000]
 gi|28851539|gb|AAO54616.1| type I restriction-modification system, S subunit [Pseudomonas 
syringae pv. tomato str. DC3000]
Length=422

 Score =  113 bits (283),  Expect = 9e-24, Method: Compositional matrix adjust.
 Identities = 49/91 (54%), Positives = 64/91 (71%), Gaps = 0/91 (0%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           MS  W+ + FG+ + +QRGHDLP  +R  G+VP++GSFG+TG HDT    GPGV IGRSG
Sbjct  1   MSSEWREITFGDFVAIQRGHDLPDQNRKLGSVPILGSFGITGYHDTVKAKGPGVTIGRSG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
           A+ G A +     WPL+T L+V DFKGN P+
Sbjct  61  ASFGVAAYTDQDYWPLNTALYVTDFKGNHPK  91


>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 =  100 bits (249),  Expect = 8e-20, Method: Compositional matrix adjust.
 Identities = 45/87 (52%), Positives = 60/87 (69%), Gaps = 0/87 (0%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           W+     +++ L+RG DLPA  R +G  PV+GS GV+G HD     GPG+ +GRSG++IG
Sbjct  4   WRRSSLADLVRLRRGFDLPAPERRAGCFPVVGSAGVSGWHDRGPIAGPGITLGRSGSSIG  63

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
           T T+V    WPL+T LFV DF+GNDPR
Sbjct  64  TVTYVPSDYWPLNTVLFVEDFQGNDPR  90


>gi|294339299|emb|CAZ87655.1| Putative Restriction modification system protein [Thiomonas sp. 
3As]
Length=407

 Score = 99.8 bits (247),  Expect = 1e-19, Method: Compositional matrix adjust.
 Identities = 47/91 (52%), Positives = 57/91 (63%), Gaps = 0/91 (0%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   W+    GE + LQRGHDL +  +  G VPV+GS G  G HD A   GPGV IGRSG
Sbjct  1   MRSDWRQSNLGEFVRLQRGHDLTSLEQRPGNVPVMGSAGPNGTHDVARATGPGVVIGRSG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
           A+IG   F +   WP +TCL+V DF GN+PR
Sbjct  61  ASIGRVHFSSSDYWPHNTCLYVTDFCGNNPR  91


>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 = 99.0 bits (245),  Expect = 2e-19, Method: Composition-based stats.
 Identities = 42/87 (49%), Positives = 59/87 (68%), Gaps = 0/87 (0%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           W+   F  ++ LQRG DLP + R  G +P++GS GVTG H+ +   GPG+ +GRSGA+IG
Sbjct  7   WRDTIFENLVLLQRGIDLPQSQRRVGDIPILGSSGVTGYHNESKVAGPGITVGRSGASIG  66

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
             T+     WPL+T L+V+DF+GN PR
Sbjct  67  VVTYSDIDFWPLNTALYVKDFRGNHPR  93


>gi|220908458|ref|YP_002483769.1| restriction modification system DNA specificity domain-containing 
protein [Cyanothece sp. PCC 7425]
 gi|219865069|gb|ACL45408.1| restriction modification system DNA specificity domain protein 
[Cyanothece sp. PCC 7425]
Length=388

 Score = 97.8 bits (242),  Expect = 5e-19, Method: Composition-based stats.
 Identities = 48/90 (54%), Positives = 58/90 (65%), Gaps = 1/90 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            + + W  L  GEV+ LQRG DLP  SR  G +P+I S GVT  H++A+  GPGV  GR G
Sbjct  18   IPENWADLLLGEVITLQRGFDLPNRSRRKGDIPIISSSGVTDTHNSASALGPGVITGRYG  77

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
              IG   FV G  WPL+T LFV +FKGNDP
Sbjct  78   -TIGEVFFVEGDYWPLNTTLFVSNFKGNDP  106


>gi|198284497|ref|YP_002220818.1| restriction modification system DNA specificity protein [Acidithiobacillus 
ferrooxidans ATCC 53993]
 gi|218667678|ref|YP_002427161.1| type I restriction-modification system, S subunit [Acidithiobacillus 
ferrooxidans ATCC 23270]
 gi|198249018|gb|ACH84611.1| restriction modification system DNA specificity domain [Acidithiobacillus 
ferrooxidans ATCC 53993]
 gi|218519891|gb|ACK80477.1| type I restriction-modification system, S subunit [Acidithiobacillus 
ferrooxidans ATCC 23270]
Length=418

 Score = 97.1 bits (240),  Expect = 8e-19, Method: Composition-based stats.
 Identities = 46/91 (51%), Positives = 58/91 (64%), Gaps = 1/91 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M++ WK    G+V+EL+RG+DLP   R SG VP++ S GVT  H  A   GPGV  GR G
Sbjct  1   MANEWKECSLGDVIELKRGYDLPQKERLSGDVPLVSSSGVTDTHAKAMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
             +G   +V    WPL+T L+V DFKGNDPR
Sbjct  61  -TLGQVFYVRQNFWPLNTTLYVYDFKGNDPR  90


>gi|291296825|ref|YP_003508223.1| restriction modification system DNA specificity domain-containing 
protein [Meiothermus ruber DSM 1279]
 gi|290471784|gb|ADD29203.1| restriction modification system DNA specificity domain protein 
[Meiothermus ruber DSM 1279]
Length=419

 Score = 95.9 bits (237),  Expect = 2e-18, Method: Compositional matrix adjust.
 Identities = 47/91 (52%), Positives = 56/91 (62%), Gaps = 1/91 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   WK    GEV+EL+RG+DLP   R  G VP++ S GVT  H  A   GPGV  GR G
Sbjct  1   MGGEWKECALGEVIELKRGYDLPQQDRRPGYVPIVSSSGVTDYHAEAMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
             +G   +V    WPL+T L+VRDFKGNDPR
Sbjct  61  -TLGEVFYVEQDFWPLNTTLYVRDFKGNDPR  90


>gi|339834363|gb|EGQ62133.1| type I restriction-modification system, S subunit [Acidithiobacillus 
sp. GGI-221]
Length=355

 Score = 95.9 bits (237),  Expect = 2e-18, Method: Composition-based stats.
 Identities = 46/91 (51%), Positives = 58/91 (64%), Gaps = 1/91 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            M++ WK    G+V+EL+RG+DLP   R SG VP++ S GVT  H  A   GPGV  GR G
Sbjct  18   MANEWKECSLGDVIELKRGYDLPQKERLSGDVPLVSSSGVTDTHAKAMVKGPGVVTGRYG  77

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
              +G   +V    WPL+T L+V DFKGNDPR
Sbjct  78   -TLGQVFYVRQNFWPLNTTLYVYDFKGNDPR  107


>gi|2408222|gb|AAB70708.1| HsdS [Klebsiella pneumoniae]
Length=439

 Score = 95.1 bits (235),  Expect = 3e-18, Method: Composition-based stats.
 Identities = 42/90 (47%), Positives = 54/90 (60%), Gaps = 1/90 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   WK    G+ +EL+RG+DLP ++R  G++P+I S G T  HD     GPGV  GR G
Sbjct  1   MGSKWKECELGDFIELKRGYDLPKSTRNEGSIPIISSSGFTDFHDKPMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             IG   +     WPL+T L+V DFKGNDP
Sbjct  61  -TIGEVFYSEEDFWPLNTTLYVVDFKGNDP  89


>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 = 94.7 bits (234),  Expect = 4e-18, Method: Compositional matrix adjust.
 Identities = 46/91 (51%), Positives = 57/91 (63%), Gaps = 1/91 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M+  WK    G+V+EL+RG+DLP   R  G VP++ S GVT  H  A   GPGV  GR G
Sbjct  1   MAGEWKECSLGDVIELKRGYDLPQKDRLPGDVPLVSSSGVTDTHAKAMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
             +G   +V    WPL+T L+VRDFKGNDPR
Sbjct  61  -TLGQVFYVEQDFWPLNTTLYVRDFKGNDPR  90


>gi|121583503|ref|YP_973929.1| restriction modification system DNA specificity subunit [Polaromonas 
naphthalenivorans CJ2]
 gi|120596753|gb|ABM40187.1| restriction modification system DNA specificity domain [Polaromonas 
naphthalenivorans CJ2]
Length=415

 Score = 94.0 bits (232),  Expect = 6e-18, Method: Compositional matrix adjust.
 Identities = 45/91 (50%), Positives = 57/91 (63%), Gaps = 1/91 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           MS  W+  + G+ +EL+RG+DLP A R SG  P++ S GV+  H      GPGV  GR G
Sbjct  1   MSSEWQFGKLGDFIELKRGYDLPQAKRTSGPFPLVSSSGVSDCHSVPMVRGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
             IG   FV    WPL+T L+VRDFKGNDP+
Sbjct  61  -TIGQVYFVEDDFWPLNTTLYVRDFKGNDPK  90


>gi|331650404|ref|ZP_08351476.1| type I restriction-modification system, S subunit, EcoA family 
[Escherichia coli M605]
 gi|331040798|gb|EGI12956.1| type I restriction-modification system, S subunit, EcoA family 
[Escherichia coli M605]
Length=440

 Score = 93.6 bits (231),  Expect = 1e-17, Method: Composition-based stats.
 Identities = 42/90 (47%), Positives = 52/90 (58%), Gaps = 1/90 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           MS  W     GE + L+RG+DLP + R  G +PVI S G +G HD     GPGV  GR G
Sbjct  1   MSSEWINTTLGEFITLKRGYDLPKSKRNDGNIPVISSSGYSGTHDVPMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             IG   +V    WP++T L+V DFKGN P
Sbjct  61  -TIGEVFYVVDDFWPINTTLYVSDFKGNSP  89


>gi|49658898|emb|CAF28524.1| putative HsdS-like DNA methylase [Yersinia pseudotuberculosis]
Length=449

 Score = 92.8 bits (229),  Expect = 2e-17, Method: Composition-based stats.
 Identities = 43/90 (48%), Positives = 53/90 (59%), Gaps = 1/90 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   W  +  GE L L+RG+DLP + R SG +P+I S G +G HD     GPGV  GR G
Sbjct  1   MGSEWLDITLGEFLNLKRGYDLPKSKRNSGNIPIISSSGFSGNHDKPMVYGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             IG   +V    WPL+T L+V DFKGN P
Sbjct  61  -TIGEVFYVNESYWPLNTTLYVDDFKGNSP  89


>gi|322691669|ref|YP_004221239.1| hypothetical protein BLLJ_1480 [Bifidobacterium longum subsp. 
longum JCM 1217]
 gi|320456525|dbj|BAJ67147.1| conserved hypothetical protein [Bifidobacterium longum subsp. 
longum JCM 1217]
Length=363

 Score = 91.3 bits (225),  Expect = 4e-17, Method: Composition-based stats.
 Identities = 42/87 (49%), Positives = 54/87 (63%), Gaps = 1/87 (1%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           W+    GE+ +L+RG DLP + R  G VPV  S G+TG + TAA +GP V  GR G  IG
Sbjct  3   WRETTLGEITDLKRGFDLPKSQRLQGDVPVYSSSGITGSNSTAAVEGPCVITGRYG-TIG  61

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
              F  GP WPL+T L+  +F GN+PR
Sbjct  62  EVFFSGGPCWPLNTALYSTEFNGNNPR  88


>gi|86146743|ref|ZP_01065063.1| type I restriction-modification system, S subunit, EcoA family 
protein [Vibrio sp. MED222]
 gi|85835393|gb|EAQ53531.1| type I restriction-modification system, S subunit, EcoA family 
protein [Vibrio sp. MED222]
Length=417

 Score = 90.9 bits (224),  Expect = 6e-17, Method: Composition-based stats.
 Identities = 42/90 (47%), Positives = 53/90 (59%), Gaps = 1/90 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   W     G+V+EL+RG+DLP   R  G VPVI S G +G H+      PGV  GR G
Sbjct  1   MGSEWIQSELGDVIELKRGYDLPKTKRIDGNVPVISSSGHSGFHNEVKIKSPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             IG   ++    WPL+T L+V+DFKGNDP
Sbjct  61  -TIGQVFYIEEDFWPLNTTLYVKDFKGNDP  89


>gi|344198531|ref|YP_004782857.1| restriction modification system DNA specificity domain-containing 
protein [Acidithiobacillus ferrivorans SS3]
 gi|343773975|gb|AEM46531.1| restriction modification system DNA specificity domain-containing 
protein [Acidithiobacillus ferrivorans SS3]
Length=360

 Score = 89.4 bits (220),  Expect = 2e-16, Method: Compositional matrix adjust.
 Identities = 44/90 (49%), Positives = 57/90 (64%), Gaps = 1/90 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            +  GW+  R  +VL LQRG DLP ASR  GTVP+I + GVTG H+ A    PGV  GRSG
Sbjct  181  IPQGWEVGRLDDVLVLQRGFDLPKASRTEGTVPIIAATGVTGFHNEAKVKAPGVVTGRSG  240

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             +IG   +V    WPL+T L+V++F   +P
Sbjct  241  -SIGEVLYVQEDFWPLNTSLWVKEFPKGEP  269


 Score = 76.3 bits (186),  Expect = 1e-12, Method: Compositional matrix adjust.
 Identities = 37/84 (45%), Positives = 51/84 (61%), Gaps = 1/84 (1%)

Query  8   LRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTAT  67
           +RF + + LQRG DL       G  PVIGS  + G H+    + PGV  GRSG ++G   
Sbjct  1   MRFSDFVTLQRGFDLTKTEMRDGPYPVIGSTSIIGHHNEFKANPPGVVTGRSG-SLGVVQ  59

Query  68  FVAGPIWPLDTCLFVRDFKGNDPR  91
           +VA   WP +T L+V+DF+GNDP+
Sbjct  60  YVAERYWPHNTSLWVKDFRGNDPK  83


>gi|288937352|ref|YP_003441411.1| restriction modification system DNA specificity domain-containing 
protein [Klebsiella variicola At-22]
 gi|288892061|gb|ADC60379.1| restriction modification system DNA specificity domain protein 
[Klebsiella variicola At-22]
Length=430

 Score = 89.4 bits (220),  Expect = 2e-16, Method: Composition-based stats.
 Identities = 40/89 (45%), Positives = 52/89 (59%), Gaps = 1/89 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           M   W     G+ +EL+RG+DLP ++R  G++P+I S G T  HD     GPGV  GR G
Sbjct  1   MGSKWIECELGDFIELKRGYDLPKSTRNEGSIPIISSSGFTDFHDKPMVKGPGVVTGRYG  60

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGND  89
             IG   +     WPL+T L+V DFKGND
Sbjct  61  -TIGEVFYSEEDFWPLNTTLYVVDFKGND  88


>gi|253687261|ref|YP_003016451.1| restriction modification system DNA specificity domain protein 
[Pectobacterium carotovorum subsp. carotovorum PC1]
 gi|251753839|gb|ACT11915.1| restriction modification system DNA specificity domain protein 
[Pectobacterium carotovorum subsp. carotovorum PC1]
Length=413

 Score = 88.6 bits (218),  Expect = 3e-16, Method: Composition-based stats.
 Identities = 41/87 (48%), Positives = 52/87 (60%), Gaps = 1/87 (1%)

Query  4    GWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
            GW   + G+VL LQRG DLP   R  G +P+I S G +G H+ A    PG+  GR G  I
Sbjct  21   GWLQCKLGDVLTLQRGFDLPQRLRKEGNIPIISSSGESGWHNNAIVSPPGIVTGRYG-TI  79

Query  64   GTATFVAGPIWPLDTCLFVRDFKGNDP  90
            G   F+  P WPL+T L+VR+FKG  P
Sbjct  80   GEVFFIDKPFWPLNTTLYVREFKGITP  106


>gi|219851731|ref|YP_002466163.1| restriction modification system DNA specificity domain protein 
[Methanosphaerula palustris E1-9c]
 gi|219545990|gb|ACL16440.1| restriction modification system DNA specificity domain protein 
[Methanosphaerula palustris E1-9c]
Length=205

 Score = 86.7 bits (213),  Expect = 9e-16, Method: Compositional matrix adjust.
 Identities = 39/89 (44%), Positives = 56/89 (63%), Gaps = 1/89 (1%)

Query  3    DGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAA  62
            + W+  RF E + LQRG+D+    +  G VPV+ S G+   H+ +  +GPGV IGR G  
Sbjct  31   ERWREGRFDEFILLQRGYDITKDEQHDGIVPVVSSSGIKSFHNESRANGPGVVIGRKG-T  89

Query  63   IGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
            +G   +V  P WP DT L+V+DFKGN+P+
Sbjct  90   LGKVFYVDCPYWPHDTSLWVKDFKGNNPK  118


>gi|253563523|ref|ZP_04840980.1| restriction modification system DNA specificity subunit [Bacteroides 
sp. 3_2_5]
 gi|251947299|gb|EES87581.1| restriction modification system DNA specificity subunit [Bacteroides 
sp. 3_2_5]
Length=151

 Score = 85.5 bits (210),  Expect = 2e-15, Method: Compositional matrix adjust.
 Identities = 38/83 (46%), Positives = 54/83 (66%), Gaps = 1/83 (1%)

Query  9   RFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATF  68
           + GE+LELQRG+DLP++    G + V GS GV G H+ A  + P + +GRSG ++G   +
Sbjct  7   KLGEILELQRGYDLPSSQMKKGDILVAGSNGVIGYHNEARSNHPCITVGRSG-SVGKVHY  65

Query  69  VAGPIWPLDTCLFVRDFKGNDPR  91
                W  +T LFV+DFKGNDP+
Sbjct  66  YEQATWAHNTALFVKDFKGNDPK  88


>gi|32455520|ref|NP_862272.1| HsdS' [Lactobacillus sakei]
 gi|24461247|gb|AAN61994.1|AF438419_4 HsdS' [Lactobacillus sakei]
Length=151

 Score = 84.7 bits (208),  Expect = 4e-15, Method: Compositional matrix adjust.
 Identities = 43/82 (53%), Positives = 52/82 (64%), Gaps = 2/82 (2%)

Query  10  FGEVLELQRGHDLPAASRGS-GTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATF  68
           F + +EL RG+DLP A R + G VPVI + G  G HD A   GPGV  GRSG  +G+  F
Sbjct  3   FSDFIELHRGYDLPKAKRNADGLVPVIAAGGYDGFHDAAKVQGPGVITGRSG-TLGSVYF  61

Query  69  VAGPIWPLDTCLFVRDFKGNDP  90
                WPL+T LFV +FKGNDP
Sbjct  62  NTTDFWPLNTTLFVSNFKGNDP  83


>gi|301162175|emb|CBW21720.1| putative type I DNA restriction-modification [Bacteroides fragilis 
638R]
Length=399

 Score = 84.7 bits (208),  Expect = 4e-15, Method: Composition-based stats.
 Identities = 38/83 (46%), Positives = 54/83 (66%), Gaps = 1/83 (1%)

Query  9   RFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATF  68
           + GE+LELQRG+DLP++    G + V GS GV G H+ A  + P + +GRSG ++G   +
Sbjct  7   KLGEILELQRGYDLPSSQMKKGDILVAGSNGVIGYHNEARSNHPCITVGRSG-SVGKVHY  65

Query  69  VAGPIWPLDTCLFVRDFKGNDPR  91
                W  +T LFV+DFKGNDP+
Sbjct  66  YEQATWAHNTALFVKDFKGNDPK  88


>gi|160945577|ref|ZP_02092803.1| hypothetical protein FAEPRAM212_03106 [Faecalibacterium prausnitzii 
M21/2]
 gi|158443308|gb|EDP20313.1| hypothetical protein FAEPRAM212_03106 [Faecalibacterium prausnitzii 
M21/2]
Length=156

 Score = 84.0 bits (206),  Expect = 6e-15, Method: Compositional matrix adjust.
 Identities = 38/87 (44%), Positives = 50/87 (58%), Gaps = 1/87 (1%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           WK    GE + L+RG+DLP   R +G +P+  S GVTG H+ A  + PGV  GR G  IG
Sbjct  4   WKIDELGEFVTLKRGYDLPQQKRKNGEIPIFSSSGVTGTHNEAMVEAPGVITGRYG-TIG  62

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
              F     WPL+T L+V +F  ND +
Sbjct  63  EVFFAETSFWPLNTTLYVENFHENDEK  89


>gi|126726006|ref|ZP_01741848.1| putative type I restriction-modification system, S subunit [Rhodobacterales 
bacterium HTCC2150]
 gi|126705210|gb|EBA04301.1| putative type I restriction-modification system, S subunit [Rhodobacterales 
bacterium HTCC2150]
Length=371

 Score = 81.6 bits (200),  Expect = 3e-14, Method: Composition-based stats.
 Identities = 40/86 (47%), Positives = 51/86 (60%), Gaps = 1/86 (1%)

Query  3   DGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAA  62
           +GW   R GEV E QRG DLP + R  G +P+I S G +G H TA  +  G+  GR G +
Sbjct  4   EGWGECRLGEVTEFQRGFDLPKSQRQVGEIPIISSAGYSGWHSTAKVERAGIVTGRYG-S  62

Query  63  IGTATFVAGPIWPLDTCLFVRDFKGN  88
           IG    V    WPL+T L+V+DF GN
Sbjct  63  IGDVFLVYEDHWPLNTTLWVKDFHGN  88


>gi|242278888|ref|YP_002991017.1| restriction modification system DNA specificity domain protein 
[Desulfovibrio salexigens DSM 2638]
 gi|242121782|gb|ACS79478.1| restriction modification system DNA specificity domain protein 
[Desulfovibrio salexigens DSM 2638]
Length=387

 Score = 81.6 bits (200),  Expect = 4e-14, Method: Composition-based stats.
 Identities = 39/91 (43%), Positives = 51/91 (57%), Gaps = 1/91 (1%)

Query  1   MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
           +  GW     GE   LQRG DLP   R  G  P+I S G    H+     GPGV  GRSG
Sbjct  3   LPKGWDKKTLGESCTLQRGFDLPKRLRVKGEHPLISSSGCIDSHNEPKVAGPGVVTGRSG  62

Query  61  AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
            +IG+  ++    WPL+T L+V+++ GNDPR
Sbjct  63  -SIGSLFYIEDDFWPLNTTLYVKNYFGNDPR  92


>gi|194364815|ref|YP_002027425.1| restriction modification system DNA specificity domain [Stenotrophomonas 
maltophilia R551-3]
 gi|194347619|gb|ACF50742.1| restriction modification system DNA specificity domain [Stenotrophomonas 
maltophilia R551-3]
Length=364

 Score = 81.3 bits (199),  Expect = 4e-14, Method: Compositional matrix adjust.
 Identities = 39/88 (45%), Positives = 49/88 (56%), Gaps = 1/88 (1%)

Query  4   GWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
           GW  +    +  L+RG+DLP   R  G VP+  + G  G HD    +GPGV  GRSG  I
Sbjct  9   GWPLVPLKNIATLKRGYDLPVGMRNKGEVPIYAANGQNGSHDEVKINGPGVITGRSG-TI  67

Query  64  GTATFVAGPIWPLDTCLFVRDFKGNDPR  91
           G   +  G  WPL+T L+V DF GN PR
Sbjct  68  GRVHYCEGGFWPLNTALYVMDFHGNHPR  95


>gi|260767611|ref|ZP_05876547.1| type I restriction-modification system specificity subunit S 
[Vibrio furnissii CIP 102972]
 gi|260617511|gb|EEX42694.1| type I restriction-modification system specificity subunit S 
[Vibrio furnissii CIP 102972]
Length=374

 Score = 81.3 bits (199),  Expect = 5e-14, Method: Composition-based stats.
 Identities = 43/87 (50%), Positives = 61/87 (71%), Gaps = 1/87 (1%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           W   + G++L L+RG+DLP ++R SG+VPV+ S G+TG H+TA  +GP V  GR G  +G
Sbjct  3   WVECQLGDILTLKRGYDLPHSARKSGSVPVVSSSGITGYHNTAKVEGPAVVTGRYG-TLG  61

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
              +V G  WPL+T L+V+DFKGN P+
Sbjct  62  EVYYVEGECWPLNTSLYVQDFKGNRPK  88


>gi|159904437|ref|YP_001548099.1| restriction modification system DNA specificity subunit [Methanococcus 
maripaludis C6]
 gi|159885930|gb|ABX00867.1| restriction modification system DNA specificity domain [Methanococcus 
maripaludis C6]
Length=397

 Score = 81.3 bits (199),  Expect = 5e-14, Method: Composition-based stats.
 Identities = 35/87 (41%), Positives = 52/87 (60%), Gaps = 1/87 (1%)

Query  5    WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
            W+    GE++  QRGHDLP  SR +G  PV+ S G+ G H+    +  G+ IGRSG  +G
Sbjct  17   WEVKEIGELVTFQRGHDLPVNSRKNGIYPVVASNGIVGYHNEYKVENEGLTIGRSG-NLG  75

Query  65   TATFVAGPIWPLDTCLFVRDFKGNDPR  91
               +V+   WPL+T L+V+ F  + P+
Sbjct  76   EPFYVSTSFWPLNTTLYVKKFHNSHPK  102


>gi|344340139|ref|ZP_08771065.1| restriction modification system DNA specificity domain protein 
[Thiocapsa marina 5811]
 gi|343799797|gb|EGV17745.1| restriction modification system DNA specificity domain protein 
[Thiocapsa marina 5811]
Length=386

 Score = 80.9 bits (198),  Expect = 5e-14, Method: Composition-based stats.
 Identities = 37/87 (43%), Positives = 52/87 (60%), Gaps = 1/87 (1%)

Query  5   WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
           W+     +++  QRG+D+  A + SG +PVI S G+T  H+ A   GPGV IGR G  +G
Sbjct  11  WRDATLQDLVFFQRGYDITKAQQHSGNIPVISSSGITSFHNEAKAQGPGVVIGRKG-TLG  69

Query  65  TATFVAGPIWPLDTCLFVRDFKGNDPR  91
           +  + A   WP DT L+ RD +GN PR
Sbjct  70  SVHYSAQDYWPHDTTLWSRDLRGNIPR  96


>gi|119357508|ref|YP_912152.1| restriction modification system DNA specificity subunit [Chlorobium 
phaeobacteroides DSM 266]
 gi|119354857|gb|ABL65728.1| restriction modification system DNA specificity domain [Chlorobium 
phaeobacteroides DSM 266]
Length=414

 Score = 80.5 bits (197),  Expect = 9e-14, Method: Composition-based stats.
 Identities = 40/87 (46%), Positives = 50/87 (58%), Gaps = 1/87 (1%)

Query  4    GWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
            GW+  R  +VL LQRG DLP A R  GTVP+  + GVTG H  A    P V  GRSG  I
Sbjct  242  GWEAGRLDDVLVLQRGFDLPKAKRMEGTVPIYAATGVTGFHCEAKVKAPCVVTGRSG-TI  300

Query  64   GTATFVAGPIWPLDTCLFVRDFKGNDP  90
            G   +V    WPL+T L+ + F  ++P
Sbjct  301  GDVIYVQEDFWPLNTSLWAKGFPKSEP  327


>gi|21226532|ref|NP_632454.1| type I restriction-modification system specificity subunit [Methanosarcina 
mazei Go1]
 gi|20904802|gb|AAM30126.1| type I restriction-modification system specificity subunit [Methanosarcina 
mazei Go1]
Length=439

 Score = 80.1 bits (196),  Expect = 9e-14, Method: Composition-based stats.
 Identities = 41/86 (48%), Positives = 54/86 (63%), Gaps = 1/86 (1%)

Query  6   KTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGT  65
           KT+   E + LQRG DLPA+ R +G VPV+ S G+ G HD A  + PGV I     +IG 
Sbjct  4   KTIALNEFIILQRGFDLPASDRKNGNVPVVASTGIGGYHDEAKVEAPGVVI-GRSGSIGG  62

Query  66  ATFVAGPIWPLDTCLFVRDFKGNDPR  91
             F+    WPL+T L+V+DFKG+ PR
Sbjct  63  GQFITENFWPLNTTLWVKDFKGHHPR  88


>gi|343920509|gb|EGV31240.1| restriction modification system DNA specificity domain protein 
[Thiorhodococcus drewsii AZ1]
Length=409

 Score = 79.7 bits (195),  Expect = 1e-13, Method: Composition-based stats.
 Identities = 39/87 (45%), Positives = 52/87 (60%), Gaps = 1/87 (1%)

Query  5    WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
            W     G+ ++LQRG D+    R  G  PVI S G+   HDTA   GPGV +GR G+A G
Sbjct  235  WDCCGLGQKIKLQRGVDITKDKRIEGPFPVISSGGIDYFHDTALCAGPGVLVGRKGSA-G  293

Query  65   TATFVAGPIWPLDTCLFVRDFKGNDPR  91
               +V    WP DT L+VR+F+GN+P+
Sbjct  294  KLHYVETDFWPHDTTLYVREFRGNEPK  320


>gi|284108343|ref|ZP_06386407.1| Restriction modification system DNA specificity domain [Candidatus 
Poribacteria sp. WGA-A3]
 gi|283829904|gb|EFC34190.1| Restriction modification system DNA specificity domain [Candidatus 
Poribacteria sp. WGA-A3]
Length=393

 Score = 78.2 bits (191),  Expect = 4e-13, Method: Composition-based stats.
 Identities = 37/79 (47%), Positives = 46/79 (59%), Gaps = 1/79 (1%)

Query  4   GWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
           GW+T R G+VL+L  G  LP  SR  G +PV GS GV G H  A  D PG+ +GR G+A 
Sbjct  3   GWQTKRLGDVLQLNYGKSLPVKSRVEGPIPVYGSNGVVGSHTEAIVDAPGLIVGRKGSA-  61

Query  64  GTATFVAGPIWPLDTCLFV  82
           G      GP  P+DT  +V
Sbjct  62  GQVHLSRGPFCPIDTTFYV  80


>gi|257051191|ref|YP_003129024.1| restriction modification system DNA specificity domain protein 
[Halorhabdus utahensis DSM 12940]
 gi|256689954|gb|ACV10291.1| restriction modification system DNA specificity domain protein 
[Halorhabdus utahensis DSM 12940]
Length=442

 Score = 78.2 bits (191),  Expect = 4e-13, Method: Composition-based stats.
 Identities = 36/80 (45%), Positives = 49/80 (62%), Gaps = 1/80 (1%)

Query  3    DGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAA  62
            + W  +R GE+L L+ G +LP+ SR SGTVPV GS G    H  AA + PG+ +GR G +
Sbjct  38   ESWNLVRLGEILTLEYGDNLPSDSRESGTVPVFGSNGQVDTHSEAAVEKPGIILGRKG-S  96

Query  63   IGTATFVAGPIWPLDTCLFV  82
            IG   F   P WP+DT  ++
Sbjct  97   IGEIDFSDRPFWPIDTTYYI  116


>gi|329913308|ref|ZP_08275914.1| Type I restriction-modification system, specificity subunit S 
[Oxalobacteraceae bacterium IMCC9480]
 gi|327545395|gb|EGF30614.1| Type I restriction-modification system, specificity subunit S 
[Oxalobacteraceae bacterium IMCC9480]
Length=517

 Score = 78.2 bits (191),  Expect = 4e-13, Method: Compositional matrix adjust.
 Identities = 37/86 (44%), Positives = 50/86 (59%), Gaps = 1/86 (1%)

Query  5    WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
            W  +  G+V+ LQRG DLP+  R  G VP++ S GV+  +      GPGV  GR G  IG
Sbjct  112  WAEVPLGDVITLQRGFDLPSQKRKPGKVPIVSSSGVSDYNSEVGVKGPGVVTGRYG-TIG  170

Query  65   TATFVAGPIWPLDTCLFVRDFKGNDP  90
                +    WPL+T L+V++F GNDP
Sbjct  171  QVFLIKEDFWPLNTTLWVKNFHGNDP  196


>gi|313140399|ref|ZP_07802592.1| type I restriction-modification system [Bifidobacterium bifidum 
NCIMB 41171]
 gi|313132909|gb|EFR50526.1| type I restriction-modification system [Bifidobacterium bifidum 
NCIMB 41171]
Length=167

 Score = 77.4 bits (189),  Expect = 6e-13, Method: Compositional matrix adjust.
 Identities = 39/91 (43%), Positives = 50/91 (55%), Gaps = 1/91 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
             +D W+  + GEV  LQRG DLP      G  PV+ S G+ G H      GPGV  GRSG
Sbjct  16   FTDPWEQRKLGEVAPLQRGFDLPVNQMTPGPYPVVMSNGIGGWHSKYMVKGPGVVTGRSG  75

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
              IG+  ++    WP +T L+V  F GN+PR
Sbjct  76   -TIGSLHYIEQNFWPHNTSLWVTSFNGNEPR  105


>gi|148262630|ref|YP_001229336.1| restriction modification system DNA specificity subunit [Geobacter 
uraniireducens Rf4]
 gi|146396130|gb|ABQ24763.1| restriction modification system DNA specificity domain [Geobacter 
uraniireducens Rf4]
Length=385

 Score = 75.9 bits (185),  Expect = 2e-12, Method: Composition-based stats.
 Identities = 33/83 (40%), Positives = 49/83 (60%), Gaps = 1/83 (1%)

Query  9   RFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATF  68
           R G+++  +RG+DLP+  R  G  P++ S G++G H      G G+  GR G  +G   +
Sbjct  8   RLGDIVNFKRGYDLPSYKRKEGPYPIVSSSGISGYHAEYKAKGEGLITGRYG-TLGEMYY  66

Query  69  VAGPIWPLDTCLFVRDFKGNDPR  91
           V G  WP +T L+V DFKGN P+
Sbjct  67  VNGKYWPHNTALYVTDFKGNYPK  89


>gi|315648621|ref|ZP_07901718.1| restriction modification system DNA specificity domain protein 
[Paenibacillus vortex V453]
 gi|315276000|gb|EFU39348.1| restriction modification system DNA specificity domain protein 
[Paenibacillus vortex V453]
Length=389

 Score = 75.9 bits (185),  Expect = 2e-12, Method: Composition-based stats.
 Identities = 38/90 (43%), Positives = 50/90 (56%), Gaps = 1/90 (1%)

Query  2    SDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGA  61
            +D W+  +  ++  LQRG DLP +   +G+ PVI S G+   H       PGV  GRSG 
Sbjct  15   TDAWEQRKVIDIAPLQRGFDLPVSEMEAGSYPVIMSNGIGAYHSKYKAKAPGVVTGRSG-  73

Query  62   AIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
             IG  TFV    WP +T L+V DFK ND +
Sbjct  74   TIGNLTFVEVDYWPHNTALWVTDFKRNDAK  103


>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 = 74.3 bits (181),  Expect = 6e-12, Method: Composition-based stats.
 Identities = 34/80 (43%), Positives = 49/80 (62%), Gaps = 1/80 (1%)

Query  10  FGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATFV  69
            G+++E QRG+DLP +    G  PV  S G+ G H+    +GPG+ IGRSG  +G    +
Sbjct  9   LGDLVEFQRGYDLPKSKFVEGKYPVQSSNGILGYHNEYKVEGPGITIGRSG-TVGNPHLI  67

Query  70  AGPIWPLDTCLFVRDFKGND  89
               +P +T LFV++FKGND
Sbjct  68  RENFFPHNTSLFVKEFKGND  87


>gi|323965458|gb|EGB60913.1| type I restriction modification DNA specificity domain-containing 
protein [Escherichia coli M863]
 gi|327250265|gb|EGE61984.1| type I restriction modification DNA specificity domain protein 
[Escherichia coli STEC_7v]
Length=438

 Score = 73.9 bits (180),  Expect = 8e-12, Method: Composition-based stats.
 Identities = 34/87 (40%), Positives = 48/87 (56%), Gaps = 1/87 (1%)

Query  4    GWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
            GW+  +  ++L LQRG DLP+  R  G  P+I + G  G HD A    PG+  GRSG  +
Sbjct  256  GWEIGKLQDLLILQRGFDLPSTQRNIGLHPIIAASGYNGTHDIAMVKAPGIVTGRSG-VL  314

Query  64   GTATFVAGPIWPLDTCLFVRDFKGNDP  90
            G    +    WPL+T L+V++ K   P
Sbjct  315  GNVFLILEDFWPLNTTLWVKELKHATP  341


>gi|75909705|ref|YP_324001.1| restriction modification system DNA specificity subunit [Anabaena 
variabilis ATCC 29413]
 gi|75703430|gb|ABA23106.1| Restriction modification system DNA specificity domain protein 
[Anabaena variabilis ATCC 29413]
Length=405

 Score = 73.2 bits (178),  Expect = 1e-11, Method: Composition-based stats.
 Identities = 35/90 (39%), Positives = 49/90 (55%), Gaps = 1/90 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            + + W+  RF E   LQRG+DLP      G  PVI S  +   H+    + P +  GRSG
Sbjct  217  IPESWERKRFDEFCLLQRGYDLPDTQVIPGQYPVIASTSIKTYHNQFKVNPPVITTGRSG  276

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
             ++G   F+    WPL+T LFV++F GN P
Sbjct  277  -SLGIILFINSQAWPLNTTLFVKNFYGNSP  305


>gi|295136493|ref|YP_003587169.1| restriction modification system DNA specificity subunit [Zunongwangia 
profunda SM-A87]
 gi|294984508|gb|ADF54973.1| restriction modification system DNA specificity subunit [Zunongwangia 
profunda SM-A87]
Length=405

 Score = 71.6 bits (174),  Expect = 3e-11, Method: Composition-based stats.
 Identities = 31/90 (35%), Positives = 50/90 (56%), Gaps = 1/90 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            + +GW      ++  LQRG DLP   R  G VP+  + G +  H+    + PGV  GRSG
Sbjct  220  IPEGWVVANLEDLFVLQRGFDLPKKKRIEGNVPIYAASGKSTYHNEYKVEAPGVTTGRSG  279

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDP  90
              +G   FV+   WPL+T L++++++ + P
Sbjct  280  -VLGNVYFVSEDFWPLNTSLWIKEYRSSTP  308


>gi|339730590|emb|CCC41940.1| type I site-specific deoxyribonuclease subunit rmeS [Haloquadratum 
walsbyi C23]
Length=424

 Score = 71.6 bits (174),  Expect = 4e-11, Method: Composition-based stats.
 Identities = 33/82 (41%), Positives = 45/82 (55%), Gaps = 1/82 (1%)

Query  3    DGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAA  62
            + W+ LR G+VL  + G +LPA +R  G +PV GS G    H  A  D  G+ IGR G +
Sbjct  42   EDWELLRLGDVLSFEYGDNLPADTREDGEIPVYGSNGQVDTHTEANVDTAGIIIGRKG-S  100

Query  63   IGTATFVAGPIWPLDTCLFVRD  84
            IG   F   P W +DT  ++ D
Sbjct  101  IGKGEFSINPFWAIDTTYYISD  122


>gi|334123285|ref|ZP_08497313.1| type I restriction system specificity protein [Enterobacter hormaechei 
ATCC 49162]
 gi|333390907|gb|EGK62032.1| type I restriction system specificity protein [Enterobacter hormaechei 
ATCC 49162]
Length=421

 Score = 71.2 bits (173),  Expect = 4e-11, Method: Composition-based stats.
 Identities = 34/86 (40%), Positives = 49/86 (57%), Gaps = 1/86 (1%)

Query  5    WKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIG  64
            W+  +  ++L LQRG DLP+++R  G  P+I + G  G H+ A    PGV  GRSG  +G
Sbjct  240  WELGKLQDLLVLQRGFDLPSSARKDGEFPLIAASGPNGTHNVAMAKAPGVITGRSG-VLG  298

Query  65   TATFVAGPIWPLDTCLFVRDFKGNDP  90
                     WPL+T L+V++FK   P
Sbjct  299  KVYLTLEDYWPLNTTLWVKEFKRATP  324


>gi|262369880|ref|ZP_06063207.1| sty SBLI [Acinetobacter johnsonii SH046]
 gi|262314919|gb|EEY95959.1| sty SBLI [Acinetobacter johnsonii SH046]
Length=434

 Score = 71.2 bits (173),  Expect = 4e-11, Method: Composition-based stats.
 Identities = 37/80 (47%), Positives = 49/80 (62%), Gaps = 1/80 (1%)

Query  12  EVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATFVAG  71
           E + LQRG DLP + R SG +PV+ S GV G H+    D PGV I     +IG   ++  
Sbjct  11  EFILLQRGFDLPQSDRISGDIPVVASTGVAGFHNEYKVDAPGVVI-GRSGSIGGGQYIKE  69

Query  72  PIWPLDTCLFVRDFKGNDPR  91
             WPL+T L+V+DFKG+D R
Sbjct  70  KFWPLNTTLWVKDFKGHDAR  89


>gi|227546694|ref|ZP_03976743.1| type I restriction enzyme HindVIIP specificity protein [Bifidobacterium 
longum subsp. infantis ATCC 55813]
 gi|227212840|gb|EEI80719.1| type I restriction enzyme HindVIIP specificity protein [Bifidobacterium 
longum subsp. infantis ATCC 55813]
Length=200

 Score = 70.5 bits (171),  Expect = 8e-11, Method: Compositional matrix adjust.
 Identities = 37/91 (41%), Positives = 47/91 (52%), Gaps = 1/91 (1%)

Query  1    MSDGWKTLRFGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
             +D W+  +   V  LQRG DLPA     G  PV+ S G+   H+     GPGV  GRSG
Sbjct  16   FTDPWEQRKLIMVAPLQRGFDLPAEKIIPGVYPVMMSNGIGAYHNEYKVKGPGVVTGRSG  75

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGNDPR  91
              IG   +V    WP +T L+V  F GN P+
Sbjct  76   -TIGNLQYVESAFWPHNTTLWVTKFYGNHPK  105


>gi|301055838|ref|YP_003794049.1| type I restriction modification enzyme protein S [Bacillus cereus 
biovar anthracis str. CI]
 gi|300378007|gb|ADK06911.1| type I restriction modification enzyme protein S [Bacillus cereus 
biovar anthracis str. CI]
Length=369

 Score = 70.1 bits (170),  Expect = 9e-11, Method: Composition-based stats.
 Identities = 36/86 (42%), Positives = 50/86 (59%), Gaps = 2/86 (2%)

Query  5    WKTLRFGEVLELQRGHDLPAASRGS-GTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAI  63
            W+T R   ++ LQRG+DLP  SR   G V + GS GV G+H+ A   G G+  GRSG +I
Sbjct  206  WETTRLDNIVVLQRGYDLPIKSRNEMGEVEIWGSNGVVGVHNEAKIIGGGIVTGRSG-SI  264

Query  64   GTATFVAGPIWPLDTCLFVRDFKGND  89
            G   +     W L+T LF ++  GN+
Sbjct  265  GNVYYTYKDFWALNTTLFSKETYGNN  290


>gi|328947972|ref|YP_004365309.1| restriction modification system DNA specificity domain protein 
[Treponema succinifaciens DSM 2489]
 gi|328448296|gb|AEB14012.1| restriction modification system DNA specificity domain protein 
[Treponema succinifaciens DSM 2489]
Length=337

 Score = 69.7 bits (169),  Expect = 1e-10, Method: Composition-based stats.
 Identities = 38/89 (43%), Positives = 50/89 (57%), Gaps = 3/89 (3%)

Query  2    SDGWKTLRFGEVLELQRGHDLPAASRGS-GTVPVIGSFGVTGMHDTAAYDGPGVAIGRSG  60
            S  + T +  +V+ +QRG+DLP   R S G +PV GS G+ G H+ A  D  G+  GRSG
Sbjct  171  SKNFPTKKVIDVVTMQRGYDLPVQDRDSKGKIPVFGSNGILGNHNLAKMD-KGIITGRSG  229

Query  61   AAIGTATFVAGPIWPLDTCLFVRDFKGND  89
              IG       P WPL+T LF  D  GN+
Sbjct  230  T-IGEVYMCETPFWPLNTTLFSNDTHGNN  257


>gi|261867039|ref|YP_003254961.1| restriction modification system DNA specificity subunit [Aggregatibacter 
actinomycetemcomitans D11S-1]
 gi|261412371|gb|ACX81742.1| restriction modification system DNA specificity subunit [Aggregatibacter 
actinomycetemcomitans D11S-1]
Length=317

 Score = 69.7 bits (169),  Expect = 1e-10, Method: Composition-based stats.
 Identities = 32/80 (40%), Positives = 47/80 (59%), Gaps = 1/80 (1%)

Query  10  FGEVLELQRGHDLPAASRGSGTVPVIGSFGVTGMHDTAAYDGPGVAIGRSGAAIGTATFV  69
            G+++E QRG+DLP  +   G  PV  S G+ G H+      PG+ IGRSG  +G    +
Sbjct  9   LGDLVEFQRGYDLPKDAFVKGEYPVQSSNGILGYHNEYKVKAPGITIGRSG-TVGIPHLI  67

Query  70  AGPIWPLDTCLFVRDFKGND  89
               +P +T L+V+DFKGN+
Sbjct  68  TKNFFPHNTALYVKDFKGNN  87



Lambda     K      H
   0.320    0.141    0.452 

Gapped
Lambda     K      H
   0.267   0.0410    0.140 

Effective search space used: 128725229700


  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