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

Query= DDB0218043
         (1139 letters)

Database: /home/dicty1/resource/WorkingDBs//nr-clean
           2,329,665 sequences; 788,375,511 total letters

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

                                                                   Score     E
Sequences producing significant alignments:                        (bits)  Value

gi|57506070|ref|ZP_00371993.1| probable type I restriction enzym...   105  8e-21
gi|551948|gb|AAA98302.1| unknown [Geobacillus stearothermophilus]     103  4e-20
gi|21229078|ref|NP_635000.1| type I restriction-modification sys...   101  1e-19
gi|48895667|ref|ZP_00328651.1| COG0610: Type I site-specific res...   100  4e-19
gi|53718588|ref|YP_107574.1| hypothetical protein BPSL0945 [Burk...    99  1e-18
gi|56459753|ref|YP_155034.1| Type I site-specific restriction-mo...    98  2e-18
gi|53800416|ref|ZP_00359771.1| COG0610: Type I site-specific res...    98  2e-18
gi|15839310|ref|NP_299998.1| type I restriction-modification sys...    97  3e-18
gi|50086398|ref|YP_047908.1| putative type I restriction-modific...    97  4e-18
gi|49243380|emb|CAG41797.1| putative type I restriction enzyme p...    97  4e-18
gi|38505824|ref|NP_942394.1| type I site-specific deoxyribonucle...    96  6e-18
gi|52549654|gb|AAU83503.1| type I site-specific restriction-modi...    94  2e-17
gi|53690938|ref|ZP_00346524.1| COG0610: Type I site-specific res...    94  2e-17
gi|48839712|ref|ZP_00296642.1| COG0610: Type I site-specific res...    94  2e-17
gi|9948809|gb|AAG06120.1| hypothetical protein PA2732 [Pseudomon...    94  3e-17
gi|48868351|ref|ZP_00321696.1| COG0610: Type I site-specific res...    93  4e-17
gi|9656287|gb|AAF94914.1| type I restriction enzyme HsdR, putati...    92  7e-17
gi|41690596|ref|ZP_00147128.1| COG0610: Type I site-specific res...    92  9e-17
gi|52007207|ref|ZP_00334585.1| COG0610: Type I site-specific res...    92  9e-17
gi|21109200|gb|AAM37743.1| type I restriction-modification syste...    92  1e-16
gi|23473873|ref|ZP_00129168.1| COG0610: Type I site-specific res...    92  1e-16
gi|53756220|gb|AAU90511.1| type I restriction-modification syste...    91  2e-16
gi|34762431|ref|ZP_00143431.1| TYPE I RESTRICTION-MODIFICATION S...    91  2e-16
gi|52857360|ref|ZP_00341546.1| COG0610: Type I site-specific res...    91  3e-16
gi|28199930|ref|NP_780244.1| type I restriction-modification sys...    90  3e-16
gi|42632163|ref|ZP_00157701.1| COG0610: Type I site-specific res...    90  3e-16
gi|42629287|ref|ZP_00154835.1| COG0610: Type I site-specific res...    90  3e-16
gi|2411488|gb|AAB70710.1| HsdR [Klebsiella pneumoniae] >gi|74672...    90  3e-16
gi|56750492|ref|YP_171193.1| hypothetical protein syc0483_c [Syn...    90  5e-16
gi|52426226|ref|YP_089363.1| HsdR protein [Mannheimia succinicip...    90  5e-16
gi|51598164|ref|YP_072355.1| possible type I restriction enzyme ...    89  6e-16
gi|48243623|gb|AAT40773.1| putative type I site-specific restric...    88  1e-15
gi|23491885|dbj|BAC16860.1| conserved hypothetical protein [Cory...    88  1e-15
gi|48853355|ref|ZP_00307528.1| COG0610: Type I site-specific res...    87  3e-15
gi|1679868|emb|CAA68059.1| Sty SBLI [Salmonella enterica]              87  4e-15
gi|16124876|ref|NP_419440.1| hypothetical protein CC0623 [Caulob...    86  7e-15
gi|54308080|ref|YP_129100.1| putative type I restriction enzyme ...    82  9e-14
gi|45522101|ref|ZP_00173617.1| COG0610: Type I site-specific res...    82  9e-14
gi|56900920|gb|AAW31753.1| type I restriction enzyme HsdR [Vibri...    81  2e-13
gi|53692041|ref|ZP_00347095.1| COG0610: Type I site-specific res...    76  7e-12
gi|58038318|ref|YP_190287.1| Putative type I restriction enzyme ...    74  2e-11
gi|49257051|dbj|BAD24840.1| type I restriction-modification syst...    74  3e-11
gi|19919302|gb|AAM08232.1| putative type I restriction enzyme [L...    72  1e-10
gi|46141215|ref|ZP_00145477.2| COG0610: Type I site-specific res...    70  5e-10
gi|42525030|ref|NP_970410.1| type I restriction-modification sys...    69  6e-10
gi|15795175|ref|NP_247089.2| Type I restriction-modification enz...    69  1e-09
gi|2826244|gb|AAB98104.1| type I restriction-modification enzyme...    69  1e-09
gi|2129335|pir||D64315 type I restriction enzyme homolog - Metha...    69  1e-09
gi|53688022|ref|ZP_00109780.2| COG0610: Type I site-specific res...    68  1e-09
gi|6968968|emb|CAB73965.1| putative type I restriction enzyme R ...    67  2e-09
gi|19881241|gb|AAM00850.1| HsdR [Campylobacter jejuni]                 67  2e-09
gi|19881235|gb|AAM00845.1| HsdR [Campylobacter jejuni]                 67  2e-09
gi|19881208|gb|AAM00823.1| HsdR [Campylobacter jejuni]                 67  2e-09
gi|19881202|gb|AAM00818.1| HsdR [Campylobacter jejuni]                 67  2e-09
gi|56419917|ref|YP_147235.1| type I restriction-modification sys...    66  5e-09
gi|15669399|ref|NP_248209.1| type I restriction-modification enz...    66  5e-09
gi|2129336|pir||E64451 type I restriction enzyme homolog - Metha...    66  5e-09
gi|46142301|ref|ZP_00204237.1| COG0610: Type I site-specific res...    66  7e-09
gi|53758784|gb|AAU93075.1| type I restriction-modification syste...    65  9e-09
gi|21228306|ref|NP_634228.1| hypothetical protein MM2204 [Methan...    65  9e-09
gi|2622038|gb|AAB85438.1| type I restriction enzyme [Methanother...    65  2e-08
gi|15789431|ref|NP_279255.1| RmeR [Halobacterium salinarum NRC-1...    64  3e-08
gi|46580113|ref|YP_010921.1| type I restriction-modification enz...    64  3e-08
gi|48853287|ref|ZP_00307467.1| COG0610: Type I site-specific res...    62  8e-08
gi|52010778|ref|ZP_00338137.1| COG0610: Type I site-specific res...    62  8e-08
gi|21232336|ref|NP_638253.1| HsdR protein [Xanthomonas campestri...    62  1e-07
gi|24461250|gb|AAN61997.1| HsdR [Lactobacillus sakei] >gi|324555...    62  1e-07
gi|52853569|ref|ZP_00341000.1| COG0610: Type I site-specific res...    62  1e-07
gi|17133741|dbj|BAB76303.1| type I site-specific deoxyribonuclea...    61  2e-07
gi|15839328|ref|NP_300016.1| type I restriction-modification sys...    61  2e-07
gi|58583075|ref|YP_202091.1| type I restriction-modification sys...    61  2e-07
gi|48782405|ref|ZP_00278934.1| COG0610: Type I site-specific res...    60  3e-07
gi|27380123|ref|NP_771652.1| putative type I restriction enzyme ...    60  3e-07
gi|54030757|ref|ZP_00362893.1| COG0610: Type I site-specific res...    60  3e-07
gi|10717104|gb|AAG22017.1| putative HsdR [Streptococcus thermoph...    60  3e-07
gi|48477148|ref|YP_022854.1| type I restriction-modification sys...    59  7e-07
gi|2314575|gb|AAD08445.1| type I restriction enzyme R protein (h...    59  7e-07
gi|28868304|ref|NP_790923.1| type I restriction-modification enz...    59  9e-07
gi|15603406|ref|NP_246480.1| HsdR [Pasteurella multocida subsp. ...    59  9e-07
gi|11499298|ref|NP_070536.1| type I restriction-modification enz...    59  9e-07
gi|47574853|ref|ZP_00244888.1| COG0610: Type I site-specific res...    59  9e-07
gi|47572323|ref|ZP_00242368.1| COG0610: Type I site-specific res...    59  9e-07
gi|49658900|emb|CAF28526.1| hsdr-like Type I restriction enzyme ...    58  2e-06
gi|52082593|ref|YP_081384.1| putative Type I site-specific deoxy...    57  2e-06
gi|15612489|ref|NP_224142.1| TYPE I RESTRICTION ENZYME (RESTRICT...    57  2e-06
gi|30248399|ref|NP_840469.1| Restriction enzymes type I helicase...    57  3e-06
gi|30250446|ref|NP_842516.1| Restriction enzymes type I helicase...    57  4e-06
gi|53757541|gb|AAU91832.1| type I restriction-modification syste...    57  4e-06
gi|2313568|gb|AAD07526.1| type I restriction enzyme R protein (h...    57  4e-06
gi|46156519|ref|ZP_00132562.2| COG0610: Type I site-specific res...    57  4e-06
gi|16273199|ref|NP_439437.1| type I restriction enzyme [Haemophi...    56  6e-06
gi|53795635|ref|ZP_00356694.1| COG0610: Type I site-specific res...    56  6e-06
gi|9789462|gb|AAF98314.1| HsdR [Lactococcus lactis subsp. lactis...    56  7e-06
gi|21226533|ref|NP_632455.1| type I restriction-modification sys...    56  7e-06
gi|2689698|gb|AAB91415.1| restriction subunit [Lactococcus lacti...    56  7e-06
gi|2865242|gb|AAC15896.1| type IC restriction subunit [Lactococc...    56  7e-06
gi|23113394|ref|ZP_00098774.1| COG0610: Type I site-specific res...    56  7e-06
gi|7467215|pir||T30818 hsdR protein - Klebsiella pneumoniae >gi|...    56  7e-06
gi|51893046|ref|YP_075737.1| type I restriction-modification sys...    55  9e-06
gi|28377764|ref|NP_784656.1| type I site-specific deoxyribonucle...    55  9e-06
gi|21228835|ref|NP_634757.1| hypothetical protein MM2733 [Methan...    55  9e-06
gi|6137147|gb|AAF04356.1| type IC restriction subunit [Streptoco...    55  9e-06
gi|23467874|ref|ZP_00123451.1| COG0610: Type I site-specific res...    55  1e-05
gi|46318360|ref|ZP_00218822.1| COG0610: Type I site-specific res...    55  2e-05
gi|42631821|ref|ZP_00157359.1| COG0610: Type I site-specific res...    55  2e-05
gi|1363459|pir||S49394 HsdR1 protein - Mycoplasma pulmonis >gi|4...    54  2e-05
gi|15964350|ref|NP_384703.1| PUTATIVE TYPE I RESTRICTION ENZYME ...    54  2e-05
gi|10954527|ref|NP_044166.1| type I restriction enzyme subunit R...    54  2e-05
gi|53732697|ref|ZP_00154929.2| COG0610: Type I site-specific res...    54  2e-05
gi|53729670|ref|ZP_00149994.2| COG0610: Type I site-specific res...    54  2e-05
gi|60492049|emb|CAH06811.1| putative type I restriction enzyme [...    54  3e-05
gi|48825428|ref|ZP_00286678.1| COG0610: Type I site-specific res...    54  4e-05
gi|15828903|ref|NP_326263.1| RESTRICTION MODIFICATION ENZYME SUB...    53  5e-05
gi|15829148|ref|NP_326508.1| RESTRICTION-MODIFICATION ENZYME SUB...    53  5e-05
gi|48839096|ref|ZP_00296031.1| COG0610: Type I site-specific res...    53  5e-05
gi|32033689|ref|ZP_00133984.1| COG0610: Type I site-specific res...    53  5e-05
gi|24371976|ref|NP_716018.1| type I restriction-modification sys...    53  6e-05
gi|28867251|ref|NP_789870.1| type I site-specific deoxyribonucle...    53  6e-05
gi|50590811|ref|ZP_00332158.1| COG0610: Type I site-specific res...    52  8e-05
gi|59713722|ref|YP_206497.1| type I restriction-modification sys...    52  8e-05
gi|48868625|ref|ZP_00321936.1| COG0610: Type I site-specific res...    52  1e-04
gi|59801121|ref|YP_207833.1| putative type I restriction enzyme ...    52  1e-04
gi|2581809|gb|AAC25971.1| endonuclease (R) subunit homolog [Myco...    52  1e-04
gi|30022537|ref|NP_834168.1| Type I restriction-modification sys...    52  1e-04
gi|17988794|ref|NP_541427.1| TYPE I RESTRICTION-MODIFICATION SYS...    52  1e-04
gi|23464209|gb|AAN34017.1| type I restriction-modification enzym...    52  1e-04
gi|46142217|ref|ZP_00147940.2| COG0610: Type I site-specific res...    52  1e-04
gi|23494172|dbj|BAC19139.1| putative type I restriction-modifica...    52  1e-04
gi|52079178|ref|YP_077969.1| HsdRI [Bacillus licheniformis ATCC ...    51  2e-04
gi|21228394|ref|NP_634316.1| type I restriction-modification sys...    51  2e-04
gi|19881310|gb|AAM00900.1| HsdR [Campylobacter jejuni] >gi|19698...    51  2e-04
gi|53731382|ref|ZP_00148539.2| COG0610: Type I site-specific res...    51  2e-04
gi|57865914|ref|YP_190017.1| type I restriction-modification sys...    50  3e-04
gi|37678997|ref|NP_933606.1| type I site-specific restriction-mo...    50  4e-04
gi|15672632|ref|NP_266806.1| type I restriction enzyme R protein...    50  5e-04
gi|24377275|gb|AAN58610.1| putative type I restriction-modificat...    50  5e-04
gi|22416338|emb|CAC87149.1| restriction-modification enzyme type...    50  5e-04
gi|13357652|ref|NP_077926.1| type I restriction enzyme R protein...    49  7e-04
gi|45522310|ref|ZP_00173825.1| COG0610: Type I site-specific res...    49  7e-04
gi|60492761|emb|CAH07535.1| putative type I restriction enzyme R...    49  0.001
gi|45510657|ref|ZP_00162989.1| COG0610: Type I site-specific res...    49  0.001
gi|1685100|gb|AAC44668.1| HSDR >gi|2126351|pir||JC5216 type I si...    49  0.001
gi|49243514|emb|CAG41938.1| putative type I restriction enzyme [...    48  0.002
gi|14245962|dbj|BAB56357.1| probable type I restriction enzyme r...    48  0.002
gi|57651196|ref|YP_185079.1| type I restriction-modification enz...    48  0.002
gi|49482437|ref|YP_039661.1| putative type I restriction enzyme ...    48  0.002
gi|34557504|ref|NP_907319.1| HSDR PROTEIN PROBABLE TYPE I RESTRI...    48  0.002
gi|42525882|ref|NP_970980.1| type I restriction-modification sys...    48  0.002
gi|15828549|ref|NP_325909.1| RESTRICTION-MODIFICATION ENZYME SUB...    47  0.003
gi|38234847|ref|NP_940614.1| Putative type I restriction/modific...    47  0.003
gi|49237263|ref|ZP_00331318.1| COG0610: Type I site-specific res...    47  0.004
gi|37680388|ref|NP_934997.1| putative type I site-specific restr...    47  0.004
gi|15611851|ref|NP_223502.1| TYPE I RESTRICTION ENZYME (RESTRICT...    47  0.004
gi|1674181|gb|AAB96137.1| Type I Restriction Enzyme (fragment) [...    47  0.004
gi|47459118|ref|YP_015980.1| type I restriction enzyme r protein...    46  0.006
gi|21229939|ref|NP_635856.1| type I restriction-modification sys...    46  0.006
gi|55822674|ref|YP_141115.1| type I restriction-modification sys...    46  0.008
gi|22995100|ref|ZP_00039582.1| COG0610: Type I site-specific res...    46  0.008
gi|3057061|gb|AAC38345.1| HsdR [Lactococcus lactis]                    46  0.008
gi|48832815|ref|ZP_00289843.1| COG0610: Type I site-specific res...    45  0.010
gi|22299773|ref|NP_683020.1| type I site-specific deoxyribonucle...    45  0.010
gi|15839314|ref|NP_300002.1| type I restriction-modification sys...    45  0.013
gi|56476904|ref|YP_158493.1| Type I restriction-modification sys...    45  0.013
gi|53765022|ref|ZP_00161201.2| COG0610: Type I site-specific res...    45  0.013
gi|28896554|ref|NP_802904.1| putative type I site-specific deoxy...    44  0.022
gi|28199934|ref|NP_780248.1| type I restriction-modification sys...    44  0.022
gi|5457726|emb|CAB49217.1| hsdR type I restriction-modification ...    44  0.022
gi|46120969|ref|ZP_00173785.2| COG0610: Type I site-specific res...    44  0.022
gi|56808771|ref|ZP_00366487.1| COG0610: Type I site-specific res...    44  0.029
gi|13622935|gb|AAK34612.1| putative type I site-specific deoxyri...    44  0.029
gi|48857009|ref|ZP_00311166.1| COG0610: Type I site-specific res...    44  0.037
gi|22996439|ref|ZP_00040695.1| COG0610: Type I site-specific res...    44  0.037
gi|45656817|ref|YP_000903.1| type I restriction enzyme [Leptospi...    44  0.037
gi|28262878|gb|EAA26382.1| hypothetical type I restriction enzym...    44  0.037
gi|22995174|ref|ZP_00039655.1| COG0610: Type I site-specific res...    44  0.037
gi|23492094|dbj|BAC17069.1| conserved hypothetical protein [Cory...    44  0.037
gi|22997649|ref|ZP_00041877.1| COG0610: Type I site-specific res...    43  0.049
gi|19746825|ref|NP_607961.1| putative type I site-specific deoxy...    43  0.049
gi|50914975|ref|YP_060947.1| Type I restriction-modification sys...    43  0.049
gi|51594886|ref|YP_069077.1| putative type I restriction enzyme,...    43  0.063
gi|4467404|emb|CAA31543.1| unnamed protein product [Escherichia ...    42  0.083
gi|20090965|ref|NP_617040.1| type I site-specific deoxyribonucle...    42  0.083
gi|21674690|ref|NP_662755.1| type I restriction system endonucle...    42  0.11
gi|15896079|ref|NP_349428.1| N-terminal HKD family nuclease fuse...    42  0.11
gi|32476945|ref|NP_869939.1| HsdR [Rhodopirellula baltica SH 1] ...    42  0.11
gi|29349943|ref|NP_813446.1| Type I restriction enzyme EcoR124II...    42  0.11
gi|31541736|gb|AAP57035.1| HsdR [Mycoplasma gallisepticum R] >gi...    42  0.14
gi|16799596|ref|NP_469864.1| hypothetical protein lin0521 [Liste...    42  0.14
gi|37528152|ref|NP_931497.1| Type I site-specific deoxyribonucle...    42  0.14
gi|38505927|ref|NP_942545.1| type I restriction-modification sys...    42  0.14
gi|15892796|ref|NP_360510.1| similarity to type I restriction en...    42  0.14
gi|46143838|ref|ZP_00133968.2| COG0610: Type I site-specific res...    42  0.14
gi|2313977|gb|AAD07892.1| type I restriction enzyme R protein (h...    41  0.18
gi|7226069|gb|AAF41246.1| type I restriction enzyme EcoR124II R ...    41  0.24
gi|56679311|gb|AAV95977.1| type I restriction-modification syste...    41  0.24
gi|34763948|ref|ZP_00144844.1| TYPE I RESTRICTION-MODIFICATION S...    41  0.24
gi|25532078|pir||AE1870 type I site-specific deoxyribonuclease [...    41  0.24
gi|57505321|ref|ZP_00371250.1| type I restriction enzyme R prote...    40  0.32
gi|48868205|ref|ZP_00321573.1| COG0610: Type I site-specific res...    40  0.32
gi|59800850|ref|YP_207562.1| putative type I site-specific deoxy...    40  0.32
gi|42630520|ref|ZP_00156059.1| COG0610: Type I site-specific res...    40  0.32
gi|46128931|ref|ZP_00154665.2| COG0610: Type I site-specific res...    40  0.32
gi|41720060|ref|ZP_00148899.1| COG4096: Type I site-specific res...    40  0.32
gi|48824781|ref|ZP_00286118.1| COG0610: Type I site-specific res...    40  0.41
gi|21227763|ref|NP_633685.1| type I restriction-modification sys...    40  0.54
gi|48839274|ref|ZP_00296207.1| COG4096: Type I site-specific res...    40  0.54
gi|32262969|gb|AAP78015.1| hypothetical protein HH1418 [Helicoba...    39  0.70
gi|52549588|gb|AAU83437.1| restriction enzymes type I helicase [...    39  0.92
gi|23482959|gb|EAA18785.1| inositol hexakisphosphate kinase [Pla...    39  0.92

>gi|57506070|ref|ZP_00371993.1| probable type I restriction enzyme
           HsdR VC1765 [Campylobacter upsaliensis RM3195]
 gb|EAL52469.1| probable type I restriction enzyme HsdR VC1765 [Campylobacter
           upsaliensis RM3195]
          Length = 980

 Score =  105 bits (262), Expect = 8e-21
 Identities = 80/209 (38%), Positives = 114/209 (54%), Gaps = 32/209 (15%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGR----KGLLIYKKYLDQLISQLPIEE 674
           I+ KA  +I H      K   +  K++AM+V S R    K  L +KKYL +        +
Sbjct: 547 IEKKATIMIEHFFQNIYKKIGE--KAKAMIVTSSRENAVKYYLFFKKYLQENY------Q 598

Query: 675 QFDSIAAFSPFIQEGKSMVESDPMINGK-YAPYGTDKNRGIIKALNDD--KQFNIRLLIV 731
           Q+ ++ AFS            D  I+GK Y+  G +   G  +A+  D  K+ N R LIV
Sbjct: 599 QYQALVAFS-----------GDKEIDGKNYSEVGLN---GFSEAMLKDEFKKDNYRFLIV 644

Query: 732 ADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIF 791
           A K QTGFDEP L  MYVDK+L+G NAVQTLSRL+RI ++KE T ++DF NT +++   F
Sbjct: 645 AQKYQTGFDEPLLHTMYVDKRLEGINAVQTLSRLNRICQNKEDTCVLDFENTHEEIRESF 704

Query: 792 NLYKHSTSLRD---LQDTQTLISSLIMAG 817
           + +   T L++   + +  TL S L   G
Sbjct: 705 STFYGQTYLKEPTNIDNIYTLKSELFDYG 733


 Score =  103 bits (257), Expect = 3e-20
 Identities = 86/291 (29%), Positives = 128/291 (43%), Gaps = 77/291 (26%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQ+D +  +LKD        + +      YLIQH+AGSGKS +IA L + L  L   
Sbjct: 270 PRFHQFDVVRKLLKD--------AKEKGVGQRYLIQHSAGSGKSNSIAWLAHNLASLHRV 321

Query: 331 NKHK----FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNR 386
             HK    FD+I+++ DR  LD Q+ + V  F  +  +      + SKQL+  +  GK +
Sbjct: 322 ENHKEKKIFDSILVVTDRKVLDRQIQENVKSFTENRNLVE-AITEGSKQLKAALKEGK-K 379

Query: 387 VVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKS 446
           +++TT+QKF  +                                        D     ++
Sbjct: 380 IIVTTIQKFPYIA---------------------------------------DEITHLQN 400

Query: 447 NSVAIISDEAHRSHGKSTTEKLHEFLSG--------------------NTRQSSKITYFS 486
            + AII DEAH S      E++ + +S                     N +     +YF+
Sbjct: 401 KTFAIIIDEAHSSQSGKNAEEMGKAISNKDDKETSEMDLEEELIKIIENKKFQKNASYFA 460

Query: 487 FTCTPTTQCLEMFG--IASKGDK--IPFYTYSMHNALEDKLILDVIKQFST 533
           FT TP  + LEMFG      G+K  IPF+ YSM  A+E+  ILDV+K + T
Sbjct: 461 FTATPKPKTLEMFGNPCEINGEKKFIPFHLYSMKQAIEEGFILDVLKGYIT 511


>gi|551948|gb|AAA98302.1| unknown [Geobacillus stearothermophilus]
          Length = 405

 Score =  103 bits (256), Expect = 4e-20
 Identities = 80/277 (28%), Positives = 129/277 (45%), Gaps = 20/277 (7%)

Query: 262 VQNSFKDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLV 321
           V+ S K   PR HQ D +  +++D+         +  +  NYLIQH+AGSGKS +IA L 
Sbjct: 18  VKTSSKLIFPRYHQLDVVRKLVEDVR--------QKGSGENYLIQHSAGSGKSNSIAWLA 69

Query: 322 YFLYKLEI-NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK-SSKQLQNE 379
           Y L  L   +N+  F +++++ DR+ LD QL   +  F  + G+      K +SK L++ 
Sbjct: 70  YHLASLHNEDNESIFTSVIVVTDRTVLDRQLQQTISSFDHTTGLVETIDDKKTSKDLRDA 129

Query: 380 ILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDD 439
           I +GK R++ITT+QKF  +                          +         + E+ 
Sbjct: 130 INNGK-RIIITTLQKFPVIYEEVEVNKGSRFAVIVDEAHSSQTGKSAKKLKAALADTEE- 187

Query: 440 NQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKITYFSFTCTPTTQCLEMF 499
                     A +  E          E + E L+    ++  +++F+FT TP  + LEMF
Sbjct: 188 -----ALREYAELEAEIEAEQLDFEDEIVQELLTHGRHKN--LSFFAFTATPKEKTLEMF 240

Query: 500 GIAS-KGDKIPFYTYSMHNALEDKLILDVIKQFSTLR 535
           G     G   PF+ YSM  A+E+  ILDV++ + T +
Sbjct: 241 GTKQPDGTFKPFHIYSMRQAIEEGFILDVLQNYMTYK 277


>gi|21229078|ref|NP_635000.1| type I restriction-modification system
           restriction subunit [Methanosarcina mazei Go1]
 gb|AAM32672.1| type I restriction-modification system restriction subunit
           [Methanosarcina mazei Goe1]
          Length = 997

 Score =  101 bits (252), Expect = 1e-19
 Identities = 85/286 (29%), Positives = 128/286 (44%), Gaps = 73/286 (25%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQ D +  +L       H +   +  RY  L+QH+AGSGKS +IA L + L  L+ +
Sbjct: 292 PRYHQLDVVRKLLA------HAEWHGAGMRY--LVQHSAGSGKSNSIAWLAHQLISLKKD 343

Query: 331 NKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVIT 390
           +K  FD+I+++ DR  LD Q+ D + +F            + S  L+  I SGK +++I+
Sbjct: 344 DKEIFDSIIVITDRRVLDKQIRDTIKQFAQVGSTVG--HAERSGDLRKFIESGK-KIIIS 400

Query: 391 TMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSVA 450
           T+QKF  +L                                      D+     +    A
Sbjct: 401 TVQKFPFIL--------------------------------------DEIGNEHRGRKFA 422

Query: 451 IISDEAHRSHGKSTTEKLHEFLSGNTRQSS--------------------KITYFSFTCT 490
           II DEAH S G  T   + E LS +  Q                        +YF+FT T
Sbjct: 423 IIIDEAHSSQGGRTAAAMSEALSDSKDQKDPEDTINDALQKRMEARKMLPNASYFAFTAT 482

Query: 491 PTTQCLEMFGIASK-GDKI---PFYTYSMHNALEDKLILDVIKQFS 532
           P  + LE+FG   K GDKI   PF+ Y+M  A+++  ILDV++ ++
Sbjct: 483 PKNKTLELFGEPYKDGDKIKHRPFHNYTMKQAIQEGFILDVLRYYT 528


 Score = 80.5 bits (197), Expect = 3e-13
 Identities = 64/192 (33%), Positives = 93/192 (48%), Gaps = 15/192 (7%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSG-RKGLLIYKKYLDQLISQLPIEEQFD 677
           I+ KA+ I++H        +    ++RAM+V +G  + +  Y    D L   L +   + 
Sbjct: 566 IRLKAEIIVDHFQEQVIALNKIGGQARAMVVTNGIERAIQYYHAIRDYL---LELNSPYK 622

Query: 678 SIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           +I AFS   + G   V ++  +NG  +    DK           ++   R LI ADK QT
Sbjct: 623 AIVAFSGEHEYGGVKV-TEATLNGFPSSQIEDKI----------EEEPYRFLICADKFQT 671

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           G+DEP L  MYVDK L G  AVQTLSRL+R    K    ++DF+N    +   F+ Y  +
Sbjct: 672 GYDEPLLHTMYVDKVLSGIKAVQTLSRLNRAHPKKHDVFVLDFMNDTDTIQDAFSDYYRT 731

Query: 798 TSLRDLQDTQTL 809
           T L D  D   L
Sbjct: 732 TVLSDETDPNKL 743


>gi|48895667|ref|ZP_00328651.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Trichodesmium erythraeum IMS101]
          Length = 991

 Score = 99.8 bits (247), Expect = 4e-19
 Identities = 74/204 (36%), Positives = 108/204 (52%), Gaps = 18/204 (8%)

Query: 607 QYEDETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQL 666
           QY +E   DQ  I+ K++ +I H L           K++AM+V +  K  + YKK  D+ 
Sbjct: 547 QYVEE---DQESIRKKSEVMIEHFLSKVIAQGKINGKAKAMVVSNSIKSAIYYKKAFDKY 603

Query: 667 ISQLPIEEQFDSIAAFSPFIQ-EGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFN 725
           + +   +  + +I AFS   + +GK   E++  +NG  +   T+K        ND K   
Sbjct: 604 LREK--KSDYQTIVAFSGSKEIDGKK--ENESSMNGFSSSKITEK-------FNDSKY-- 650

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
            R LIVA+K QTGFDEP L  MYVDK L    AVQTLSRL+R  + K  T ++DFVN+  
Sbjct: 651 -RFLIVANKYQTGFDEPLLHTMYVDKVLSDVKAVQTLSRLNRSCEGKTDTFVLDFVNSAD 709

Query: 786 DMDFIFNLYKHSTSLRDLQDTQTL 809
           ++   F  Y  +T L +  D+  L
Sbjct: 710 EIQRAFEPYYKTTILSEETDSDRL 733


 Score = 94.4 bits (233), Expect = 2e-17
 Identities = 78/287 (27%), Positives = 129/287 (44%), Gaps = 72/287 (25%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLE-- 328
           PR HQ D +  +L         K     NRY  LIQH+AGSGKS +I  L + L +L+  
Sbjct: 285 PRYHQLDLVKQLLTSA------KKHGVGNRY--LIQHSAGSGKSNSITWLSHQLVELKNI 336

Query: 329 INNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVV 388
              ++ FD+++++ DR  LD Q+ + + +F   + +      K+SK+L++  L  K +++
Sbjct: 337 TEKENIFDSVLVVTDRKILDKQIRENIQQFAQEDKVVE--ATKNSKKLKS-ALENKRKII 393

Query: 389 ITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNS 448
           ITT+QKF  ++                                           S   + 
Sbjct: 394 ITTVQKFPYVVKEIQ---------------------------------------SLSDHK 414

Query: 449 VAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKIT------------------YFSFTCT 490
            AII DEAH S    +   + E LS    +  + T                  Y++FT T
Sbjct: 415 FAIIIDEAHSSQTGKSAASMSESLSKKDSEVEETTEDKIIRIIESQKLCPNANYYAFTAT 474

Query: 491 PTTQCLEMFGIAS--KGDKIPFYTYSMHNALEDKLILDVIKQFSTLR 535
           P  + LE+FG+ +   G   PF++YSM  A+E+  IL+V++ ++T +
Sbjct: 475 PKNKTLELFGVKNPEDGKFYPFHSYSMKQAIEEGFILNVLQHYTTYK 521


>gi|53718588|ref|YP_107574.1| hypothetical protein BPSL0945
           [Burkholderia pseudomallei K96243]
 emb|CAH34941.1| conserved hypothetical protein [Burkholderia pseudomallei K96243]
          Length = 1037

 Score = 98.6 bits (244), Expect = 1e-18
 Identities = 80/293 (27%), Positives = 126/293 (42%), Gaps = 71/293 (24%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQ D    ++ D+++          +   YLIQH+AGSGK+ +IA   +FL  L  +
Sbjct: 279 PRYHQLDATRRLVADVLA--------QGSGQRYLIQHSAGSGKTNSIAWTAHFLADLHDD 330

Query: 331 NKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNG-ITNFCRPKSSKQLQ-NEILSGKNRV 387
           N HK FD++++++DR+ LDAQL + +  F  + G +    +   SK  Q  + L    ++
Sbjct: 331 NHHKLFDSVLVVSDRTVLDAQLQEAIFDFERTTGVVATITKEHGSKSAQLGQALKDGKKI 390

Query: 388 VITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSN 447
           ++ T+Q F   L                                             +  
Sbjct: 391 IVCTIQTFPFALQAVQELAAT------------------------------------EGK 414

Query: 448 SVAIISDEAHRSHGKSTTEKLHEFLSG----------------------NTRQSSK-ITY 484
             A+I+DEAH S       KL + LS                         R  +K +TY
Sbjct: 415 RFAVIADEAHCSQTGEAAAKLKQLLSAEEWADLQDGGEIDTEALLAAQMEARTGAKGLTY 474

Query: 485 FSFTCTPTTQCLEMFG-IASKGDKIPFYTYSMHNALEDKLILDVIKQFSTLRL 536
            +FT TP  + LE+FG     G   PF+ YSM  A+E+  ILDV+K +++ +L
Sbjct: 475 VAFTATPKQKTLELFGRPGPDGLPEPFHVYSMRQAIEEGFILDVLKNYTSYKL 527


 Score = 82.8 bits (203), Expect = 6e-14
 Identities = 71/224 (31%), Positives = 103/224 (45%), Gaps = 34/224 (15%)

Query: 599 NGIDFDDDQYEDETIQDQML---------IKTKAQFIINHLLYTKNKHDNKLFKSRAMLV 649
           NG ++D+ Q E       ++         I  K Q ++ H  Y +N       K++AM+V
Sbjct: 534 NGQEYDEKQVERSAAMKGIMQWVRLHPYNIAQKVQIVVEH--YRENVQPLLDGKAKAMVV 591

Query: 650 CSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMV--------ESDPMING 701
            + RK  + ++K +   I +         + AFS  I + +S          E +P +NG
Sbjct: 592 VASRKEAVRWQKAIRAYIEKR--NYPLGVLVAFSGEIDDPESYPAPVTEASKELNPGLNG 649

Query: 702 KYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQT 761
                     R I  A  + +     LL+VA+K QTGFD+P L  MYVDK L G  AVQT
Sbjct: 650 ----------RDIRDAFAEPEY---HLLLVANKFQTGFDQPLLCGMYVDKMLGGIQAVQT 696

Query: 762 LSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQD 805
           LSRL+R    K+ T I+DFVN   D+   F  Y  +  L    D
Sbjct: 697 LSRLNRAHPGKDTTYILDFVNEAGDILKAFKTYYETAELEATTD 740


>gi|56459753|ref|YP_155034.1| Type I site-specific
           restriction-modification system, R subunit (helicase)
           [Idiomarina loihiensis L2TR]
 gb|AAV81485.1| Type I site-specific restriction-modification system, R subunit
           (helicase) [Idiomarina loihiensis L2TR]
          Length = 1058

 Score = 97.8 bits (242), Expect = 2e-18
 Identities = 83/295 (28%), Positives = 132/295 (44%), Gaps = 76/295 (25%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQWD    V+  LI     +   S+    YLIQH+AGSGKS +IA   + L  L  +
Sbjct: 289 PRYHQWD----VVNQLIDTTRREGPGSK----YLIQHSAGSGKSNSIAWTAHQLSSLYAD 340

Query: 331 NKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNG----ITNFCRPKSSKQLQNEILSGKN 385
           +  K F++IV++ DR+ LD+QL + + +F  ++G    IT     KS  +   E L  + 
Sbjct: 341 DGSKLFNSIVVVTDRTVLDSQLQNTIRQFEHADGVVRAITRDVGNKSKSEQLAEALQEQV 400

Query: 386 RVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFK 445
           R+++ T+Q F  L                                     +  D   + +
Sbjct: 401 RIIVVTIQTFPALY------------------------------------DAIDKYPALQ 424

Query: 446 SNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQ-------------------SSKITYFS 486
           + + AII+DEAH S   ++  KL   L  +  +                   + +I+Y++
Sbjct: 425 AGNYAIIADEAHSSQTGASASKLKALLGSDRPEGDDISAEELLDAAVQARMPNERISYYA 484

Query: 487 FTCTPTTQCLEMFG--------IASKGDKIPFYTYSMHNALEDKLILDVIKQFST 533
           FT TP  + LE+FG        I+       F+ YSM  A+E+  ILDV+K ++T
Sbjct: 485 FTATPKGKTLELFGRPDDPHSPISDDNKPKAFHVYSMRQAIEEGFILDVLKNYTT 539


 Score = 89.7 bits (221), Expect = 5e-16
 Identities = 69/212 (32%), Positives = 106/212 (49%), Gaps = 20/212 (9%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I  + + I+ H  + +N       +++AM+V   R+  + Y+  +   I    IE+ +  
Sbjct: 576 INQRVEIIVEH--FRENVRHLLNGQAKAMVVTGSRQEAVRYQLAMTSYI----IEKGYQD 629

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGII-----KALNDDKQFNIRLLIVAD 733
           + A   F  E    V +D  I      Y    N G+      KAL+ D  FN+  +IVA+
Sbjct: 630 VHALVAFSGE----VPADDYITESVTEYDKKLNPGLNGRKHDKALDSD-DFNV--MIVAN 682

Query: 734 KLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNL 793
           K QTGFD+P L  MYVDKKL+G   VQTLSRL+R    K+ T ++DF N+E+D+   F  
Sbjct: 683 KYQTGFDQPKLCAMYVDKKLQGVECVQTLSRLNRTFPGKDQTFVLDFFNSEEDILESFKP 742

Query: 794 YKHSTSLRDLQDTQTLISSLIMAGRSIQKVFE 825
           Y  +  L  + D   L    IM+    + ++E
Sbjct: 743 YYTTAELDSVTDPHILFD--IMSKLDEEHIYE 772


>gi|53800416|ref|ZP_00359771.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Xylella fastidiosa Dixon]
          Length = 1023

 Score = 97.8 bits (242), Expect = 2e-18
 Identities = 80/300 (26%), Positives = 133/300 (43%), Gaps = 84/300 (28%)

Query: 267 KDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYK 326
           K   PR HQ  T+ ++L+   +    K         YLIQH+AGSGKS TIA L + L +
Sbjct: 290 KQIFPRFHQLRTVRALLRRTRTDGVGK--------RYLIQHSAGSGKSNTIAWLAHQLVE 341

Query: 327 LEINN---KHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSG 383
           L   +   + +FD+I+++ DR  LD Q+ D++  +   +    F    ++++L+  +  G
Sbjct: 342 LRRKDDPLRVQFDSIIVITDRRTLDKQIADIIKGY--DHVAAIFGHSDNAQELREYLRRG 399

Query: 384 KNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFS 443
           K ++++TT+QKF  +L                                       D    
Sbjct: 400 K-KIIVTTVQKFPFIL---------------------------------------DELGD 419

Query: 444 FKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------------------------ 479
                 A++ DEAH S G  TT ++HE L G   +                         
Sbjct: 420 LSGKRFALLIDEAHSSQGGKTTARMHEALGGKVAEEEFEEDSTQDAVNAEIERRIASRKL 479

Query: 480 -SKITYFSFTCTPTTQCLEMFGIAS-KGDKIPF-----YTYSMHNALEDKLILDVIKQFS 532
            +  +YF+FT TP  + LE+FG  +  GDK+ F      TY+   A+++K ILDV+++++
Sbjct: 480 LANASYFAFTATPKNKTLELFGERTLVGDKVQFRSPEELTYTTKQAIQEKFILDVVEKYT 539


 Score = 80.5 bits (197), Expect = 3e-13
 Identities = 61/181 (33%), Positives = 94/181 (51%), Gaps = 15/181 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           ++RAM+VC+G    + Y + +   ++Q  I+  + +
Sbjct: 577 IRRKAEIMVDHFIACVAGKQIIGGQARAMIVCNGIARAIDYWREVSDYLTQ--IKSPYKA 634

Query: 679 IAAFS-PFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           I A+S  F   G+   E+D  +NG         ++ I + L  D     R LIVA+K  T
Sbjct: 635 IVAYSGSFEIGGQKKTEAD--LNGF-------PSKKIPENLKKDPY---RFLIVANKFVT 682

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           GFDEP L  MYVDK L G  AVQTLSRL+R    K  T ++DF +  + +   F  Y  +
Sbjct: 683 GFDEPLLHTMYVDKPLAGVLAVQTLSRLNRAHPQKRDTFVLDFADNAEAVKAAFQDYYRA 742

Query: 798 T 798
           T
Sbjct: 743 T 743


>gi|15839310|ref|NP_299998.1| type I restriction-modification system
           endonuclease [Xylella fastidiosa 9a5c]
 gb|AAF85757.1| type I restriction-modification system endonuclease [Xylella
           fastidiosa 9a5c]
 pir||A82516 type I restriction-modification system endonuclease XF2721
           [imported] - Xylella fastidiosa (strain 9a5c)
          Length = 1025

 Score = 97.1 bits (240), Expect = 3e-18
 Identities = 78/300 (26%), Positives = 127/300 (42%), Gaps = 84/300 (28%)

Query: 267 KDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYK 326
           K   PR HQ  T+ ++L+        ++        YLIQH+AGSGKS TIA L + L +
Sbjct: 292 KQIFPRFHQLRTVRALLR--------RAHSDGVGRRYLIQHSAGSGKSNTIAWLAHQLVE 343

Query: 327 LEINN---KHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSG 383
           L   +     +FD+I+++ DR  LD Q+ D +  +   + +        + Q   E L  
Sbjct: 344 LRRKDDPMTPQFDSIIVITDRRALDTQIADTIKGY---DHVAAIFGHSDNAQELREYLRR 400

Query: 384 KNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFS 443
             ++++TT+QKF  +L                                       D    
Sbjct: 401 DKKIIVTTVQKFPFIL---------------------------------------DELGD 421

Query: 444 FKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------------------------ 479
                 A++ DEAH S G  TT ++HE L G   +                         
Sbjct: 422 LSGKRFALLIDEAHSSQGGKTTARMHEALGGKVAEEEFEEDRTQDAVNAEIERRIASRKL 481

Query: 480 -SKITYFSFTCTPTTQCLEMFGIASK-GDKIPF-----YTYSMHNALEDKLILDVIKQFS 532
            +  +YF+FT TP  + LE+FG  ++ GDK+ F      TY+   A+++K ILDV+++++
Sbjct: 482 LANASYFAFTATPKNKTLELFGEKTRVGDKVQFRSPEELTYTTKQAIQEKFILDVVEKYT 541


 Score = 80.1 bits (196), Expect = 4e-13
 Identities = 61/181 (33%), Positives = 92/181 (50%), Gaps = 15/181 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           ++RAM+VC+G    + Y + +   ++Q  I+  + +
Sbjct: 579 IRRKAEIMVDHFIAEVVGKQKIGGQARAMIVCNGIARAIDYWREVSDYLTQ--IKSPYRA 636

Query: 679 IAAFS-PFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           I A+S  F   G+   E+D  +NG         ++ I   L  D     R LIVA+K  T
Sbjct: 637 IVAYSGSFEIGGQKKTEAD--LNGF-------PSKDIPDTLKKDPY---RFLIVANKFVT 684

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           GFDEP L  MYVDK L G  AVQTLSRL+R    K  T ++DF +  + +   F  Y   
Sbjct: 685 GFDEPLLHTMYVDKPLAGVLAVQTLSRLNRAHPQKHDTFVLDFADNAEAVKAAFQDYYRV 744

Query: 798 T 798
           T
Sbjct: 745 T 745


>gi|50086398|ref|YP_047908.1| putative type I
           restriction-modification system (HsdR) [Acinetobacter
           sp. ADP1]
 emb|CAG70086.1| putative type I restriction-modification system (HsdR)
           [Acinetobacter sp. ADP1]
          Length = 1020

 Score = 96.7 bits (239), Expect = 4e-18
 Identities = 74/271 (27%), Positives = 129/271 (47%), Gaps = 18/271 (6%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL-EI 329
           PR HQ D +  + K      H K+  + + +NYLIQH+AGSGKS TIA L + L  +   
Sbjct: 282 PRYHQLDAVRKLTK------HSKA--NGSGHNYLIQHSAGSGKSNTIAWLAHQLASMHNA 333

Query: 330 NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVI 389
           ++K  F++I+++ DR  LD QL + + +F   +G+      +++ QL N + S    ++I
Sbjct: 334 DDKKIFNSIIVITDRIVLDRQLQETIAQFEHKDGVVQKI-DENTLQLTNALASNV-PIII 391

Query: 390 TTMQKFSQLLPN-----KSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSF 444
           TT+QKF  ++ +     K                     +++          +  N+   
Sbjct: 392 TTIQKFPYIMQSIRTQAKKGIKVDLSTEGKNFAVIVDEAHSSQSGETAMELRKILNKDGI 451

Query: 445 KSNSVAIISDEAHRSHGKSTTE-KLHEFLSGNTR-QSSKITYFSFTCTPTTQCLEMFGIA 502
           +S       D+      + T + K H F   + R +   +++F+FT TP  + L  F  A
Sbjct: 452 ESAIAQEFLDDGDEDESELTDDVKKHLFAEASKRSRQPNLSFFAFTATPKWKTLACFDEA 511

Query: 503 SKGDKIPFYTYSMHNALEDKLILDVIKQFST 533
               + PF+ YSM  A+++  ILDV+  ++T
Sbjct: 512 GDNGEAPFHHYSMKQAIQEGFILDVLANYTT 542


 Score = 88.6 bits (218), Expect = 1e-15
 Identities = 63/171 (36%), Positives = 94/171 (54%), Gaps = 20/171 (11%)

Query: 643 KSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVE----SDPM 698
           +++AM+V + R+  + YK   D+ I     E+++  I +   F   GK  V+    ++P 
Sbjct: 601 RAKAMVVTNSREQAVRYKLAFDEYIK----EKKYTGIKSLVAF--SGKLEVDGTEYTEPK 654

Query: 699 INGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANA 758
           +NG       D+         DD Q    +L+VA+K QTGFD+P L  M+VDKKL G  A
Sbjct: 655 MNGFKETELPDQFD------TDDYQ----VLLVAEKYQTGFDQPLLHTMFVDKKLSGIQA 704

Query: 759 VQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDTQTL 809
           VQTLSRL+R +  KE T ++DFVNT +D+   F  +   T L D+ + Q L
Sbjct: 705 VQTLSRLNRCAYGKEDTFVLDFVNTLEDIYKAFKPFYEQTKLGDIPNEQKL 755


>gi|49243380|emb|CAG41797.1| putative type I restriction enzyme
           protein [Staphylococcus aureus subsp. aureus MSSA476]
 ref|YP_042158.1| putative type I restriction enzyme protein [Staphylococcus aureus
           subsp. aureus MSSA476]
          Length = 1012

 Score = 96.7 bits (239), Expect = 4e-18
 Identities = 80/296 (27%), Positives = 125/296 (42%), Gaps = 80/296 (27%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQ D +  ++ D    QH    K     NYLIQH+AGSGKS +IA L Y L  L  N
Sbjct: 275 PRYHQLDVVKKLISDT---QHHGPGK-----NYLIQHSAGSGKSNSIAWLAYHLASLHSN 326

Query: 331 NKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVI 389
                + +++++ DR+ LD QL + ++ F  + G         + Q     ++ + +++I
Sbjct: 327 EDQPIYSSVIVVTDRTVLDRQLQNTLMSFEHTVGQVETIGDNKTSQDLKTAINDQKKIII 386

Query: 390 TTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSV 449
           TT+QKF  +                                     +E DN       S 
Sbjct: 387 TTLQKFPVIY------------------------------------DEIDNT---AGRSF 407

Query: 450 AIISDEAHRSHGKSTTEKLHEFLSG-------------------------------NTRQ 478
           AI+ DEAH S   S+ +KL   L+                                +  Q
Sbjct: 408 AILVDEAHSSQTGSSAQKLKHALADKEASLKEYQEFEESVEDQSLDDQDQIVQTLLSQGQ 467

Query: 479 SSKITYFSFTCTPTTQCLEMFGI-ASKGDKIPFYTYSMHNALEDKLILDVIKQFST 533
              +++F+FT TP  + +EMFG   + G   PF+ YSM  A+E+  IL+V++ + T
Sbjct: 468 HDNLSFFAFTATPKEKTIEMFGTKQADGSYHPFHIYSMRQAIEEGFILNVLQNYMT 523


 Score = 89.7 bits (221), Expect = 5e-16
 Identities = 65/195 (33%), Positives = 107/195 (54%), Gaps = 13/195 (6%)

Query: 643 KSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDPMINGK 702
           + +AM+V + R   + Y +   + I +    +  D + AFS  +++   +  ++P +N  
Sbjct: 582 RGKAMIVTASRLHAVRYMQEFKKYIKEKGYTDM-DVLVAFSGTVEDD-GIEYTEPKMNLT 639

Query: 703 YAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTL 762
           ++     +N+  +K     + FN  LLIVA+K QTGFDEP L  M+VDK+L+G  AVQTL
Sbjct: 640 HSGERIQENQ--LKEAFHSEDFN--LLIVAEKYQTGFDEPLLHTMFVDKRLRGVKAVQTL 695

Query: 763 SRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRD------LQDTQTLISSL-IM 815
           SRL+R    K  T I+DFVNT +++   F  +  +TSL++      + DTQ+ +    I 
Sbjct: 696 SRLNRTMPGKTDTFILDFVNTAEEIYEAFQPFYEATSLKEEININLIYDTQSKLRKFHIY 755

Query: 816 AGRSIQKVFEKEKRA 830
               IQ V +  ++A
Sbjct: 756 NDDDIQTVMKLNRQA 770


>gi|38505824|ref|NP_942394.1| type I site-specific deoxyribonuclease
           chain R [Synechocystis sp. PCC 6803]
 dbj|BAD02057.1| type I site-specific deoxyribonuclease chain R [Synechocystis sp.
           PCC 6803]
          Length = 989

 Score = 95.9 bits (237), Expect = 6e-18
 Identities = 83/299 (27%), Positives = 132/299 (43%), Gaps = 78/299 (26%)

Query: 263 QNSFKDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVY 322
           +   K   PR HQ    ++V++ L++              YLIQH+AGSGKS +IA L +
Sbjct: 266 RKKLKQVFPRYHQ----LTVVRQLLATAKQDGVGKR----YLIQHSAGSGKSNSIAWLAH 317

Query: 323 FLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILS 382
            L  LE + K  FD+I+++ DR  LD Q+ + + +F   +        + S  L+  I  
Sbjct: 318 QLVGLEQDGKALFDSIIVVTDRQVLDRQIRNTIKQFAQVSATVG--HAERSGDLRQFIKD 375

Query: 383 GKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQF 442
           GK +V+ITT+QKF  +L                                    NE  ++ 
Sbjct: 376 GK-KVIITTVQKFPFIL------------------------------------NEIGDE- 397

Query: 443 SFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNT------------------------RQ 478
             + +  AI+ DEAH S G  TT  ++  L G                          + 
Sbjct: 398 -HRQSKFAIVIDEAHSSQGGKTTAAMNRVLEGTASYNVANEEEEETTEDKINRILEGRKM 456

Query: 479 SSKITYFSFTCTPTTQCLEMFG-IASKGD----KIPFYTYSMHNALEDKLILDVIKQFS 532
            +  +YF+FT TP  + LE+FG  A + D      PF++Y+M  A+++  ILDV+K ++
Sbjct: 457 VTNASYFAFTATPKNKTLEIFGQPAPQPDGTVKHYPFHSYTMKQAIQEGFILDVLKSYT 515


 Score = 82.4 bits (202), Expect = 7e-14
 Identities = 63/195 (32%), Positives = 96/195 (48%), Gaps = 15/195 (7%)

Query: 616 QMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSG-RKGLLIYKKYLDQLISQLPIEE 674
           Q  I+ K + +++H       H     ++RAM++ +G  + +  +  + D L  +   + 
Sbjct: 550 QHAIREKTEIMVDHFHTQVINHGKIGGQARAMVITNGIGQAIQYFYAFKDYLRER---KS 606

Query: 675 QFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADK 734
            + +I AFS     G   V ++  +NG  +   TDK       + +D     R LIVADK
Sbjct: 607 PYQAIVAFSGEYDYGGQKV-TEATLNGFPSSQITDK-------IEEDPY---RFLIVADK 655

Query: 735 LQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLY 794
            QTG+D+P L  MYVDK L G  AVQTLSRL+R    K  T ++DF N    +   F+ Y
Sbjct: 656 YQTGYDQPLLHTMYVDKALSGIKAVQTLSRLNRAHPQKYDTFVLDFYNDSGTIQSSFDPY 715

Query: 795 KHSTSLRDLQDTQTL 809
             +T L +  D   L
Sbjct: 716 YRTTILSNETDPNKL 730


>gi|52549654|gb|AAU83503.1| type I site-specific
           restriction-modification system R subunit and related
           helicases [uncultured archaeon GZfos29E12]
          Length = 889

 Score = 94.4 bits (233), Expect = 2e-17
 Identities = 70/214 (32%), Positives = 110/214 (50%), Gaps = 29/214 (13%)

Query: 604 DDDQYEDETIQ---------DQMLIKTKAQFIINHLLYTKNKHDNKLF-KSRAMLVCSGR 653
           DD +Y+++  +           + I  K   ++NH +   N   NK+  K++AM+V   R
Sbjct: 430 DDPEYDEKKAKRLLRTFVEKHPVAIARKTDIMLNHFM---NSTINKINGKAKAMVVTRSR 486

Query: 654 KGLLIYKKYLDQLISQ--LPIEEQFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKN 711
              ++YKK  D+LI +   PI+    ++ AF+        +V+ D     + +       
Sbjct: 487 LHAVLYKKAFDKLIKENNYPIK----TLVAFT-------GVVKHDEQEYTENSRNDLPSQ 535

Query: 712 RGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKD 771
           + I KA  +D+    ++LIVA+K QTGFD+P L  MYVDK L G  AVQTL RL+RI  +
Sbjct: 536 KSIRKAFEEDQY---KILIVANKFQTGFDQPLLHTMYVDKMLNGITAVQTLCRLNRIHTN 592

Query: 772 KEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQD 805
           K+ T I+DF N  + +   F  Y  +T L +  D
Sbjct: 593 KKDTLILDFANKTEVIQKAFQPYYETTFLEEATD 626


 Score = 90.5 bits (223), Expect = 3e-16
 Identities = 77/294 (26%), Positives = 121/294 (40%), Gaps = 80/294 (27%)

Query: 267 KDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLY- 325
           K   PR HQ D    +L D           ++   NYLI+H+ GSGK+ TIA L + L  
Sbjct: 173 KPIFPRYHQLDCANLLLAD-----------AKPGKNYLIEHSNGSGKTKTIAWLAHGLIN 221

Query: 326 KLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFC-RPKSSKQLQNEILSGK 384
           K +  +   +D +++++DR  +D QL +         G      + K+S  L+  + +G 
Sbjct: 222 KFDSLDNRVYDMVIVVSDRRVIDTQLQNQAKSIEKVKGTVEVIDKKKTSDDLRKALETGS 281

Query: 385 NRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSF 444
           N +V+TT+QKF  +L                                       +     
Sbjct: 282 N-IVVTTLQKFPYIL---------------------------------------EEVKDL 301

Query: 445 KSNSVAIISDEAHRSHGKSTTEKLHEFLSGNT--------------------------RQ 478
                A+I DEAH S   +   K+ + LS N+                          R 
Sbjct: 302 PKRKYAVIIDEAHSSQTGTMARKMKQVLSTNSLEEAEALDIEEMDEVDEELLREIESYRN 361

Query: 479 SSKITYFSFTCTPTTQCLEMFGIASK-GDKIPFYTYSMHNALEDKLILDVIKQF 531
              I++F+FT TP  + LEMFG+  + G   PF+TY+M  A+E+  I DV+K +
Sbjct: 362 LKNISFFAFTATPKNKTLEMFGMLDEYGHYYPFHTYTMKQAIEEGFIFDVLKNY 415


>gi|53690938|ref|ZP_00346524.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Desulfovibrio desulfuricans G20]
          Length = 810

 Score = 94.4 bits (233), Expect = 2e-17
 Identities = 79/297 (26%), Positives = 130/297 (43%), Gaps = 80/297 (26%)

Query: 266 FKDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLY 325
           +K   PR HQ   +  +L +        + +S     YLIQH+AGSGKS +IA L + L 
Sbjct: 269 YKQIFPRYHQLKVVRMLLAN--------AAESGIGRRYLIQHSAGSGKSNSIAWLAHQLV 320

Query: 326 KLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKN 385
            LE  +K  FD+++++ DR  LD Q+ D + +F   +        + S  L+  + +GK 
Sbjct: 321 GLEHESKALFDSVIVVTDRRVLDKQIRDTIKQFAQVSATVG--HAEHSGDLRKFLKAGK- 377

Query: 386 RVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFK 445
           +++ITT+QKF  +L                                      D+     +
Sbjct: 378 KIIITTVQKFPFIL--------------------------------------DEIGDEHR 399

Query: 446 SNSVAIISDEAHRSHGKSTT-------------------------EKLHEFLSGNTRQSS 480
            +  AII DEAH S G  TT                         +K+++ + G  +  +
Sbjct: 400 KSKFAIIIDEAHSSQGGKTTAAMNRVLEETAPYGGSDDEEEETVEDKINKIMEGR-KMVT 458

Query: 481 KITYFSFTCTPTTQCLEMFGIASKG-----DKIPFYTYSMHNALEDKLILDVIKQFS 532
             +YF+FT TP  + LE+FG             PF++Y+M  A+++  ILDV+K ++
Sbjct: 459 NASYFAFTATPKNKTLEIFGEPDPQPDGTVKHHPFHSYTMKQAIQEGFILDVLKNYT 515


>gi|48839712|ref|ZP_00296642.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Methanosarcina barkeri str.
           fusaro]
          Length = 965

 Score = 94.0 bits (232), Expect = 2e-17
 Identities = 75/268 (27%), Positives = 119/268 (43%), Gaps = 22/268 (8%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL-EI 329
           PR HQ  ++    K+L++  H K   +   YNYLIQH+AGSGKS TIA L   L  L + 
Sbjct: 246 PRYHQLYSV----KNLVT--HAKYHGTG--YNYLIQHSAGSGKSNTIAWLCLQLCGLHDT 297

Query: 330 NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVI 389
            NK  FD+I+++ DR  LD QL + + +F    G     +   +K L   +  GK+ ++I
Sbjct: 298 ENKRVFDSIIVVTDRRVLDRQLQETITQFEHVKGTVATIKKDKAKSLTAALQEGKD-III 356

Query: 390 TTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSV 449
            T+Q F                              +           D  +    S  +
Sbjct: 357 CTLQTFP-----------FAVNAISEMSGKHFAVVVDEAHSSQSGEGADSVKKILGSADL 405

Query: 450 AIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKITYFSFTCTPTTQCLEMFGIAS-KGDKI 508
                E      +    +  E    +  +   +++F+FT TP T+ LE+FG     G   
Sbjct: 406 ESDEQEEEPEDDEDLVNQRVEAELKHKGRLPNVSFFAFTATPKTKTLELFGTQQPDGSYE 465

Query: 509 PFYTYSMHNALEDKLILDVIKQFSTLRL 536
           PF  Y+M  A+E+K ILDV++ ++T ++
Sbjct: 466 PFSLYTMKQAIEEKFILDVLENYTTFKV 493


 Score = 92.8 bits (229), Expect = 5e-17
 Identities = 63/168 (37%), Positives = 91/168 (53%), Gaps = 13/168 (7%)

Query: 643 KSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVE-SDPMING 701
           +++AMLV   R   + YK   D+ I     E  F ++ AFS  +Q+G++  E ++ ++NG
Sbjct: 549 QAKAMLVTKSRASAVKYKLAFDRYIKTH--EYPFKALVAFSGTVQDGETGREFTEAIMNG 606

Query: 702 KYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQT 761
               +   +   + K      +   R+LIVA+K QTGFD+P L  MYVDK+L G NAVQT
Sbjct: 607 ----FSESQTAEVFK------REEYRILIVANKFQTGFDQPLLYAMYVDKRLGGVNAVQT 656

Query: 762 LSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDTQTL 809
           LSRL+R   DKE T ++DF N    +   F  Y   T L +  D   L
Sbjct: 657 LSRLNRKYPDKEDTMVLDFENEADLIQRSFQPYYEKTLLTEATDPNKL 704


>gi|9948809|gb|AAG06120.1| hypothetical protein PA2732 [Pseudomonas
           aeruginosa PAO1]
 pir||C83304 hypothetical protein PA2732 [imported] - Pseudomonas aeruginosa
           (strain PAO1)
 ref|NP_251422.1| hypothetical protein PA2732 [Pseudomonas aeruginosa PAO1]
          Length = 1146

 Score = 93.6 bits (231), Expect = 3e-17
 Identities = 86/308 (27%), Positives = 136/308 (43%), Gaps = 84/308 (27%)

Query: 261 RVQNSFKDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACL 320
           R +   K   PR HQ  T+ ++L+   S    K         YLIQH+AGSGKS TIA L
Sbjct: 287 RKKKKRKPVFPRFHQLRTVRALLRRTRSDGVGK--------RYLIQHSAGSGKSNTIAWL 338

Query: 321 VYFLYKLEI---NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQ 377
            + L +L         +FD+I+++ DR  LD Q+   +  +   +  + F    ++++L+
Sbjct: 339 AHQLVELRSVADPLMPQFDSIIVITDRRALDTQIARTIKGY--DHVASIFGHSDNAQELR 396

Query: 378 NEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
           + +  G+ ++++TT+QKF  +L                                     E
Sbjct: 397 DYLRRGR-KIIVTTVQKFPFIL-------------------------------------E 418

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNT--------------------- 476
           +    S    + A++ DEAH S G  TT ++HE LSG T                     
Sbjct: 419 EIGDLS--GRTFALLIDEAHSSQGGKTTARMHEALSGKTDEEAFEEDATQDAVNTEIEKR 476

Query: 477 ----RQSSKITYFSFTCTPTTQCLEMFGIAS-KGDKIPF-----YTYSMHNALEDKLILD 526
               R  S  +YF+FT TP  + LE+FG  +   DK+ F      TY+   A+++  ILD
Sbjct: 477 IQSRRMLSNASYFAFTATPKNRTLELFGEKTVVDDKVEFRSPEELTYTTKQAIQEGFILD 536

Query: 527 VIKQFSTL 534
           VI  ++TL
Sbjct: 537 VIAHYTTL 544


 Score = 83.2 bits (204), Expect = 4e-14
 Identities = 62/181 (34%), Positives = 92/181 (50%), Gaps = 15/181 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           K+RAM+VC+G    + Y + +   + Q  I+  + +
Sbjct: 580 IRRKAEIMVDHFVAQVMGAKKIGGKARAMIVCNGIARAIDYWREVSDYLEQ--IKSPYKA 637

Query: 679 IAAFS-PFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           I A+S  F   G+   E+D  +N   +    D+ R          Q   R LIVA+K  T
Sbjct: 638 IVAYSGEFEIGGEKKTEAD--LNDFPSKEIPDRLR----------QDPYRFLIVANKFVT 685

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           GFDEP L  MYVDKKL G  AVQTLSRL+R    K  T ++DF +  + +   F  Y  +
Sbjct: 686 GFDEPLLHTMYVDKKLAGVQAVQTLSRLNRAHPQKHDTFVLDFADNAEAVKKAFQDYYRA 745

Query: 798 T 798
           T
Sbjct: 746 T 746


>gi|48868351|ref|ZP_00321696.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Haemophilus influenzae 86-028NP]
          Length = 646

 Score = 93.2 bits (230), Expect = 4e-17
 Identities = 90/312 (28%), Positives = 132/312 (41%), Gaps = 88/312 (28%)

Query: 261 RVQNSFKD-------FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGK 313
           R+  S KD       F PR HQ D +  ++ D+      K         YLIQH+AGSGK
Sbjct: 248 RLDGSTKDPLDKRTLFFPRYHQLDVVRRLIADVSEHGVGK--------RYLIQHSAGSGK 299

Query: 314 SLTIACLVYFLYKLEINNKHK----------FDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           S +I  L Y L +    N+            FD+++++ DR  LD QL D +  F   + 
Sbjct: 300 SNSITWLAYQLIEAYPRNEKAANGREADRPIFDSVIVVTDRRLLDKQLRDNIKDF---SE 356

Query: 364 ITNFCRPK-SSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXX 422
           + N   P  SS +L+  +  GK +++ITT+QKF  ++   +                   
Sbjct: 357 VKNIVAPALSSAELRQSLEQGK-KIIITTIQKFPFIVDGIA------------------- 396

Query: 423 XYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSG-------- 474
                        +  D QF       A+I DEAH S   S  + ++  +          
Sbjct: 397 -------------DLGDKQF-------AVIIDEAHSSQSGSAHDNMNRAIGKTEDLDAED 436

Query: 475 ------NTRQSSKI----TYFSFTCTPTTQCLEMFGIASKGDKI-PFYTYSMHNALEDKL 523
                  T QS K+    +YF+FT TP    LE FG      K  PF+ YSM  A+E+  
Sbjct: 437 VQDLILQTMQSRKMHGNASYFAFTATPKNSTLEKFGEKQADGKFKPFHLYSMKQAIEEGF 496

Query: 524 ILDVIKQFSTLR 535
           ILDVI  ++T +
Sbjct: 497 ILDVIANYTTYK 508


>gi|9656287|gb|AAF94914.1| type I restriction enzyme HsdR, putative
           [Vibrio cholerae O1 biovar eltor str. N16961]
 ref|NP_231400.1| type I restriction enzyme HsdR, putative [Vibrio cholerae O1 biovar
           eltor str. N16961]
 pir||E82160 probable type I restriction enzyme HsdR VC1765 [imported] - Vibrio
           cholerae (strain N16961 serogroup O1)
          Length = 1022

 Score = 92.4 bits (228), Expect = 7e-17
 Identities = 83/293 (28%), Positives = 124/293 (41%), Gaps = 77/293 (26%)

Query: 269 FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL- 327
           F PR HQ D +V  L D  S   +          YLIQH+AGSGKS +I    Y L +  
Sbjct: 290 FFPRYHQMD-VVRRLVDHCSVNGVGQ-------TYLIQHSAGSGKSNSITWAAYQLIETY 341

Query: 328 ---------EINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK-SSKQLQ 377
                           FD+++++ DR  LD QL D + +F   + + N   P   S +L+
Sbjct: 342 PISDDLPGSRGKEMPLFDSVIVVTDRRLLDKQLRDNIKEF---SEVKNIVAPAFKSSELK 398

Query: 378 NEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
           + + +GK +++ITT+QKF  ++                                      
Sbjct: 399 SALENGK-KIIITTIQKFPYIV-------------------------------------- 419

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFL--------------SGNTRQSSKIT 483
            D          A+I DEAH S      +KL+E +              + N +  S  +
Sbjct: 420 -DGIADLSDRRFAVIIDEAHSSQDGHNQDKLNEAMGFVSEDVLDKALQSAKNRKMRSNAS 478

Query: 484 YFSFTCTPTTQCLEMFGIA-SKGDKIPFYTYSMHNALEDKLILDVIKQFSTLR 535
           YF+FT TP    LE FG   + G  +PF+ YSM  A+E+  ILDVI  ++T +
Sbjct: 479 YFAFTATPKNTTLEKFGQRQADGTYVPFHLYSMKQAIEEGFILDVIANYTTYK 531


 Score = 71.2 bits (173), Expect = 2e-10
 Identities = 63/208 (30%), Positives = 95/208 (45%), Gaps = 20/208 (9%)

Query: 616 QMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQ 675
           Q  I TKA+ ++ H +           K + M+V    +  + Y + L + ++++     
Sbjct: 563 QETIDTKAEIMLEHFIKHVVNGKKLKGKGKGMVVTQNIESAIRYYRALTRQLNKMG--NP 620

Query: 676 FDSIAAFSPFIQ-EGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNI-------- 726
           F    AFS   + +G    E+D  ING   P G  K+   +     +    I        
Sbjct: 621 FKVAIAFSGSKEVDGIEYTEAD--ING--FPEGDTKDYFDVNYKRKEPDSPIPKHVDQDA 676

Query: 727 -RLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFV 781
            RLL+VA+K  TGFD+P L  MYVDKKL     VQ LSRL+R      K  E   ++DF 
Sbjct: 677 YRLLVVANKYLTGFDQPKLCAMYVDKKLASVLCVQALSRLNRSAPKYGKKTEDLFVLDFF 736

Query: 782 NTEKDMDFIFNLYKHSTSLRDLQDTQTL 809
           N+  D+   F+ +  ST+L +  D   L
Sbjct: 737 NSVDDIKTAFDPFYTSTTLSEATDVNVL 764


>gi|41690596|ref|ZP_00147128.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Psychrobacter sp. 273-4]
          Length = 1102

 Score = 92.0 bits (227), Expect = 9e-17
 Identities = 79/307 (25%), Positives = 129/307 (41%), Gaps = 85/307 (27%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQWD +  ++   ++    +         YLIQH+AGSGKS +IA   + L  L  +
Sbjct: 316 PRYHQWDVVTKLVNAAVTEGAGQK--------YLIQHSAGSGKSNSIAWTAHQLSTLYND 367

Query: 331 NKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK----SSKQLQNEILSGKN 385
              K F +++++ DR+ LD QL + + +F  ++G+      K    S  +   E L    
Sbjct: 368 QGDKQFHSVIVITDRTILDDQLQETIYQFEHADGVVGKINRKEGDGSKSEKLAEALVNSQ 427

Query: 386 RVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFK 445
            ++I T+Q F  +L                                       ++    K
Sbjct: 428 PIIIVTIQTFPFVL------------------------------------KAIEDSSVLK 451

Query: 446 SNSVAIISDEAHRSHGKSTTEKLHEFL-----SGNT-----------------------R 477
           S   AII+DEAH S   ST  +L E L       NT                       R
Sbjct: 452 SRRYAIIADEAHSSQTGSTARQLKEVLVSGDIESNTDDDKPLSSEDSLDFILDATLAARR 511

Query: 478 QSSKITYFSFTCTPTTQCLEMFG------IASKGDKIP--FYTYSMHNALEDKLILDVIK 529
            S  ++Y++FT TP  + +E+FG      + +  D +P  ++ YSM  A+E+  ILDV+K
Sbjct: 512 SSPNLSYYAFTATPKPKTIELFGRLPNPDLPASKDNLPAAYHVYSMRQAIEEGFILDVLK 571

Query: 530 QFSTLRL 536
            ++  ++
Sbjct: 572 NYTNYKV 578


 Score = 92.0 bits (227), Expect = 9e-17
 Identities = 72/236 (30%), Positives = 114/236 (47%), Gaps = 22/236 (9%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I  K + II H  +  N       +++AM+V   RK  + YK   D  I+        ++
Sbjct: 613 IAQKVKIIIEH--FNDNIKGLLGGQAKAMVVTGSRKEAVRYKIAFDNYIADYGYRN-INA 669

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTD---KNRGIIKALNDDKQFNIRLLIVADKL 735
           + AFS  +    +  +S  +I  K+  +  +   K R + KA + D   + ++++VA+K 
Sbjct: 670 MVAFSGEVTFNDNDPDSGALIGEKFTEHNMNIGLKGRDLRKAFDSD---DYQVMLVANKF 726

Query: 736 QTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDK--EATAIIDFVNTEKDMDFIFNL 793
           QTGFD+P L  MYVDKKL G + VQTLSRL+R  K K    T ++DF N  +D+   F  
Sbjct: 727 QTGFDQPKLCAMYVDKKLTGVDCVQTLSRLNRTYKGKAESGTFVLDFFNDPEDILEAFQP 786

Query: 794 YKHSTSLRDLQDTQ---TLISSLIMAG--------RSIQKVFEKEKRAGLITPFCK 838
           Y  +  L D+ D      L   L ++G        + ++  + K K    I+  CK
Sbjct: 787 YYETAHLADVSDPDQVFDLYEKLRVSGIFTKNEVEQFVEAFYSKNKSNAAISNICK 842


>gi|52007207|ref|ZP_00334585.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Thiobacillus denitrificans ATCC
           25259]
          Length = 1009

 Score = 92.0 bits (227), Expect = 9e-17
 Identities = 70/262 (26%), Positives = 112/262 (42%), Gaps = 58/262 (22%)

Query: 301 YNYLIQHAAGSGKSLTIACLVYFLYKLEINNKHK-FDTIVIMNDRSQLDAQLGDVVIKFL 359
           +NYL++H+AGSGKS TI  L + L  L   + HK FD+++++ DR  LD QL D + +F 
Sbjct: 300 HNYLVEHSAGSGKSNTIGWLAHRLASLHNESDHKVFDSVIVITDRRVLDRQLQDTIYQFE 359

Query: 360 TSNGITNFCRPKSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXX 419
              G+      + S+QL  E L+    ++I+T+QKF                        
Sbjct: 360 HRQGVVQKI-DEDSRQLA-EALTNSVPIIISTLQKFP----------------------- 394

Query: 420 XXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNT--- 476
               +             +       +   A+I DEAH S    T   L E L G+    
Sbjct: 395 ----FVARQLARMAEERREGADGILPTRRCAVIIDEAHSSQSGETATDLKEVLGGHVLRE 450

Query: 477 ------------------------RQSSKITYFSFTCTPTTQCLEMFGIASKGDK-IPFY 511
                                    + + +++F+FT TP  + L +FG    G    PF+
Sbjct: 451 EAARYAVDEGVDNLDELYRSMAKRGRQTNLSFFAFTATPKHKTLAVFGRDEAGKAGQPFH 510

Query: 512 TYSMHNALEDKLILDVIKQFST 533
            Y+M  A+E+  I+DV+K ++T
Sbjct: 511 RYTMSQAIEEGFIMDVLKNYTT 532


 Score = 78.6 bits (192), Expect = 1e-12
 Identities = 59/190 (31%), Positives = 92/190 (48%), Gaps = 16/190 (8%)

Query: 619 IKTKAQFIINHL-LYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFD 677
           I  K + ++ H    T++K   K   ++AM+V   R   + YK   D  I     E+ + 
Sbjct: 569 IAQKTEVMVEHFNAVTRHKVGGK---AKAMVVTGSRLEAVRYKLAFDAYIR----EKGYS 621

Query: 678 SIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKN--RGIIKALNDDKQFNIRLLIVADKL 735
            I     F  E    V  D +    +   G +     G++    D  ++ +  L+VA+K 
Sbjct: 622 WIKTLVAFSGE----VVDDELPGKTWTEVGMNGGIKEGVLPEKFDSSEYQV--LLVAEKY 675

Query: 736 QTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYK 795
           QTG+D+P L  MYVDK+L G  AVQTLSRL+R+   KE T ++DFVN  +++   F  Y 
Sbjct: 676 QTGYDQPLLHTMYVDKRLAGIQAVQTLSRLNRMHPLKEDTFVLDFVNDRQEIKEAFKQYH 735

Query: 796 HSTSLRDLQD 805
               + +  D
Sbjct: 736 EGAEMGEEAD 745


>gi|21109200|gb|AAM37743.1| type I restriction-modification system
           endonuclease [Xanthomonas axonopodis pv. citri str. 306]
 ref|NP_643207.1| type I restriction-modification system endonuclease [Xanthomonas
           axonopodis pv. citri str. 306]
          Length = 1084

 Score = 91.7 bits (226), Expect = 1e-16
 Identities = 60/199 (30%), Positives = 97/199 (48%), Gaps = 16/199 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           +  K +FI+ H  +  N        ++AM+V S R   + YKK  D  I   P  + + +
Sbjct: 578 VTQKVEFIVEH--FAANVAPLLGGTAKAMIVTSSRAAAVRYKKAFDAYIRANPQHKPYQA 635

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFN-----------IR 727
           + AFS  +   + M  SD  +  K  P+    +    ++  + +Q+             R
Sbjct: 636 LVAFSGSLTGNEVMHASDEQL--KNDPFAAQDDEQFTESTMNPRQYGTDLRKVFDRPEFR 693

Query: 728 LLIVADKLQTGFDEPSLSIMYVDKKLKG-ANAVQTLSRLSRISKDKEATAIIDFVNTEKD 786
           ++IVA+K QTGFD+P L  MY+DKK+      VQTLSRL+R +  K+ T ++DFVN  + 
Sbjct: 694 VMIVANKFQTGFDQPKLCAMYLDKKIANPVEIVQTLSRLNRTAPGKDTTFVVDFVNDPQS 753

Query: 787 MDFIFNLYKHSTSLRDLQD 805
           +   F +Y     +   QD
Sbjct: 754 ILDAFAMYDAGAKIEKAQD 772


 Score = 70.9 bits (172), Expect = 2e-10
 Identities = 72/280 (25%), Positives = 116/280 (40%), Gaps = 48/280 (17%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQW  +  ++ D  +    +S        YL +H+AGSGK+ TIA   + L +L  +
Sbjct: 289 PRFHQWKAVTRMIDDAKAKGPGQS--------YLCEHSAGSGKTSTIAWTAHDLIRLRTD 340

Query: 331 -NKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK-----SSKQLQNEILSGK 384
              + F +++++ DR+ LD QL   V +      +      K      SK+L  E L   
Sbjct: 341 AGVNLFHSVIVVTDRTVLDKQLQYAVQQIDHQKDLIAAIGEKDSGVSKSKELA-EALERG 399

Query: 385 NRVVITTMQKF----SQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDN 440
             +++ T+Q F     Q+L N                        N        +     
Sbjct: 400 TPIIVVTIQTFPYAMEQILTNGKLAD------------------RNFAVIIDEAHTSQTG 441

Query: 441 QFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKITYFSFTCTPTTQCLEMFG 500
             + K  +   + D A    G +  E L +  S   R ++ ++YF+FT TP    L +FG
Sbjct: 442 TAASKLQATLALRD-AKAMEGLTVEELLEKIQSSRVRPNN-VSYFAFTATPKHATLTLFG 499

Query: 501 I-------ASKGDKIP--FYTYSMHNALEDKLILDVIKQF 531
                   ASK D +P  F  Y M  A+E+  ILDV++ +
Sbjct: 500 TRPDPGKPASKEDNLPVSFDKYPMRQAIEEGFILDVLQGY 539


>gi|23473873|ref|ZP_00129168.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Desulfovibrio desulfuricans G20]
          Length = 1085

 Score = 91.7 bits (226), Expect = 1e-16
 Identities = 68/201 (33%), Positives = 103/201 (50%), Gaps = 20/201 (9%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           +  K QFI+ H  +TKN       K++AM+V S R   + YKK  D+ I Q        S
Sbjct: 584 VTQKVQFIVEH--FTKNVAHRLDGKAKAMVVTSSRAAAIRYKKAFDRYIEQHSEYGFIHS 641

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAP--YGTDKNRGIIK-ALNDD----------KQFN 725
           + AFS  +  GK ++  D   +G++    +  D+N    + ++N D           +  
Sbjct: 642 LVAFSGKMT-GKQVMHQD---DGEFKDDVFIVDENEEFTEQSMNPDIHGQDLRFAFDRPE 697

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKG-ANAVQTLSRLSRISKDKEATAIIDFVNTE 784
            R+++VADK QTGFD+P L  MYVDKK+      VQT SRL+R +  K+   IIDFVN  
Sbjct: 698 YRVMLVADKFQTGFDQPKLVAMYVDKKIANHVEIVQTFSRLNRTAPGKDEVFIIDFVNDP 757

Query: 785 KDMDFIFNLYKHSTSLRDLQD 805
           +++   F  Y     + ++QD
Sbjct: 758 ENVRQAFTTYDKGAHIDEVQD 778


 Score = 87.4 bits (215), Expect = 2e-15
 Identities = 88/295 (29%), Positives = 124/295 (41%), Gaps = 79/295 (26%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL-EI 329
           PR HQW  +  +L D        + ++     YL  H+AGSGK+ TI+   + L KL E 
Sbjct: 296 PRFHQWTAVNQILAD--------ARQNGAGMTYLCDHSAGSGKTSTISWTAHDLVKLRED 347

Query: 330 NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNG-ITNFCRPKS----SKQLQNEILSGK 384
           N    F++++I+ DR+ LD QL D V +     G I    R KS    SKQL   +LSG 
Sbjct: 348 NGDAVFNSVIIVTDRNVLDGQLQDAVKQIDHQFGVIAAIDRQKSSKSKSKQLSEALLSG- 406

Query: 385 NRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSF 444
             +V+ T+Q F                            YA                 + 
Sbjct: 407 TPIVVVTIQTFP---------------------------YAMEAIITDK---------AL 430

Query: 445 KSNSVAIISDEAHRSHGKSTTEKLHEFL--SGNTRQSS------------------KITY 484
           K  + A+I DEAH S   +T  KL   L  SG  + S+                   I+Y
Sbjct: 431 KGKNFAVIIDEAHNSQTGATAAKLQAALGMSGQGKMSTMTVDELLEQLQKSRARPKNISY 490

Query: 485 FSFTCTPTTQCLEMFG------IASKGDKIP--FYTYSMHNALEDKLILDVIKQF 531
           F+FT TP    L +FG        +  D  P  F+ Y+M  A+E+K ILDV+K +
Sbjct: 491 FAFTGTPKHSTLMLFGRPTDPSQPASDDNPPQAFHLYTMRQAIEEKFILDVLKGY 545


>gi|53756220|gb|AAU90511.1| type I restriction-modification system,
           R subunit, putative [Methylococcus capsulatus str. Bath]
 ref|YP_112808.1| type I restriction-modification system, R subunit, putative
           [Methylococcus capsulatus str. Bath]
          Length = 1002

 Score = 91.3 bits (225), Expect = 2e-16
 Identities = 79/304 (25%), Positives = 132/304 (42%), Gaps = 82/304 (26%)

Query: 261 RVQNSFKDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACL 320
           R + + K   PR HQ  T+ ++L+        +S        YLIQH+AGSGKS TIA L
Sbjct: 266 RKRKTRKQIFPRYHQLRTVRALLR--------RSQADGVGKRYLIQHSAGSGKSNTIAWL 317

Query: 321 VYFLYKLEI---NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQ 377
            + L +L+      + +FD+++++ DR  LD Q+   +  +   +  + +   +S+++L+
Sbjct: 318 AHQLVELKTAADAGQAQFDSVIVITDRRALDTQIARTIRSY--DHVASIYGHSESAEELR 375

Query: 378 NEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
             +  GK ++++TT+QKF  +L                                      
Sbjct: 376 TFLRRGK-KIIVTTVQKFPFIL-------------------------------------- 396

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------------------ 479
            D          A++ DEAH S G  TT K+H  LSG   +                   
Sbjct: 397 -DELGDLGDKKFALLIDEAHSSQGGKTTAKMHLALSGQAAEGGEDEEEESVEDKVNALIE 455

Query: 480 -----SKITYFSFTCTPTTQCLEMFGIAS-KGDKIPF-----YTYSMHNALEDKLILDVI 528
                +  +Y++FT TP T+ LE+FG     GD + F      TY+   A+++  ILDVI
Sbjct: 456 SRKMLANASYYAFTATPKTKTLELFGERQVVGDTVQFRSPEELTYTTKQAIQEGFILDVI 515

Query: 529 KQFS 532
             ++
Sbjct: 516 ANYT 519


 Score = 79.7 bits (195), Expect = 5e-13
 Identities = 58/180 (32%), Positives = 92/180 (50%), Gaps = 13/180 (7%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           K+RAM+VC+G    + Y + +   + +  I+  + +
Sbjct: 557 IRRKAEIMVDHFIEQVIGAKKIGGKARAMIVCNGIARAIDYFREVSDYLRE--IKSPYKA 614

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTG 738
           I A+S   + G  + +++  +NG         ++ I   L  D     R LIVA+K  TG
Sbjct: 615 IVAYSGDFEVG-GVKKTEADLNGF-------PSKDIPAKLRQDPY---RFLIVANKFVTG 663

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHST 798
           FDEP L  MYVDK L G  AVQTLSRL+R    K  T ++DF +  + +   F  Y  +T
Sbjct: 664 FDEPLLHTMYVDKPLAGVLAVQTLSRLNRAHPQKADTFVLDFADNAEAVKAAFQEYYRAT 723


>gi|34762431|ref|ZP_00143431.1| TYPE I RESTRICTION-MODIFICATION
           SYSTEM RESTRICTION SUBUNIT [Fusobacterium nucleatum
           subsp. vincentii ATCC 49256]
 gb|EAA24967.1| TYPE I RESTRICTION-MODIFICATION SYSTEM RESTRICTION SUBUNIT
           [Fusobacterium nucleatum subsp. vincentii ATCC 49256]
          Length = 980

 Score = 91.3 bits (225), Expect = 2e-16
 Identities = 75/275 (27%), Positives = 121/275 (43%), Gaps = 53/275 (19%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL--- 327
           PR HQ+  +  +++D++        K+    NYLIQH+AGSGKS +I  L Y L ++   
Sbjct: 281 PRYHQYRVVTKLVEDVM--------KNGVGKNYLIQHSAGSGKSNSITWLAYRLVEVARK 332

Query: 328 ----EINNKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILS 382
               +I+   K F++++++ DR  LD Q+ D + KF+    +       SS  L+N +L 
Sbjct: 333 VQNGDISKFEKIFNSVIVVTDRLNLDKQIDDNIRKFIDVRSVVG--HASSSTDLKN-LLV 389

Query: 383 GKNRVVITTMQKFSQLLPN-----KSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
            + +++I T+QKF  LL       K                       N         N 
Sbjct: 390 NEKKIIIATIQKFPYLLEKIGTELKGKNFAIIIDEAHSSQSGRAAASLNMAVSATIDIN- 448

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKITYFSFTCTPTTQCLE 497
           DDN F  +                    +KL+E +    +     ++F+FT TP  + L+
Sbjct: 449 DDNDFEIE--------------------DKLNELIEAR-KMPENASFFAFTATPKAKTLQ 487

Query: 498 MFGIASKGDKIPFYTYSMHNALEDKLILDVIKQFS 532
           MFG         F  YSM  A+E+  ILDV+K ++
Sbjct: 488 MFGNV-------FDLYSMKQAIEEGFILDVLKNYT 515


 Score = 85.1 bits (209), Expect = 1e-14
 Identities = 69/221 (31%), Positives = 110/221 (49%), Gaps = 19/221 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ ++ H L+  N       +++AM+V       + Y   ++ ++ +      +++
Sbjct: 553 IREKAEVMVEHFLH--NTATKIAGQAKAMVVTQSILSAIEYYHCINDMLKKS--NTGYEA 608

Query: 679 IAAFSPFIQ-EGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           + AFS   +  GK++ E+   +NG      T K           K+   + LIVADK QT
Sbjct: 609 LIAFSGEKEYNGKTVTEAS--LNGFPDTETTTKF----------KETKYKFLIVADKYQT 656

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           G+DEP L  MYVDK L    AVQTLSRL+R +K+K  T IIDF N  + +   F  Y   
Sbjct: 657 GYDEPLLHTMYVDKVLNDVKAVQTLSRLNRCTKNKIDTCIIDFANKPEHIQEAFQPYYKE 716

Query: 798 TSLRDLQDTQTLISSLIMAGRSIQKVFEKEKRAGLITPFCK 838
           T L    D   L +  +++    + V+EK++   L+  F K
Sbjct: 717 TKLDGEVDPNKLFN--LLSVLDSKYVYEKDEVDELVELFLK 755


>gi|52857360|ref|ZP_00341546.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Xylella fastidiosa Ann-1]
          Length = 1027

 Score = 90.5 bits (223), Expect = 3e-16
 Identities = 78/304 (25%), Positives = 123/304 (39%), Gaps = 88/304 (28%)

Query: 267 KDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYK 326
           K   PR HQ  T+    + L+   H           YLIQH+AGSGKS TIA L + L +
Sbjct: 290 KQIFPRFHQLRTV----RALLGRAHSDGVGGR----YLIQHSAGSGKSNTIAWLAHQLVE 341

Query: 327 LEINN---KHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSG 383
           L   +     +FD+I+++ DR  LD Q+ D +  +   + +        + Q   E L  
Sbjct: 342 LRRKDDPMTPQFDSIIVITDRRALDTQIADTIKGY---DHVAAIFGHSDNAQELREYLRR 398

Query: 384 KNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFS 443
             ++++TT+QKF  +L                                       D    
Sbjct: 399 DKKIIVTTVQKFPFIL---------------------------------------DELGD 419

Query: 444 FKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------------------------ 479
                 A++ DEAH S G   T ++HE L G   +                         
Sbjct: 420 LSGKRFALLIDEAHSSQGGKITARMHETLGGKVAEEEFEEEFEEDRTQDAVNAEIERRIA 479

Query: 480 -----SKITYFSFTCTPTTQCLEMFGIAS-KGDKIPF-----YTYSMHNALEDKLILDVI 528
                +  +YF+FT TP  + LE+FG  +  GDK+ F      TY+   A+++K ILDV+
Sbjct: 480 SRKLLANASYFAFTATPKNKTLELFGEKTLVGDKVQFRSPEELTYTTKQAIQEKFILDVV 539

Query: 529 KQFS 532
           ++++
Sbjct: 540 EKYT 543


 Score = 81.3 bits (199), Expect = 2e-13
 Identities = 61/181 (33%), Positives = 93/181 (50%), Gaps = 15/181 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           ++RAM+VC+G    + Y + +   ++Q  I+  + +
Sbjct: 581 IRRKAEIMVDHFIAEVVGKQKIGGQARAMIVCNGIARAIDYWREVSDYLTQ--IKSPYRA 638

Query: 679 IAAFS-PFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           I A+S  F   G+   E+D  +NG         ++ I   L  D     R LIVA+K  T
Sbjct: 639 IVAYSGSFEIGGQKKTEAD--LNGF-------PSKDIPDTLKKDPY---RFLIVANKFVT 686

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           GFDEP L  MYVDK L G  AVQTLSRL+R    K  T ++DF +  + +   F  Y  +
Sbjct: 687 GFDEPLLHTMYVDKPLAGVLAVQTLSRLNRAHPQKRDTFVLDFADNAEAVKAAFQDYYRA 746

Query: 798 T 798
           T
Sbjct: 747 T 747


>gi|28199930|ref|NP_780244.1| type I restriction-modification system
           endonuclease [Xylella fastidiosa Temecula1]
 gb|AAO29893.1| type I restriction-modification system endonuclease [Xylella
           fastidiosa Temecula1]
          Length = 1031

 Score = 90.1 bits (222), Expect = 3e-16
 Identities = 79/308 (25%), Positives = 124/308 (39%), Gaps = 92/308 (29%)

Query: 267 KDFKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYK 326
           K   PR HQ  T+    + L+   H           YLIQH+AGSGKS TIA L + L +
Sbjct: 290 KQIFPRFHQLCTV----RALLGRAHSDGVGGR----YLIQHSAGSGKSNTIAWLAHQLVE 341

Query: 327 LEINN---KHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSG 383
           L   +     +FD+I+++ DR  LD Q+ D +  +   + +        + Q   E L G
Sbjct: 342 LRRKDDPMTPQFDSIIVITDRRALDTQIADTIKGY---DHVAAIFGHSDNAQELCEYLRG 398

Query: 384 KNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFS 443
             ++++TT+QKF  +L                                       D    
Sbjct: 399 GKKIIVTTVQKFPFIL---------------------------------------DQLGD 419

Query: 444 FKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------------------------ 479
                 A++ DEAH S G   T ++HE L G   +                         
Sbjct: 420 LSGKRFALLIDEAHSSQGGKITARMHETLGGKVAEEEFEEEFEEEFEEDRTQDAVNAEIE 479

Query: 480 ---------SKITYFSFTCTPTTQCLEMFGIAS-KGDKIPF-----YTYSMHNALEDKLI 524
                    +  +YF+FT TP  + LE+FG  +  GDK+ F      TY+   A+++K I
Sbjct: 480 RRIASRKLLANASYFAFTATPKNKTLELFGEKTLVGDKVQFRSPEELTYTTKQAIQEKFI 539

Query: 525 LDVIKQFS 532
           LDV+++++
Sbjct: 540 LDVVEKYT 547


 Score = 81.3 bits (199), Expect = 2e-13
 Identities = 61/181 (33%), Positives = 93/181 (50%), Gaps = 15/181 (8%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+ KA+ +++H +           ++RAM+VC+G    + Y + +   ++Q  I+  + +
Sbjct: 585 IRRKAEIMVDHFIAEVVGKQKIGGQARAMIVCNGIARAIDYWREVSDYLTQ--IKSPYRA 642

Query: 679 IAAFS-PFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQT 737
           I A+S  F   G+   E+D  +NG         ++ I   L  D     R LIVA+K  T
Sbjct: 643 IVAYSGSFEIGGQKKTEAD--LNGF-------PSKDIPDTLKKDPY---RFLIVANKFVT 690

Query: 738 GFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           GFDEP L  MYVDK L G  AVQTLSRL+R    K  T ++DF +  + +   F  Y  +
Sbjct: 691 GFDEPLLHTMYVDKPLAGVLAVQTLSRLNRAHPQKRDTFVLDFADNAEAVKAAFQDYYRA 750

Query: 798 T 798
           T
Sbjct: 751 T 751


>gi|42632163|ref|ZP_00157701.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Haemophilus influenzae R2866]
          Length = 1006

 Score = 90.1 bits (222), Expect = 3e-16
 Identities = 91/313 (29%), Positives = 135/313 (43%), Gaps = 90/313 (28%)

Query: 261 RVQNSFKD-------FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGK 313
           R+  S KD       F PR HQ + +  ++ D +S Q +          YLIQH+AGSGK
Sbjct: 273 RLDGSTKDPLEKRSLFFPRYHQLEVVRRLIAD-VSEQGVGK-------RYLIQHSAGSGK 324

Query: 314 SLTIACLVYFLYKLEINNKHK----------FDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           S +I  L Y L +    N+            FD+++++ DR  LD QL D +  F   + 
Sbjct: 325 SNSITWLAYQLIEAYPCNEKAANGREADRPIFDSVIVVTDRRLLDKQLRDNIKDF---SE 381

Query: 364 ITNFCRPK-SSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXX 422
           + N   P  SS +L+  +  GK +++ITT+QKF  ++   +                   
Sbjct: 382 VKNIVAPALSSAELRQSLEQGK-KIIITTIQKFPFIVDGIA------------------- 421

Query: 423 XYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTR----- 477
                        +  D QF       A+I DEAH S   S  + ++  + G T      
Sbjct: 422 -------------DLGDKQF-------AVIIDEAHSSQSGSAHDNMNRAI-GKTEDLDAE 460

Query: 478 ----------QSSKI----TYFSFTCTPTTQCLEMFGIASKGDKI-PFYTYSMHNALEDK 522
                     QS K+    +YF+FT TP    LE FG      K  PF+ YSM  A+E+ 
Sbjct: 461 DVQDLILQAMQSRKMRGNASYFAFTATPKNSTLEKFGEKQADGKFKPFHLYSMKQAIEEG 520

Query: 523 LILDVIKQFSTLR 535
            ILDVI  ++T +
Sbjct: 521 FILDVIANYTTYK 533


 Score = 75.9 bits (185), Expect = 7e-12
 Identities = 64/201 (31%), Positives = 98/201 (47%), Gaps = 19/201 (9%)

Query: 614 QDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIE 673
           + Q  I TKA+ +++H +           K++ M+V    +  + Y + L +L++     
Sbjct: 563 RSQQTIDTKAEIMLDHFIRQIFNRKKLKGKAKGMVVTQNIETAIRYFQALKRLLAGRG-- 620

Query: 674 EQFDSIAAFSPF-IQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVA 732
             F    AFS   + +G    E++  +NG +A   T       K   D  ++  RLL+VA
Sbjct: 621 NPFKIAIAFSGSKVVDGVEYTEAE--MNG-FAESET-------KEYFDQDEY--RLLVVA 668

Query: 733 DKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFVNTEKDMD 788
           +K  TGFD+P L  MYVDKKL G   VQ LSRL+R    + K  E   ++DF N+ +D+ 
Sbjct: 669 NKYLTGFDQPKLCAMYVDKKLSGVLCVQALSRLNRSANKLGKRTEDLFVLDFFNSVEDIQ 728

Query: 789 FIFNLYKHSTSLRDLQDTQTL 809
             F  +  STSL    D   L
Sbjct: 729 QAFEPFYTSTSLSQATDVNVL 749


>gi|42629287|ref|ZP_00154835.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Haemophilus influenzae R2846]
          Length = 1007

 Score = 90.1 bits (222), Expect = 3e-16
 Identities = 91/313 (29%), Positives = 135/313 (43%), Gaps = 90/313 (28%)

Query: 261 RVQNSFKD-------FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGK 313
           R+  S KD       F PR HQ + +  ++ D +S Q +          YLIQH+AGSGK
Sbjct: 274 RLDGSTKDPLEKRSLFFPRYHQLEVVRRLIAD-VSEQGVGK-------RYLIQHSAGSGK 325

Query: 314 SLTIACLVYFLYKLEINNKHK----------FDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           S +I  L Y L +    N+            FD+++++ DR  LD QL D +  F   + 
Sbjct: 326 SNSITWLAYQLIEAYPCNEKAANGREADRPIFDSVIVVTDRRLLDKQLRDNIKDF---SE 382

Query: 364 ITNFCRPK-SSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXX 422
           + N   P  SS +L+  +  GK +++ITT+QKF  ++   +                   
Sbjct: 383 VKNIVAPALSSAELRQSLEQGK-KIIITTIQKFPFIVDGIA------------------- 422

Query: 423 XYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTR----- 477
                        +  D QF       A+I DEAH S   S  + ++  + G T      
Sbjct: 423 -------------DLGDKQF-------AVIIDEAHSSQSGSAHDNMNRAI-GKTEDLDAE 461

Query: 478 ----------QSSKI----TYFSFTCTPTTQCLEMFGIASKGDKI-PFYTYSMHNALEDK 522
                     QS K+    +YF+FT TP    LE FG      K  PF+ YSM  A+E+ 
Sbjct: 462 DVQDLILQAMQSRKMRGNASYFAFTATPKNSTLEKFGEKQADGKFKPFHLYSMKQAIEEG 521

Query: 523 LILDVIKQFSTLR 535
            ILDVI  ++T +
Sbjct: 522 FILDVIANYTTYK 534


 Score = 77.0 bits (188), Expect = 3e-12
 Identities = 65/201 (32%), Positives = 98/201 (48%), Gaps = 19/201 (9%)

Query: 614 QDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIE 673
           + Q  I TKA+ +++H +           K++ M+V    +  + Y + L  L++     
Sbjct: 564 RSQQTIDTKAEIMLDHFIRQIFSRKKLKGKAKGMVVTQNIETAIRYFQALKHLLAGRG-- 621

Query: 674 EQFDSIAAFSPF-IQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVA 732
             F    AFS   + +G    E++  +NG +A   T       K   D  ++  RLL+VA
Sbjct: 622 NPFKIAIAFSGSKVVDGVEYTEAE--MNG-FAESET-------KEYFDQDEY--RLLVVA 669

Query: 733 DKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFVNTEKDMD 788
           +K  TGFD+P L  MYVDKKL G   VQ LSRL+R    +SK  E   ++DF N+ +D+ 
Sbjct: 670 NKYLTGFDQPKLCAMYVDKKLSGVLCVQALSRLNRSANKLSKRTEDLFVLDFFNSVEDIQ 729

Query: 789 FIFNLYKHSTSLRDLQDTQTL 809
             F  +  STSL    D   L
Sbjct: 730 QAFEPFYTSTSLSQATDVNVL 750


>gi|2411488|gb|AAB70710.1| HsdR [Klebsiella pneumoniae]
 pir||T30329 hsdR protein - Klebsiella pneumoniae
          Length = 1102

 Score = 90.1 bits (222), Expect = 3e-16
 Identities = 121/532 (22%), Positives = 201/532 (37%), Gaps = 133/532 (25%)

Query: 308 AAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNF 367
           A G GKS T+  L   L K     + +  TIV++ DR+ LD QL   +       G    
Sbjct: 327 ATGCGKSYTMQFLTRLLMK---KRRVRSPTIVLITDRTDLDDQLSAQMCNAKNYIGDDTV 383

Query: 368 CRPKSSKQLQNEILSGKNR--VVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYA 425
             P +S+    E L+G+N   V +TT+ KF+                             
Sbjct: 384 V-PVTSRDDLREKLAGRNSGGVFLTTIHKFT----------------------------- 413

Query: 426 NSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------ 479
                      ED    S +SN + I SDEAHRS      + + +  SG  R++      
Sbjct: 414 -----------EDTELLSGRSNIICI-SDEAHRSQVNLDQKVIVDKESGRVRKTYGFAKY 461

Query: 480 -----SKITYFSFTCTPTTQCLEMFGIASKGDKIPFYTYSMHNALEDKLILDVIKQ---- 530
                 + TY  FT TP    +E+FG           +Y+M  +++D++ + ++ +    
Sbjct: 462 LYDSLPQATYVGFTGTPIDATMEVFGNIVD-------SYTMTESVQDEITVRIVYEGRAA 514

Query: 531 ---FSTLRLNEVNVDSFXXXXXXXXXXXXXXXXXXXXXXXXYNEEAXXXXXXXXXXXXXX 587
                T +L EV                             Y EE               
Sbjct: 515 KVILDTRKLEEVE---------------------------KYYEECANAGTNEWQIDESK 547

Query: 588 XXXXXXXXCFLNGIDFDDDQYEDETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAM 647
                     +N I  DDD+ E             A+    H  Y K  ++    K +AM
Sbjct: 548 KASAT-----MNAILGDDDRLE-----------ALAEDFAKH--YEKRVNEGSTVKGKAM 589

Query: 648 LVCSGRK-GLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDP----MINGK 702
            VC+ R+     Y++  D     + +++  D +       QE K +  S+     M  GK
Sbjct: 590 FVCASREIDWDFYRQLKDFRPEWVEVKQAPDDVELTE---QEEKELPPSEMVKMVMTRGK 646

Query: 703 ------YAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGA 756
                 Y   G+ + R  +     + + N ++ IV D   TGFD P L  +Y+DK L+  
Sbjct: 647 DDHAKLYELLGSKEYRKELDKQFKNAKSNFKIAIVVDMWLTGFDVPELDTIYIDKPLQKH 706

Query: 757 NAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKH--STSLRDLQDT 806
           N +QT+SR++R  + K    ++D++  ++ M+    +Y    +T+  D+Q +
Sbjct: 707 NLIQTISRVNRKLEGKSKGLVVDYIGIKRQMNQALAMYSRIDATNFEDIQQS 758


>gi|56750492|ref|YP_171193.1| hypothetical protein syc0483_c
           [Synechococcus elongatus PCC 6301]
 dbj|BAD78673.1| hypothetical protein [Synechococcus elongatus PCC 6301]
 ref|ZP_00164194.1| COG0610: Type I site-specific restriction-modification system, R
           (restriction) subunit and related helicases
           [Synechococcus elongatus PCC 7942]
          Length = 1039

 Score = 89.7 bits (221), Expect = 5e-16
 Identities = 71/272 (26%), Positives = 117/272 (42%), Gaps = 74/272 (27%)

Query: 301 YNYLIQHAAGSGKSLTIACLVYFLYKL-EINNKHKFDTIVIMNDRSQLDAQLGDVVIKFL 359
           + YL+QH+AGSGK+ +IA   +F  +L + +N+  FD++++++DR+ +D QL + +  F 
Sbjct: 305 HKYLVQHSAGSGKTNSIAWTAHFFSELHDADNRKVFDSVIVVSDRNVIDTQLQEAIESFE 364

Query: 360 TSNG----ITNFCRPKSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXX 415
              G    IT     KSSK    E L G  ++V+ T+Q F   L                
Sbjct: 365 RHKGVVAAITRDDGSKSSKLA--EALKGDKKIVVCTIQTFPFALQAVRELAAT------- 415

Query: 416 XXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLS-- 473
                                        +    A+I+DEAH S       KL + LS  
Sbjct: 416 -----------------------------EGKRFAVIADEAHSSQTGEAASKLKQLLSPA 446

Query: 474 ---------------------GNTRQSSKITYFSFTCTPTTQCLEMFGIASKGDK----- 507
                                 +  + S ITY +FT TP  + LE+FG     ++     
Sbjct: 447 ELADLNDGGEVDLEAVLAAQMSDRARESGITYVAFTATPKAKTLELFGRRPDPNRPAGPD 506

Query: 508 ---IPFYTYSMHNALEDKLILDVIKQFSTLRL 536
              +PF+ YSM  A+E+  ILDV++ +++ ++
Sbjct: 507 NLPMPFHVYSMRQAIEEGFILDVLQNYTSYKV 538


 Score = 79.7 bits (195), Expect = 5e-13
 Identities = 55/172 (31%), Positives = 89/172 (50%), Gaps = 16/172 (9%)

Query: 643 KSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKS----MVESDPM 698
           K++AM+V   RK  + +K  +D+ I+    + +  ++ AFS  +++  S      E+ P 
Sbjct: 596 KAKAMVVVGSRKEAVRWKLAIDRYIAAK--QYRLGTLVAFSGEVRDRDSGPDPFTETSPN 653

Query: 699 INGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANA 758
           +N           +G I+      ++ I  LIVA+K QTGFD+P L  MY+D++L G  A
Sbjct: 654 LNPNL--------QGDIRETFKSDRYQI--LIVANKFQTGFDQPLLCGMYIDRRLAGIQA 703

Query: 759 VQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDTQTLI 810
           VQTLSRL+R    K+ T I+DF N   ++   F  Y  +  L    D   ++
Sbjct: 704 VQTLSRLNRCHPGKDTTYIVDFSNDPAEILAAFKTYYTTAELASATDPNLIL 755


>gi|52426226|ref|YP_089363.1| HsdR protein [Mannheimia
           succiniciproducens MBEL55E]
 gb|AAU38778.1| HsdR protein [Mannheimia succiniciproducens MBEL55E]
          Length = 1020

 Score = 89.7 bits (221), Expect = 5e-16
 Identities = 65/197 (32%), Positives = 103/197 (51%), Gaps = 21/197 (10%)

Query: 618 LIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFD 677
           +I  K + I+ H   T         +++AM+V +GR+  + YK   D+ I +        
Sbjct: 572 VIAQKVEIIVEHFRTTTMHKIGG--RAKAMVVTNGREHAVRYKLAFDEYIKEKGYTG-IK 628

Query: 678 SIAAFSPFI--QEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQF---NIRLLIVA 732
           S+ AFS  I  +E      ++ ++NG             I+ ++  +QF   + ++L+VA
Sbjct: 629 SLVAFSGGITLKEAPEKEYTEALMNG-------------IREVDLPEQFASEHYQVLLVA 675

Query: 733 DKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFN 792
           +K QTGFD+P L  M+VDKKL G  AVQTLSRL+R +K K  T ++DFVN  +D+   F 
Sbjct: 676 EKYQTGFDQPLLHTMFVDKKLSGIQAVQTLSRLNRCAKGKTDTFVLDFVNQPEDIYKAFK 735

Query: 793 LYKHSTSLRDLQDTQTL 809
            +   T L D+   + L
Sbjct: 736 PFYEVTELGDIPSNEKL 752


 Score = 88.6 bits (218), Expect = 1e-15
 Identities = 74/270 (27%), Positives = 124/270 (45%), Gaps = 17/270 (6%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLEIN 330
           PR HQ D +    + LI+  H +   +    NYLIQH+AGSGKS TIA L + L  L   
Sbjct: 277 PRFHQLDAV----RQLIA--HSREHGAGR--NYLIQHSAGSGKSNTIAWLAHQLSSLHNR 328

Query: 331 NKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVI 389
           +  K F++++++ DR  LD QL   + +F   +G+       +S+QL   I S    ++I
Sbjct: 329 DDQKIFNSVIVVTDRVVLDRQLQATISQFEHKDGVVQKIE-HNSQQLAVAIASD-TPIII 386

Query: 390 TTMQKFSQLLP-----NKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSF 444
           TT+QKF  ++       +S                    +++          +  N+   
Sbjct: 387 TTIQKFPFVMAALARKQESGINVAISTEGKQFAVIVDEAHSSQSGEAAMELRKVLNKDGI 446

Query: 445 KSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSS-KITYFSFTCTPTTQCLEMFGIAS 503
           ++  +A   D+     G S   K   F+    RQ    +++F+FT TP  +   +F    
Sbjct: 447 EAAVMAEFLDDDDDETGLSDEAKKQLFIEAAKRQRQPNLSFFAFTATPKWKTKALFDEPG 506

Query: 504 KGDKIPFYTYSMHNALEDKLILDVIKQFST 533
                PF+ Y+M  A+E+  ILDV++ ++T
Sbjct: 507 ADGNTPFHHYTMKQAIEEGFILDVLENYAT 536


>gi|51598164|ref|YP_072355.1| possible type I restriction enzyme
           (restriction subunit) [Yersinia pseudotuberculosis IP
           32953]
 emb|CAH23117.1| possible type I restriction enzyme (restriction subunit) [Yersinia
           pseudotuberculosis IP 32953]
          Length = 1092

 Score = 89.4 bits (220), Expect = 6e-16
 Identities = 123/513 (23%), Positives = 198/513 (37%), Gaps = 121/513 (23%)

Query: 299 NRYNYLIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKF 358
           +R   +I H  GSGKSLT   +V+    L ++   K   I+++ DR  L+ QL    I  
Sbjct: 339 SRNGGVIWHTTGSGKSLT---MVFLSKALILSETLKQCRIIVVTDRVDLENQLSRTFI-- 393

Query: 359 LTSNGITNFCRPK------SSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXX 412
             S G     R K      S ++L  +I  G  R++ + +QKF                 
Sbjct: 394 --STGELAGKRDKADAIATSGRRLAEQIGRGTERILFSLIQKF----------------- 434

Query: 413 XXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFL 472
                        N+        N         S+ + ++ DE HRSHG  + +++   L
Sbjct: 435 -------------NTATRLPECVNT--------SSDIIVLIDEGHRSHGGESHQRMKLAL 473

Query: 473 SGNTRQSSKITYFSFTCTPTTQCLEMFGIASKGDKIPFYTYSMHNALEDKLILDVIKQFS 532
                      + +FT TP    L+     +K   I  + Y+M  A+EDK +  ++ +  
Sbjct: 474 PN-------AAFVAFTGTPL---LKDDKTTNKFGPI-VHAYTMQRAVEDKTVTPLLYEE- 521

Query: 533 TLRLNEVNVDSFXXXXXXXXXXXXXXXXXXXXXXXXYNEEAXXXXXXXXXXXXXXXXXXX 592
             R  E+NV                            NE A                   
Sbjct: 522 --RPPELNV----------------------------NERAIDSWFDRITEGLTDKQKAD 551

Query: 593 XXXCFLNGIDFDDDQYEDETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSG 652
               F N         + +       I   A  I NH  + KN  D      +  L C  
Sbjct: 552 LKRKFAN---------KGQVYNADDRIHLIALDIANH--FVKNIDDGL----KGQLACDS 596

Query: 653 RKGLLIYKKYLDQLISQLPIEEQFDS-IAAFSPFIQEGKSMVESD--PMINGKYAP-YGT 708
           +   + YKKYLD+          F+S +   SP  +EG + V+    P +   +    G+
Sbjct: 597 KAAAIKYKKYLDEA-------GLFESAVVISSPDTREGNTEVDESGIPDVTKWWQDNVGS 649

Query: 709 DKNRGIIKALND--DKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLS 766
              +   + + D  DK  +++LLIV DKL TGFDEP  +++Y+DK LKG N +Q ++R++
Sbjct: 650 QDEQVYTRNVIDRFDKDDSLKLLIVVDKLLTGFDEPKNAVLYIDKPLKGHNLIQAIARVN 709

Query: 767 RISKDKEATAIIDFVNTEKDMDFIFNLYKHSTS 799
           R+   K+   +ID+    K++D     Y+   S
Sbjct: 710 RLHPKKKFGLLIDYRGILKELDTTIEDYQKLAS 742


>gi|48243623|gb|AAT40773.1| putative type I site-specific
           restriction-modification enzyme HsdR [Haemophilus
           influenzae]
          Length = 567

 Score = 88.2 bits (217), Expect = 1e-15
 Identities = 78/263 (29%), Positives = 115/263 (43%), Gaps = 73/263 (27%)

Query: 303 YLIQHAAGSGKSLTIACLVYFLYKLEINNKHK----------FDTIVIMNDRSQLDAQLG 352
           YLIQH+AGSGKS +I  L Y L +    N+            FD+++++ DR  LD QL 
Sbjct: 5   YLIQHSAGSGKSNSITWLAYQLIEAYPRNEKAANGREADRPIFDSVIVVTDRRLLDKQLR 64

Query: 353 DVVIKFLTSNGITNFCRPK-SSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXX 411
           D +  F   + + N   P  SS +L+  +  GK +++ITT+QKF  ++   +        
Sbjct: 65  DNIKDF---SEVKNIVAPALSSAELRQSLEQGK-KIIITTIQKFPFIVDGIA-------- 112

Query: 412 XXXXXXXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEF 471
                                   +  D QF       A+I DEAH S   S  + ++  
Sbjct: 113 ------------------------DLGDKQF-------AVIIDEAHSSQSGSAHDNMNRA 141

Query: 472 LSG--------------NTRQSSKI----TYFSFTCTPTTQCLEMFGIASKGDKI-PFYT 512
           +                 T QS K+    +YF+FT TP    LE FG      K  PF+ 
Sbjct: 142 IGKTEDLDAEDVQDLILQTMQSRKMHGNASYFAFTATPKNSTLEKFGEKQADGKFKPFHL 201

Query: 513 YSMHNALEDKLILDVIKQFSTLR 535
           YSM  A+E+  ILDVI  ++T +
Sbjct: 202 YSMKQAIEEGFILDVIANYTTYK 224


 Score = 80.9 bits (198), Expect = 2e-13
 Identities = 66/201 (32%), Positives = 99/201 (48%), Gaps = 19/201 (9%)

Query: 614 QDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIE 673
           + Q  I TKA+ +++H +Y          K++ M+V    +  + Y + L  L++     
Sbjct: 254 RSQQTIDTKAEIMLDHFIYQVFNRKKLKGKAKGMVVTQNIETAIRYFQALKHLLAGRG-- 311

Query: 674 EQFDSIAAFSPF-IQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVA 732
             F    AFS   + +G    E++  +NG +A   T       K   D  ++  RLL+VA
Sbjct: 312 NPFKIAIAFSGSKVVDGVEYTEAE--MNG-FAESET-------KEYFDQDEY--RLLVVA 359

Query: 733 DKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFVNTEKDMD 788
           +K  TGFD+P L  MYVDKKL G   VQ LSRL+R    +SK  E   ++DF N+ +D+ 
Sbjct: 360 NKYLTGFDQPKLCAMYVDKKLSGVLCVQALSRLNRSANKLSKRTEDLFVLDFFNSVEDIQ 419

Query: 789 FIFNLYKHSTSLRDLQDTQTL 809
             F  +  STSL    D   L
Sbjct: 420 QAFEPFYTSTSLSQATDVNVL 440


>gi|23491885|dbj|BAC16860.1| conserved hypothetical protein
           [Corynebacterium efficiens YS-314]
 ref|NP_736660.1| hypothetical protein CE0050 [Corynebacterium efficiens YS-314]
          Length = 1064

 Score = 88.2 bits (217), Expect = 1e-15
 Identities = 80/300 (26%), Positives = 124/300 (40%), Gaps = 78/300 (26%)

Query: 271 PRTHQWDTI--VSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLE 328
           PR HQW  +  ++    +  P H          NYLIQH+AGSGK+ +IA   + L  L 
Sbjct: 290 PRFHQWRAVTKIAAAAAVEGPGH----------NYLIQHSAGSGKTDSIAWTAHRLASLH 339

Query: 329 INNKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTSNG----ITNFCRPKSSKQLQNEILSG 383
                K FDT+ ++ DR  LD QL D V +  T  G    I +      + +L   + SG
Sbjct: 340 DKAGEKVFDTVFVIADRQVLDRQLQDAVQQLETVAGTFQAIDSGGDGSKTSRLTQVLTSG 399

Query: 384 KNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDDNQFS 443
             ++V  T+Q F   L                                     +  NQ  
Sbjct: 400 TAKIVGVTLQTFPHAL---------------------------------AALKDPANQDK 426

Query: 444 FKSNSVAIISDEAHRSHGKSTTEKLHEFL------------------------SGNTRQS 479
                 A+I+DEAH S   S+ + + E L                        + ++   
Sbjct: 427 LAGRRFAVIADEAHSSQTGSSAKAVREMLYLTREIPEFDLEEPGADQDALAAMAAHSDAD 486

Query: 480 SKITYFSFTCTPTTQCLEMFG---IASKGDKIPFYTYSMHNALEDKLILDVIKQFSTLRL 536
            +++YF+FT TP  + LE+FG   +A++ +  PF  Y M  A+E+  ILDV+K ++T  +
Sbjct: 487 KRLSYFAFTATPKAKTLELFGRRNLATQ-ELEPFDLYPMKQAIEEGFILDVLKNYTTYEM 545


 Score = 75.9 bits (185), Expect = 7e-12
 Identities = 60/200 (30%), Positives = 101/200 (50%), Gaps = 14/200 (7%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQ--LPIEEQF 676
           + +K   I+ H   T   H     +++AM+V + R   + Y +  ++ I++  LP++   
Sbjct: 598 VASKVDVILEHYRSTVQPHLGG--RAKAMVVTASRAAAVRYARAFEKAITERNLPLQ--- 652

Query: 677 DSIAAFSPFIQEGK--SMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADK 734
            ++ AFS  + +    S+  +   +  + +     K R +      D Q    +L+VA+K
Sbjct: 653 -TLVAFSGEVPDPDVASLPGTAQKMVTEMSMNPALKGRDLASVFAQDGQ---NILVVANK 708

Query: 735 LQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLY 794
            QTGFD+P L  MYVDK+L G  AVQTLSRL+R +  K  T ++DFVN  + +   F  Y
Sbjct: 709 YQTGFDQPLLVGMYVDKQLSGIAAVQTLSRLNRRAPGKTDTYVLDFVNDAEQILESFKTY 768

Query: 795 KHSTSLRDLQDTQTLISSLI 814
             +  +    D   L++ LI
Sbjct: 769 YEAAHIETESDPD-LVADLI 787


>gi|48853355|ref|ZP_00307528.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Cytophaga hutchinsonii]
          Length = 1034

 Score = 87.0 bits (214), Expect = 3e-15
 Identities = 78/299 (26%), Positives = 128/299 (42%), Gaps = 85/299 (28%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL--E 328
           PR HQ D +  ++ D        + K      YLIQH+AGSGKS +I    Y L +   E
Sbjct: 299 PRYHQLDVVRKLVDD--------ASKKGVGQTYLIQHSAGSGKSNSITWAAYQLIETYPE 350

Query: 329 INN--------KHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK-SSKQLQNE 379
            +N           FD+++++ DR  LD Q+ + + +F   + + N   P  SSK+L+  
Sbjct: 351 SDNLPGSKGVSNPLFDSVIVVTDRRLLDKQIRENIKEF---SEVKNIVAPAYSSKELKES 407

Query: 380 ILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDD 439
           + SGK +++ITT+QKF  ++                                       D
Sbjct: 408 LESGK-KIIITTIQKFPFIV---------------------------------------D 427

Query: 440 NQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFL---SGNTRQSSKI-------------- 482
                     A+I DEAH S   ++   +++ +   +G+  +  +I              
Sbjct: 428 GIADLSDKRFAVIIDEAHSSQSGTSAGTMNQVMGNAAGDATEDEEIDPQDKILEAMRARK 487

Query: 483 -----TYFSFTCTPTTQCLEMFGIASKGDKI-PFYTYSMHNALEDKLILDVIKQFSTLR 535
                +Y +FT TP    LE FG+  +  K  PF+ YSM  A+E+  ILDV+  ++T +
Sbjct: 488 MRGNASYLAFTATPKNNTLERFGVKQEDGKYKPFHLYSMKQAIEEGFILDVLANYTTYK 546


 Score = 80.1 bits (196), Expect = 4e-13
 Identities = 67/208 (32%), Positives = 109/208 (52%), Gaps = 19/208 (9%)

Query: 614 QDQMLIKTKAQFIINHLLYTKNKHDNKLF-KSRAMLVCSGRKGLLIYKKYLDQLISQLPI 672
           Q++  I TKA  +++H + TK  +  +L  K++AM+     +  + Y   L +L+ +   
Sbjct: 576 QNKQTIATKADIMLDHFM-TKLVNTKRLKGKAKAMVATQSIESAINYYFALKKLLDEQG- 633

Query: 673 EEQFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDK-------NRGIIKALNDDKQFN 725
              F    AFS  +++ K +  ++  +N   A   T K       +  I + +N D+   
Sbjct: 634 -NPFKIAIAFSG-VKKIKGIEYTEDSLNDFPAHLDTAKPSDPGYISDKIARFVNMDEY-- 689

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFV 781
            R+L+VA+K  TG+D+P L+ MYVDKKL+G  AVQ LSRL+R    + K  E   I+DF 
Sbjct: 690 -RMLVVANKYLTGYDQPKLTAMYVDKKLQGVIAVQALSRLNRSADKLGKKTEDLFILDFY 748

Query: 782 NTEKDMDFIFNLYKHSTSLRDLQDTQTL 809
           NT +D+   F+ +  +TSL    D   L
Sbjct: 749 NTTEDIKKAFDPFYTATSLSKATDINVL 776


>gi|1679868|emb|CAA68059.1| Sty SBLI [Salmonella enterica]
          Length = 1088

 Score = 86.7 bits (213), Expect = 4e-15
 Identities = 117/525 (22%), Positives = 203/525 (38%), Gaps = 119/525 (22%)

Query: 308 AAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNF 367
           A G GKS T+  L   L K   + +    TIV++ DR+ LD QL   +       G    
Sbjct: 327 ATGCGKSYTMQFLTRLLMK---SVEFASPTIVLITDRTDLDDQLSAQMCNAKNYIGDDTV 383

Query: 368 CRPKSSKQLQNEILSGKNR--VVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYA 425
             P +S+    E L+G+N   V +TT+ KF+                             
Sbjct: 384 V-PVTSRDDLREKLAGRNSGGVFLTTIHKFT----------------------------- 413

Query: 426 NSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS------ 479
                      ED    S +SN + I SDEAHRS      + + +  SG  R++      
Sbjct: 414 -----------EDTELLSERSNIICI-SDEAHRSQVNLDQKVIVDKESGRVRKTYGFAKY 461

Query: 480 -----SKITYFSFTCTPTTQCLEMFGIASKGDKIPFYTYSMHNALEDKLILDVIKQFSTL 534
                 + TY  FT TP    +++FG+          +Y+M  +++D++ + ++ +    
Sbjct: 462 LHDSLPRATYVGFTGTPIDATMDVFGVIVD-------SYTMTESVQDEITVRIVYEG--- 511

Query: 535 RLNEVNVDSFXXXXXXXXXXXXXXXXXXXXXXXXYNEEAXXXXXXXXXXXXXXXXXXXXX 594
           R  +V +D+                         Y EE                      
Sbjct: 512 RAAKVILDA-----------------GKLEEVEKYYEECANAGTNEWQIDESKKASAT-- 552

Query: 595 XCFLNGIDFDDDQYEDETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRK 654
              +N I  D+D+           +K  A+    H  Y K   +    K +AM VC+ R+
Sbjct: 553 ---MNAILGDEDR-----------LKALAEDFAKH--YEKRVAEGSTVKGKAMFVCASRE 596

Query: 655 -GLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDP----MINGK------Y 703
                Y++  D   +   +++  + +       QE K +  S+     M  GK      Y
Sbjct: 597 IAWDFYRQLKDFRPAWFEVKQAPEGVELTG---QEEKELPPSEMVKMVMTRGKDDDAKLY 653

Query: 704 APYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLS 763
              G+ + R  +     + + N ++ IV D   TGFD P L  +Y+DK L+  N +QT+S
Sbjct: 654 DLLGSKEYRKELDKQFKNAKSNFKIAIVVDMWLTGFDVPELDTIYIDKPLQKHNLIQTIS 713

Query: 764 RLSRISKDKEATAIIDFVNTEKDMDFIFNLYKH--STSLRDLQDT 806
           R++R  + K    ++D++  ++ M+    +Y    +T+  D+Q +
Sbjct: 714 RVNRKMEGKSKGLVVDYIGIKRQMNQALAMYSRIDATNFEDIQQS 758


>gi|16124876|ref|NP_419440.1| hypothetical protein CC0623
           [Caulobacter crescentus CB15]
 gb|AAK22608.1| conserved hypothetical protein [Caulobacter crescentus CB15]
 pir||D87326 conserved hypothetical protein CC0623 [imported] - Caulobacter
           crescentus
          Length = 964

 Score = 85.9 bits (211), Expect = 7e-15
 Identities = 70/263 (26%), Positives = 114/263 (42%), Gaps = 66/263 (25%)

Query: 303 YLIQHAAGSGKSLTIACLVYFLYKLEINNKHK-FDTIVIMNDRSQLDAQLGDVVIKFLTS 361
           +LIQH+AGSGK+ +IA   + L +L     +K F ++++++DR+ +D QL + +  F  +
Sbjct: 301 FLIQHSAGSGKTNSIAWSAHILSELHDAADNKLFSSVIVVSDRNVIDTQLQEALFDFQRT 360

Query: 362 NGITNFCRPKS---SKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXX 418
            G+    + +S   S QL   + +GK +V++ T+Q F   L                   
Sbjct: 361 AGVVATIKSESGSKSGQLAEALAAGK-KVIVCTIQTFPFALEAVR--------------- 404

Query: 419 XXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLS----- 473
                                 Q + +  + A+I+DEAH S       KL + LS     
Sbjct: 405 ---------------------EQAATQGKAFAVIADEAHSSQTGEAAGKLKQMLSDAEIA 443

Query: 474 -------------------GNTRQSSKITYFSFTCTPTTQCLEMFGIA-SKGDKIPFYTY 513
                              G       ITY +FT TP  + L++FG   + G   PF+ Y
Sbjct: 444 ELADGGEIDFETQLTANMAGRAASDKSITYVAFTATPKAKTLQLFGRPDADGLPAPFHVY 503

Query: 514 SMHNALEDKLILDVIKQFSTLRL 536
           SM  A+E+  ILDV+K ++   L
Sbjct: 504 SMRQAIEEGFILDVLKNYTPYSL 526


 Score = 79.3 bits (194), Expect = 6e-13
 Identities = 65/222 (29%), Positives = 105/222 (47%), Gaps = 32/222 (14%)

Query: 599 NGIDFDDDQYEDETIQDQML---------IKTKAQFIINHLLYTKNKHDNKLFKSRAMLV 649
           +G ++D+ + E +T    ++         I  K Q ++ H  Y +N       K++AM+V
Sbjct: 533 DGQEWDEREVERDTALKGIMRWVRLHPYNIAQKVQIVVEH--YRENVQPLLDGKAKAMVV 590

Query: 650 CSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESD------PMINGKY 703
              R+  + +   ++  I ++       ++ AFS  + + +    S+      P +NG  
Sbjct: 591 LGSRQEAVRWSVAINAYIRKMGYP--LGTLVAFSGEVNDDELGAVSETSAALNPGLNG-- 646

Query: 704 APYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLS 763
                   R I  A    K  N   L+VA+K QTGFD+P L  MYVD++L G  AVQTLS
Sbjct: 647 --------RDIRDAF---KAENFHFLLVANKFQTGFDQPLLCGMYVDRRLAGIQAVQTLS 695

Query: 764 RLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQD 805
           RL+R    K+ T I+DF N+ +D+   F  Y  +  L    D
Sbjct: 696 RLNRAHPGKDTTYILDFCNSAEDILTAFRTYYETAELEAATD 737


>gi|54308080|ref|YP_129100.1| putative type I restriction enzyme
           HsdR [Photobacterium profundum SS9]
 emb|CAG19298.1| putative type I restriction enzyme HsdR [Photobacterium profundum]
          Length = 1016

 Score = 82.0 bits (201), Expect = 9e-14
 Identities = 78/296 (26%), Positives = 125/296 (41%), Gaps = 81/296 (27%)

Query: 269 FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKLE 328
           F PR  Q D +  ++ D        + K+     YLIQH+AGSGKS +I    Y L +  
Sbjct: 299 FFPRYQQMDVVRKLIAD--------AEKNGVGQTYLIQHSAGSGKSNSITWAAYQLIEAY 350

Query: 329 INNKHK----------FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRP-KSSKQLQ 377
             N+            FD+++++ DR  LD QL + + +F   + + N   P   S +L+
Sbjct: 351 PGNETAAGGKSLEQPLFDSVIVVTDRRILDKQLRENIQQF---SEVKNIIVPANKSSELK 407

Query: 378 NEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
           + +  GK +++ITT+QKF  ++   S                                + 
Sbjct: 408 SALEQGK-KIIITTIQKFPFIIDGIS--------------------------------DL 434

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSS----------------- 480
            D +F       A+I DEAH S   S  + ++  + GN  + +                 
Sbjct: 435 SDKRF-------AVIIDEAHSSQSGSAHDNMNRAM-GNVEEEADAQDKILEAMKSRKMRG 486

Query: 481 KITYFSFTCTPTTQCLEMFGIAS-KGDKIPFYTYSMHNALEDKLILDVIKQFSTLR 535
             +Y +FT TP    LE FG     G   PF+ YSM  A+E+  ILDV+  ++T +
Sbjct: 487 NASYLAFTATPKNTTLEKFGQKQPDGSYKPFHLYSMKQAIEEDFILDVLASYNTYK 542


 Score = 77.4 bits (189), Expect = 2e-12
 Identities = 68/204 (33%), Positives = 98/204 (47%), Gaps = 25/204 (12%)

Query: 614 QDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRA--MLVCSGRKGLLIYKKYLDQLISQLP 671
           + Q  I  KA+ +++H +   +  + KL K +A  M+V    +  + Y K + + + +  
Sbjct: 572 RQQQTINIKAEIMLDHFI--PHVVNPKLLKGKAKGMVVTQNIETAIRYYKAICRQLKERG 629

Query: 672 IEEQFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDK--QFNIRLL 729
               F  I AFS     G+  V+           Y   +  G       DK      RLL
Sbjct: 630 --NPFKVILAFS-----GEKTVDG--------IEYTESEVNGFADTETKDKFDTDEYRLL 674

Query: 730 IVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSR----ISKDKEATAIIDFVNTEK 785
           IVA+K  TGFD+P LS MYVDKKL+G  AVQ LSRL+R    + K  E   I+DF N+  
Sbjct: 675 IVANKYLTGFDQPKLSAMYVDKKLQGVLAVQALSRLNRSAQKLDKRTEDLFILDFFNSTD 734

Query: 786 DMDFIFNLYKHSTSLRDLQDTQTL 809
           D+   F+ +  STSL +  D   L
Sbjct: 735 DIKKSFDDFYTSTSLSEATDVNVL 758


>gi|45522101|ref|ZP_00173617.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Methylobacillus flagellatus KT]
          Length = 1022

 Score = 82.0 bits (201), Expect = 9e-14
 Identities = 72/271 (26%), Positives = 119/271 (43%), Gaps = 17/271 (6%)

Query: 271 PRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFLYKL-EI 329
           PR HQ D +  ++       H +   S    NYLIQH+AGSGKS +IA L + L  L + 
Sbjct: 287 PRYHQLDAVRKMVA------HAREQGSGR--NYLIQHSAGSGKSNSIAWLAHRLASLHDA 338

Query: 330 NNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPKSSKQLQNEILSGKNRVVI 389
            ++  F ++V++ DR  LD QL + + +F    G+      + ++QL   +  G   ++I
Sbjct: 339 QDQKVFHSVVVVTDRRVLDQQLQNTIYQFEHKTGVVEKI-DEDTQQLARALSQG-TPIII 396

Query: 390 TTMQKF---SQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNEDD--NQFSF 444
           TT+QKF   SQ L                         A+S              N+   
Sbjct: 397 TTIQKFPFISQALSTLERKGAGVKIDTAGRRFAVIVDEAHSSQSGETATALRGMLNRDGI 456

Query: 445 KSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKITYFSFTCTPTTQCLEMFGIASK 504
           ++   A +SDE            L        RQ + +++F+FT TP  +   +F     
Sbjct: 457 EAAIAAQMSDEEDDDLSAEARAALLRDALKRARQPN-LSFFAFTATPKFKTKVLFDEPGP 515

Query: 505 GDKIPFYTYSMHNALEDKLILDVIKQFSTLR 535
               PF+ Y+M  A+E+  I+DV++ ++T +
Sbjct: 516 SGTSPFHEYTMRQAIEEGFIMDVLQNYTTYK 546


 Score = 81.3 bits (199), Expect = 2e-13
 Identities = 56/167 (33%), Positives = 83/167 (49%), Gaps = 8/167 (4%)

Query: 643 KSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDPMINGK 702
           +++AM+V   R   + YK   DQ I          S+ AFS  +++       DP     
Sbjct: 603 RAKAMVVTGSRLAAVKYKLAFDQYIKANGYTG-IRSLVAFSGTVED-----PDDP--GAS 654

Query: 703 YAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTL 762
           Y     +      +     ++ + R+L+VA+K QTGFD+P L  MYV KKL G  AVQTL
Sbjct: 655 YTEVSMNDGLAESELPETFERDDYRVLLVAEKYQTGFDQPLLQTMYVVKKLAGVQAVQTL 714

Query: 763 SRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDTQTL 809
           SRL+R++  K  T ++DF N E D+   F  Y  +T + +  D   L
Sbjct: 715 SRLNRMAPGKARTFVLDFANEEDDIFKAFKPYYEATPIGENADPHRL 761


>gi|56900920|gb|AAW31753.1| type I restriction enzyme HsdR [Vibrio
           cholerae]
          Length = 339

 Score = 81.3 bits (199), Expect = 2e-13
 Identities = 79/279 (28%), Positives = 120/279 (42%), Gaps = 77/279 (27%)

Query: 269 FKPRTHQWDTIVSVLKDLISPQHMKSFKSENRYNYLIQHAAGSGKSLTIACLVYFL---Y 325
           F PR HQ D +V  L D  S   +          YLIQH+AGSGKS +I    Y L   Y
Sbjct: 110 FFPRYHQMD-VVRRLVDHCSVNGVGQ-------TYLIQHSAGSGKSNSITWAAYQLIETY 161

Query: 326 KLEIN---NKHK----FDTIVIMNDRSQLDAQLGDVVIKFLTSNGITNFCRPK-SSKQLQ 377
            +  +   ++ K    FD+++++ DR  LD QL D + +F   + + N   P   S +L+
Sbjct: 162 PISDDLPGSRGKEVPLFDSVIVVTDRRLLDKQLRDNIKEF---SEVKNIVAPAFKSSELK 218

Query: 378 NEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXXYANSXXXXXXXNNE 437
           + + +GK +++ITT+QKF  ++                                      
Sbjct: 219 SALENGK-KIIITTIQKFPYIV-------------------------------------- 239

Query: 438 DDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFL--------------SGNTRQSSKIT 483
            D          A+I DEAH S      +KL+E +              + N +  S  +
Sbjct: 240 -DGIADLSDRRFAVIIDEAHSSQDGHNQDKLNEAMGFVSEDVLDKALQSAKNRKMRSNAS 298

Query: 484 YFSFTCTPTTQCLEMFG-IASKGDKIPFYTYSMHNALED 521
           YF+FT TP    LE FG   + G  +PF+ YSM  A+E+
Sbjct: 299 YFAFTATPKNTTLEKFGQRQADGTYVPFHLYSMKQAIEE 337


>gi|53692041|ref|ZP_00347095.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Desulfovibrio desulfuricans G20]
          Length = 984

 Score = 75.9 bits (185), Expect = 7e-12
 Identities = 111/491 (22%), Positives = 185/491 (37%), Gaps = 127/491 (25%)

Query: 304 LIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           LI H  GSGKSL    +V+   KL ++ K    T++I+ DR  LD Q    +     +  
Sbjct: 303 LIWHFQGSGKSLL---MVFAAQKLRMHRKLGNPTVMIVVDRIDLDTQ----ITATFNAAD 355

Query: 364 ITNFCRPKSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXX 423
           I N     + ++LQ  + +   +++ITT+ KF                            
Sbjct: 356 IPNMIGAATRQELQGLLAADTRKIIITTIHKFG--------------------------- 388

Query: 424 YANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQSSKIT 483
                        E D + + +SN + ++ DEAHR+       K+ + L           
Sbjct: 389 -------------EADGRLNERSNIIVMV-DEAHRTQEGDLGRKMRDALPN-------AF 427

Query: 484 YFSFTCTPTTQ----CLEMFGIASKGDKIPFYTYSMHNALEDKLILDVIKQFSTLRLNEV 539
            F  T TP  +        FG A + ++     YS  +++ DK  L +  +   ++L+ +
Sbjct: 428 LFGLTGTPINKRDRNTFWAFG-ADEDEQGYMSRYSFQDSIRDKATLPLHFEAVDVKLH-I 485

Query: 540 NVDSFXXXXXXXXXXXXXXXXXXXXXXXXYNEEAXXXXXXXXXXXXXXXXXXXXXXCFLN 599
           N D+                           +EA                          
Sbjct: 486 NKDAI--------------------------DEAYSQMTDELS----------------- 502

Query: 600 GIDFDDDQYEDETIQDQMLIKTKAQF--IINHLL-YTKNKHDNKLFKSRAMLVCSGRKGL 656
             + D D       +  +LIK  A+   I  H++ + + K +   FK  A +V   R+  
Sbjct: 503 --ELDRDDLAKRAAKMAVLIKAPARVNAICQHIVKHFQEKVEPNGFK--AQVVTFDRECC 558

Query: 657 LIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIK 716
           ++YKK +D+L+            +A     Q GKS    D     K A    +K     +
Sbjct: 559 VLYKKAMDELVGP--------EASAIVMHTQGGKS----DEYAEWKLAKDEEEKLLDRFR 606

Query: 717 ALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATA 776
             ND  +F    LIV  KL TGFD P L +MY+DK +K  N +Q + R +R+   K    
Sbjct: 607 DPNDPLKF----LIVTSKLLTGFDAPILQVMYLDKPMKDHNLLQAICRTNRVYSGKTHGL 662

Query: 777 IIDFVNTEKDM 787
           ++D++    D+
Sbjct: 663 VVDYLGIFDDV 673


>gi|58038318|ref|YP_190287.1| Putative type I restriction enzyme
           [Gluconobacter oxydans 621H]
 gb|AAW59631.1| Putative type I restriction enzyme [Gluconobacter oxydans 621H]
          Length = 1098

 Score = 74.3 bits (181), Expect = 2e-11
 Identities = 47/152 (30%), Positives = 76/152 (49%), Gaps = 17/152 (11%)

Query: 645 RAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSP--------FIQEGKSMVESD 696
           + ML C  +   + YK+YLD++         F+S    SP         + E K     D
Sbjct: 586 KGMLACDSKLSAIRYKQYLDEV-------GLFESAIVMSPPDTREGHTEVDERKQPEIQD 638

Query: 697 PMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGA 756
              N   +    D  RG+I+    D   ++RL+IV DKL TGFDEP  +++Y+DK LK  
Sbjct: 639 WWKNNVGSQSEDDYTRGVIERFEKDP--DLRLIIVVDKLLTGFDEPKNAVLYIDKPLKQH 696

Query: 757 NAVQTLSRLSRISKDKEATAIIDFVNTEKDMD 788
           N +Q ++R++R+   K+   +ID+     ++D
Sbjct: 697 NLLQAIARVNRLHDQKKHGLLIDYRGILAELD 728


>gi|49257051|dbj|BAD24840.1| type I restriction-modification system
           endonuclease homologue [Staphylococcus aureus]
          Length = 1039

 Score = 73.9 bits (180), Expect = 3e-11
 Identities = 54/203 (26%), Positives = 91/203 (44%), Gaps = 13/203 (6%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           ++  A+ II H      +   +  K + M+V   R+   I     D++I   P     D 
Sbjct: 541 VEALAKDIIQHF-----ETRQQAMKGKGMIVTMSRR---IAVDLYDEIIRLKPEWHSDDD 592

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTG 738
                  +  G S   SDP    ++   G  K R +++    D    ++L+IV D   TG
Sbjct: 593 DKGVIKVVMTGSS---SDPTSFQRHI--GPKKRRNVLEKRMKDVNDELQLVIVRDMWLTG 647

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHST 798
           FD PS+  MY+DK +KG N +Q ++R++R+ KDK    I+D+V   + +      Y  S 
Sbjct: 648 FDVPSMHTMYIDKPMKGHNLMQAIARVNRVFKDKPGGLIVDYVGIAESLKEALKEYTESD 707

Query: 799 SLRDLQDTQTLISSLIMAGRSIQ 821
             +   DT   +  +++    IQ
Sbjct: 708 QAQTAIDTDKAVELMLLKYDVIQ 730


 Score = 42.4 bits (98), Expect = 0.083
 Identities = 36/99 (36%), Positives = 53/99 (53%), Gaps = 15/99 (15%)

Query: 304 LIQHAAGSGKSLTIACLVYFLYKL--EINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTS 361
           +I H  GSGKSLT   +V+F  KL   +NN     T+V++ DR+ LD QL    +K    
Sbjct: 298 VIWHTQGSGKSLT---MVFFSGKLIQMLNNP----TLVVVTDRNDLDNQLYSTFVKSKGR 350

Query: 362 NGITNFCR--PKSS---KQLQNEILSGKNRVVITTMQKF 395
           +G     R  PK +   K+L++ +      ++ TTMQKF
Sbjct: 351 SG-KGLLRQTPKQAETRKELKSLLSVESGGIIFTTMQKF 388


>gi|19919302|gb|AAM08232.1| putative type I restriction enzyme
           [Legionella pneumophila]
          Length = 1090

 Score = 71.6 bits (174), Expect = 1e-10
 Identities = 57/210 (27%), Positives = 101/210 (47%), Gaps = 36/210 (17%)

Query: 645 RAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSP-FIQEGKSMVES-------- 695
           +  L C  +   + YKKYLD +         F+S    SP   +EG + V+         
Sbjct: 585 KGQLACDSKLSAIRYKKYLDDV-------GLFESAVVMSPPDTREGNTEVDESKIPEVTK 637

Query: 696 ---DPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKK 752
              D + N +   Y     R +I+   +D+  +++++IV DKL TGFDEP  +++Y+DK+
Sbjct: 638 WWKDNVGNQEEREY----TRQVIERFANDE--DLKIIIVVDKLLTGFDEPRNTVLYIDKQ 691

Query: 753 LKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTS-------LRDLQD 805
           LK  N +Q ++R++R+ + K+   +ID+    K++D     Y+   S       + DL  
Sbjct: 692 LKEHNLIQAIARVNRLHEKKKFGLLIDYRGILKELDTTIAKYQDLASRTQQGYDITDLDG 751

Query: 806 TQTLISS----LIMAGRSIQKVFEKEKRAG 831
               +SS    L +   ++ K+F + K  G
Sbjct: 752 LYHQMSSEYKRLPLLYDALWKIFNEVKNKG 781


 Score = 51.2 bits (121), Expect = 2e-04
 Identities = 59/247 (23%), Positives = 101/247 (40%), Gaps = 60/247 (24%)

Query: 304 LIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           +I H  GSGKSLT   +V+    L +++  K   IV++ DR  L+ QL +        +G
Sbjct: 340 VIWHTTGSGKSLT---MVFLSRALILHDSLKQCRIVVVTDRVDLERQLSNTFHSGGEFSG 396

Query: 364 ITNFCRPK--SSKQLQNEILSGKNRVVITTMQKFSQL--LPNKSXXXXXXXXXXXXXXXX 419
            T F   K  S ++L  +I  G  RV+ + +QKF+    LP                   
Sbjct: 397 KTEFSNAKATSGRRLAEQISKGTERVIFSLIQKFNTAAKLP------------------- 437

Query: 420 XXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQS 479
                                Q    S  + ++ DE HRS G     ++ + L       
Sbjct: 438 ---------------------QCKNLSRDIVVLVDEGHRSQGGENYIRMKQAL------- 469

Query: 480 SKITYFSFTCTPTTQCLEMFGIASKGDKIPFYTYSMHNALEDKLILDVI--KQFSTLRLN 537
               + +FT TP    L+     +K   I  + Y++  A+EDK++  ++  ++   L +N
Sbjct: 470 PNAAFIAFTGTP---LLKEDKTTNKFGSI-IHAYTIQRAVEDKMVTPLLYEERRPELMVN 525

Query: 538 EVNVDSF 544
           +  +DS+
Sbjct: 526 DRAIDSW 532


>gi|46141215|ref|ZP_00145477.2| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Psychrobacter sp. 273-4]
          Length = 1060

 Score = 69.7 bits (169), Expect = 5e-10
 Identities = 58/207 (28%), Positives = 100/207 (48%), Gaps = 31/207 (14%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+  A+  + H  Y     +    K +AM VCS R+    Y+ Y D +     +  ++D 
Sbjct: 538 IRAIAEDFVQH--YEARVAEGTTQKGKAMFVCSSREAA--YQLYKDII----ELHPEWDE 589

Query: 679 IAAFSPFIQEGKSMVESD-----PMINGKYA-PYGTDKNRGIIKALND-------DKQF- 724
           I A     +EG ++ ESD     P+   K       D  +G+   L +       D+QF 
Sbjct: 590 IRA----CEEGSTLSESDKKKVMPIERVKMIMTRNKDDEKGMWDNLGNKDYRKELDRQFK 645

Query: 725 ----NIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDF 780
               N ++ IV D   TGFD P L  +Y+DK L+  + +QT+SR++R  + KE+  ++D+
Sbjct: 646 NGKSNFKIAIVVDMWLTGFDVPFLDTIYIDKPLQKHSLIQTISRVNRKFEGKESGLVVDY 705

Query: 781 VNTEKDMDFIFNLYKHSTSLRDLQDTQ 807
           +  +K M+   + Y     +++L+D Q
Sbjct: 706 IGIKKQMNLALSHYT-GGQIQNLEDIQ 731


>gi|42525030|ref|NP_970410.1| type I restriction-modification system
           restriction subunit [Bdellovibrio bacteriovorus HD100]
 emb|CAE81064.1| type I restriction-modification system restriction subunit
           [Bdellovibrio bacteriovorus HD100]
          Length = 1075

 Score = 69.3 bits (168), Expect = 6e-10
 Identities = 61/211 (28%), Positives = 99/211 (46%), Gaps = 22/211 (10%)

Query: 609 EDETIQDQMLIKTKAQFIINHLLYTKN-KHDNKLFKS--RAMLVCSGRKGLLIYKKYLDQ 665
           E+E ++    I+  A  I NH  Y +N + + +  K   +  L  + +   L YK+YLD+
Sbjct: 555 EEEILKADDRIRAIALDIKNH--YRENFRQEGEPTKDFMKGQLATNTKGEALAYKRYLDE 612

Query: 666 LISQLPIEEQFDSIAAFSPFIQEGKSMVESD--PMI----NGKYAPYGTDKNRGIIKALN 719
              +        S+    P ++EG   +E D  P +    +     YG+ K        N
Sbjct: 613 FGIRT-------SLVISPPEMREGSKSIEDDDRPEVVKFWDECMKRYGSAKKYQDTIIDN 665

Query: 720 DDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIID 779
             K     +LIV DKL TGFD P  +++Y+DKKLK  N +Q ++R++R+   K    +ID
Sbjct: 666 FKKDDVPEILIVVDKLLTGFDAPRNAVLYIDKKLKEHNVLQAIARVNRLYPRKSEGLVID 725

Query: 780 FVNTEKDM----DFIFNLYKHSTSLRDLQDT 806
           +     DM    DF   +  +   L D++DT
Sbjct: 726 YRGIFDDMTQAIDFYRKMEGNDIDLEDIKDT 756


>gi|15795175|ref|NP_247089.2| Type I restriction-modification
           enzyme, helicase subunit [Methanocaldococcus jannaschii
           DSM 2661]
          Length = 1163

 Score = 68.6 bits (166), Expect = 1e-09
 Identities = 58/210 (27%), Positives = 101/210 (47%), Gaps = 18/210 (8%)

Query: 623 AQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQL--PIEEQFDSIA 680
           A++I++ +     + D + FK +AM+V   R G + +KK LD+ + +      E++  + 
Sbjct: 638 AKYIVDRI-----EEDTENFKFKAMVVAVNRLGCVRFKKALDKYLKEKFGDEAEKWAEVV 692

Query: 681 AFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKA--LNDDKQFNIRLLIVADKLQTG 738
                 +E K ++E    +  +      ++   II+   LN +   N ++LIV D L TG
Sbjct: 693 MTYHHNEEEKEIIEYMKKLKKERNSNDFNEINQIIREEFLNSE---NPKILIVTDMLLTG 749

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHST 798
           FD P L +MY+DK L G   +Q ++R +R   DKE   I+D V   K +     LY    
Sbjct: 750 FDAPRLKVMYLDKPLYGHRLLQAIARTNRPYPDKEFGLIVDSVGLFKVLTETMALYNMLA 809

Query: 799 SLRDLQDTQTLISSLIMAGRSIQKVFEKEK 828
                +D +  + S      SI ++F++ K
Sbjct: 810 EEEIREDFKNNLIS------SIDEIFQEFK 833


>gi|2826244|gb|AAB98104.1| type I restriction-modification enzyme 2,
           R subunit [Methanocaldococcus jannaschii DSM 2661]
 sp|Q57588|Y124_METJA Hypothetical protein MJ0124
          Length = 1075

 Score = 68.6 bits (166), Expect = 1e-09
 Identities = 58/210 (27%), Positives = 101/210 (47%), Gaps = 18/210 (8%)

Query: 623 AQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQL--PIEEQFDSIA 680
           A++I++ +     + D + FK +AM+V   R G + +KK LD+ + +      E++  + 
Sbjct: 550 AKYIVDRI-----EEDTENFKFKAMVVAVNRLGCVRFKKALDKYLKEKFGDEAEKWAEVV 604

Query: 681 AFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKA--LNDDKQFNIRLLIVADKLQTG 738
                 +E K ++E    +  +      ++   II+   LN +   N ++LIV D L TG
Sbjct: 605 MTYHHNEEEKEIIEYMKKLKKERNSNDFNEINQIIREEFLNSE---NPKILIVTDMLLTG 661

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHST 798
           FD P L +MY+DK L G   +Q ++R +R   DKE   I+D V   K +     LY    
Sbjct: 662 FDAPRLKVMYLDKPLYGHRLLQAIARTNRPYPDKEFGLIVDSVGLFKVLTETMALYNMLA 721

Query: 799 SLRDLQDTQTLISSLIMAGRSIQKVFEKEK 828
                +D +  + S      SI ++F++ K
Sbjct: 722 EEEIREDFKNNLIS------SIDEIFQEFK 745


>gi|2129335|pir||D64315 type I restriction enzyme homolog -
           Methanococcus jannaschii
          Length = 1163

 Score = 68.6 bits (166), Expect = 1e-09
 Identities = 58/210 (27%), Positives = 101/210 (47%), Gaps = 18/210 (8%)

Query: 623 AQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQL--PIEEQFDSIA 680
           A++I++ +     + D + FK +AM+V   R G + +KK LD+ + +      E++  + 
Sbjct: 638 AKYIVDRI-----EEDTENFKFKAMVVAVNRLGCVRFKKALDKYLKEKFGDEAEKWAEVV 692

Query: 681 AFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKA--LNDDKQFNIRLLIVADKLQTG 738
                 +E K ++E    +  +      ++   II+   LN +   N ++LIV D L TG
Sbjct: 693 MTYHHNEEEKEIIEYMKKLKKERNSNDFNEINQIIREEFLNSE---NPKILIVTDMLLTG 749

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHST 798
           FD P L +MY+DK L G   +Q ++R +R   DKE   I+D V   K +     LY    
Sbjct: 750 FDAPRLKVMYLDKPLYGHRLLQAIARTNRPYPDKEFGLIVDSVGLFKVLTETMALYNMLA 809

Query: 799 SLRDLQDTQTLISSLIMAGRSIQKVFEKEK 828
                +D +  + S      SI ++F++ K
Sbjct: 810 EEEIREDFKNNLIS------SIDEIFQEFK 833


>gi|53688022|ref|ZP_00109780.2| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Nostoc punctiforme PCC 73102]
          Length = 1078

 Score = 68.2 bits (165), Expect = 1e-09
 Identities = 62/250 (24%), Positives = 116/250 (45%), Gaps = 34/250 (13%)

Query: 603 FDDDQYED----ETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLI 658
           F++  Y+D    +  +++  I   AQ+I +  L  K+ +       +AML     K  + 
Sbjct: 615 FEEPNYKDADKKQIYENRDRINAIAQYIAD--LLVKDVYRQIRGTGKAMLAVYSIKAAIA 672

Query: 659 YKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNR-GIIKA 717
           YK+ + Q  ++L               +++ K    +D  I+  Y+    ++N  G+   
Sbjct: 673 YKQAVTQHFNEL---------------VKQPKYAKYTDAPIHVVYSSNQDEQNATGLNNG 717

Query: 718 LNDDK---QFNIR---LLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKD 771
           L ++K    F ++   L+IV  KLQTGFDE  L  +++DK++KG  A+Q +SR+ R +K 
Sbjct: 718 LTEEKVLESFALKKNGLIIVVAKLQTGFDEKRLHTLFLDKEIKGIAAIQAISRVDRTAKY 777

Query: 772 KEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDTQTLISSLIMAGRSIQKVFEKEKRAG 831
           K    I+DF     ++  I + ++H        D            R ++ +F + K++G
Sbjct: 778 KNDCKIVDFSYNNVNVQNIKDAFEH------FSDVVVSDFDPFSDKRVLEILFAELKKSG 831

Query: 832 LITPFCKAIL 841
           + T F +  L
Sbjct: 832 VYTEFFETFL 841


 Score = 59.7 bits (143), Expect = 5e-07
 Identities = 61/249 (24%), Positives = 96/249 (38%), Gaps = 53/249 (21%)

Query: 297 SENRYNYLIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVI 356
           ++N Y+ L+Q+AAG GKS  I      L  L  N ++ +D I+I+ DR QL +Q+     
Sbjct: 383 NKNVYSLLLQYAAGFGKSNIIGWTALQLKDLRRNGQYVYDKIMIVVDRLQLRSQIDS--- 439

Query: 357 KFLTSNGITNFCRPKSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXX 416
           K L  N         ++K+   + L    R+VI  +QKF  +                  
Sbjct: 440 KMLNMNIDNRMYIEANNKKTFQDALQSDTRLVIVNLQKFGAV------------------ 481

Query: 417 XXXXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSG-- 474
                               + D         +  + DE HRS+     E++        
Sbjct: 482 ----------------REMLDADTLQKLAQLRIVFLIDEIHRSNSGDQHEEMVNIFDELQ 525

Query: 475 ---------NTRQSSKITYFSFTCTPTTQCLEMFGIAS---KGDK--IPFYTYSMHNALE 520
                    N  ++ K     FT TP    L  FG  S   + +K  +PF +Y+M  A+E
Sbjct: 526 TPFDQQTQYNQTRTKKNLIVGFTATPDDHTLARFGEFSGYAESEKLWVPFDSYTMKEAIE 585

Query: 521 DKLILDVIK 529
           D  IL+ +K
Sbjct: 586 DGFILNPLK 594


>gi|6968968|emb|CAB73965.1| putative type I restriction enzyme R
           protein [Campylobacter jejuni subsp. jejuni NCTC 11168]
 gb|AAM00885.1| HsdR [Campylobacter jejuni]
 pir||C81302 probable type I site-specific deoxyribonuclease (EC 3.1.21.3) chain
           R Cj1549c [imported] - Campylobacter jejuni  (strain
           NCTC 11168)
 ref|NP_282680.1| putative type I restriction enzyme R protein [Campylobacter jejuni
           subsp. jejuni NCTC 11168]
          Length = 1031

 Score = 67.4 bits (163), Expect = 2e-09
 Identities = 34/106 (32%), Positives = 59/106 (55%), Gaps = 5/106 (4%)

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
           ++L+IV D   TGFD PSL  +Y+DK +KG N +Q ++R++R+ KDK +  ++D++    
Sbjct: 630 LKLVIVVDMWLTGFDVPSLHTLYIDKPMKGHNLMQAIARVNRVYKDKPSGLVVDYLGIAS 689

Query: 786 DMDFIFNLYKHSTSLRDLQDTQ-----TLISSLIMAGRSIQKVFEK 826
           D+    + Y  S    D+  TQ      +I  L +  + +   +EK
Sbjct: 690 DLKKALSFYSESGGKGDIITTQEDAVKVMIEKLEIVEQMLPNDYEK 735


>gi|19881241|gb|AAM00850.1| HsdR [Campylobacter jejuni]
          Length = 1031

 Score = 67.4 bits (163), Expect = 2e-09
 Identities = 34/106 (32%), Positives = 59/106 (55%), Gaps = 5/106 (4%)

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
           ++L+IV D   TGFD PSL  +Y+DK +KG N +Q ++R++R+ KDK +  ++D++    
Sbjct: 630 LKLVIVVDMWLTGFDVPSLHTLYIDKPMKGHNLMQAIARVNRVYKDKPSGLVVDYLGIAS 689

Query: 786 DMDFIFNLYKHSTSLRDLQDTQ-----TLISSLIMAGRSIQKVFEK 826
           D+    + Y  S    D+  TQ      +I  L +  + +   +EK
Sbjct: 690 DLKKALSFYSESGGKGDIITTQEDAVKVMIEKLEIVEQMLPNDYEK 735


>gi|19881235|gb|AAM00845.1| HsdR [Campylobacter jejuni]
          Length = 1031

 Score = 67.4 bits (163), Expect = 2e-09
 Identities = 34/106 (32%), Positives = 59/106 (55%), Gaps = 5/106 (4%)

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
           ++L+IV D   TGFD PSL  +Y+DK +KG N +Q ++R++R+ KDK +  ++D++    
Sbjct: 630 LKLVIVVDMWLTGFDVPSLHTLYIDKPMKGHNLMQAIARVNRVYKDKPSGLVVDYLGIAS 689

Query: 786 DMDFIFNLYKHSTSLRDLQDTQ-----TLISSLIMAGRSIQKVFEK 826
           D+    + Y  S    D+  TQ      +I  L +  + +   +EK
Sbjct: 690 DLKKALSFYSESGGKGDIITTQEDAVKVMIEKLEIVEQMLPNDYEK 735


>gi|19881208|gb|AAM00823.1| HsdR [Campylobacter jejuni]
          Length = 1031

 Score = 67.4 bits (163), Expect = 2e-09
 Identities = 34/106 (32%), Positives = 59/106 (55%), Gaps = 5/106 (4%)

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
           ++L+IV D   TGFD PSL  +Y+DK +KG N +Q ++R++R+ KDK +  ++D++    
Sbjct: 630 LKLVIVVDMWLTGFDVPSLHTLYIDKPMKGHNLMQAIARVNRVYKDKPSGLVVDYLGIAS 689

Query: 786 DMDFIFNLYKHSTSLRDLQDTQ-----TLISSLIMAGRSIQKVFEK 826
           D+    + Y  S    D+  TQ      +I  L +  + +   +EK
Sbjct: 690 DLKKALSFYSESGGKGDIITTQEDAVKVMIEKLEIVEQMLPNDYEK 735


>gi|19881202|gb|AAM00818.1| HsdR [Campylobacter jejuni]
          Length = 1031

 Score = 67.4 bits (163), Expect = 2e-09
 Identities = 34/106 (32%), Positives = 59/106 (55%), Gaps = 5/106 (4%)

Query: 726 IRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEK 785
           ++L+IV D   TGFD PSL  +Y+DK +KG N +Q ++R++R+ KDK +  ++D++    
Sbjct: 630 LKLVIVVDMWLTGFDVPSLHTLYIDKPMKGHNLMQAIARVNRVYKDKPSGLVVDYLGIAS 689

Query: 786 DMDFIFNLYKHSTSLRDLQDTQ-----TLISSLIMAGRSIQKVFEK 826
           D+    + Y  S    D+  TQ      +I  L +  + +   +EK
Sbjct: 690 DLKKALSFYSESGGKGDIITTQEDAVKVMIEKLEIVEQMLPNDYEK 735


>gi|56419917|ref|YP_147235.1| type I restriction-modification system
           specificity protein [Geobacillus kaustophilus HTA426]
 dbj|BAD75667.1| type I restriction-modification system specificity protein
           [Geobacillus kaustophilus HTA426]
          Length = 1021

 Score = 66.2 bits (160), Expect = 5e-09
 Identities = 47/179 (26%), Positives = 81/179 (44%), Gaps = 13/179 (7%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           ++  A+ I++H     N         +AM+VC  R+  +    YL ++    P+E     
Sbjct: 522 LRRLAKDIVSHFNQVANPFQ------KAMIVCMSRRIAVDLYNYLREIPECPPVEVIMTG 575

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTG 738
             +  P     K   E  P    KYA   + K +  +KA   D +  ++ +IV D   TG
Sbjct: 576 DVSKDP-----KEWREVQP--GSKYAHIKSKKEQEEVKAKLRDPEDPLKFVIVVDMWLTG 628

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
            D   LS +YVDK +KG N +Q ++R++R+   KE   ++D++     +    + Y  S
Sbjct: 629 MDAKPLSYLYVDKPMKGHNLMQAIARVNRVFPGKEGGVVVDYIGIATSLKEATHKYTQS 687


 Score = 40.0 bits (92), Expect = 0.41
 Identities = 50/232 (21%), Positives = 90/232 (38%), Gaps = 64/232 (27%)

Query: 304 LIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           ++ H  GSGKSL+   +++F   L  + + +  +IVI  DR+ LD QL D  ++    + 
Sbjct: 293 VMYHTTGSGKSLS---MLFFAGILSRHPEMENPSIVIQVDRNDLDEQLHDTFVQ--GKSY 347

Query: 364 ITNFCRPKSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXXXXXXX 423
           I +    ++++QL+  +     +++ +T++KF +                          
Sbjct: 348 IGHVHHAENTEQLRTLLRGEAGQIIFSTIEKFRK-------------------------- 381

Query: 424 YANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRSHGKSTTEKLHEFLSGNTRQ----S 479
                        E  +    +  ++ +I+DEAHR+           F  G   Q     
Sbjct: 382 ----------KEGETRHPVLSERRNIVVIADEAHRTQA--------GFAGGFAAQLRYAL 423

Query: 480 SKITYFSFTCTPT----TQCLEMFGIASKGDKIPFYTYSMHNALEDKLILDV 527
              +Y  FT TP         E+FG          + Y M  A+EDK  L +
Sbjct: 424 PNASYIGFTGTPVDFVDNNTEEIFGHI-------IHRYDMAQAVEDKATLPI 468


>gi|15669399|ref|NP_248209.1| type I restriction-modification enzyme
           1, R subunit [Methanocaldococcus jannaschii DSM 2661]
 gb|AAB99215.1| type I restriction-modification enzyme 1, R subunit
           [Methanocaldococcus jannaschii DSM 2661]
 sp|Q58611|YC14_METJA Hypothetical protein MJ1214
          Length = 1018

 Score = 66.2 bits (160), Expect = 5e-09
 Identities = 51/198 (25%), Positives = 94/198 (46%), Gaps = 12/198 (6%)

Query: 611 ETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQL 670
           + ++ +  IK  +++I  H+L      D + FK +AM+V   RK  +++KKYLD+ + + 
Sbjct: 522 DLLKSESSIKEASKYIAEHIL-----EDTENFKFKAMVVAQDRKSCILFKKYLDEYLKEK 576

Query: 671 PIEEQFDSIAAFSPFIQEGKSMVES-DPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLL 729
                 +       +I      +E+    I  KY     + N+   +   + +  N ++L
Sbjct: 577 IKNYNENWTQVVITYIHNDDVEIENYKKEIEKKYGKNVDELNKKWTEDFINKE--NPKIL 634

Query: 730 IVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIID----FVNTEK 785
           IV  KL TGFD P L  +Y+ + LK    +Q  +R +R +K+K+   I+D     +   K
Sbjct: 635 IVNKKLLTGFDAPILKTIYIHQFLKDYLLLQASARANRPAKNKKYGLIVDLTGILIENYK 694

Query: 786 DMDFIFNLYKHSTSLRDL 803
                +NLY+     +D+
Sbjct: 695 KAIENYNLYRDEAINKDI 712


>gi|2129336|pir||E64451 type I restriction enzyme homolog -
           Methanococcus jannaschii
          Length = 1018

 Score = 66.2 bits (160), Expect = 5e-09
 Identities = 51/198 (25%), Positives = 94/198 (46%), Gaps = 12/198 (6%)

Query: 611 ETIQDQMLIKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQL 670
           + ++ +  IK  +++I  H+L      D + FK +AM+V   RK  +++KKYLD+ + + 
Sbjct: 522 DLLKSESSIKEASKYIAEHIL-----EDTENFKFKAMVVAQDRKSCILFKKYLDEYLKEK 576

Query: 671 PIEEQFDSIAAFSPFIQEGKSMVES-DPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLL 729
                 +       +I      +E+    I  KY     + N+   +   + +  N ++L
Sbjct: 577 IKNYNENWTQVVITYIHNDDVEIENYKKEIEKKYGKNVDELNKKWTEDFINKE--NPKIL 634

Query: 730 IVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIID----FVNTEK 785
           IV  KL TGFD P L  +Y+ + LK    +Q  +R +R +K+K+   I+D     +   K
Sbjct: 635 IVNKKLLTGFDAPILKTIYIHQFLKDYLLLQASARANRPAKNKKYGLIVDLTGILIENYK 694

Query: 786 DMDFIFNLYKHSTSLRDL 803
                +NLY+     +D+
Sbjct: 695 KAIENYNLYRDEAINKDI 712


>gi|46142301|ref|ZP_00204237.1| COG0610: Type I site-specific
           restriction-modification system, R (restriction) subunit
           and related helicases [Methanococcoides burtonii DSM
           6242]
          Length = 1074

 Score = 65.9 bits (159), Expect = 7e-09
 Identities = 47/179 (26%), Positives = 81/179 (44%), Gaps = 13/179 (7%)

Query: 619 IKTKAQFIINHLLYTKNKHDNKLFKSRAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDS 678
           I+  AQ I+ H    K     ++F+ + M+V   R+   I     D++I   P     D 
Sbjct: 576 IRQIAQDIVTHFEQRK-----EVFEGKGMIVAMSRR---IAADLYDEIIKIKPQWHNDDF 627

Query: 679 IAAFSPFIQEGKSMVESDPMINGKYAPYGTDKNRGIIKALNDDKQFNIRLLIVADKLQTG 738
                  +    S     P I+  +    T + R  +     D +  ++L+IV D   TG
Sbjct: 628 KKGEIKIVMTSAS--SDGPKISKHHT---TKEQRRALAERMKDPEDELKLVIVRDMWLTG 682

Query: 739 FDEPSLSIMYVDKKLKGANAVQTLSRLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHS 797
           FD P L  MY+DK +KG   +Q ++R++R+ KDK  + ++D++    D+    + Y  S
Sbjct: 683 FDVPCLHTMYIDKPMKGHTLMQAIARVNRVYKDKPGSLVVDYLGIASDLKEALSFYSDS 741


>gi|53758784|gb|AAU93075.1| type I restriction-modification system,
           R subunit [Methylococcus capsulatus str. Bath]
 ref|YP_113334.1| type I restriction-modification system, R subunit [Methylococcus
           capsulatus str. Bath]
          Length = 1210

 Score = 65.5 bits (158), Expect = 9e-09
 Identities = 46/163 (28%), Positives = 80/163 (48%), Gaps = 14/163 (8%)

Query: 645 RAMLVCSGRKGLLIYKKYLDQLISQLPIEEQFDSIAAFSPFIQEGKSMVESDPMINGKYA 704
           +AM+VC  R+   I     D+++   P     D  A     +  G +   SDP    ++ 
Sbjct: 623 KAMIVCMSRR---ICVALYDEIVRLRPDWHSADDKAGAIKIVMTGAA---SDPP---EWQ 673

Query: 705 PY-GTDKNRGIIKALNDDKQFNIRLLIVADKLQTGFDEPSLSIMYVDKKLKGANAVQTLS 763
           P+ G    R ++     D    ++L+IV D   TGFD P +  MYVDK ++G   +Q ++
Sbjct: 674 PHIGNKARRDLLAKRARDPNDPLKLVIVRDMWLTGFDAPCMHTMYVDKPMRGHGLMQAIA 733

Query: 764 RLSRISKDKEATAIIDFVNTEKDMDFIFNLYKHSTSLRDLQDT 806
           R++R+ +DK A  I+D++   +++      Y    S RD ++T
Sbjct: 734 RVNRVFRDKPAGLIVDYIGIAQNLKSALAQY----SPRDRENT 772


 Score = 41.6 bits (96), Expect = 0.14
 Identities = 44/161 (27%), Positives = 63/161 (38%), Gaps = 17/161 (10%)

Query: 304 LIQHAAGSGKSLTIACLVYFLYKLEINNKHKFDTIVIMNDRSQLDAQLGDVVIKFLTSNG 363
           +I H  GSGKSL +A     L K  +       T+V++ DR+ LD QL      F T + 
Sbjct: 327 VIWHTQGSGKSLLMAFYAGMLVKHPLLEN---PTLVVITDRNDLDDQL------FATFSM 377

Query: 364 ITNFCRP-----KSSKQLQNEILSGKNRVVITTMQKFSQLLPNKSXXXXXXXXXXXXXXX 418
             +  R      +S + LQ  +      V+ TT+QKF  L P  S               
Sbjct: 378 CRDLIRQTPVQAESREHLQALLNRASGGVIFTTLQKFGPLSPGPS---PNGRGEDTGPAP 434

Query: 419 XXXXXYANSXXXXXXXNNEDDNQFSFKSNSVAIISDEAHRS 459
                +A+           D         +V +I+DEAHRS
Sbjct: 435 NVREEHASPSPSGRGVGVRDTPPPLTLRRNVVVIADEAHRS 475


>gi|21228306|ref|NP_634228.1| hypothetical protein MM2204