BLASTP 2.2.17 [Aug-26-2007]


Reference:
Altschul, Stephen F., 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:
Schäffer, Alejandro A., L. Aravind, Thomas L. Madden, 
Schaffer, Alejandro A., 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.

Query= TDE2740 
         (162 letters)

Database: nr 
           5,470,121 sequences; 1,894,087,724 total letters

Searching..................................................done



                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

gi|42528240|ref|NP_973338.1|  type I restriction-modificatio...   328   8e-89
gi|23452768|gb|AAN33154.1|  putative type I specificity subu...   154   2e-36
gi|40467|emb|CAA35604.1|  HsdS polypeptide, part of CfrA fam...   116   3e-25
gi|20807979|ref|NP_623150.1|  Restriction endonuclease S sub...   114   1e-24
gi|68550910|ref|ZP_00590345.1|  Restriction modification sys...   105   1e-21
gi|19881243|gb|AAM00852.1|AF486551_3  HsdS [Campylobacter je...   101   1e-20
gi|15792859|ref|NP_282682.1|  putative type I restriction en...   101   1e-20
gi|126661170|ref|ZP_01732247.1|  type I restriction-modifica...    98   2e-19
gi|84489266|ref|YP_447498.1|  hypothetical protein Msp_0455 ...    93   4e-18
gi|15669726|ref|NP_248539.1|  type I restriction-modificatio...    91   2e-17
gi|150388684|ref|YP_001318733.1|  restriction modification s...    91   2e-17
gi|134097471|ref|YP_001103132.1|  putative restriction modif...    87   3e-16
gi|121541757|ref|ZP_01673502.1|  Restriction endonuclease S ...    86   1e-15
gi|21228395|ref|NP_634317.1|  type I restriction-modificatio...    85   1e-15
gi|52082594|ref|YP_081385.1|  HsdS [Bacillus licheniformis A...    85   1e-15
gi|91772548|ref|YP_565240.1|  restriction modification syste...    85   2e-15
gi|42779918|ref|NP_977165.1|  type I restriction-modificatio...    83   4e-15
gi|148544100|ref|YP_001271470.1|  restriction modification s...    82   1e-14
gi|31983512|ref|NP_858124.1|  putative type i restriction en...    82   1e-14
gi|55820774|ref|YP_139216.1|  type I restriction-modificatio...    81   2e-14
gi|89210922|ref|ZP_01189304.1|  Restriction modification sys...    80   3e-14
gi|111222733|ref|YP_713527.1|  Type I restriction modificati...    79   9e-14
gi|77415025|ref|ZP_00791100.1|  Type I restriction modificat...    79   1e-13
gi|119491619|ref|ZP_01623491.1|  hypothetical protein L8106_...    77   2e-13
gi|146319105|ref|YP_001198817.1|  type I restriction-modific...    76   7e-13
gi|154173663|ref|YP_001408732.1|  type I restriction-modific...    75   9e-13
gi|156864174|gb|EDO57605.1|  hypothetical protein CLOL250_01...    75   9e-13
gi|156864459|gb|EDO57890.1|  hypothetical protein CLOL250_01...    75   9e-13
gi|146321309|ref|YP_001201020.1|  type I restriction-modific...    75   1e-12
gi|60681333|ref|YP_211477.1|  putative type IC restriction-m...    75   1e-12
gi|71904273|ref|YP_281076.1|  type I restriction-modificatio...    75   2e-12
gi|19746826|ref|NP_607962.1|  putative specificity determina...    74   2e-12
gi|139474405|ref|YP_001129121.1|  type I restriction-modific...    74   3e-12
gi|94989255|ref|YP_597356.1|  type I restriction-modificatio...    74   3e-12
gi|15675718|ref|NP_269892.1|  putative type I site-specific ...    74   4e-12
gi|21911179|ref|NP_665447.1|  putative type I site-specific ...    74   4e-12
gi|94993143|ref|YP_601242.1|  Type I restriction-modificatio...    73   4e-12
gi|91213998|ref|YP_543984.1|  HsdS, type I site-specific deo...    72   9e-12
gi|126090307|ref|YP_001041762.1|  restriction modification s...    71   2e-11
gi|71024881|ref|YP_263290.1|  HsdS [Lactococcus lactis subsp...    70   3e-11
gi|32455448|ref|NP_862562.1|  HsdS [Lactococcus lactis] >gi|...    70   3e-11
gi|153808175|ref|ZP_01960843.1|  hypothetical protein BACCAC...    69   8e-11
gi|15678961|ref|NP_276078.1|  type I restriction modificatio...    67   3e-10
gi|2129238|pir||B64316  restriction modification system S ch...    67   4e-10
gi|126007743|ref|ZP_01708319.1|  hypothetical protein Faci_0...    66   6e-10
gi|99078515|ref|YP_611773.1|  restriction modification syste...    65   1e-09
gi|15669898|ref|NP_247095.1|  Type I restriction-modificatio...    65   1e-09
gi|15611482|ref|NP_223133.1|  TYPE I RESTRICTION ENZYME (SPE...    65   2e-09
gi|21673510|ref|NP_661575.1|  type I restriction system spec...    64   4e-09
gi|149369906|ref|ZP_01889757.1|  type I restriction-modifica...    63   4e-09
gi|26554275|ref|NP_758209.1|  type I restriction-modificatio...    63   4e-09
gi|154174911|ref|YP_001408734.1|  putative type I restrictio...    63   5e-09
gi|17232094|ref|NP_488642.1|  type I site-specific deoxyribo...    62   8e-09
gi|15645090|ref|NP_207260.1|  type I restriction enzyme S pr...    62   1e-08
gi|154245043|ref|YP_001416001.1|  restriction modification s...    61   2e-08
gi|145622248|ref|ZP_01778209.1|  hypothetical protein PmobDR...    60   3e-08
gi|126208116|ref|YP_001053341.1|  putative type I restrictio...    60   4e-08
gi|75909474|ref|YP_323770.1|  Restriction modification syste...    60   4e-08
gi|2826248|gb|AAB98112.1|  type I restriction-modification e...    60   6e-08
gi|120601902|ref|YP_966302.1|  restriction modification syst...    59   7e-08
gi|21228301|ref|NP_634223.1|  type I restriction-modificatio...    59   8e-08
gi|32266920|ref|NP_860952.1|  hypothetical protein HH1421 [H...    59   9e-08
gi|67919739|ref|ZP_00513304.1|  Restriction modification sys...    59   1e-07
gi|23452743|gb|AAN33142.1|  putative type I specificity subu...    59   1e-07
gi|156866399|gb|EDO59771.1|  hypothetical protein CLOLEP_038...    58   1e-07
gi|94267246|ref|ZP_01290822.1|  Restriction modification sys...    57   3e-07
gi|26554274|ref|NP_758208.1|  type I restriction-modificatio...    57   4e-07
gi|94263107|ref|ZP_01286925.1|  Restriction modification sys...    57   4e-07
gi|10954528|ref|NP_044167.1|  type I restriction enzyme subu...    57   5e-07
gi|51893047|ref|YP_075738.1|  type I restriction-modificatio...    56   6e-07
gi|23452748|gb|AAN33144.1|  putative type I specificity subu...    55   9e-07
gi|149115916|ref|ZP_01842652.1|  restriction modification sy...    55   1e-06
gi|26991425|ref|NP_746850.1|  type I restriction-modificatio...    55   1e-06
gi|94266803|ref|ZP_01290467.1|  Restriction modification sys...    55   1e-06
gi|150401945|ref|YP_001329239.1|  restriction modification s...    55   1e-06
gi|119356723|ref|YP_911367.1|  N-6 DNA methylase [Chlorobium...    55   2e-06
gi|42528243|ref|NP_973341.1|  type I restriction-modificatio...    55   2e-06
gi|89098144|ref|ZP_01171029.1|  type I restriction modificat...    54   2e-06
gi|118743998|ref|ZP_01591997.1|  restriction modification sy...    54   2e-06
gi|148544101|ref|YP_001271471.1|  restriction modification s...    54   2e-06
gi|150006638|ref|YP_001301382.1|  type I restriction enzyme ...    54   3e-06
gi|124520753|ref|ZP_01695856.1|  restriction modification sy...    54   4e-06
gi|121608003|ref|YP_995810.1|  restriction modification syst...    54   4e-06
gi|152991445|ref|YP_001357167.1|  type I restriction-modific...    53   6e-06
>gi|42528240|ref|NP_973338.1| type I restriction-modification system, S subunit, truncation
           [Treponema denticola ATCC 35405]
 gi|41819510|gb|AAS13257.1| type I restriction-modification system, S subunit, truncation
           [Treponema denticola ATCC 35405]
          Length = 162

 Score =  328 bits (840), Expect = 8e-89,   Method: Composition-based stats.
 Identities = 162/162 (100%), Positives = 162/162 (100%)

Query: 1   MWCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKN 60
           MWCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKN
Sbjct: 1   MWCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKN 60

Query: 61  AILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTF 120
           AILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTF
Sbjct: 61  AILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTF 120

Query: 121 QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
           QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL
Sbjct: 121 QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
>gi|23452768|gb|AAN33154.1| putative type I specificity subunit HsdS [Campylobacter jejuni]
 gi|23452773|gb|AAN33157.1| putative type I specificity subunit HsdS [Campylobacter jejuni]
 gi|23452787|gb|AAN33165.1| putative type I specificity subunit HsdS [Campylobacter jejuni]
          Length = 403

 Score =  154 bits (388), Expect = 2e-36,   Method: Composition-based stats.
 Identities = 77/156 (49%), Positives = 108/156 (69%), Gaps = 4/156 (2%)

Query: 4   RLGEICSITMGQSPESSFI---SNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKN 60
           RL E+C + MGQSP  + I     N D +EFHQGKI F++KYI ++N  T ++ KIA KN
Sbjct: 12  RLEEVCEVVMGQSPNGNCIFDKDKNKDLIEFHQGKIAFSDKYIDESNFVTSDVKKIAKKN 71

Query: 61  AILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTF 120
           +++LCVRAPVG VNIT ++I IGRGLCS+    +I + F F++L   K  F   STG+TF
Sbjct: 72  SVVLCVRAPVGEVNITTKDIAIGRGLCSL-NGVKINNNFLFFYLLTLKKYFNDNSTGSTF 130

Query: 121 QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
           +AI++ +IK   IPLPP  EQ+RIV  ++  F+++D
Sbjct: 131 KAINVRVIKETKIPLPPLKEQERIVGILDFAFSKID 166
>gi|40467|emb|CAA35604.1| HsdS polypeptide, part of CfrA family [Citrobacter freundii]
          Length = 578

 Score =  116 bits (291), Expect = 3e-25,   Method: Composition-based stats.
 Identities = 64/157 (40%), Positives = 92/157 (58%), Gaps = 1/157 (0%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNA 61
           W RLGE   I MGQS  S + + + +G+ F QGK  F +KY   A  +  + TK+A KN 
Sbjct: 107 WVRLGEAFYIEMGQSXSSQYYNQSEEGIPFFQGKADFGKKY-PTARYWCTSPTKLAQKND 165

Query: 62  ILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQ 121
           +LL VRAPVG  N++    CIGRGL ++        E+  Y L+  +   E+ +TGTTF 
Sbjct: 166 VLLSVRAPVGPTNLSPYHCCIGRGLAAIRCLSDAPHEYLLYILKASQRRLEELATGTTFV 225

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           A+S   I+ +L+P+PP +EQ RIV  I+ L +  D +
Sbjct: 226 AVSKTDIEPLLMPIPPLNEQIRIVDTIDRLMSLCDQL 262
>gi|20807979|ref|NP_623150.1| Restriction endonuclease S subunits [Thermoanaerobacter
           tengcongensis MB4]
 gi|20516553|gb|AAM24754.1| Restriction endonuclease S subunits [Thermoanaerobacter
           tengcongensis MB4]
          Length = 398

 Score =  114 bits (286), Expect = 1e-24,   Method: Composition-based stats.
 Identities = 57/157 (36%), Positives = 91/157 (57%), Gaps = 1/157 (0%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNA 61
           W RLGE+C I MGQSP SS  +   +G+ F QGK  F + +      +     K+A    
Sbjct: 216 WVRLGEVCDIIMGQSPPSSTYNFEGNGLPFFQGKADFGDLH-PTPRIWCSAPQKVARPGD 274

Query: 62  ILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQ 121
           +L+ VRAPVG  N+     CIGRGL ++ P+  ++  +  Y+L   +    +   G+TF 
Sbjct: 275 VLISVRAPVGSTNVANLACCIGRGLAALRPRDSLERFWLLYYLHYLEPELSKMGAGSTFN 334

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           AI+ + ++N+ IPLPP  EQ+RIVA ++++  Q+ ++
Sbjct: 335 AITKKDLQNVFIPLPPLEEQRRIVAYLDQIQQQVAAL 371
>gi|68550910|ref|ZP_00590345.1| Restriction modification system DNA specificity domain [Pelodictyon
           phaeoclathratiforme BU-1]
 gi|68242209|gb|EAN24435.1| Restriction modification system DNA specificity domain [Pelodictyon
           phaeoclathratiforme BU-1]
          Length = 392

 Score =  105 bits (261), Expect = 1e-21,   Method: Composition-based stats.
 Identities = 57/150 (38%), Positives = 88/150 (58%), Gaps = 3/150 (2%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L ++ +I  GQSPESS  ++ +DG+ F QGK  F +K+ +          K A    IL+
Sbjct: 6   LQQVTTIIAGQSPESSTYNSIADGLPFFQGKADFQDKFPKVRIWCNSAKRKEADPGDILM 65

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAIS 124
            VRAPVG VNI  ++  IGRGL ++ P   + + F +Y+L+C +       TG+TFQAI+
Sbjct: 66  SVRAPVGSVNICNQKCIIGRGLSAIRPDANLNNYFLYYYLKCNEKNVASLGTGSTFQAIT 125

Query: 125 IELIKNILIPLPPSSEQKR---IVAKIEEL 151
              +K + +PLPP  +Q R   +++K+E L
Sbjct: 126 QTTLKRLDVPLPPLDDQIRSATLLSKVENL 155
>gi|19881243|gb|AAM00852.1|AF486551_3 HsdS [Campylobacter jejuni]
          Length = 380

 Score =  101 bits (252), Expect = 1e-20,   Method: Composition-based stats.
 Identities = 63/159 (39%), Positives = 86/159 (54%), Gaps = 9/159 (5%)

Query: 3   CRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAI 62
           C+LG+I  ITMGQSP+S F + ++ GM F QG   F  KY    + Y     KIA K  I
Sbjct: 7   CKLGDIAEITMGQSPKSEFYNFDNIGMPFLQGNRTFGRKY-PYFDTYCTEYKKIAKKGEI 65

Query: 63  LLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQA 122
           L  VRAPVG +N    +ICIGRGLCS+  K   ++EF +Y L   +       +G+ F +
Sbjct: 66  LFSVRAPVGDINFANNDICIGRGLCSMNAK-NGENEFLYYLLHNLRSVIINNESGSVFGS 124

Query: 123 ISIELIKNILIPLPPSSEQKRIVA-------KIEELFAQ 154
           ++   ++ I I LP   EQ++I         KIE L  Q
Sbjct: 125 VNKNDLQTIEILLPLLEEQRQIATILSSIDDKIELLHEQ 163
>gi|15792859|ref|NP_282682.1| putative type I restriction enzyme S protein [Campylobacter jejuni
           subsp. jejuni NCTC 11168]
 gi|86149261|ref|ZP_01067492.1| HsdS [Campylobacter jejuni subsp. jejuni CF93-6]
 gi|88596768|ref|ZP_01100005.1| HsdS [Campylobacter jejuni subsp. jejuni 84-25]
 gi|121612440|ref|YP_001001191.1| HsdS [Campylobacter jejuni subsp. jejuni 81-176]
 gi|19881204|gb|AAM00820.1|AF486544_3 HsdS [Campylobacter jejuni]
 gi|19881210|gb|AAM00825.1|AF486545_3 HsdS [Campylobacter jejuni]
 gi|19881237|gb|AAM00847.1|AF486550_3 HsdS [Campylobacter jejuni]
 gi|19881285|gb|AAM00887.1|AF486558_3 HsdS [Campylobacter jejuni]
 gi|19881289|gb|AAM00890.1|AF486559_1 HsdS [Campylobacter jejuni]
 gi|19881291|gb|AAM00891.1|AF486560_1 HsdS [Campylobacter jejuni]
 gi|19881293|gb|AAM00892.1|AF486561_1 HsdS [Campylobacter jejuni]
 gi|19881295|gb|AAM00893.1|AF486562_1 HsdS [Campylobacter jejuni]
 gi|19881297|gb|AAM00894.1|AF486563_1 HsdS [Campylobacter jejuni]
 gi|19881301|gb|AAM00896.1|AF486565_1 HsdS [Campylobacter jejuni]
 gi|19881303|gb|AAM00897.1|AF486566_1 HsdS [Campylobacter jejuni]
 gi|19881306|gb|AAM00898.1|AF486568_1 HsdS [Campylobacter jejuni]
 gi|19881308|gb|AAM00899.1|AF486569_1 HsdS [Campylobacter jejuni]
 gi|85840043|gb|EAQ57301.1| HsdS [Campylobacter jejuni subsp. jejuni CF93-6]
 gi|87249550|gb|EAQ72509.1| HsdS [Campylobacter jejuni subsp. jejuni 81-176]
 gi|88191609|gb|EAQ95581.1| HsdS [Campylobacter jejuni subsp. jejuni 84-25]
 gi|112360850|emb|CAL35651.1| putative type I restriction enzyme S protein [Campylobacter jejuni
           subsp. jejuni NCTC 11168]
          Length = 380

 Score =  101 bits (251), Expect = 1e-20,   Method: Composition-based stats.
 Identities = 63/159 (39%), Positives = 86/159 (54%), Gaps = 9/159 (5%)

Query: 3   CRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAI 62
           C+LG+I  ITMGQSP+S F + ++ GM F QG   F  KY    + Y     KIA K  I
Sbjct: 7   CKLGDIAEITMGQSPKSEFYNFDNIGMPFLQGNRTFGRKY-PYFDTYCTEYKKIAKKGEI 65

Query: 63  LLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQA 122
           L  VRAPVG +N    +ICIGRGLCS+  K   ++EF +Y L   +       +G+ F +
Sbjct: 66  LFSVRAPVGDINFANNDICIGRGLCSMNAK-NGENEFLYYLLHNLRSVIINNESGSVFGS 124

Query: 123 ISIELIKNILIPLPPSSEQKRIVA-------KIEELFAQ 154
           ++   ++ I I LP   EQ++I         KIE L  Q
Sbjct: 125 VNKNDLQTIEILLPLLEEQRQIATILSSIDDKIELLHEQ 163
>gi|126661170|ref|ZP_01732247.1| type I restriction-modification enzyme, S subunit, putative
           [Cyanothece sp. CCY0110]
 gi|126617543|gb|EAZ88335.1| type I restriction-modification enzyme, S subunit, putative
           [Cyanothece sp. CCY0110]
          Length = 191

 Score = 98.2 bits (243), Expect = 2e-19,   Method: Composition-based stats.
 Identities = 56/143 (39%), Positives = 81/143 (56%), Gaps = 1/143 (0%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +LG+I  I  GQSP+S F + N  G+ F+QGKI F   Y+++   +T  ITK + K+ IL
Sbjct: 49  KLGKITEIISGQSPQSKFYNKNQQGLPFYQGKIEFGNMYLKEPKTWTTQITKESIKDDIL 108

Query: 64  LCVRAPVGVVNITE-REICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQA 122
           + VRAPVG +NI    +ICIGRGL ++  K       + Y+           + G  F +
Sbjct: 109 MSVRAPVGSLNINRFDKICIGRGLAAIRSKAENVFIKYIYYFLLFNPELIVGTEGLIFSS 168

Query: 123 ISIELIKNILIPLPPSSEQKRIV 145
           IS + I  I IPLPP   Q++I+
Sbjct: 169 ISRDQISKISIPLPPKEVQEQII 191
>gi|84489266|ref|YP_447498.1| hypothetical protein Msp_0455 [Methanosphaera stadtmanae DSM 3091]
 gi|84372585|gb|ABC56855.1| conserved hypothetical protein [Methanosphaera stadtmanae DSM 3091]
          Length = 180

 Score = 93.2 bits (230), Expect = 4e-18,   Method: Composition-based stats.
 Identities = 57/145 (39%), Positives = 87/145 (60%), Gaps = 5/145 (3%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +L +I  ITMGQSP S + + N +     QG       Y+     YT  ITKIA K  IL
Sbjct: 23  KLSQINKITMGQSPSSKYYTKNQNDTILVQGNQDIANNYV-IPRIYTSKITKIAKKGEIL 81

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQ--CQKDTFEQKSTGTTFQ 121
           L VRAPVG + IT+ ++CIGRG+CS+ P   I + F F++L     K+ + +   G+TF+
Sbjct: 82  LTVRAPVGDIVITQYDVCIGRGVCSIKPS--ISTGFMFFYLAKLNSKNQWNKYIQGSTFE 139

Query: 122 AISIELIKNILIPLPPSSEQKRIVA 146
           +I+ + IK++ I +P SS+++  +A
Sbjct: 140 SINSKDIKSMKIKIPKSSKEQEKIA 164
>gi|15669726|ref|NP_248539.1| type I restriction-modification enzyme, S subunit, putative
           [Methanocaldococcus jannaschii DSM 2661]
 gi|2496187|sp|Q58926|Y1531_METJA Uncharacterized protein MJ1531
 gi|1592162|gb|AAB99552.1| type I restriction-modification enzyme, S subunit, putative
           [Methanocaldococcus jannaschii DSM 2661]
          Length = 425

 Score = 91.3 bits (225), Expect = 2e-17,   Method: Composition-based stats.
 Identities = 56/142 (39%), Positives = 78/142 (54%), Gaps = 2/142 (1%)

Query: 11  ITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPV 70
           I MGQSP SS  +   +G+ F QGK  F   Y      YT    K+     IL+ VRAPV
Sbjct: 251 IIMGQSPPSSSYNKEGEGVPFLQGKAEFGNIYPNPVL-YTNKPLKVVDDEDILISVRAPV 309

Query: 71  GVVNITEREICIGRGLCSVYP-KYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIK 129
           G VNI   ++CIGRGL  +   K ++ + F FY+L   K   E    G  F+AI+ + ++
Sbjct: 310 GDVNIAPFKLCIGRGLAGIKSNKEKVDNFFVFYYLSYIKPKIEYLGGGAVFKAITKKDLE 369

Query: 130 NILIPLPPSSEQKRIVAKIEEL 151
           +I IPLPP  EQK I  +++ +
Sbjct: 370 SIKIPLPPLEEQKAIAKRLKAI 391
>gi|150388684|ref|YP_001318733.1| restriction modification system DNA specificity domain
           [Alkaliphilus metalliredigens QYMF]
 gi|149948546|gb|ABR47074.1| restriction modification system DNA specificity domain
           [Alkaliphilus metalliredigens QYMF]
          Length = 467

 Score = 91.3 bits (225), Expect = 2e-17,   Method: Composition-based stats.
 Identities = 53/177 (29%), Positives = 89/177 (50%), Gaps = 21/177 (11%)

Query: 1   MWCRLGEICSITMGQSPESSFISNNSDG---------------MEFHQGKIHFTEKYIQK 45
           +W RLG + +I  G +P S  I    +G               +    GK + TE  ++K
Sbjct: 33  VWTRLGNVTTIIGGGTPPSRVIEYYENGSIPWISPVDLSGYTDIYISHGKKNITELGLKK 92

Query: 46  ANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQ 105
           ++       ++ P+N +LL  RAP+G V I + E+C  +G  S  P      ++ +++L+
Sbjct: 93  SS------ARLLPENTVLLSSRAPIGYVAIADNELCTNQGFKSFLPSPCYLPKYLYFYLK 146

Query: 106 CQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
             K   E  ++GTTF  +S      +  PLPP +EQ+RIV +IE LF +L+   A++
Sbjct: 147 SSKKLLEAYASGTTFLELSGRKAAIVEFPLPPLAEQQRIVDRIESLFEKLNQAKALI 203
>gi|134097471|ref|YP_001103132.1| putative restriction modification system DNA specificity domain
           [Saccharopolyspora erythraea NRRL 2338]
 gi|133910094|emb|CAM00207.1| putative restriction modification system DNA specificity domain
           [Saccharopolyspora erythraea NRRL 2338]
          Length = 411

 Score = 87.0 bits (214), Expect = 3e-16,   Method: Composition-based stats.
 Identities = 54/158 (34%), Positives = 86/158 (54%), Gaps = 4/158 (2%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKY-IQKANNYTFNITKIAPKNAIL 63
           L  +  I  GQ+P S+ +S+  DG+ F QG   F   + + +    T   +KIA    +L
Sbjct: 6   LKYVAYIRPGQAPPSTEVSDLIDGLPFLQGNAEFQAAHPVPRLQCDT--ASKIAKCGDVL 63

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAI 123
           L VRAPVG +NI +RE  IGRGLCSV       + F ++WL    +  +  STGTT++A+
Sbjct: 64  LSVRAPVGALNIADREYGIGRGLCSV-SATGCDARFLWWWLHSAGERLDAVSTGTTYRAV 122

Query: 124 SIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAV 161
           + E +  +  P     EQ+RI   ++    ++D ++A+
Sbjct: 123 TGEDVGMLPFPRVSLEEQRRIADFLDAETTRIDKLSAL 160
>gi|121541757|ref|ZP_01673502.1| Restriction endonuclease S subunits-like [Candidatus Desulfococcus
           oleovorans Hxd3]
 gi|121518303|gb|EAX55166.1| Restriction endonuclease S subunits-like [Candidatus Desulfococcus
           oleovorans Hxd3]
          Length = 385

 Score = 85.5 bits (210), Expect = 1e-15,   Method: Composition-based stats.
 Identities = 51/162 (31%), Positives = 80/162 (49%), Gaps = 11/162 (6%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGK-----IHFTEKYIQKANNYTFNITKIAPK 59
           L ++C + MGQ+P+    + N+ G     G      IH + K       +T + TK    
Sbjct: 14  LNQLCLVVMGQAPKGDTYNENTLGTPLIAGAADLGLIHPSPK------KWTTSPTKTGKA 67

Query: 60  NAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTT 119
             I+LCVRA +G +N  + + C GRG+C +        EF ++WL   KD       G T
Sbjct: 68  GDIILCVRATIGDLNWADSKYCYGRGVCGLRIIEGHDPEFLWFWLMACKDHLLSLGRGAT 127

Query: 120 FQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAV 161
           F+ IS   I N+ +P     EQ+RIVA+I+E   +++ I  +
Sbjct: 128 FKQISKTDIANLPVPALAVDEQRRIVARIKECMERVEEIEGL 169
>gi|21228395|ref|NP_634317.1| type I restriction-modification system specificity subunit
           [Methanosarcina mazei Go1]
 gi|20906868|gb|AAM31989.1| type I restriction-modification system specificity subunit
           [Methanosarcina mazei Go1]
          Length = 406

 Score = 85.1 bits (209), Expect = 1e-15,   Method: Composition-based stats.
 Identities = 51/156 (32%), Positives = 91/156 (58%), Gaps = 8/156 (5%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIH-FTEKYIQKANNYTFNITKIAPKNAIL 63
           LG++  ITMGQSP S+  ++  +G+   QG    F  K I++   +T  IT+      I+
Sbjct: 231 LGDVAEITMGQSPSSANYNSKGEGLPLIQGNADIFNRKTIKRV--FTTEITRRGKCGDII 288

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAI 123
           + VRAPVG V+  E +IC+GRG+C++    R  + F ++ L  ++ T+ + S G+TF ++
Sbjct: 289 MSVRAPVGEVSRAEFDICLGRGVCAI----RYSNNFLYHTLISKESTWAKLSKGSTFDSV 344

Query: 124 SIELIKNILIPLPP-SSEQKRIVAKIEELFAQLDSI 158
           +   +K   I LP  S+EQ+ I   + ++ A++ ++
Sbjct: 345 NSADVKAFDIELPTDSAEQEAIATILSDMDAEITAL 380
>gi|52082594|ref|YP_081385.1| HsdS [Bacillus licheniformis ATCC 14580]
 gi|52787992|ref|YP_093821.1| hypothetical protein BLi04316 [Bacillus licheniformis ATCC 14580]
 gi|52005805|gb|AAU25747.1| HsdS [Bacillus licheniformis ATCC 14580]
 gi|52350494|gb|AAU43128.1| putative protein [Bacillus licheniformis DSM 13]
          Length = 387

 Score = 85.1 bits (209), Expect = 1e-15,   Method: Composition-based stats.
 Identities = 51/144 (35%), Positives = 81/144 (56%), Gaps = 3/144 (2%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L +I  ITMGQSP  +  ++  DG+    G   FT KY      +T   TK+     ILL
Sbjct: 14  LSDIADITMGQSPPGNSYNDIKDGIGLINGPTEFTNKY-PVVKQWTSKPTKLCKAGDILL 72

Query: 65  CVR-APVGVVNITEREICIGRGLCSVYPKY-RIQSEFWFYWLQCQKDTFEQKSTGTTFQA 122
           CVR +  G +NI + E CIGRG+ S+  K  + ++ F +Y L  + +   QK+ G+TF  
Sbjct: 73  CVRGSSTGRMNIADDEYCIGRGVASIRAKKDKAETSFIYYTLNYKVNQLLQKTAGSTFPN 132

Query: 123 ISIELIKNILIPLPPSSEQKRIVA 146
           +S   IK++++ +P  +EQ++I +
Sbjct: 133 LSSNEIKDMIVGIPLFAEQQKIAS 156
>gi|91772548|ref|YP_565240.1| restriction modification system DNA specificity subunit
           [Methanococcoides burtonii DSM 6242]
 gi|91711563|gb|ABE51490.1| restriction modification system DNA specificity subunit
           [Methanococcoides burtonii DSM 6242]
          Length = 511

 Score = 84.7 bits (208), Expect = 2e-15,   Method: Composition-based stats.
 Identities = 42/131 (32%), Positives = 76/131 (58%), Gaps = 6/131 (4%)

Query: 38  FTEKYIQKANNYTFNI------TKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYP 91
           ++EKYI K      ++       ++ PK ++L   RAP+G + I   E+C  +G  ++ P
Sbjct: 55  YSEKYIYKGRKSITHLGLKNSSARLIPKGSVLFSSRAPIGYIAIAGNELCTNQGFKTLIP 114

Query: 92  KYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEEL 151
              +  +F +Y+L+  K   E +++GTTF+ +S +    + + +PP  EQ+ IV+KIE+L
Sbjct: 115 SEALNRDFLYYYLKSIKQLAEGRASGTTFKELSGKAFAELPLCVPPLPEQRAIVSKIEQL 174

Query: 152 FAQLDSITAVL 162
           F++LD+  A L
Sbjct: 175 FSELDNGIANL 185
>gi|42779918|ref|NP_977165.1| type I restriction-modification enzyme, S subunit, putative
           [Bacillus cereus ATCC 10987]
 gi|42735836|gb|AAS39773.1| type I restriction-modification enzyme, S subunit, putative
           [Bacillus cereus ATCC 10987]
          Length = 476

 Score = 83.2 bits (204), Expect = 4e-15,   Method: Composition-based stats.
 Identities = 50/177 (28%), Positives = 85/177 (48%), Gaps = 21/177 (11%)

Query: 1   MWCRLGEICSITMGQSPESSFISNNSDG---------------MEFHQGKIHFTEKYIQK 45
           +W   G I  +  G +P+S       DG               M   +GK + TE  + K
Sbjct: 26  IWTWTGAIAEVISGGTPKSKVEEYYKDGTISWITPADLSGYQDMYISKGKRNITELGLNK 85

Query: 46  ANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQ 105
           ++       K+ P N +LL  RAP+G V I  +++C  +G  S  P      ++ +++L+
Sbjct: 86  SS------AKMLPINTVLLSSRAPIGYVAIAAKDLCTNQGFKSFAPSNAYYPKYLYWYLK 139

Query: 106 CQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
             K   E  ++G+TF+ +S    K I IPLPP +EQKR+  K+E L  +++    ++
Sbjct: 140 FSKYYMESMASGSTFKELSSNKSKEIPIPLPPINEQKRVSEKVERLLNKVEEAKTLI 196
>gi|148544100|ref|YP_001271470.1| restriction modification system DNA specificity domain
           [Lactobacillus reuteri F275]
 gi|148531134|gb|ABQ83133.1| restriction modification system DNA specificity domain
           [Lactobacillus reuteri F275]
          Length = 195

 Score = 81.6 bits (200), Expect = 1e-14,   Method: Composition-based stats.
 Identities = 47/154 (30%), Positives = 86/154 (55%), Gaps = 2/154 (1%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +L ++  I MGQSP+S F + N +G  F QG   F E Y Q  + +T +  + A    IL
Sbjct: 7   QLKDVADIVMGQSPKSVFYNTNGNGTPFLQGVRTFGENYPQ-IDTWTTSYNRKAKSGEIL 65

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAI 123
             VRAPVG VN   +++ +GRGL ++  K     E+ +Y  +      +  +TGT F +I
Sbjct: 66  FSVRAPVGRVNWANQDLAVGRGLAALKIKSGYSKEYLYYLFKKIGGQLDSLATGTVFTSI 125

Query: 124 SIELIKNILIPLPPS-SEQKRIVAKIEELFAQLD 156
           + + ++ I + +P + S+Q++I   ++++  +++
Sbjct: 126 NKKELEAIELKIPVNLSDQEKIADYLQKIDQEIE 159
>gi|31983512|ref|NP_858124.1| putative type i restriction enzyme hsds subunit [Lactobacillus
           delbrueckii subsp. lactis]
 gi|18077746|emb|CAD13349.1| putative Type I restriction enzyme hsdS subunit [Lactobacillus
           delbrueckii subsp. lactis]
          Length = 396

 Score = 81.6 bits (200), Expect = 1e-14,   Method: Composition-based stats.
 Identities = 51/144 (35%), Positives = 78/144 (54%), Gaps = 8/144 (5%)

Query: 3   CRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAI 62
           C+LG++  ITMGQSP S   ++N       QG     +  +     +T  ITKIA K  +
Sbjct: 23  CKLGDVAKITMGQSPNSKNYTDNPKDHILVQGNADMKDGQVH-PRIWTTEITKIADKGDL 81

Query: 63  LLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTF 120
           +L VRAPVG +  T  ++ IGRG+ ++       +EF F  L+  K    + + STG+TF
Sbjct: 82  ILSVRAPVGDIGKTSYDVVIGRGVAAIK-----GNEFIFQLLKRMKTVGYWTKYSTGSTF 136

Query: 121 QAISIELIKNILIPLPPSSEQKRI 144
           ++I+   I N +I LP   EQ ++
Sbjct: 137 ESINSLEINNAVINLPKEHEQNKV 160
>gi|55820774|ref|YP_139216.1| type I restriction-modification system specificty subunit
           [Streptococcus thermophilus LMG 18311]
 gi|55822677|ref|YP_141118.1| type I restriction-modification system specificty subunit
           [Streptococcus thermophilus CNRZ1066]
 gi|55736759|gb|AAV60401.1| type I restriction-modification system specificty subunit
           [Streptococcus thermophilus LMG 18311]
 gi|55738662|gb|AAV62303.1| type I restriction-modification system specificty subunit
           [Streptococcus thermophilus CNRZ1066]
          Length = 406

 Score = 80.9 bits (198), Expect = 2e-14,   Method: Composition-based stats.
 Identities = 56/156 (35%), Positives = 87/156 (55%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L E+ +ITMGQSP+S   ++N +     QG     +K +     +T  +TK+A    I+L
Sbjct: 237 LSEVTNITMGQSPKSENYTDNPNDYILVQGNADIKDKQVV-PRLWTTEVTKMAEIGDIIL 295

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG +  T+  + IGRG+ ++       ++F FY L+  K T  + + STG+TF++
Sbjct: 296 TVRAPVGDIGKTDYNVVIGRGVAAIK-----GNDFIFYTLEKMKMTGFWNKFSTGSTFES 350

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   IK  +I +P   EQK I A     F+ LD++
Sbjct: 351 ISSNDIKEAIIQIPTLEEQKAIGA----YFSNLDNL 382
>gi|89210922|ref|ZP_01189304.1| Restriction modification system DNA specificity domain
           [Halothermothrix orenii H 168]
 gi|89159468|gb|EAR79134.1| Restriction modification system DNA specificity domain
           [Halothermothrix orenii H 168]
          Length = 565

 Score = 80.5 bits (197), Expect = 3e-14,   Method: Composition-based stats.
 Identities = 51/173 (29%), Positives = 82/173 (47%), Gaps = 15/173 (8%)

Query: 2   WCRLGEICSITMGQSPES---SFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIA- 57
           W RLG I  I  G +P++   ++  N +              KYI +      NITK+  
Sbjct: 88  WVRLGNIGRIVGGGTPKTKVHAYWENGNIAWLTPADLNGLKSKYISRGRR---NITKLGL 144

Query: 58  --------PKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKD 109
                   PK ++L   RAP+G V I + ++   +G  S  P     +++ +Y+L     
Sbjct: 145 QNSSAKLLPKGSVLFSSRAPIGYVAIAQNDLATNQGFKSCVPYIMDMNQYIYYFLMYDAK 204

Query: 110 TFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
                ++GTTF+ +S + + N + PLPP +EQKRIV K++EL    D +   L
Sbjct: 205 RINDNASGTTFKEVSGKEVANFIFPLPPLNEQKRIVNKLDELMTFCDQLEVSL 257
>gi|111222733|ref|YP_713527.1| Type I restriction modification enzyme protein S [Frankia alni
           ACN14a]
 gi|111150265|emb|CAJ61962.1| Type I restriction modification enzyme protein S [Frankia alni
           ACN14a]
          Length = 399

 Score = 79.0 bits (193), Expect = 9e-14,   Method: Composition-based stats.
 Identities = 52/173 (30%), Positives = 81/173 (46%), Gaps = 29/173 (16%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHF----------------TEKYIQKANN 48
           LGE C I  G +P+++         E+  G+I +                T + I +A  
Sbjct: 9   LGEFCEIISGATPKTA-------SEEYWGGEIPWATPRDLGSLNSKFLASTSRAITEAGL 61

Query: 49  YTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKY-RIQSEFWFYWLQCQ 107
            +   T + P  ++LL  RAP+G V I  R +   +G  S+ P   R    + ++WL+CQ
Sbjct: 62  RSC-ATHVLPAGSVLLTSRAPIGSVAINARPMATNQGFKSLVPDTSRALPGYLYHWLRCQ 120

Query: 108 KDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITA 160
           +   +    G TF+ +S      I +PLPP SEQKR    IE++  Q D+I A
Sbjct: 121 RSRLQSLGNGATFKELSKSATARIAVPLPPLSEQKR----IEQMLDQADTIRA 169
>gi|77415025|ref|ZP_00791100.1| Type I restriction modification DNA specificity domain protein
           [Streptococcus agalactiae 515]
 gi|77158925|gb|EAO70161.1| Type I restriction modification DNA specificity domain protein
           [Streptococcus agalactiae 515]
          Length = 385

 Score = 78.6 bits (192), Expect = 1e-13,   Method: Composition-based stats.
 Identities = 55/170 (32%), Positives = 88/170 (51%), Gaps = 22/170 (12%)

Query: 2   WCRLGEICSITMGQSPESS----------FISNNSDGME-----FHQGKIHFTEKYIQKA 46
           W +L  I SIT G +P+SS          +I+    G +     F +     TE  +QK+
Sbjct: 90  WVKLRNIGSITSGGTPKSSEPSYYGGNITWITPADMGKQQNNKFFAKSSKKITELGLQKS 149

Query: 47  NNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQC 106
           +       ++  KN+I+   RAP+G +NI   +    +G  S+ P   +   F ++ LQ 
Sbjct: 150 S------AQLISKNSIVYSSRAPIGHINIVTEDYTTNQGCKSITP-LLVDLIFLYWLLQF 202

Query: 107 QKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
           +      +S+GTTF+ IS     + L+PLPP +EQKRIVA+IE+  A++D
Sbjct: 203 RTKDIILRSSGTTFKEISASGFGDTLLPLPPLAEQKRIVAQIEKALAKVD 252
>gi|119491619|ref|ZP_01623491.1| hypothetical protein L8106_03529 [Lyngbya sp. PCC 8106]
 gi|119453348|gb|EAW34512.1| hypothetical protein L8106_03529 [Lyngbya sp. PCC 8106]
          Length = 433

 Score = 77.4 bits (189), Expect = 2e-13,   Method: Composition-based stats.
 Identities = 50/158 (31%), Positives = 82/158 (51%), Gaps = 6/158 (3%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTF--NITKIAPKNA 61
           +L  I  + MGQSP+S+  +    G+ F QG   F    I   N      +  K A K+ 
Sbjct: 29  KLKFIADLIMGQSPDSTDYNYEEIGVPFLQGTAEFG---IINPNPRLSCESAKKYARKDD 85

Query: 62  ILLCVRAPVGVVNITEREICIGRGLCSVYPKYRI-QSEFWFYWLQCQKDTFEQKSTGTTF 120
           +LL VRAPVG +N+ ++   IGRGLC++ PK  +    F  Y+L+  K      +TG+ +
Sbjct: 86  LLLSVRAPVGEINVADQVYGIGRGLCAIRPKINVFNKTFTRYFLEIGKVELVSGATGSIY 145

Query: 121 QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
            A+++  + N+    PP  EQK I   ++    ++D++
Sbjct: 146 DAVTVNQVANLQCLTPPLKEQKLIATFLDRETTRIDTL 183
>gi|146319105|ref|YP_001198817.1| type I restriction-modification system, S subunit [Streptococcus
           suis 05ZYH33]
 gi|145689911|gb|ABP90417.1| type I restriction-modification system, S subunit [Streptococcus
           suis 05ZYH33]
          Length = 522

 Score = 75.9 bits (185), Expect = 7e-13,   Method: Composition-based stats.
 Identities = 54/170 (31%), Positives = 84/170 (49%), Gaps = 22/170 (12%)

Query: 2   WCRLGEICSITMGQSPESS----------FISNNSDGME-----FHQGKIHFTEKYIQKA 46
           W RL  +  IT G +P+SS          +I+    G +     F       TE  +QK+
Sbjct: 90  WVRLRNLGVITSGGTPKSSESTYYDGNITWITPADMGKQQNNKLFATSSKKITELGVQKS 149

Query: 47  NNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQC 106
           +       ++  KN+I+   RAP+G +NI   +    +G  SV P   +  +F ++ LQ 
Sbjct: 150 S------AQLISKNSIVYSSRAPIGHINIVNYDFTTNQGCKSVTP-ILVNLDFMYWILQF 202

Query: 107 QKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
           +      +S+GTTF+ IS     + L+PLPP +EQKRIVA IE    Q++
Sbjct: 203 RTKDIILRSSGTTFKEISASGFGDTLLPLPPLAEQKRIVAHIERALEQVE 252
>gi|154173663|ref|YP_001408732.1| type I restriction-modification system S subunit [Campylobacter
           curvus 525.92]
 gi|153792995|gb|EAT99440.2| type I restriction-modification system S subunit [Campylobacter
           curvus 525.92]
          Length = 323

 Score = 75.5 bits (184), Expect = 9e-13,   Method: Composition-based stats.
 Identities = 61/169 (36%), Positives = 82/169 (48%), Gaps = 20/169 (11%)

Query: 2   WCRLGEICSITMG----QSPESSFISNNSDGMEF-------HQGKIHFTEKYIQKANNYT 50
           W RLGEIC I  G    Q+ + +  +N  DG  +         G I +        N   
Sbjct: 158 WVRLGEICQIYTGDSINQTQKLTKYTNLEDGRCYIATKDVDFDGSIDYENGVKIPFNESR 217

Query: 51  FNITKIAPKNAILLCVRAPVG--VVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQK 108
           F   KIAPKN++LLCV        +   + ++C G  LC   P   I+ +F +Y+LQ Q 
Sbjct: 218 F---KIAPKNSVLLCVEGGSAGKKIGYLDCDVCFGNKLCCFNP-LLIEPKFIYYYLQSQI 273

Query: 109 --DTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQL 155
              +F QK +G     IS+  IK I+I +PP  EQKRIV KIE L   L
Sbjct: 274 FIYSFMQKMSGI-ISGISLNSIKTIVIAIPPLPEQKRIVEKIELLLPLL 321
>gi|156864174|gb|EDO57605.1| hypothetical protein CLOL250_01695 [Clostridium sp. L2-50]
          Length = 372

 Score = 75.5 bits (184), Expect = 9e-13,   Method: Composition-based stats.
 Identities = 50/150 (33%), Positives = 79/150 (52%), Gaps = 8/150 (5%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L +I  IT+GQSP  S  + +  G  F QG+  F  ++      YT    ++A  N  L+
Sbjct: 190 LSDIADITIGQSPSGSSYNEDGTGTIFFQGRAEFGFRF-PSVRLYTTEPKRMARSNDTLM 248

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQS--EFWFYWLQCQKDTFEQKSTGTTFQA 122
            VRAPVG +N+   + CIGRGL +++ K   QS   +  + L+ Q D F  +  GT F +
Sbjct: 249 SVRAPVGDLNVAHTDCCIGRGLAAIHSKSNHQSFVLYTMFSLKKQLDVFNGE--GTVFGS 306

Query: 123 ISIELIKNILIPLPPS---SEQKRIVAKIE 149
           I+   + ++ I +P      E +RIVA ++
Sbjct: 307 INRNSLNDMPILIPSDDILDEFERIVAPMD 336
>gi|156864459|gb|EDO57890.1| hypothetical protein CLOL250_01430 [Clostridium sp. L2-50]
          Length = 215

 Score = 75.5 bits (184), Expect = 9e-13,   Method: Composition-based stats.
 Identities = 50/150 (33%), Positives = 79/150 (52%), Gaps = 8/150 (5%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L +I  IT+GQSP  S  + +  G  F QG+  F  ++      YT    ++A  N  L+
Sbjct: 33  LSDIADITIGQSPSGSSYNEDGTGTIFFQGRAEFGFRF-PSVRLYTTEPKRMARSNDTLM 91

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQS--EFWFYWLQCQKDTFEQKSTGTTFQA 122
            VRAPVG +N+   + CIGRGL +++ K   QS   +  + L+ Q D F  +  GT F +
Sbjct: 92  SVRAPVGDLNVAHTDCCIGRGLAAIHSKSNHQSFVLYTMFSLKKQLDVFNGE--GTVFGS 149

Query: 123 ISIELIKNILIPLPPS---SEQKRIVAKIE 149
           I+   + ++ I +P      E +RIVA ++
Sbjct: 150 INRNSLNDMPILIPSDDILDEFERIVAPMD 179
>gi|146321309|ref|YP_001201020.1| type I restriction-modification system, S subunit [Streptococcus
           suis 98HAH33]
 gi|145692115|gb|ABP92620.1| type I restriction-modification system, S subunit [Streptococcus
           suis 98HAH33]
          Length = 284

 Score = 75.1 bits (183), Expect = 1e-12,   Method: Composition-based stats.
 Identities = 54/170 (31%), Positives = 84/170 (49%), Gaps = 22/170 (12%)

Query: 2   WCRLGEICSITMGQSPESS----------FISNNSDGME-----FHQGKIHFTEKYIQKA 46
           W RL  +  IT G +P+SS          +I+    G +     F       TE  +QK+
Sbjct: 90  WVRLRNLGVITSGGTPKSSESTYYDGNITWITPADMGKQQNDKVFATSSKKITELGVQKS 149

Query: 47  NNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQC 106
           +       ++  KN+I+   RAP+G +NI   +    +G  SV P   +  +F ++ LQ 
Sbjct: 150 S------AQLISKNSIVYSSRAPIGHINIVNYDFTTNQGCKSVTP-ILVNLDFMYWILQF 202

Query: 107 QKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
           +      +S+GTTF+ IS     + L+PLPP +EQKRIVA IE    Q++
Sbjct: 203 RTKDIILRSSGTTFKEISASGFGDTLLPLPPLAEQKRIVAHIERALEQVE 252
>gi|60681333|ref|YP_211477.1| putative type IC restriction-modification system specificity
           subunit, partial [Bacteroides fragilis NCTC 9343]
 gi|60492767|emb|CAH07541.1| putative type IC restriction-modification system specificity
           subunit, partial [Bacteroides fragilis NCTC 9343]
          Length = 376

 Score = 75.1 bits (183), Expect = 1e-12,   Method: Composition-based stats.
 Identities = 54/156 (34%), Positives = 79/156 (50%), Gaps = 6/156 (3%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNA 61
           +  +GE+C + MGQSP S   +   +G+   QG +   E  +     +T +ITK      
Sbjct: 202 YTSIGEMCDVVMGQSPSSVAYNYTKNGLPLIQGNLDIFEG-VTSPRMWTSDITKQCDIGD 260

Query: 62  ILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQ 121
           I+L VRAPVG V  +    C+GRG+C++  K    SE+ + +L   K  +     G+TF 
Sbjct: 261 IILTVRAPVGDVAKSNMIACVGRGVCAIKVKESGCSEYVYQYLLYFKAKWGSIEQGSTFS 320

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDS 157
           AIS   I NI IP+     ++ IVA    L A  DS
Sbjct: 321 AISRNDILNINIPV---ITKRLIVA--SHLLALFDS 351
>gi|71904273|ref|YP_281076.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS6180]
 gi|71803368|gb|AAX72721.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS6180]
          Length = 395

 Score = 74.7 bits (182), Expect = 2e-12,   Method: Composition-based stats.
 Identities = 60/170 (35%), Positives = 84/170 (49%), Gaps = 18/170 (10%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 226 LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 284

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 285 SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 339

Query: 123 ISIELIKNILIPLPPSSEQKRI----------VAKIEELFAQLDSITAVL 162
           IS   IK   I +P  SEQ+ I          +A+ EE   +L ++   L
Sbjct: 340 ISSSDIKYAKIQIPSLSEQEAIGNFFQTLDQQIAQSEEKLTELKALKQTL 389
>gi|19746826|ref|NP_607962.1| putative specificity determinant HsdS [Streptococcus pyogenes
           MGAS8232]
 gi|19749064|gb|AAL98461.1| putative specificity determinant HsdS [Streptococcus pyogenes
           MGAS8232]
          Length = 395

 Score = 74.3 bits (181), Expect = 2e-12,   Method: Composition-based stats.
 Identities = 60/170 (35%), Positives = 84/170 (49%), Gaps = 18/170 (10%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 226 LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 284

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 285 SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 339

Query: 123 ISIELIKNILIPLPPSSEQKRI----------VAKIEELFAQLDSITAVL 162
           IS   IK   I +P  SEQ+ I          +A+ EE   +L ++   L
Sbjct: 340 ISSSDIKYAKIQIPSLSEQEAIGNFFQTLDQQIAQSEEKLTELKALKQTL 389
>gi|139474405|ref|YP_001129121.1| type I restriction-modification system S protein [Streptococcus
           pyogenes str. Manfredo]
 gi|134272652|emb|CAM30919.1| type I restriction-modification system S protein [Streptococcus
           pyogenes str. Manfredo]
          Length = 391

 Score = 73.9 bits (180), Expect = 3e-12,   Method: Composition-based stats.
 Identities = 60/170 (35%), Positives = 84/170 (49%), Gaps = 18/170 (10%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 222 LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 280

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 281 SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 335

Query: 123 ISIELIKNILIPLPPSSEQKRI----------VAKIEELFAQLDSITAVL 162
           IS   IK   I +P  SEQ+ I          +A+ EE   +L ++   L
Sbjct: 336 ISSSDIKYAKIQIPSLSEQEAIGNFFQTLDQQIAQSEEKLTELKALKQTL 385
>gi|94989255|ref|YP_597356.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS9429]
 gi|94542763|gb|ABF32812.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS9429]
          Length = 399

 Score = 73.9 bits (180), Expect = 3e-12,   Method: Composition-based stats.
 Identities = 59/156 (37%), Positives = 79/156 (50%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       YI     +T  ITK A K  I+L
Sbjct: 23  LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYI-FPRVWTTQITKQADKGDIIL 81

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 82  SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 136

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   IK   I +P   EQ+     I ELF  +D +
Sbjct: 137 ISSSDIKYAKIQIPSLPEQE----AIGELFQTVDQL 168
>gi|15675718|ref|NP_269892.1| putative type I site-specific deoxyribonuclease [Streptococcus
           pyogenes M1 GAS]
 gi|71911435|ref|YP_282985.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS5005]
 gi|13622936|gb|AAK34613.1| putative type I site-specific deoxyribonuclease [Streptococcus
           pyogenes M1 GAS]
 gi|71854217|gb|AAZ52240.1| type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS5005]
          Length = 399

 Score = 73.6 bits (179), Expect = 4e-12,   Method: Composition-based stats.
 Identities = 58/156 (37%), Positives = 79/156 (50%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 23  LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 81

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 82  SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 136

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   IK   I +P   EQ+     I ELF  +D +
Sbjct: 137 ISSSDIKYAKIQIPSLPEQE----AIGELFQMVDQL 168
>gi|21911179|ref|NP_665447.1| putative type I site-specific deoxyribonuclease hsdS subunit
           [Streptococcus pyogenes MGAS315]
 gi|28896555|ref|NP_802905.1| putative type I site-specific deoxyribonuclease [Streptococcus
           pyogenes SSI-1]
 gi|94995073|ref|YP_603171.1| Type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS10750]
 gi|21905391|gb|AAM80250.1| putative type I site-specific deoxyribonuclease hsdS subunit
           [Streptococcus pyogenes MGAS315]
 gi|28811809|dbj|BAC64738.1| putative type I site-specific deoxyribonuclease [Streptococcus
           pyogenes SSI-1]
 gi|94548581|gb|ABF38627.1| Type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS10750]
          Length = 391

 Score = 73.6 bits (179), Expect = 4e-12,   Method: Composition-based stats.
 Identities = 58/156 (37%), Positives = 79/156 (50%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 23  LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 81

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 82  SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 136

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   IK   I +P   EQ+     I ELF  +D +
Sbjct: 137 ISSSNIKYAKIQIPSLPEQE----AIGELFQTVDQL 168
>gi|94993143|ref|YP_601242.1| Type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS2096]
 gi|94546651|gb|ABF36698.1| Type I restriction-modification system specificity subunit
           [Streptococcus pyogenes MGAS2096]
          Length = 399

 Score = 73.2 bits (178), Expect = 4e-12,   Method: Composition-based stats.
 Identities = 58/156 (37%), Positives = 79/156 (50%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG+I  ITMGQSP S   + N       QG       Y+     +T  ITK A K  I+L
Sbjct: 23  LGDIVQITMGQSPSSQNYTTNPSDYILVQGNADIKNGYV-FPRVWTTQITKQADKGDIIL 81

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQA 122
            VRAPVG V  T   + IGRG+ ++       +EF F  L+  K+   +++ STG+TF +
Sbjct: 82  SVRAPVGDVGKTNYHVIIGRGVAAIK-----GNEFIFQILKYLKEIGYWKRISTGSTFDS 136

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   IK   I +P   EQ+     I ELF  +D +
Sbjct: 137 ISSSDIKYAKIQIPSLPEQE----AIGELFQTVDQL 168
>gi|91213998|ref|YP_543984.1| HsdS, type I site-specific deoxyribonuclease [Escherichia coli
           UTI89]
 gi|117626660|ref|YP_859983.1| HsdS, type I site-specific deoxyribonuclease [Escherichia coli APEC
           O1]
 gi|91075572|gb|ABE10453.1| HsdS, type I site-specific deoxyribonuclease [Escherichia coli
           UTI89]
 gi|115515784|gb|ABJ03859.1| HsdS, type I site-specific deoxyribonuclease [Escherichia coli APEC
           O1]
          Length = 455

 Score = 72.0 bits (175), Expect = 9e-12,   Method: Composition-based stats.
 Identities = 53/169 (31%), Positives = 80/169 (47%), Gaps = 11/169 (6%)

Query: 5   LGEICSITMGQSPESSFISN---NSDGMEFHQGK--IHFTEKYIQKANN------YTFNI 53
           +G+I  +  G +P+S    N   + +G+ +        + EKYI           Y+   
Sbjct: 15  IGDIADVISGGTPKSGVAENFAPSGEGVAWLTPADLSGYKEKYISHGARDLTTLGYSSCS 74

Query: 54  TKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ 113
            K+ PK  IL   RAP+G V I   EI   +G  S      I  ++ +Y+L+  +   E+
Sbjct: 75  AKLMPKGTILFSSRAPIGYVAIAANEIATNQGFKSFAFPSDIFPDYAYYFLRNIRHIAEE 134

Query: 114 KSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
             TGTTF+ IS    K +   L P +EQK I  K++ L AQ+DS  A L
Sbjct: 135 MGTGTTFKEISGSSAKTLPFVLVPFAEQKIIAEKLDTLLAQVDSTKARL 183
>gi|126090307|ref|YP_001041762.1| restriction modification system DNA specificity domain [Shewanella
           baltica OS155]
 gi|125999938|gb|ABN64007.1| restriction modification system DNA specificity domain [Shewanella
           baltica OS155]
          Length = 349

 Score = 70.9 bits (172), Expect = 2e-11,   Method: Composition-based stats.
 Identities = 52/156 (33%), Positives = 79/156 (50%), Gaps = 12/156 (7%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LG +  ITMGQSP     +NNS+     QG       ++     +T  +TK A + A++ 
Sbjct: 180 LGSVTQITMGQSPSGENYTNNSNDFILVQGNADLKNGFVV-PRVWTSEVTKTATQGALIF 238

Query: 65  CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWL-QCQKDTFEQK-STGTTFQA 122
            VRAPVG V  T  ++ +GRG+ ++       +EF F  L + + D +  K S G+TF A
Sbjct: 239 SVRAPVGEVGKTNYDVVLGRGVAAIN-----ANEFIFQQLKKLKSDNYWHKVSAGSTFDA 293

Query: 123 ISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           IS   + + LI +   SEQ      I   F +LD++
Sbjct: 294 ISSTELDSTLIWVSSDSEQ----TAIGNYFQKLDTL 325
>gi|71024881|ref|YP_263290.1| HsdS [Lactococcus lactis subsp. cremoris]
 gi|70067198|dbj|BAE06236.1| HsdS [Lactococcus lactis subsp. cremoris]
          Length = 388

 Score = 70.5 bits (171), Expect = 3e-11,   Method: Composition-based stats.
 Identities = 52/160 (32%), Positives = 81/160 (50%), Gaps = 12/160 (7%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +LG+   I MGQSP S   +++ +     QG        +     +T  +TK A K+ ++
Sbjct: 32  KLGDEVRIVMGQSPNSENYTDDPNDYILVQGNADMKNGRVL-PRVWTTQVTKQAEKDDLI 90

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQ 121
           L VRAPVG +  T  ++ IGRG+ ++       +EF F  L   K    + + STG+TF+
Sbjct: 91  LSVRAPVGDIGKTAYDVVIGRGVAAIK-----GNEFIFQNLGKMKSDGYWTRYSTGSTFE 145

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAV 161
           +I+   IK  +I +P   EQ     KI   F QLD+  A+
Sbjct: 146 SINSTDIKEAIISVPAIEEQD----KIGSFFKQLDNTIAL 181
>gi|32455448|ref|NP_862562.1| HsdS [Lactococcus lactis]
 gi|14251229|gb|AAK57812.1|U35629_2 HsdS [Lactococcus lactis]
          Length = 387

 Score = 70.5 bits (171), Expect = 3e-11,   Method: Composition-based stats.
 Identities = 52/160 (32%), Positives = 81/160 (50%), Gaps = 12/160 (7%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +LG+   I MGQSP S   +++ +     QG        +     +T  +TK A K+ ++
Sbjct: 32  KLGDEVRIVMGQSPNSENYTDDPNDYILVQGNADMKNGRVL-PRVWTTQVTKQAEKDDLI 90

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDT--FEQKSTGTTFQ 121
           L VRAPVG +  T  ++ IGRG+ ++       +EF F  L   K    + + STG+TF+
Sbjct: 91  LSVRAPVGDIGKTAYDVVIGRGVAAIK-----GNEFIFQNLGKMKSDGYWTRYSTGSTFE 145

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAV 161
           +I+   IK  +I +P   EQ     KI   F QLD+  A+
Sbjct: 146 SINSTDIKEAIISVPAIEEQD----KIGSFFKQLDNTIAL 181
>gi|153808175|ref|ZP_01960843.1| hypothetical protein BACCAC_02461 [Bacteroides caccae ATCC 43185]
 gi|149129078|gb|EDM20294.1| hypothetical protein BACCAC_02461 [Bacteroides caccae ATCC 43185]
          Length = 473

 Score = 68.9 bits (167), Expect = 8e-11,   Method: Composition-based stats.
 Identities = 57/170 (33%), Positives = 78/170 (45%), Gaps = 12/170 (7%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDG-------MEFHQGKI-HFTEKYIQKANNYTFNI 53
           W  LGEI     G +P S    +  +G        + + G I H  E   +KA N T   
Sbjct: 307 WVTLGEIGKWQSGSTP-SRLNKDYYNGDIPWLKTGDLNDGYITHIPEYITEKALNET--S 363

Query: 54  TKIAPKNAILLCVR-APVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFE 112
            K+ P  +IL+ +  A +G + I        +  C+      I+ EF FY+L   +  F 
Sbjct: 364 VKLNPTGSILMAMYGATIGKLGILTYPATTNQACCACEIYTGIEKEFLFYFLLSHRADFI 423

Query: 113 QKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
           +   G     IS E I N  IPLPPS EQKRIV  + ++FAQLD I   L
Sbjct: 424 KLGGGGAQPNISKEKIINTYIPLPPSEEQKRIVNAVNDVFAQLDVIMESL 473
>gi|15678961|ref|NP_276078.1| type I restriction modification system, subunit S
           [Methanothermobacter thermautotrophicus str. Delta H]
 gi|2622039|gb|AAB85439.1| type I restriction modification system, subunit S
           [Methanothermobacter thermautotrophicus str. Delta H]
          Length = 407

 Score = 67.4 bits (163), Expect = 3e-10,   Method: Composition-based stats.
 Identities = 35/135 (25%), Positives = 71/135 (52%), Gaps = 6/135 (4%)

Query: 24  NNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIG 83
           +N +G+     +   T K +++ N        + PK +IL+  RAPVG V +   E    
Sbjct: 57  SNLNGIMIADSERKVTVKALKRTN------LNLIPKLSILISTRAPVGYVALNTVECVFN 110

Query: 84  RGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKR 143
           +G  ++ PK  + + ++ Y+L   K   +  S G+TF+ ++ + ++ I +P+PP  EQKR
Sbjct: 111 QGCKALVPKSHVDTRYFAYYLLINKKRLQDLSGGSTFKELNKKTLEKIYLPVPPLEEQKR 170

Query: 144 IVAKIEELFAQLDSI 158
           I   ++++   ++ +
Sbjct: 171 ISEILQDVDGAIEKV 185
>gi|2129238|pir||B64316 restriction modification system S chain homolog - Methanococcus
           jannaschii
          Length = 425

 Score = 66.6 bits (161), Expect = 4e-10,   Method: Composition-based stats.
 Identities = 37/95 (38%), Positives = 58/95 (61%), Gaps = 3/95 (3%)

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR--IQSEFWFYWLQCQKDTFEQ 113
           + PK +I++  RAPVG V +   E    +G   + PK    + +EF+ Y+L+ +K+  E 
Sbjct: 294 LIPKGSIIISTRAPVGYVAVLTVESTFNQGCKGLVPKNNDSVNTEFYAYYLKFKKNLLEN 353

Query: 114 KSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKI 148
            S G+TF+ +S  +++N  IPLPP  EQK+I AKI
Sbjct: 354 LSGGSTFKELSKSMLENFKIPLPPLEEQKQI-AKI 387

 Score = 60.1 bits (144), Expect = 4e-08,   Method: Composition-based stats.
 Identities = 37/120 (30%), Positives = 65/120 (54%), Gaps = 10/120 (8%)

Query: 35  KIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR 94
           KI  TE+ +  +N +      I PKN++L  +   +G   I + E+   + +  + PK  
Sbjct: 70  KIKITEEGLNNSNAW------IVPKNSVLFAMYGSIGETAINKIEVATNQAILGIIPKDN 123

Query: 95  I-QSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQK---RIVAKIEE 150
           I +SEF +Y L   K+ + +    TT + ++ +++K+  IPLPP  EQK   +I+ KI+E
Sbjct: 124 ILESEFLYYILAKNKNYYSKLGMQTTQKNLNAQIVKSFKIPLPPLEEQKQIAKILTKIDE 183
>gi|126007743|ref|ZP_01708319.1| hypothetical protein Faci_03001417 [Ferroplasma acidarmanus fer1]
          Length = 420

 Score = 66.2 bits (160), Expect = 6e-10,   Method: Composition-based stats.
 Identities = 52/171 (30%), Positives = 80/171 (46%), Gaps = 28/171 (16%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHF-TEKYIQKANNYTF------------ 51
           LGEI     G +P +           +H G I + T K + K NN               
Sbjct: 228 LGEIFKTITGTTPSTKVKD------YWHGGTIEWLTPKDLNKLNNTITLPPSERKVTEKA 281

Query: 52  ---NITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSE--FWFYWLQC 106
              N   I P+N+IL+  RAPVG V I   +I   +G   + P  R  S   F+ Y+L+ 
Sbjct: 282 LKENNLNILPENSILISTRAPVGYVGINNTKITFNQGCKGLVPLNRDVSFPFFYAYYLKS 341

Query: 107 QKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDS 157
           +       STG+TF+ +S E + ++++PLPP  EQ+    KI E+ + +D+
Sbjct: 342 KTTFLNSLSTGSTFKELSKEGLDDVVVPLPPLPEQQ----KIGEILSTVDN 388
>gi|99078515|ref|YP_611773.1| restriction modification system DNA specificity domain
           [Silicibacter sp. TM1040]
 gi|99035653|gb|ABF62511.1| type I restriction-modification system; S subunit [Silicibacter sp.
           TM1040]
          Length = 387

 Score = 65.5 bits (158), Expect = 1e-09,   Method: Composition-based stats.
 Identities = 41/153 (26%), Positives = 75/153 (49%), Gaps = 9/153 (5%)

Query: 5   LGEICSITMGQSPESSFISNNSDG-------MEFHQGKIHFTEKYIQKANNYTFNITKIA 57
           LGE+  I  G +P+   + +  DG        +F    +  T   I +A     + T++ 
Sbjct: 6   LGELVEIRGGGTPDKK-VPDYWDGDIPWASVKDFKSTSLASTIDRITQAG-VANSATQVI 63

Query: 58  PKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTG 117
           P   I++  R  VG   I E ++ I + L ++ P  RI  ++  + L     T E ++TG
Sbjct: 64  PAGNIIVPTRMAVGKAAINEIDLAINQDLKALIPSQRIDRQYLLHALLANAKTLEDQATG 123

Query: 118 TTFQAISIELIKNILIPLPPSSEQKRIVAKIEE 150
            T + I ++ ++++ IPLPP  EQ+RI   +++
Sbjct: 124 ATVKGIKLDALRSLQIPLPPLQEQRRIAGILDQ 156
>gi|15669898|ref|NP_247095.1| Type I restriction-modification enzyme, S subunit
           [Methanocaldococcus jannaschii DSM 2661]
 gi|115502461|sp|Q57594|T1S1_METJA Type I restriction enzyme MjaXIP specificity protein (S protein)
           (S.MjaXIP)
 gi|61680619|pdb|1YF2|A Chain A, Three-Dimensional Structure Of Dna Sequence Specificity
           (S) Subunit Of A Type I Restriction-Modification Enzyme
           And Its Functional Implications
 gi|61680620|pdb|1YF2|B Chain B, Three-Dimensional Structure Of Dna Sequence Specificity
           (S) Subunit Of A Type I Restriction-Modification Enzyme
           And Its Functional Implications
          Length = 425

 Score = 65.1 bits (157), Expect = 1e-09,   Method: Composition-based stats.
 Identities = 36/95 (37%), Positives = 58/95 (61%), Gaps = 3/95 (3%)

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR--IQSEFWFYWLQCQKDTFEQ 113
           + PK +I++  RAPVG V +   E    +G   ++ K    + +EF+ Y+L+ +K+  E 
Sbjct: 294 LIPKGSIIISTRAPVGYVAVLTVESTFNQGCKGLFQKNNDSVNTEFYAYYLKFKKNLLEN 353

Query: 114 KSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKI 148
            S G+TF+ +S  +++N  IPLPP  EQK+I AKI
Sbjct: 354 LSGGSTFKELSKSMLENFKIPLPPLEEQKQI-AKI 387

 Score = 60.1 bits (144), Expect = 4e-08,   Method: Composition-based stats.
 Identities = 37/120 (30%), Positives = 65/120 (54%), Gaps = 10/120 (8%)

Query: 35  KIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR 94
           KI  TE+ +  +N +      I PKN++L  +   +G   I + E+   + +  + PK  
Sbjct: 70  KIKITEEGLNNSNAW------IVPKNSVLFAMYGSIGETAINKIEVATNQAILGIIPKDN 123

Query: 95  I-QSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQK---RIVAKIEE 150
           I +SEF +Y L   K+ + +    TT + ++ +++K+  IPLPP  EQK   +I+ KI+E
Sbjct: 124 ILESEFLYYILAKNKNYYSKLGMQTTQKNLNAQIVKSFKIPLPPLEEQKQIAKILTKIDE 183
>gi|15611482|ref|NP_223133.1| TYPE I RESTRICTION ENZYME (SPECIFICITY SUBUNIT) [Helicobacter
           pylori J99]
 gi|4154940|gb|AAD05986.1| TYPE I RESTRICTION ENZYME (SPECIFICITY SUBUNIT) [Helicobacter
           pylori J99]
          Length = 409

 Score = 64.7 bits (156), Expect = 2e-09,   Method: Composition-based stats.
 Identities = 32/99 (32%), Positives = 53/99 (53%)

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKS 115
           + PK+AIL   RAP+G V I E+ +C  +G  S+ P  +I  EF +Y L+  KD      
Sbjct: 68  LLPKHAILFSSRAPIGYVAIAEKRLCTNQGFKSIIPNKKIYFEFLYYLLKYYKDNISNIG 127

Query: 116 TGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQ 154
            GTTF+ +S   +    + +PP+  +++ +A    +  Q
Sbjct: 128 GGTTFKEVSGATLGLFQVKIPPTYYEQQKIAHTLSILDQ 166
>gi|21673510|ref|NP_661575.1| type I restriction system specificity protein [Chlorobium tepidum
           TLS]
 gi|21646618|gb|AAM71917.1| type I restriction system specificity protein [Chlorobium tepidum
           TLS]
          Length = 444

 Score = 63.5 bits (153), Expect = 4e-09,   Method: Composition-based stats.
 Identities = 43/142 (30%), Positives = 71/142 (50%), Gaps = 3/142 (2%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LGE C I MGQSP S   +    G+    G   F   +   A  +T ++ K A    IL 
Sbjct: 12  LGEECEIVMGQSPPSETCNTVGIGIPLLNGPTEFGPHHPSPA-QFTTDVRKRAIPGDILF 70

Query: 65  CVR-APVGVVNITEREICIGRGLCSVYPKYRIQSE-FWFYWLQCQKDTFEQKSTGTTFQA 122
           CVR +  G +N  ++E  IGRG+ ++  K++ + + F    ++C       ++TG+TF  
Sbjct: 71  CVRGSTTGRMNWADQEYAIGRGIAAIRHKFKPELQPFVRAVIECYLPELLAQATGSTFPN 130

Query: 123 ISIELIKNILIPLPPSSEQKRI 144
           +S + + N+  P   + EQ+ I
Sbjct: 131 VSAQQLSNLKWPELAADEQRAI 152
>gi|149369906|ref|ZP_01889757.1| type I restriction-modification system specificity determinant
           protein [unidentified eubacterium SCB49]
 gi|149356397|gb|EDM44953.1| type I restriction-modification system specificity determinant
           protein [unidentified eubacterium SCB49]
          Length = 415

 Score = 63.2 bits (152), Expect = 4e-09,   Method: Composition-based stats.
 Identities = 36/116 (31%), Positives = 58/116 (50%), Gaps = 6/116 (5%)

Query: 22  ISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREIC 81
           +SNN       +G+   TE+ I+K++       K+ P   +LL  RAP+G + I    + 
Sbjct: 235 LSNNKGKKYITRGEYDVTEQGIKKSS------LKLMPSGTVLLSSRAPIGYLAIARETVT 288

Query: 82  IGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPP 137
             +G  S  PK    SEF +Y ++ +    E +S G+TF+ +S   +K I I  PP
Sbjct: 289 TNQGFKSFEPKSYFTSEFLYYQIKNKIPLIEARSGGSTFKEVSASTLKTIKIITPP 344
>gi|26554275|ref|NP_758209.1| type I restriction-modification system S subunit [Mycoplasma
           penetrans HF-2]
 gi|26454284|dbj|BAC44613.1| type I restriction-modification system S subunit [Mycoplasma
           penetrans HF-2]
          Length = 519

 Score = 63.2 bits (152), Expect = 4e-09,   Method: Composition-based stats.
 Identities = 45/168 (26%), Positives = 87/168 (51%), Gaps = 13/168 (7%)

Query: 2   WCRLGEICSITMGQSPESSFISN-NSDGMEF-------HQGKIHFTEKYIQKANNYTFNI 53
           W +L  IC + +G++P+ S  +N NS+ + +         GKI  T++YI     +    
Sbjct: 353 WIKLSNICELILGKTPKRSINTNWNSNDINWVTISDMKDLGKIFSTKEYITN-EAFKNEF 411

Query: 54  TKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQ---SEFWFYWLQCQKDT 110
           T+I+ K ++L+  +  +G  +I E +      + ++ P Y       +F FY L     +
Sbjct: 412 TRISKKESLLMSFKLTIGRTSILEIDAVHNEAIVTINPYYDKDYAIRDFLFYTLGTFV-S 470

Query: 111 FEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           F +K++      I+ E + N+L+ LPP +EQ+RI+  I ++ + ++SI
Sbjct: 471 FIEKTSAIKGSTINKEKMINMLVSLPPINEQRRIIKSISKIHSLINSI 518
>gi|154174911|ref|YP_001408734.1| putative type I restriction-modification system, S subunit
           [Campylobacter curvus 525.92]
 gi|153793147|gb|EAT99394.2| putative type I restriction-modification system, S subunit
           [Campylobacter curvus 525.92]
          Length = 528

 Score = 63.2 bits (152), Expect = 5e-09,   Method: Composition-based stats.
 Identities = 52/172 (30%), Positives = 82/172 (47%), Gaps = 17/172 (9%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEF-------HQGKIHFTEKYIQKANNYTFNIT 54
           W RLG +     G +P +       D M F       + G  +  EK  +K       + 
Sbjct: 362 WVRLGSMGVAQTGSTPSTQVRDFYGDYMPFIKPADITNSGIDYNNEKLSKKGTE----VG 417

Query: 55  KIAPKNAILL-CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ 113
           ++A K +IL+ C+   +G     +R +   + + S+ P +    +F FY+L      FEQ
Sbjct: 418 RVAEKGSILMVCIGGSLGKCYFNDRIVSFNQQINSLTPFFS-SYKFIFYYL-LSSHFFEQ 475

Query: 114 ---KSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
              ++TGT    ++    ++ILIPLPP  EQKRIV KIEEL   +D + + L
Sbjct: 476 LQDRATGTATPIVNKTSWESILIPLPPLPEQKRIVTKIEELLKFVDILQSSL 527
>gi|17232094|ref|NP_488642.1| type I site-specific deoxyribonuclease chain S [Nostoc sp. PCC
           7120]
 gi|17133739|dbj|BAB76301.1| type I site-specific deoxyribonuclease chain S [Nostoc sp. PCC
           7120]
          Length = 390

 Score = 62.4 bits (150), Expect = 8e-09,   Method: Composition-based stats.
 Identities = 49/166 (29%), Positives = 82/166 (49%), Gaps = 10/166 (6%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           LGEI  I MGQSP     +NN  G+    G   F ++       +T +  KIA    +L 
Sbjct: 9   LGEIADIIMGQSPTGETCNNNGQGLPLLNGPTEFGDRN-PLPTQFTIDPKKIAEAGDLLF 67

Query: 65  CVR-APVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYW-LQCQKDTFEQKSTGTTFQA 122
           CVR +  G +N  +++  IGRG+ S+  K  I  + +    ++ +  +    +TG+TF  
Sbjct: 68  CVRGSTTGRMNWADQKYAIGRGIASIRAKDGILFQPYIRAIIEKELKSLLAVATGSTFPN 127

Query: 123 ISIELIKNILIPLPPS------SEQKRIV-AKIEELFAQLDSITAV 161
           IS + + N+++ LP        S   RI+  KI  L +Q +++ A+
Sbjct: 128 ISKDHLLNLIVQLPSKNIKIYISNLARILDEKIYNLRSQNETLEAI 173
>gi|15645090|ref|NP_207260.1| type I restriction enzyme S protein (hsdS) [Helicobacter pylori
           26695]
 gi|2313566|gb|AAD07524.1| type I restriction enzyme S protein (hsdS) [Helicobacter pylori
           26695]
          Length = 365

 Score = 62.0 bits (149), Expect = 1e-08,   Method: Composition-based stats.
 Identities = 31/99 (31%), Positives = 53/99 (53%)

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKS 115
           + PK+AIL   RAP+G V I E+ +C  +G  S+ P  +I  EF +Y L+  K+ F    
Sbjct: 68  LLPKHAILFSSRAPIGYVAIAEKRLCTNQGFKSIIPNKKIYFEFLYYLLKYHKNNFINMG 127

Query: 116 TGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQ 154
            GTT + I    +    + +PP+  +++ +A+   +  Q
Sbjct: 128 EGTTIKGIYNIALGLFKVKIPPTYYEQQKIARTLSILDQ 166
>gi|154245043|ref|YP_001416001.1| restriction modification system DNA specificity domain
           [Xanthobacter autotrophicus Py2]
 gi|154159128|gb|ABS66344.1| restriction modification system DNA specificity domain
           [Xanthobacter autotrophicus Py2]
          Length = 450

 Score = 61.2 bits (147), Expect = 2e-08,   Method: Composition-based stats.
 Identities = 48/163 (29%), Positives = 76/163 (46%), Gaps = 9/163 (5%)

Query: 1   MWCRLGEICSITMGQSPESSFISN-NSDGMEFHQGK--IHFTEKYIQKA------NNYTF 51
           +W   GE+  I  G +P +   +N    G+ +        + E YI +         Y  
Sbjct: 13  LWASFGEVADIVGGGTPPTGDEANFTKQGVPWLTPADLTGYRETYISRGRRDLSEKGYRE 72

Query: 52  NITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTF 111
           +  ++ PK  +L   RAPVG   I    +   +G  S   K  I  E+  ++L    +  
Sbjct: 73  SAARLLPKGTVLFSSRAPVGYCAIASENVSTNQGFKSFILKGDISPEYVRHYLLGSTEYA 132

Query: 112 EQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQ 154
           E K++GTTF+ +S      + +PLPP  EQ+RIVAKI+ L A+
Sbjct: 133 ESKASGTTFKELSGSRATELALPLPPLPEQRRIVAKIDSLTAK 175
>gi|145622248|ref|ZP_01778209.1| hypothetical protein PmobDRAFT_1688 [Petrotoga mobilis SJ95]
 gi|144947397|gb|EDJ82428.1| hypothetical protein PmobDRAFT_1688 [Petrotoga mobilis SJ95]
          Length = 121

 Score = 60.5 bits (145), Expect = 3e-08,   Method: Composition-based stats.
 Identities = 35/81 (43%), Positives = 46/81 (56%), Gaps = 1/81 (1%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILL 64
           L E+  I MGQSP+S   S   +G+ F+QG+  F  KY    + Y  +  KIA K  +LL
Sbjct: 24  LQEVRQINMGQSPKSETYSEVGNGLPFYQGRRDFGRKY-PTVSVYCTDPKKIAIKGDVLL 82

Query: 65  CVRAPVGVVNITEREICIGRG 85
            VRAPVG +N    + CIG G
Sbjct: 83  TVRAPVGEINEAIEDCCIGLG 103
>gi|126208116|ref|YP_001053341.1| putative type I restriction system specificity protein
           [Actinobacillus pleuropneumoniae L20]
 gi|126096908|gb|ABN73736.1| putative type I restriction system specificity protein
           [Actinobacillus pleuropneumoniae L20]
          Length = 413

 Score = 60.1 bits (144), Expect = 4e-08,   Method: Composition-based stats.
 Identities = 42/131 (32%), Positives = 65/131 (49%), Gaps = 2/131 (1%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +LG+I  I MGQSP SS ++    G     G   FT  Y      YT    K A K  +L
Sbjct: 3   KLGDIADIVMGQSPSSSDVNMERIGDPLLNGPTEFTSFYPSPV-QYTEKGKKFAEKGDLL 61

Query: 64  LCVR-APVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQA 122
            CVR +  G +N  +++  IGRGL ++  K    ++F    L+   +   Q +TG+TF  
Sbjct: 62  FCVRGSTTGRINFADQKYAIGRGLAAIRGKNGYPTKFIELILKDCLERILQSATGSTFPN 121

Query: 123 ISIELIKNILI 133
           +S  ++ ++ I
Sbjct: 122 VSQAMLLDLDI 132
>gi|75909474|ref|YP_323770.1| Restriction modification system DNA specificity domain [Anabaena
           variabilis ATCC 29413]
 gi|75703199|gb|ABA22875.1| Restriction modification system DNA specificity domain [Anabaena
           variabilis ATCC 29413]
          Length = 557

 Score = 60.1 bits (144), Expect = 4e-08,   Method: Composition-based stats.
 Identities = 47/157 (29%), Positives = 77/157 (49%), Gaps = 4/157 (2%)

Query: 4   RLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAIL 63
           +LG +  +  GQ       ++  +G+ +  G   F E +   A  +T  +  +A KN IL
Sbjct: 90  KLGYLIELVSGQHLGQEEQNDQGEGLPYLTGPADFGE-FNPVATRWTNTVKALAKKNDIL 148

Query: 64  LCVR-APVGVVNITERE-ICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQ 121
           + V+ A VG  NI   E   IGR L ++ P   ++ EF +  +    + F+  S G+T  
Sbjct: 149 ITVKGAGVGKANILSMEKAAIGRQLMAIRP-ILLEYEFIYLLIISSYEKFQALSIGSTVP 207

Query: 122 AISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
            +  + I +  + LPP +EQKRIV K + L +  D I
Sbjct: 208 GMGRKDILDFSLGLPPLAEQKRIVEKCDRLLSTCDEI 244
>gi|2826248|gb|AAB98112.1| type I restriction-modification enzyme 2, S subunit
           [Methanocaldococcus jannaschii DSM 2661]
          Length = 343

 Score = 59.7 bits (143), Expect = 6e-08,   Method: Composition-based stats.
 Identities = 37/120 (30%), Positives = 65/120 (54%), Gaps = 10/120 (8%)

Query: 35  KIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR 94
           KI  TE+ +  +N +      I PKN++L  +   +G   I + E+   + +  + PK  
Sbjct: 70  KIKITEEGLNNSNAW------IVPKNSVLFAMYGSIGETAINKIEVATNQAILGIIPKDN 123

Query: 95  I-QSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQK---RIVAKIEE 150
           I +SEF +Y L   K+ + +    TT + ++ +++K+  IPLPP  EQK   +I+ KI+E
Sbjct: 124 ILESEFLYYILAKNKNYYSKLGMQTTQKNLNAQIVKSFKIPLPPLEEQKQIAKILTKIDE 183
>gi|120601902|ref|YP_966302.1| restriction modification system DNA specificity domain
           [Desulfovibrio vulgaris subsp. vulgaris DP4]
 gi|120562131|gb|ABM27875.1| restriction modification system DNA specificity domain
           [Desulfovibrio vulgaris subsp. vulgaris DP4]
          Length = 595

 Score = 59.3 bits (142), Expect = 7e-08,   Method: Composition-based stats.
 Identities = 27/60 (45%), Positives = 41/60 (68%)

Query: 98  EFWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDS 157
           +F+F  LQ   D   +K++  T + +S + +  I +PLPP +EQ RIVAKIEELF++LD+
Sbjct: 122 KFFFLALQPYLDAIHEKTSSVTVKHLSSKTVNEIPLPLPPLNEQNRIVAKIEELFSELDA 181
>gi|21228301|ref|NP_634223.1| type I restriction-modification system specificity subunit
           [Methanosarcina mazei Go1]
 gi|20906763|gb|AAM31895.1| type I restriction-modification system specificity subunit
           [Methanosarcina mazei Go1]
          Length = 384

 Score = 58.9 bits (141), Expect = 8e-08,   Method: Composition-based stats.
 Identities = 43/142 (30%), Positives = 71/142 (50%), Gaps = 11/142 (7%)

Query: 12  TMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVG 71
           TMGQSP  +  +    G    QG   F  ++ ++   YT   T++A  +  L+ VRAPVG
Sbjct: 210 TMGQSPPGTSYNQEGVGKPMFQGNADFGFRFPEE-RVYTTEPTRLAYPHDTLISVRAPVG 268

Query: 72  VVNITEREICIGRGLCSVYPKYRIQSEFW---FYWLQCQKDTFEQ-KSTGTTFQAISIEL 127
             N+ + E CIGRG+ +   +Y+  ++F+   ++ L+   D  ++    GT F +IS   
Sbjct: 269 AQNMAKVECCIGRGVSAF--RYKANNDFYTYTYFKLRSLMDEIKKFNDEGTVFGSISKTD 326

Query: 128 IKNILIPLPPSSEQKRIVAKIE 149
              + I +PP      I+ K E
Sbjct: 327 FLQMGIAIPPED----IIEKFE 344
>gi|32266920|ref|NP_860952.1| hypothetical protein HH1421 [Helicobacter hepaticus ATCC 51449]
 gi|32262972|gb|AAP78018.1| hypothetical protein HH_1421 [Helicobacter hepaticus ATCC 51449]
          Length = 422

 Score = 58.9 bits (141), Expect = 9e-08,   Method: Composition-based stats.
 Identities = 43/162 (26%), Positives = 79/162 (48%), Gaps = 10/162 (6%)

Query: 5   LGEICSITMGQSPESS---FISNN-----SDGMEFHQGKIHFTEKYIQKANNYTFNITKI 56
           LGE+  +  G +P ++   F S N     S  +E     I   +K+I ++     + T +
Sbjct: 234 LGEVFKVIGGGTPSTANPKFWSGNIAWITSANIENENFTI-IPKKFINQSA-IQASATNL 291

Query: 57  APKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKST 116
            PKN I++  R  +G V IT+ E C  +   ++ P   +  +F  + ++ +   F   S 
Sbjct: 292 VPKNTIIVVTRVGLGKVGITDVETCFSQDSQALLPLIDLNVKFMAFQIRNKAQNFIVSSR 351

Query: 117 GTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           GTT   I+ + +K + + +PP + Q +IV  +E  FA L+ +
Sbjct: 352 GTTINGITKDTLKKVALKIPPLATQNQIVQILESKFAHLEKL 393
>gi|67919739|ref|ZP_00513304.1| Restriction modification system DNA specificity domain:N-6 DNA
           methylase:Type I restriction-modification system, M
           subunit [Chlorobium limicola DSM 245]
 gi|67782637|gb|EAM42051.1| Restriction modification system DNA specificity domain:N-6 DNA
           methylase:Type I restriction-modification system, M
           subunit [Chlorobium limicola DSM 245]
          Length = 846

 Score = 58.5 bits (140), Expect = 1e-07,   Method: Composition-based stats.
 Identities = 46/176 (26%), Positives = 77/176 (43%), Gaps = 24/176 (13%)

Query: 5   LGEICSITMGQSPESS------------FISNNSDGMEFHQGKIHFTEKYIQKANNYTFN 52
           L EIC+   G +P S+             +S +  G E    +   T+K ++ +N     
Sbjct: 461 LDEICTFMTGGTPTSTIAEYYEGGTVPWLVSGDIHGFEIMACEKRITQKAVENSN----- 515

Query: 53  ITKIAPKNAILLCV------RAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQC 106
             K+ PK+++L+ +      R  V ++ +T           +  P  R  SEF F+ L+ 
Sbjct: 516 -AKVLPKDSVLIALNGQGKTRGTVALLRMTGATCNQSLVAITPAPPPRAISEFIFWALRS 574

Query: 107 QKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
                   +  T    ++I ++KNI IPLPP   QK +VA+IE     ++   AVL
Sbjct: 575 MYSDIRALTGDTERSGLNIPILKNIQIPLPPLEVQKEVVAEIEGYQNVINGARAVL 630
>gi|23452743|gb|AAN33142.1| putative type I specificity subunit HsdS [Campylobacter jejuni]
          Length = 414

 Score = 58.5 bits (140), Expect = 1e-07,   Method: Composition-based stats.
 Identities = 41/175 (23%), Positives = 83/175 (47%), Gaps = 20/175 (11%)

Query: 2   WCRLGEICSITMGQSPESS----FISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITK-- 55
           W  LGEIC I  G +P++     + S+ +D ++F   K ++    +     Y +   +  
Sbjct: 215 WKSLGEICEILGGGTPDTKNPIFWYSSQADEVQFE--KSYYWATLVDTKQKYLYGTKRKI 272

Query: 56  -----------IAPKNAILLCVRAPVGVVNITEREICIGRGLCS-VYPKYRIQSEFWFYW 103
                      + P N+++   RA +G ++I + E    +G  + +  +  +  EF ++ 
Sbjct: 273 TQKGLDCSNAILLPINSVIFSSRASIGEISIAKVETATNQGYKNFICDESILYYEFLYFA 332

Query: 104 LQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
           L+      E  + GTT++ +S   IK+  IPLPP  EQ++I   ++ +F +  ++
Sbjct: 333 LKHFTKEIELLAQGTTYKEVSKAKIKDFKIPLPPLKEQEQITKHLDFIFEKAKAL 387
>gi|156866399|gb|EDO59771.1| hypothetical protein CLOLEP_03821 [Clostridium leptum DSM 753]
          Length = 444

 Score = 58.2 bits (139), Expect = 1e-07,   Method: Composition-based stats.
 Identities = 43/156 (27%), Positives = 75/156 (48%), Gaps = 3/156 (1%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNA 61
           W R  +I ++  G+  + +  ++   G+ +  G  +  E  +     +  N   I+ KN 
Sbjct: 36  WVRFSKIINLISGRDAKLTDCNSLGIGIPYILGASNL-ENNVFTIERWIENPQVISLKND 94

Query: 62  ILLCVRAPVGVVNITERE-ICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTTF 120
           +LL V+  +G V + + E + I R + ++     +   F ++ +    D+F Q   G   
Sbjct: 95  VLLSVKGTIGKVYLQKEEKVNISRQIMAIRTSSTLFPRFTYWLVNNISDSFRQAGNGL-I 153

Query: 121 QAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
             IS E I    +P PP  EQ+RIV +IE LFA+LD
Sbjct: 154 PGISREDILQKEVPFPPLPEQQRIVDRIESLFAKLD 189
>gi|94267246|ref|ZP_01290822.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
 gi|93452076|gb|EAT02762.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
          Length = 578

 Score = 57.0 bits (136), Expect = 3e-07,   Method: Composition-based stats.
 Identities = 41/174 (23%), Positives = 75/174 (43%), Gaps = 17/174 (9%)

Query: 2   WCRLGEICSITMGQSPE------------SSFISNNSDGMEFHQGKIHFTEKYIQKANNY 49
           WCR  +I  +  G+S               + +    D     +    FT  Y    N  
Sbjct: 390 WCRFTDIGELARGKSKHRPRNDPALYIGGKTPLVQTGDVARADRKITTFTALY----NQA 445

Query: 50  TFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYP-KYRIQSEFWFYWLQCQK 108
               +K+     + + + A +G   I   + C    +    P   R+++E++ Y+L+  K
Sbjct: 446 GVEQSKLWKAGTLCITIAANIGDTGILGFDACFPDSVVGFTPFDDRLKNEYFEYFLRTAK 505

Query: 109 DTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
              E+ +  T  + I++E+++N+L+PLPP+ E  RIV K ++L    D   A L
Sbjct: 506 KNLEEFAPSTAQKNINLEVLQNVLVPLPPARELVRIVEKTDKLMGLCDQFKASL 559
>gi|26554274|ref|NP_758208.1| type I restriction-modification system S subunit [Mycoplasma
           penetrans HF-2]
 gi|26454283|dbj|BAC44612.1| type I restriction-modification system S subunit [Mycoplasma
           penetrans HF-2]
          Length = 415

 Score = 56.6 bits (135), Expect = 4e-07,   Method: Composition-based stats.
 Identities = 36/90 (40%), Positives = 51/90 (56%), Gaps = 6/90 (6%)

Query: 55  KIAPKNAILLCVRAPVG--VVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQ--KDT 110
           KIAPKN ILLC+        + IT +++C G  L  +   +   + + FY+LQ    K+ 
Sbjct: 312 KIAPKNKILLCIEGGSAGKKIGITSKDVCFGNKLVCI-NDFLSNNLYLFYFLQSYYFKNI 370

Query: 111 FEQKSTGTTFQAISIELIKNILIPLPPSSE 140
           F Q +TG     ISI+ +KNI+IPLPP  E
Sbjct: 371 FNQLTTGI-IGGISIQNLKNIMIPLPPKRE 399
>gi|94263107|ref|ZP_01286925.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
 gi|93456478|gb|EAT06592.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
          Length = 456

 Score = 56.6 bits (135), Expect = 4e-07,   Method: Composition-based stats.
 Identities = 38/122 (31%), Positives = 61/122 (50%), Gaps = 5/122 (4%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKY-IQKANNYTFNITKIAPKNAIL 63
           +G    + MGQSP     +    G+ F QG+  F  +Y  Q+   Y     ++A     L
Sbjct: 268 VGAKFDVIMGQSPPGQSYNEIGQGLPFFQGRRDFGFRYPTQRV--YCTEPKRLANPGDTL 325

Query: 64  LCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ-KSTGTTFQA 122
           + VRAPVG +N+   + CIGRG+ +V  K   +S F +Y ++   + F   +  GT F +
Sbjct: 326 ISVRAPVGDINMARVKCCIGRGVAAVRHKSESRS-FTYYSMRALTEQFSSYEGEGTVFGS 384

Query: 123 IS 124
           I+
Sbjct: 385 IN 386
>gi|10954528|ref|NP_044167.1| type I restriction enzyme subunit S [Methanocaldococcus jannaschii
           DSM 2661]
 gi|12229988|sp|Q60296|T1SH_METJA Putative type I restriction enzyme MjaXP specificity protein (S
           protein) (S.MjaXP)
 gi|1522674|gb|AAC37110.1| hypothetical protein MJ_ECL41 [Methanocaldococcus jannaschii DSM
           2661]
          Length = 432

 Score = 56.6 bits (135), Expect = 5e-07,   Method: Composition-based stats.
 Identities = 34/122 (27%), Positives = 63/122 (51%), Gaps = 8/122 (6%)

Query: 25  NSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGR 84
           N + +   +G+ + TEK +++ +       +I PK  ILL  RAP+G V I +  +   +
Sbjct: 64  NYEYIYISRGERNITEKAVKECS------LRIFPKGTILLTSRAPIGYVAIAKNPLTTNQ 117

Query: 85  GLCSVYPKYRIQSEFWFYWLQCQ--KDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQK 142
           G  ++ PK  + SE+ +Y  + +   +  +  S G+TF  +    +K + IP P   EQ+
Sbjct: 118 GFRNIIPKDGVVSEYLYYLFKTKTMSEYLKDISGGSTFPELKGSTLKEVEIPYPSPEEQQ 177

Query: 143 RI 144
           +I
Sbjct: 178 KI 179
>gi|51893047|ref|YP_075738.1| type I restriction-modification system specificity determinant
           protein [Symbiobacterium thermophilum IAM 14863]
 gi|51856736|dbj|BAD40894.1| type I restriction-modification system specificity determinant
           protein [Symbiobacterium thermophilum IAM 14863]
          Length = 400

 Score = 56.2 bits (134), Expect = 6e-07,   Method: Composition-based stats.
 Identities = 30/105 (28%), Positives = 54/105 (51%), Gaps = 1/105 (0%)

Query: 39  TEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSE 98
           TE+ I +A     + + + P+  +LL  RAP+G + IT+    I +G  ++ P  ++  E
Sbjct: 246 TERMISEAGLKKIS-SGLLPEGTLLLSSRAPIGYLAITKIPTAINQGFIAMPPAGQLSPE 304

Query: 99  FWFYWLQCQKDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKR 143
           +  +W     DT +Q + G+TF  IS    + I + +PP+    R
Sbjct: 305 YMLFWSHYNLDTIKQHANGSTFMEISKAAFRKIKLVVPPAQLVNR 349
>gi|23452748|gb|AAN33144.1| putative type I specificity subunit HsdS [Campylobacter jejuni]
          Length = 378

 Score = 55.5 bits (132), Expect = 9e-07,   Method: Composition-based stats.
 Identities = 43/165 (26%), Positives = 79/165 (47%), Gaps = 9/165 (5%)

Query: 6   GEICSITMGQSPESSFISNNSDGM-------EFHQGKIHFTEKYIQKANNYTFNITKIAP 58
           G+I  I+ G +P  +      +G+       +  +  I  TE++I + +    +  K+  
Sbjct: 1   GDIAEISSGGTPSRNKKEYWENGIIPWVKIKDIKENFISTTEEFITE-DGLKNSSAKLFK 59

Query: 59  KNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYR-IQSEFWFYWLQCQKDTFEQKSTG 117
           K  +L  + A +G V I + +    + +  +  K   I S +  Y+L+  KD    K  G
Sbjct: 60  KGTLLYSIFATLGEVAILDIDATTNQAIAGINIKENNINSLYLMYFLKSIKDKICSKGRG 119

Query: 118 TTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
                +++ ++K I IPLPP  EQ+RIV  ++E FA++D    +L
Sbjct: 120 VAQNNLNLTILKQIQIPLPPLKEQERIVGILDESFAKIDESIKIL 164
>gi|149115916|ref|ZP_01842652.1| restriction modification system DNA specificity domain [Shewanella
           baltica OS223]
 gi|146864200|gb|EDK49637.1| restriction modification system DNA specificity domain [Shewanella
           baltica OS223]
          Length = 406

 Score = 55.5 bits (132), Expect = 1e-06,   Method: Composition-based stats.
 Identities = 33/104 (31%), Positives = 54/104 (51%), Gaps = 1/104 (0%)

Query: 55  KIAPKNAILL-CVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ 113
           K  PK++ILL  + A +G V + E E C  + +  ++ K +I  ++ +YWL   K     
Sbjct: 75  KKMPKDSILLVAIGATIGKVGLAEDESCFNQQVTGIHFKEKIHPKYAYYWLSYIKPEIIT 134

Query: 114 KSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDS 157
           KS+  T   I+   IK +    P   EQK IV K++ L  ++D+
Sbjct: 135 KSSQATLPIINQTGIKGLSFLYPEKEEQKCIVEKLDALLTRIDT 178
>gi|26991425|ref|NP_746850.1| type I restriction-modification system, S subunit [Pseudomonas
           putida KT2440]
 gi|24986497|gb|AAN70314.1|AE016672_5 type I restriction-modification system, S subunit [Pseudomonas
           putida KT2440]
          Length = 576

 Score = 55.1 bits (131), Expect = 1e-06,   Method: Composition-based stats.
 Identities = 40/131 (30%), Positives = 61/131 (46%), Gaps = 10/131 (7%)

Query: 37  HFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGVVNITEREIC---IGRGLCSVYPKY 93
           H  EK I    +  +  T+I     IL+ V   +G + I         I R +C V P  
Sbjct: 133 HKPEKSISPEVDRQYERTRIR-GGEILMGVVGSIGKLGIAPESWAGANIARAICRVVPSV 191

Query: 94  RIQSEFWFYWLQCQKDTFEQKSTGTTFQ----AISIELIKNILIPLPPSSEQKRIVAKIE 149
            +  ++  + LQ   D   ++  G T       +++ LI++   PLPP +EQ RIVAK+E
Sbjct: 192 HVSKDYIIWLLQ--SDLMRKQFLGDTRTLAQPTLNVGLIRSAAAPLPPLAEQHRIVAKVE 249

Query: 150 ELFAQLDSITA 160
           EL A  D + A
Sbjct: 250 ELMALCDRLEA 260
>gi|94266803|ref|ZP_01290467.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
 gi|93452525|gb|EAT03114.1| Restriction modification system DNA specificity domain [delta
           proteobacterium MLMS-1]
          Length = 603

 Score = 55.1 bits (131), Expect = 1e-06,   Method: Composition-based stats.
 Identities = 42/163 (25%), Positives = 71/163 (43%), Gaps = 6/163 (3%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGME-FHQGK-----IHFTEKYIQKANNYTFNITK 55
           +CRL EI +   G +P+         G+  F  G+     I   E+ I +      ++  
Sbjct: 107 YCRLNEIGTWGSGATPKRGITEYYDGGIPWFKSGELVGDFISSAEETITERALKETSVRL 166

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKS 115
             P + ++    A +G  +I +      + +C+  P   I + +   +L+  K  F    
Sbjct: 167 NLPGDVLIAMYGATIGKASILKCHATTNQAVCACTPFSGILNTYLLNFLKASKRHFTSMG 226

Query: 116 TGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
            G     IS E I  ++ PLPP +EQ RIV K++EL A  D +
Sbjct: 227 AGGAQPNISKEKIIAVVFPLPPLAEQHRIVEKVDELMALCDRL 269
>gi|150401945|ref|YP_001329239.1| restriction modification system DNA specificity domain
           [Methanococcus maripaludis C7]
 gi|150032975|gb|ABR65088.1| restriction modification system DNA specificity domain
           [Methanococcus maripaludis C7]
          Length = 432

 Score = 55.1 bits (131), Expect = 1e-06,   Method: Composition-based stats.
 Identities = 42/160 (26%), Positives = 70/160 (43%), Gaps = 9/160 (5%)

Query: 5   LGEICSITMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNIT---------K 55
           L EIC++  G +P+++      DG          T K     +     IT         K
Sbjct: 244 LDEICTVMGGGTPKTNVPEYWQDGTILWATPTDMTSKKSPVIDTTEKKITELGLKESSAK 303

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKS 115
           + PK +IL+  RA +G  +I  +EI   +G  ++    ++ + F  Y L+  KD     +
Sbjct: 304 LVPKGSILMTSRATIGYSSIAMKEISTNQGFINIICDKKVSNYFILYLLEHIKDKIIALA 363

Query: 116 TGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQL 155
            G+TF  IS    KNI + +P     ++    IEEL  ++
Sbjct: 364 NGSTFLEISKTNFKNIRVIVPDYQTMEKYNEIIEELINKI 403
>gi|119356723|ref|YP_911367.1| N-6 DNA methylase [Chlorobium phaeobacteroides DSM 266]
 gi|119354072|gb|ABL64943.1| N-6 DNA methylase [Chlorobium phaeobacteroides DSM 266]
          Length = 834

 Score = 54.7 bits (130), Expect = 2e-06,   Method: Composition-based stats.
 Identities = 42/153 (27%), Positives = 72/153 (47%), Gaps = 14/153 (9%)

Query: 12  TMGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVG 71
           T+   P + FI+      E    +   +E+ +++++       K+ P+N++++  RA +G
Sbjct: 491 TLADLPSTDFIT------EIRSTRRTISERGLRESS------AKMIPENSVIVSTRATIG 538

Query: 72  VVNITEREICIGRGLCSVY--PKYRIQSEFWFYWLQCQKDTFEQKSTGTTFQAISIELIK 129
            + I    +   +G  +V    K ++ SE+    L     T    +TG TF+ I      
Sbjct: 539 RIAINRIPMATNQGFKNVIIEDKSKVISEYVALALTKLVPTMNAWATGGTFKEIPKSRFC 598

Query: 130 NILIPLPPSSEQKRIVAKIEELFAQLDSITAVL 162
            + IPLPP   QK+IVAKIE     +D   AVL
Sbjct: 599 ELEIPLPPLEVQKKIVAKIEGYQKVIDGARAVL 631
>gi|42528243|ref|NP_973341.1| type I restriction-modification system, S subunit, truncation
           [Treponema denticola ATCC 35405]
 gi|41819513|gb|AAS13260.1| type I restriction-modification system, S subunit, truncation
           [Treponema denticola ATCC 35405]
          Length = 175

 Score = 54.7 bits (130), Expect = 2e-06,   Method: Composition-based stats.
 Identities = 41/159 (25%), Positives = 72/159 (45%), Gaps = 19/159 (11%)

Query: 2   WCRLGEICSITMGQSPESSFISNNSDGMEFHQGK--IHFTEKYIQKANNYTFNITKIAPK 59
           WC  G++  +  G+         N   +E   G+  I+ +   +  AN+Y      I PK
Sbjct: 29  WCHFGDVADVINGK---------NQSQVEDDTGEYPIYGSGGIMGYANDY------ICPK 73

Query: 60  NAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKSTGTT 119
           N  ++  +  +      E +         + P   +   + FY+  C+   F    + TT
Sbjct: 74  NCTIIGRKGSINNPIFVEEKFWNVDTAFGLAPSSIVLPRYLFYF--CKSFDFTSLDSSTT 131

Query: 120 FQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSI 158
             +++   I++IL PLPP   Q+RI+ KI+ELF+QL+ I
Sbjct: 132 LPSLTKTSIRSILFPLPPFVAQQRILDKIDELFSQLEKI 170
>gi|89098144|ref|ZP_01171029.1| type I restriction modification system, subunit S [Bacillus sp.
           NRRL B-14911]
 gi|89087001|gb|EAR66117.1| type I restriction modification system, subunit S [Bacillus sp.
           NRRL B-14911]
          Length = 435

 Score = 54.3 bits (129), Expect = 2e-06,   Method: Composition-based stats.
 Identities = 28/102 (27%), Positives = 57/102 (55%)

Query: 56  IAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQKS 115
           + P  +IL+  RA +G  +I    I   +G  S      + +E+ +Y+L+  K  F   +
Sbjct: 90  LLPAGSILMTSRATIGERSIATAPISTNQGFKSFVCHDGLSNEYMYYYLEILKQYFLLNA 149

Query: 116 TGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDS 157
           +G+TF  +S ++I+N ++ +PP  EQ++IV  +  +  Q+++
Sbjct: 150 SGSTFLEVSKQVIENQVMAIPPHKEQQKIVEVLSTVDEQIEN 191
>gi|118743998|ref|ZP_01591997.1| restriction modification system DNA specificity domain [Geobacter
           lovleyi SZ]
 gi|118683155|gb|EAV89555.1| restriction modification system DNA specificity domain [Geobacter
           lovleyi SZ]
          Length = 422

 Score = 54.3 bits (129), Expect = 2e-06,   Method: Composition-based stats.
 Identities = 30/102 (29%), Positives = 49/102 (48%)

Query: 55  KIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQK 114
           K  P   +LL  RAP+G + I   E+   +G  S  P     S F +Y L+    T  Q 
Sbjct: 293 KKYPAGTVLLSSRAPIGYMAIAREELTTNQGFKSFIPTNDYSSAFIYYTLKNSLKTIVQH 352

Query: 115 STGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLD 156
           ++G+TF+ +S  ++K + I LP S   ++    +   F + D
Sbjct: 353 ASGSTFKEVSGAVLKTVKICLPASGVVEQFTNAVAPTFKRQD 394
>gi|148544101|ref|YP_001271471.1| restriction modification system DNA specificity domain
           [Lactobacillus reuteri F275]
 gi|148531135|gb|ABQ83134.1| restriction modification system DNA specificity domain
           [Lactobacillus reuteri F275]
          Length = 211

 Score = 54.3 bits (129), Expect = 2e-06,   Method: Composition-based stats.
 Identities = 28/101 (27%), Positives = 52/101 (51%)

Query: 55  KIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQK 114
           K+ P N IL   RAP+G ++I +  +   +G  SV P      +F +  L+ +    + +
Sbjct: 89  KLLPTNTILFSSRAPIGYISIAKNNLATNQGFKSVIPNKEYSFQFIYELLKHETAAIKNE 148

Query: 115 STGTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQL 155
           + G+TF+ IS + +K  +I +P S +  +     + +F QL
Sbjct: 149 ANGSTFKEISGKKLKQHIINIPNSEDTSKFNEITKPIFKQL 189
>gi|150006638|ref|YP_001301382.1| type I restriction enzyme EcoR124II specificity protein
           [Bacteroides vulgatus ATCC 8482]
 gi|149935062|gb|ABR41760.1| type I restriction enzyme EcoR124II specificity protein
           [Bacteroides vulgatus ATCC 8482]
          Length = 447

 Score = 53.9 bits (128), Expect = 3e-06,   Method: Composition-based stats.
 Identities = 28/69 (40%), Positives = 45/69 (65%), Gaps = 4/69 (5%)

Query: 93  YRIQSEFWFYWLQCQ---KDTFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKIE 149
           Y+I +++ ++ ++ +   K  ++  S GTT   I+IE  K  ++PLPP +EQ+RIVAKIE
Sbjct: 374 YQINAKYIYHIMRSEYMLKQMYDN-SKGTTVDTITIEKAKQYILPLPPLAEQQRIVAKIE 432

Query: 150 ELFAQLDSI 158
           E F+  D I
Sbjct: 433 ETFSIFDGI 441
>gi|124520753|ref|ZP_01695856.1| restriction modification system DNA specificity domain [Bacillus
           coagulans 36D1]
 gi|124497354|gb|EAY44916.1| restriction modification system DNA specificity domain [Bacillus
           coagulans 36D1]
          Length = 427

 Score = 53.5 bits (127), Expect = 4e-06,   Method: Composition-based stats.
 Identities = 38/159 (23%), Positives = 76/159 (47%), Gaps = 20/159 (12%)

Query: 4   RLGEICSITMGQSPESSFISN-------------NSDGMEFHQGKIHFTEKYIQKANNYT 50
           ++ E+C++  G +P ++   N              S G      K   TEK ++ ++   
Sbjct: 237 KIKEVCNVVSGGTPSTNDKKNWDGNIPWCTPTDITSSGKFIRNTKQTITEKGLKNSS--- 293

Query: 51  FNITKIAPKNAILLCVRAPVGVVNITEREICIGRGLCS-VYPKYRIQSEFWFYWLQCQKD 109
                + PKN+IL+C RA +G  +I   E+   +G  S V  +  +  EF++Y L     
Sbjct: 294 ---ANLLPKNSILMCSRATIGPRSINRVEMATNQGFKSFVCNEEYLDYEFFYYLLSIYIP 350

Query: 110 TFEQKSTGTTFQAISIELIKNILIPLPPSSEQKRIVAKI 148
            F++ ++G+TF  +S + ++N  I +P   ++++ +  I
Sbjct: 351 IFKKLASGSTFLEVSKKDVENTKIFIPKDVKEQKAIGSI 389
>gi|121608003|ref|YP_995810.1| restriction modification system DNA specificity domain
           [Verminephrobacter eiseniae EF01-2]
 gi|121552643|gb|ABM56792.1| restriction modification system DNA specificity domain
           [Verminephrobacter eiseniae EF01-2]
          Length = 575

 Score = 53.5 bits (127), Expect = 4e-06,   Method: Composition-based stats.
 Identities = 42/169 (24%), Positives = 75/169 (44%), Gaps = 14/169 (8%)

Query: 2   WCRLGEICSITMGQSPESSFISN-NSDGMEF-------HQGKIHFTEKYIQKANNYTFNI 53
           W RLG +    +G++P S      + +G  +       H G++  T++ +       F  
Sbjct: 86  WVRLGALLPFRIGKTPASEDPQYWDQEGYAWVSISDMAHLGEVFDTQRKLTARGAQVFGY 145

Query: 54  TKIAPKNAILLCVRAPVGVVNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ 113
             + P   +++  +  +G +++          + S+ P   + ++F  Y L     T   
Sbjct: 146 EPL-PVGTLIMSFKLTIGKISVLRVPAYHNEAIVSLMPLCGLVTDFLKYMLPTVSKTGVS 204

Query: 114 KST--GTTFQAISIELIKNILIPLPPSSEQKRIVAKIEELFAQLDSITA 160
           K    GTT    S+    N+LI LPP+ EQ RIVA++EEL    D++ A
Sbjct: 205 KEALMGTTLNTQSLS---NLLIALPPAVEQSRIVARVEELMRLCDTLEA 250
>gi|152991445|ref|YP_001357167.1| type I restriction-modification system, S subunit [Nitratiruptor
           sp. SB155-2]
 gi|151423306|dbj|BAF70810.1| type I restriction-modification system, S subunit [Nitratiruptor
           sp. SB155-2]
          Length = 373

 Score = 53.1 bits (126), Expect = 6e-06,   Method: Composition-based stats.
 Identities = 39/129 (30%), Positives = 60/129 (46%), Gaps = 7/129 (5%)

Query: 13  MGQSPESSFISNNSDGMEFHQGKIHFTEKYIQKANNYTFNITKIAPKNAILLCVRAPVGV 72
           MG SP+ S  +    G+  +QG   F  ++ +K   +T    + A K   L+ VRAPVG 
Sbjct: 198 MGHSPKGSSFNEYGFGIPMYQGNADFGFRFPKK-RIFTTEPKRFAEKFDTLISVRAPVGE 256

Query: 73  VNITEREICIGRGLCSVYPKYRIQSEFWFYWLQCQKDTFEQ----KSTGTTFQAISIELI 128
            N+   + CIGRGL     +Y++   F+ Y     K    +       GT F +IS    
Sbjct: 257 QNMALEKCCIGRGLARF--RYKLNPNFYSYTYYKLKYLINKIKLFNDEGTVFGSISKGDF 314

Query: 129 KNILIPLPP 137
           + + I +PP
Sbjct: 315 QKLEIMIPP 323
  Database: nr
    Posted date:  Sep 17, 2007 11:41 AM
  Number of letters in database: 999,999,834
  Number of sequences in database:  2,976,859
  
  Database: /nucleus1/users/jsaw/ncbi/db/nr.01
    Posted date:  Sep 17, 2007 11:48 AM
  Number of letters in database: 894,087,890
  Number of sequences in database:  2,493,262
  
Lambda     K      H
   0.322    0.136    0.415 

Lambda     K      H
   0.267   0.0410    0.140 


Matrix: BLOSUM62
Gap Penalties: Existence: 11, Extension: 1
Number of Hits to DB: 568,573,335
Number of Sequences: 5470121
Number of extensions: 21127388
Number of successful extensions: 43255
Number of sequences better than 1.0e-05: 94
Number of HSP's better than  0.0 without gapping: 29
Number of HSP's successfully gapped in prelim test: 65
Number of HSP's that attempted gapping in prelim test: 43084
Number of HSP's gapped (non-prelim): 153
length of query: 162
length of database: 1,894,087,724
effective HSP length: 122
effective length of query: 40
effective length of database: 1,226,732,962
effective search space: 49069318480
effective search space used: 49069318480
T: 11
A: 40
X1: 16 ( 7.4 bits)
X2: 38 (14.6 bits)
X3: 64 (24.7 bits)
S1: 41 (21.9 bits)
S2: 124 (52.4 bits)