]> git.donarmstrong.com Git - biopieces.git/commitdiff
added RestrictEnz.pm
authormartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Fri, 14 Aug 2009 15:37:38 +0000 (15:37 +0000)
committermartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Fri, 14 Aug 2009 15:37:38 +0000 (15:37 +0000)
git-svn-id: http://biopieces.googlecode.com/svn/trunk@634 74ccb610-7750-0410-82ae-013aeee3265d

code_perl/Maasha/RestrictEnz.pm [new file with mode: 0644]

diff --git a/code_perl/Maasha/RestrictEnz.pm b/code_perl/Maasha/RestrictEnz.pm
new file mode 100644 (file)
index 0000000..9b1dd27
--- /dev/null
@@ -0,0 +1,1008 @@
+package Maasha::RestrictEnz;
+
+
+# Copyright (C) 2006-2007 Martin A. Hansen.
+
+# This program is free software; you can redistribute it and/or
+# modify it under the terms of the GNU General Public License
+# as published by the Free Software Foundation; either version 2
+# of the License, or (at your option) any later version.
+
+# This program is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU General Public License for more details.
+
+# You should have received a copy of the GNU General Public License
+# along with this program; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
+
+# http://www.gnu.org/copyleft/gpl.html
+
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+# This module contains routines for matching restriction enzyme cleavage sites.
+
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+use warnings;
+use strict;
+use Data::Dumper;
+
+use Exporter;
+
+use vars qw( @ISA @EXPORT @EXPORT_OK );
+
+@ISA = qw( Exporter ) ;
+
+use Inline ( C => <<'END_C', DIRECTORY => $ENV{ "BP_TMP" } );
+
+/*
+# UIPAC ambiguity codes for nucleotides:
+#
+# http://droog.gs.washington.edu/parc/images/iupac.html
+#
+#  ACGTUMRWSYKVHDBN
+# A1000011100011101
+# C0100010011011011
+# G0010001010110111
+# T0001100101101111
+# U0001100101101111
+# M1100000000000000
+# R1010000000000000
+# W1001100000000000
+# S0110000000000000
+# Y0101100000000000
+# K0011100000000000
+# V1110000000000000
+# H1101100000000000
+# D1011100000000000
+# B0111100000000000
+# N1111100000000000
+*/
+
+
+char ambi_match[16][16] = {
+    "1000011100011101",
+    "0100010011011011",
+    "0010001010110111",
+    "0001100101101111",
+    "0001100101101111",
+    "1100000000000000",
+    "1010000000000000",
+    "1001100000000000",
+    "0110000000000000",
+    "0101100000000000",
+    "0011100000000000",
+    "1110000000000000",
+    "1101100000000000",
+    "1011100000000000",
+    "0111100000000000",
+    "1111100000000000"
+};
+
+
+int hash( char c )
+{
+    switch ( toupper( c ) )
+    {
+        case 'A': return 0;
+        case 'C': return 1;
+        case 'G': return 2;
+        case 'T': return 3;
+        case 'U': return 4;
+        case 'M': return 5;
+        case 'R': return 6;
+        case 'W': return 7;
+        case 'S': return 8;
+        case 'Y': return 9;
+        case 'K': return 10;
+        case 'V': return 11;
+        case 'H': return 12;
+        case 'D': return 13;
+        case 'B': return 14;
+        case 'N': return 0;
+        default: return -1;
+    }
+}
+
+
+void scan( char *seq, char *pat, int seq_len, int pat_len )
+{
+    int i;
+
+    Inline_Stack_Vars;
+    Inline_Stack_Reset;
+
+    for ( i = 0; i < seq_len - pat_len; i++ )
+    {
+        if ( match( &seq[ i ], pat, pat_len ) ) {
+          Inline_Stack_Push( sv_2mortal( newSViv( i ) ) );
+        }
+    }
+
+    Inline_Stack_Done;
+}
+
+
+int match( char *seq1, char *seq2, int len )
+{
+    int  i = 0;
+    char c1;
+    char c2;
+
+    while ( i < len )
+    {
+        c1 = seq1[ i ];
+        c2 = seq2[ i ];
+
+        if ( ambi_match[ hash( c1 ) ][ hash( c2 ) ] == '0' ) {
+            return 0;
+        }
+
+        i++;
+    }
+
+    return 1;
+}
+
+END_C
+
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+sub re_scan
+{
+    my ( $seq,   # sequence to scan
+         $re,    # hashref with RE info
+       ) = @_; 
+
+    my ( @matches );
+
+    @matches = scan( $seq, $re->{ "pattern" }, length $seq, $re->{ "len" } );
+
+    return wantarray ? @matches: \@matches;
+}
+
+
+sub parse_re_data
+{
+    my ( @lines, $line, @fields, @re_data );
+
+    @lines = <DATA>;
+
+    chomp @lines;
+
+    foreach $line ( @lines )
+    {
+        next if $line =~ /^(#|$)/;
+
+        @fields = split " ", $line;
+
+        push @re_data, {
+            name    => $fields[ 0 ],
+            pattern => $fields[ 1 ],
+            len     => $fields[ 2 ],
+            ncuts   => $fields[ 3 ],
+            blunt   => $fields[ 4 ],
+            c1      => $fields[ 5 ],
+            c2      => $fields[ 6 ],
+            c3      => $fields[ 7 ],
+            c4      => $fields[ 8 ],
+        };
+    }
+
+    return wantarray ? @re_data : \@re_data;
+}
+
+
+# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+
+
+__DATA__
+
+# From REBASE:
+# ftp://ftp.neb.com/pub/rebase/emboss_e.908
+
+# REBASE version 908                                              emboss_e.908
+#  
+#     =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
+#     REBASE, The Restriction Enzyme Database   http://rebase.neb.com
+#     Copyright (c)  Dr. Richard J. Roberts, 2009.   All rights reserved.
+#     =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
+#  
+# Rich Roberts                                                    Jul 28 2009
+#  
+# REBASE enzyme patterns for EMBOSS (embossre.enz)
+#
+# Format:
+# 
+# name<ws>pattern<ws>len<ws>ncuts<ws>blunt<ws>c1<ws>c2<ws>c3<ws>c4
+#
+# Where:
+# name = name of enzyme
+# pattern = recognition site
+# len = length of pattern
+# ncuts = number of cuts made by enzyme
+#         Zero represents unknown
+# blunt = true if blunt end cut, false if sticky
+# c1 = First 5' cut
+# c2 = First 3' cut
+# c3 = Second 5' cut
+# c4 = Second 3' cut
+#
+# Examples:
+# AAC^TGG -> 6 2 1 3 3 0 0
+# A^ACTGG -> 6 2 0 1 5 0 0
+# AACTGG  -> 6 0 0 0 0 0 0
+# AACTGG(-5/-1) -> 6 2 0 1 5 0 0
+# (8/13)GACNNNNNNTCA(12/7) -> 12 4 0 -9 -14 24 19
+#
+# i.e. cuts are always to the right of the given
+# residue and sequences are always with reference to
+# the 5' strand.
+# Sequences are numbered ... -3 -2 -1 1 2 3 ... with
+# the first residue of the pattern at base number 1.
+#
+#
+
+AanI   TTATAA  6       2       1       3       3       0       0
+AarI   CACCTGC 7       2       0       11      15      0       0
+AasI   GACNNNNNNGTC    12      2       0       7       5       0       0
+AatI   AGGCCT  6       2       1       3       3       0       0
+AatII  GACGTC  6       2       0       5       1       0       0
+AbsI   CCTCGAGG        8       2       0       2       6       0       0
+AccI   GTMKAC  6       2       0       2       4       0       0
+AccII  CGCG    4       2       1       2       2       0       0
+AccIII TCCGGA  6       2       0       1       5       0       0
+Acc16I TGCGCA  6       2       1       3       3       0       0
+Acc36I ACCTGC  6       2       0       10      14      0       0
+Acc65I GGTACC  6       2       0       1       5       0       0
+AccB1I GGYRCC  6       2       0       1       5       0       0
+AccB7I CCANNNNNTGG     11      2       0       7       4       0       0
+AccBSI CCGCTC  6       2       1       3       3       0       0
+AceIII cagctc  6       2       0       13      17      0       0
+AciI   CCGC    4       2       0       1       3       0       0
+AclI   AACGTT  6       2       0       2       4       0       0
+AclWI  GGATC   5       2       0       9       10      0       0
+AcoI   YGGCCR  6       2       0       1       5       0       0
+AcsI   RAATTY  6       2       0       1       5       0       0
+AcuI   CTGAAG  6       2       0       22      20      0       0
+AcvI   CACGTG  6       2       1       3       3       0       0
+AcyI   GRCGYC  6       2       0       2       4       0       0
+AdeI   CACNNNGTG       9       2       0       6       3       0       0
+AfaI   GTAC    4       2       1       2       2       0       0
+AfeI   AGCGCT  6       2       1       3       3       0       0
+AfiI   CCNNNNNNNGG     11      2       0       7       4       0       0
+AflII  CTTAAG  6       2       0       1       5       0       0
+AflIII ACRYGT  6       2       0       1       5       0       0
+AgeI   ACCGGT  6       2       0       1       5       0       0
+AgsI   TTSAA   5       2       0       3       2       0       0
+AhaIII tttaaa  6       2       1       3       3       0       0
+AhdI   GACNNNNNGTC     11      2       0       6       5       0       0
+AhlI   ACTAGT  6       2       0       1       5       0       0
+AjiI   CACGTC  6       2       1       3       3       0       0
+AjnI   CCWGG   5       2       0       -1      5       0       0
+AjuI   GAANNNNNNNTTGG  14      4       0       -8      -13     25      20
+AleI   CACNNNNGTG      10      2       1       5       5       0       0
+AlfI   GCANNNNNNTGC    12      4       0       -11     -13     24      22
+AloI   GAACNNNNNNTCC   13      4       0       -8      -13     25      20
+AluI   AGCT    4       2       1       2       2       0       0
+AluBI  AGCT    4       2       1       2       2       0       0
+AlwI   GGATC   5       2       0       9       10      0       0
+Alw21I GWGCWC  6       2       0       5       1       0       0
+Alw26I GTCTC   5       2       0       6       10      0       0
+Alw44I GTGCAC  6       2       0       1       5       0       0
+AlwFI  gaaaynnnnnrtg   13      0       0       0       0       0       0
+AlwNI  CAGNNNCTG       9       2       0       6       3       0       0
+Ama87I CYCGRG  6       2       0       1       5       0       0
+Aor13HI        TCCGGA  6       2       0       1       5       0       0
+Aor51HI        AGCGCT  6       2       1       3       3       0       0
+ApaI   GGGCCC  6       2       0       5       1       0       0
+ApaBI  gcannnnntgc     11      2       0       8       3       0       0
+ApaLI  GTGCAC  6       2       0       1       5       0       0
+ApeKI  GCWGC   5       2       0       1       4       0       0
+ApoI   RAATTY  6       2       0       1       5       0       0
+ApyPI  atcgac  6       2       0       26      24      0       0
+AquIII gaggag  6       2       0       26      24      0       0
+AquIV  grggaag 7       2       0       26      24      0       0
+ArsI   GACNNNNNNTTYG   13      4       0       -9      -14     24      19
+AscI   GGCGCGCC        8       2       0       2       6       0       0
+AseI   ATTAAT  6       2       0       2       4       0       0
+Asi256I        gatc    4       2       0       1       3       0       0
+AsiGI  ACCGGT  6       2       0       1       5       0       0
+AsiSI  GCGATCGC        8       2       0       5       3       0       0
+AspI   GACNNNGTC       9       2       0       4       5       0       0
+Asp700I        GAANNNNTTC      10      2       1       5       5       0       0
+Asp718I        GGTACC  6       2       0       1       5       0       0
+AspA2I CCTAGG  6       2       0       1       5       0       0
+AspCNI gccgc   5       0       0       0       0       0       0
+AspEI  GACNNNNNGTC     11      2       0       6       5       0       0
+AspLEI GCGC    4       2       0       3       1       0       0
+AspS9I GGNCC   5       2       0       1       4       0       0
+AssI   AGTACT  6       2       1       3       3       0       0
+AsuI   ggncc   5       2       0       1       4       0       0
+AsuII  TTCGAA  6       2       0       2       4       0       0
+AsuC2I CCSGG   5       2       0       2       3       0       0
+AsuHPI GGTGA   5       2       0       13      12      0       0
+AsuNHI GCTAGC  6       2       0       1       5       0       0
+AvaI   CYCGRG  6       2       0       1       5       0       0
+AvaII  GGWCC   5       2       0       1       4       0       0
+AvaIII atgcat  6       0       0       0       0       0       0
+AviII  TGCGCA  6       2       1       3       3       0       0
+AvrII  CCTAGG  6       2       0       1       5       0       0
+AxyI   CCTNAGG 7       2       0       2       5       0       0
+BaeI   ACNNNNGTAYC     11      4       0       -11     -16     23      18
+BaeGI  GKGCMC  6       2       0       5       1       0       0
+BalI   TGGCCA  6       2       1       3       3       0       0
+BamHI  GGATCC  6       2       0       1       5       0       0
+BanI   GGYRCC  6       2       0       1       5       0       0
+BanII  GRGCYC  6       2       0       5       1       0       0
+BanIII ATCGAT  6       2       0       2       4       0       0
+BarI   GAAGNNNNNNTAC   13      4       0       -8      -13     25      20
+BasI   CCANNNNNTGG     11      2       0       7       4       0       0
+BauI   CACGAG  6       2       0       1       5       0       0
+BbeI   GGCGCC  6       2       0       5       1       0       0
+Bbr7I  gaagac  6       2       0       13      17      0       0
+BbrPI  CACGTG  6       2       1       3       3       0       0
+BbsI   GAAGAC  6       2       0       8       12      0       0
+BbuI   GCATGC  6       2       0       5       1       0       0
+BbvI   GCAGC   5       2       0       13      17      0       0
+BbvII  gaagac  6       2       0       8       12      0       0
+Bbv12I GWGCWC  6       2       0       5       1       0       0
+BbvCI  CCTCAGC 7       2       0       2       5       0       0
+BccI   CCATC   5       2       0       9       10      0       0
+Bce83I cttgag  6       2       0       22      20      0       0
+BceAI  ACGGC   5       2       0       17      19      0       0
+BcefI  acggc   5       2       0       17      18      0       0
+BcgI   CGANNNNNNTGC    12      4       0       -11     -13     24      22
+BciVI  GTATCC  6       2       0       12      11      0       0
+BclI   TGATCA  6       2       0       1       5       0       0
+BcnI   CCSGG   5       2       0       2       3       0       0
+BcuI   ACTAGT  6       2       0       1       5       0       0
+BdaI   TGANNNNNNTCA    12      4       0       -11     -13     24      22
+BetI   wccggw  6       2       0       1       5       0       0
+BfaI   CTAG    4       2       0       1       3       0       0
+BfiI   ACTGGG  6       2       0       11      10      0       0
+BfmI   CTRYAG  6       2       0       1       5       0       0
+BfoI   RGCGCY  6       2       0       5       1       0       0
+BfrI   CTTAAG  6       2       0       1       5       0       0
+BfuI   GTATCC  6       2       0       12      11      0       0
+BfuAI  ACCTGC  6       2       0       10      14      0       0
+BfuCI  GATC    4       2       0       -1      4       0       0
+BglI   GCCNNNNNGGC     11      2       0       7       4       0       0
+BglII  AGATCT  6       2       0       1       5       0       0
+BinI   ggatc   5       2       0       9       10      0       0
+BisI   GCNGC   5       2       0       2       3       0       0
+BlnI   CCTAGG  6       2       0       1       5       0       0
+BlpI   GCTNAGC 7       2       0       2       5       0       0
+BlsI   GCNGC   5       2       0       3       2       0       0
+BmcAI  AGTACT  6       2       1       3       3       0       0
+Bme18I GGWCC   5       2       0       1       4       0       0
+Bme1390I       CCNGG   5       2       0       2       3       0       0
+BmeRI  GACNNNNNGTC     11      2       0       6       5       0       0
+BmeT110I       CYCGRG  6       2       0       1       5       0       0
+BmgI   gkgccc  6       0       0       0       0       0       0
+BmgBI  CACGTC  6       2       1       3       3       0       0
+BmgT120I       GGNCC   5       2       0       2       3       0       0
+BmiI   GGNNCC  6       2       1       3       3       0       0
+BmrI   ACTGGG  6       2       0       11      10      0       0
+BmrFI  CCNGG   5       2       0       2       3       0       0
+BmsI   GCATC   5       2       0       10      14      0       0
+BmtI   GCTAGC  6       2       0       5       1       0       0
+BmuI   ACTGGG  6       2       0       11      10      0       0
+BoxI   GACNNNNGTC      10      2       1       5       5       0       0
+BpiI   GAAGAC  6       2       0       8       12      0       0
+BplI   GAGNNNNNCTC     11      4       0       -9      -14     24      19
+BpmI   CTGGAG  6       2       0       22      20      0       0
+Bpu10I CCTNAGC 7       2       0       2       5       0       0
+Bpu14I TTCGAA  6       2       0       2       4       0       0
+Bpu1102I       GCTNAGC 7       2       0       2       5       0       0
+BpuAI  GAAGAC  6       2       0       8       12      0       0
+BpuEI  CTTGAG  6       2       0       22      20      0       0
+BpuMI  CCSGG   5       2       0       2       3       0       0
+BpvUI  CGATCG  6       2       0       4       2       0       0
+BsaI   GGTCTC  6       2       0       7       11      0       0
+Bsa29I ATCGAT  6       2       0       2       4       0       0
+BsaAI  YACGTR  6       2       1       3       3       0       0
+BsaBI  GATNNNNATC      10      2       1       5       5       0       0
+BsaHI  GRCGYC  6       2       0       2       4       0       0
+BsaJI  CCNNGG  6       2       0       1       5       0       0
+BsaMI  GAATGC  6       2       0       7       5       0       0
+BsaWI  WCCGGW  6       2       0       1       5       0       0
+BsaXI  ACNNNNNCTCC     11      4       0       -10     -13     21      18
+BsbI   caacac  6       2       0       27      25      0       0
+Bsc4I  CCNNNNNNNGG     11      2       0       7       4       0       0
+BscAI  gcatc   5       2       0       9       11      0       0
+BscGI  cccgt   5       0       0       0       0       0       0
+Bse1I  ACTGG   5       2       0       6       4       0       0
+Bse8I  GATNNNNATC      10      2       1       5       5       0       0
+Bse21I CCTNAGG 7       2       0       2       5       0       0
+Bse118I        RCCGGY  6       2       0       1       5       0       0
+BseAI  TCCGGA  6       2       0       1       5       0       0
+BseBI  CCWGG   5       2       0       2       3       0       0
+BseCI  ATCGAT  6       2       0       2       4       0       0
+BseDI  CCNNGG  6       2       0       1       5       0       0
+Bse3DI GCAATG  6       2       0       8       6       0       0
+BseGI  GGATG   5       2       0       7       5       0       0
+BseJI  GATNNNNATC      10      2       1       5       5       0       0
+BseLI  CCNNNNNNNGG     11      2       0       7       4       0       0
+BseMI  GCAATG  6       2       0       8       6       0       0
+BseMII CTCAG   5       2       0       15      13      0       0
+BseNI  ACTGG   5       2       0       6       4       0       0
+BsePI  GCGCGC  6       2       0       1       5       0       0
+BseRI  GAGGAG  6       2       0       16      14      0       0
+BseSI  GKGCMC  6       2       0       5       1       0       0
+BseXI  GCAGC   5       2       0       13      17      0       0
+BseX3I CGGCCG  6       2       0       1       5       0       0
+BseYI  CCCAGC  6       2       0       1       5       0       0
+BsgI   GTGCAG  6       2       0       22      20      0       0
+Bsh1236I       CGCG    4       2       1       2       2       0       0
+Bsh1285I       CGRYCG  6       2       0       4       2       0       0
+BshFI  GGCC    4       2       1       2       2       0       0
+BshNI  GGYRCC  6       2       0       1       5       0       0
+BshTI  ACCGGT  6       2       0       1       5       0       0
+BshVI  ATCGAT  6       2       0       2       4       0       0
+BsiI   cacgag  6       2       0       1       5       0       0
+BsiEI  CGRYCG  6       2       0       4       2       0       0
+BsiHKAI        GWGCWC  6       2       0       5       1       0       0
+BsiHKCI        CYCGRG  6       2       0       1       5       0       0
+BsiSI  CCGG    4       2       0       1       3       0       0
+BsiWI  CGTACG  6       2       0       1       5       0       0
+BsiYI  ccnnnnnnngg     11      2       0       7       4       0       0
+BslI   CCNNNNNNNGG     11      2       0       7       4       0       0
+BslFI  GGGAC   5       2       0       15      19      0       0
+BsmI   GAATGC  6       2       0       7       5       0       0
+BsmAI  GTCTC   5       2       0       6       10      0       0
+BsmBI  CGTCTC  6       2       0       7       11      0       0
+BsmFI  GGGAC   5       2       0       15      19      0       0
+BsnI   GGCC    4       2       1       2       2       0       0
+Bso31I GGTCTC  6       2       0       7       11      0       0
+BsoBI  CYCGRG  6       2       0       1       5       0       0
+Bsp13I TCCGGA  6       2       0       1       5       0       0
+Bsp19I CCATGG  6       2       0       1       5       0       0
+Bsp24I gacnnnnnntgg    12      4       0       -9      -14     24      19
+Bsp68I TCGCGA  6       2       1       3       3       0       0
+Bsp119I        TTCGAA  6       2       0       2       4       0       0
+Bsp120I        GGGCCC  6       2       0       1       5       0       0
+Bsp143I        GATC    4       2       0       -1      4       0       0
+Bsp1286I       GDGCHC  6       2       0       5       1       0       0
+Bsp1407I       TGTACA  6       2       0       1       5       0       0
+Bsp1720I       GCTNAGC 7       2       0       2       5       0       0
+BspACI CCGC    4       2       0       1       3       0       0
+BspCNI CTCAG   5       2       0       14      12      0       0
+BspDI  ATCGAT  6       2       0       2       4       0       0
+BspD6I gactc   5       2       0       9       11      0       0
+BspEI  TCCGGA  6       2       0       1       5       0       0
+BspFNI CGCG    4       2       1       2       2       0       0
+BspGI  ctggac  6       0       0       0       0       0       0
+BspHI  TCATGA  6       2       0       1       5       0       0
+BspLI  GGNNCC  6       2       1       3       3       0       0
+BspLU11I       acatgt  6       2       0       1       5       0       0
+BspMI  ACCTGC  6       2       0       10      14      0       0
+BspMII tccgga  6       2       0       1       5       0       0
+BspNCI ccaga   5       0       0       0       0       0       0
+BspOI  GCTAGC  6       2       0       5       1       0       0
+BspPI  GGATC   5       2       0       9       10      0       0
+BspQI  GCTCTTC 7       2       0       8       11      0       0
+BspTI  CTTAAG  6       2       0       1       5       0       0
+BspT104I       TTCGAA  6       2       0       2       4       0       0
+BspT107I       GGYRCC  6       2       0       1       5       0       0
+BspTNI GGTCTC  6       2       0       7       11      0       0
+BspXI  ATCGAT  6       2       0       2       4       0       0
+BsrI   ACTGG   5       2       0       6       4       0       0
+BsrBI  CCGCTC  6       2       1       3       3       0       0
+BsrDI  GCAATG  6       2       0       8       6       0       0
+BsrFI  RCCGGY  6       2       0       1       5       0       0
+BsrGI  TGTACA  6       2       0       1       5       0       0
+BsrSI  ACTGG   5       2       0       6       4       0       0
+BssAI  RCCGGY  6       2       0       1       5       0       0
+BssECI CCNNGG  6       2       0       1       5       0       0
+BssHII GCGCGC  6       2       0       1       5       0       0
+BssKI  CCNGG   5       2       0       -1      5       0       0
+BssMI  GATC    4       2       0       -1      4       0       0
+BssNI  GRCGYC  6       2       0       2       4       0       0
+BssNAI GTATAC  6       2       1       3       3       0       0
+BssSI  CACGAG  6       2       0       1       5       0       0
+BssT1I CCWWGG  6       2       0       1       5       0       0
+Bst6I  CTCTTC  6       2       0       7       10      0       0
+Bst98I CTTAAG  6       2       0       1       5       0       0
+Bst1107I       GTATAC  6       2       1       3       3       0       0
+BstACI GRCGYC  6       2       0       2       4       0       0
+BstAFI CTTAAG  6       2       0       1       5       0       0
+BstAPI GCANNNNNTGC     11      2       0       7       4       0       0
+BstAUI TGTACA  6       2       0       1       5       0       0
+BstBI  TTCGAA  6       2       0       2       4       0       0
+Bst2BI CACGAG  6       2       0       1       5       0       0
+BstBAI YACGTR  6       2       1       3       3       0       0
+Bst4CI ACNGT   5       2       0       3       2       0       0
+BstC8I GCNNGC  6       2       1       3       3       0       0
+BstDEI CTNAG   5       2       0       1       4       0       0
+BstDSI CCRYGG  6       2       0       1       5       0       0
+BstEII GGTNACC 7       2       0       1       6       0       0
+BstENI CCTNNNNNAGG     11      2       0       5       6       0       0
+BstF5I GGATG   5       2       0       7       5       0       0
+BstFNI CGCG    4       2       1       2       2       0       0
+BstH2I RGCGCY  6       2       0       5       1       0       0
+BstHHI GCGC    4       2       0       3       1       0       0
+BstKTI GATC    4       2       0       3       1       0       0
+BstMAI GTCTC   5       2       0       6       10      0       0
+BstMBI GATC    4       2       0       -1      4       0       0
+BstMCI CGRYCG  6       2       0       4       2       0       0
+BstMWI GCNNNNNNNGC     11      2       0       7       4       0       0
+BstNI  CCWGG   5       2       0       2       3       0       0
+BstNSI RCATGY  6       2       0       5       1       0       0
+BstOI  CCWGG   5       2       0       2       3       0       0
+BstPI  GGTNACC 7       2       0       1       6       0       0
+BstPAI GACNNNNGTC      10      2       1       5       5       0       0
+BstSCI CCNGG   5       2       0       -1      5       0       0
+BstSFI CTRYAG  6       2       0       1       5       0       0
+BstSLI GKGCMC  6       2       0       5       1       0       0
+BstSNI TACGTA  6       2       1       3       3       0       0
+BstUI  CGCG    4       2       1       2       2       0       0
+Bst2UI CCWGG   5       2       0       2       3       0       0
+BstV1I GCAGC   5       2       0       13      17      0       0
+BstV2I GAAGAC  6       2       0       8       12      0       0
+BstXI  CCANNNNNNTGG    12      2       0       8       4       0       0
+BstX2I RGATCY  6       2       0       1       5       0       0
+BstYI  RGATCY  6       2       0       1       5       0       0
+BstZI  CGGCCG  6       2       0       1       5       0       0
+BstZ17I        GTATAC  6       2       1       3       3       0       0
+Bsu15I ATCGAT  6       2       0       2       4       0       0
+Bsu36I CCTNAGG 7       2       0       2       5       0       0
+BsuRI  GGCC    4       2       1       2       2       0       0
+BsuTUI ATCGAT  6       2       0       2       4       0       0
+BtgI   CCRYGG  6       2       0       1       5       0       0
+BtgZI  GCGATG  6       2       0       16      20      0       0
+BthCI  gcngc   5       2       0       4       1       0       0
+BtrI   CACGTC  6       2       1       3       3       0       0
+BtsI   GCAGTG  6       2       0       8       6       0       0
+BtsCI  GGATG   5       2       0       7       5       0       0
+BtuMI  TCGCGA  6       2       1       3       3       0       0
+BveI   ACCTGC  6       2       0       10      14      0       0
+Cac8I  GCNNGC  6       2       1       3       3       0       0
+CaiI   CAGNNNCTG       9       2       0       6       3       0       0
+CauII  ccsgg   5       2       0       2       3       0       0
+CciI   TCATGA  6       2       0       1       5       0       0
+CciNI  GCGGCCGC        8       2       0       2       6       0       0
+CdiI   catcg   5       2       1       4       4       0       0
+CdpI   gcggag  6       2       0       26      24      0       0
+CelII  GCTNAGC 7       2       0       2       5       0       0
+CfoI   GCGC    4       2       0       3       1       0       0
+CfrI   YGGCCR  6       2       0       1       5       0       0
+Cfr9I  CCCGGG  6       2       0       1       5       0       0
+Cfr10I RCCGGY  6       2       0       1       5       0       0
+Cfr13I GGNCC   5       2       0       1       4       0       0
+Cfr42I CCGCGG  6       2       0       4       2       0       0
+ChaI   gatc    4       2       0       4       -1      0       0
+CjeI   ccannnnnngt     11      4       0       -9      -15     26      20
+CjeNII gagnnnnngt      10      0       0       0       0       0       0
+CjePI  ccannnnnnntc    12      4       0       -8      -14     26      20
+CjuI   caynnnnnrtg     11      0       0       0       0       0       0
+CjuII  caynnnnnctc     11      0       0       0       0       0       0
+ClaI   ATCGAT  6       2       0       2       4       0       0
+CpoI   CGGWCCG 7       2       0       2       5       0       0
+CseI   GACGC   5       2       0       10      15      0       0
+CsiI   ACCWGGT 7       2       0       1       6       0       0
+CspI   CGGWCCG 7       2       0       2       5       0       0
+Csp6I  GTAC    4       2       0       1       3       0       0
+Csp45I TTCGAA  6       2       0       2       4       0       0
+CspAI  ACCGGT  6       2       0       1       5       0       0
+CspCI  CAANNNNNGTGG    12      4       0       -12     -14     24      22
+CstMI  aaggag  6       2       0       26      24      0       0
+CviAII CATG    4       2       0       1       3       0       0
+CviJI  RGCY    4       2       1       2       2       0       0
+CviKI-1        RGCY    4       2       1       2       2       0       0
+CviQI  GTAC    4       2       0       1       3       0       0
+CviRI  tgca    4       2       1       2       2       0       0
+DdeI   CTNAG   5       2       0       1       4       0       0
+DinI   GGCGCC  6       2       1       3       3       0       0
+DpnI   GATC    4       2       1       2       2       0       0
+DpnII  GATC    4       2       0       -1      4       0       0
+DraI   TTTAAA  6       2       1       3       3       0       0
+DraII  RGGNCCY 7       2       0       2       5       0       0
+DraIII CACNNNGTG       9       2       0       6       3       0       0
+DraRI  caagnac 7       2       0       27      25      0       0
+DrdI   GACNNNNNNGTC    12      2       0       7       5       0       0
+DrdII  gaacca  6       0       0       0       0       0       0
+DrdIV  tacgac  6       2       0       26      24      0       0
+DriI   GACNNNNNGTC     11      2       0       6       5       0       0
+DsaI   ccrygg  6       2       0       1       5       0       0
+DseDI  GACNNNNNNGTC    12      2       0       7       5       0       0
+EaeI   YGGCCR  6       2       0       1       5       0       0
+EagI   CGGCCG  6       2       0       1       5       0       0
+Eam1104I       CTCTTC  6       2       0       7       10      0       0
+Eam1105I       GACNNNNNGTC     11      2       0       6       5       0       0
+EarI   CTCTTC  6       2       0       7       10      0       0
+EciI   GGCGGA  6       2       0       17      15      0       0
+Ecl136II       GAGCTC  6       2       1       3       3       0       0
+EclXI  CGGCCG  6       2       0       1       5       0       0
+Eco24I GRGCYC  6       2       0       5       1       0       0
+Eco31I GGTCTC  6       2       0       7       11      0       0
+Eco32I GATATC  6       2       1       3       3       0       0
+Eco47I GGWCC   5       2       0       1       4       0       0
+Eco47III       AGCGCT  6       2       1       3       3       0       0
+Eco52I CGGCCG  6       2       0       1       5       0       0
+Eco57I CTGAAG  6       2       0       22      20      0       0
+Eco72I CACGTG  6       2       1       3       3       0       0
+Eco81I CCTNAGG 7       2       0       2       5       0       0
+Eco88I CYCGRG  6       2       0       1       5       0       0
+Eco91I GGTNACC 7       2       0       1       6       0       0
+Eco105I        TACGTA  6       2       1       3       3       0       0
+Eco130I        CCWWGG  6       2       0       1       5       0       0
+Eco147I        AGGCCT  6       2       1       3       3       0       0
+EcoHI  ccsgg   5       2       0       -1      5       0       0
+EcoICRI        GAGCTC  6       2       1       3       3       0       0
+Eco57MI        CTGRAG  6       2       0       22      20      0       0
+EcoNI  CCTNNNNNAGG     11      2       0       5       6       0       0
+EcoO65I        GGTNACC 7       2       0       1       6       0       0
+EcoO109I       RGGNCCY 7       2       0       2       5       0       0
+EcoRI  GAATTC  6       2       0       1       5       0       0
+EcoRII CCWGG   5       2       0       -1      5       0       0
+EcoRV  GATATC  6       2       1       3       3       0       0
+EcoT14I        CCWWGG  6       2       0       1       5       0       0
+EcoT22I        ATGCAT  6       2       0       5       1       0       0
+EcoT38I        GRGCYC  6       2       0       5       1       0       0
+Eco53kI        GAGCTC  6       2       1       3       3       0       0
+EgeI   GGCGCC  6       2       1       3       3       0       0
+EheI   GGCGCC  6       2       1       3       3       0       0
+ErhI   CCWWGG  6       2       0       1       5       0       0
+EsaBC3I        tcga    4       2       1       2       2       0       0
+EsaSSI gaccac  6       0       0       0       0       0       0
+EspI   gctnagc 7       2       0       2       5       0       0
+Esp3I  CGTCTC  6       2       0       7       11      0       0
+FaeI   CATG    4       2       0       4       -1      0       0
+FaiI   YATR    4       2       1       2       2       0       0
+FalI   AAGNNNNNCTT     11      4       0       -9      -14     24      19
+FaqI   GGGAC   5       2       0       15      19      0       0
+FatI   CATG    4       2       0       -1      4       0       0
+FauI   CCCGC   5       2       0       9       11      0       0
+FauNDI CATATG  6       2       0       2       4       0       0
+FbaI   TGATCA  6       2       0       1       5       0       0
+FblI   GTMKAC  6       2       0       2       4       0       0
+FinI   gggac   5       0       0       0       0       0       0
+FmuI   ggncc   5       2       0       4       1       0       0
+FnuDII cgcg    4       2       1       2       2       0       0
+Fnu4HI GCNGC   5       2       0       2       3       0       0
+FokI   GGATG   5       2       0       14      18      0       0
+FriOI  GRGCYC  6       2       0       5       1       0       0
+FseI   GGCCGGCC        8       2       0       6       2       0       0
+FspI   TGCGCA  6       2       1       3       3       0       0
+FspAI  RTGCGCAY        8       2       1       4       4       0       0
+FspBI  CTAG    4       2       0       1       3       0       0
+Fsp4HI GCNGC   5       2       0       2       3       0       0
+GdiII  cggccr  6       2       0       1       5       0       0
+GlaI   GCGC    4       2       1       2       2       0       0
+GluI   GCNGC   5       2       0       2       3       0       0
+GsaI   CCCAGC  6       2       0       5       1       0       0
+GsuI   CTGGAG  6       2       0       22      20      0       0
+HaeI   wggccw  6       2       1       3       3       0       0
+HaeII  RGCGCY  6       2       0       5       1       0       0
+HaeIII GGCC    4       2       1       2       2       0       0
+HaeIV  gaynnnnnrtc     11      4       0       -8      -14     25      20
+HapII  CCGG    4       2       0       1       3       0       0
+HgaI   GACGC   5       2       0       10      15      0       0
+HgiAI  gwgcwc  6       2       0       5       1       0       0
+HgiCI  ggyrcc  6       2       0       1       5       0       0
+HgiEII accnnnnnnggt    12      0       0       0       0       0       0
+HgiJII grgcyc  6       2       0       5       1       0       0
+HhaI   GCGC    4       2       0       3       1       0       0
+Hin1I  GRCGYC  6       2       0       2       4       0       0
+Hin1II CATG    4       2       0       4       -1      0       0
+Hin4I  GAYNNNNNVTC     11      4       0       -9      -14     24      19
+Hin4II ccttc   5       2       0       11      10      0       0
+Hin6I  GCGC    4       2       0       1       3       0       0
+HinP1I GCGC    4       2       0       1       3       0       0
+HincII GTYRAC  6       2       1       3       3       0       0
+HindII GTYRAC  6       2       1       3       3       0       0
+HindIII        AAGCTT  6       2       0       1       5       0       0
+HinfI  GANTC   5       2       0       1       4       0       0
+HpaI   GTTAAC  6       2       1       3       3       0       0
+HpaII  CCGG    4       2       0       1       3       0       0
+HphI   GGTGA   5       2       0       13      12      0       0
+Hpy8I  GTNNAC  6       2       1       3       3       0       0
+Hpy99I CGWCG   5       2       0       5       -1      0       0
+Hpy166II       GTNNAC  6       2       1       3       3       0       0
+Hpy178III      tcnnga  6       2       0       2       4       0       0
+Hpy188I        TCNGA   5       2       0       3       2       0       0
+Hpy188III      TCNNGA  6       2       0       2       4       0       0
+HpyAV  CCTTC   5       2       0       11      10      0       0
+HpyCH4III      ACNGT   5       2       0       3       2       0       0
+HpyCH4IV       ACGT    4       2       0       1       3       0       0
+HpyCH4V        TGCA    4       2       1       2       2       0       0
+HpyF3I CTNAG   5       2       0       1       4       0       0
+HpyF10VI       GCNNNNNNNGC     11      2       0       7       4       0       0
+Hsp92I GRCGYC  6       2       0       2       4       0       0
+Hsp92II        CATG    4       2       0       4       -1      0       0
+HspAI  GCGC    4       2       0       1       3       0       0
+ItaI   GCNGC   5       2       0       2       3       0       0
+KasI   GGCGCC  6       2       0       1       5       0       0
+KflI   GGGWCCC 7       2       0       2       5       0       0
+KpnI   GGTACC  6       2       0       5       1       0       0
+Kpn2I  TCCGGA  6       2       0       1       5       0       0
+KspI   CCGCGG  6       2       0       4       2       0       0
+Ksp22I TGATCA  6       2       0       1       5       0       0
+Ksp632I        ctcttc  6       2       0       7       10      0       0
+KspAI  GTTAAC  6       2       1       3       3       0       0
+Kzo9I  GATC    4       2       0       -1      4       0       0
+LguI   GCTCTTC 7       2       0       8       11      0       0
+LpnI   rgcgcy  6       2       1       3       3       0       0
+Lsp1109I       GCAGC   5       2       0       13      17      0       0
+LweI   GCATC   5       2       0       10      14      0       0
+MabI   ACCWGGT 7       2       0       1       6       0       0
+MaeI   CTAG    4       2       0       1       3       0       0
+MaeII  ACGT    4       2       0       1       3       0       0
+MaeIII GTNAC   5       2       0       -1      5       0       0
+MalI   GATC    4       2       1       2       2       0       0
+MaqI   crttgac 7       2       0       28      26      0       0
+MauBI  CGCGCGCG        8       2       0       2       6       0       0
+MbiI   CCGCTC  6       2       1       3       3       0       0
+MboI   GATC    4       2       0       -1      4       0       0
+MboII  GAAGA   5       2       0       13      12      0       0
+McaTI  gcgcgc  6       2       0       4       2       0       0
+McrI   cgrycg  6       2       0       4       2       0       0
+MfeI   CAATTG  6       2       0       1       5       0       0
+MflI   RGATCY  6       2       0       1       5       0       0
+MhlI   GDGCHC  6       2       0       5       1       0       0
+MjaIV  gtnnac  6       0       0       0       0       0       0
+MlsI   TGGCCA  6       2       1       3       3       0       0
+MluI   ACGCGT  6       2       0       1       5       0       0
+MluNI  TGGCCA  6       2       1       3       3       0       0
+MlyI   GAGTC   5       2       1       10      10      0       0
+Mly113I        GGCGCC  6       2       0       2       4       0       0
+MmeI   TCCRAC  6       2       0       26      24      0       0
+MnlI   CCTC    4       2       0       11      10      0       0
+Mph1103I       ATGCAT  6       2       0       5       1       0       0
+MreI   CGCCGGCG        8       2       0       2       6       0       0
+MroI   TCCGGA  6       2       0       1       5       0       0
+MroNI  GCCGGC  6       2       0       1       5       0       0
+MroXI  GAANNNNTTC      10      2       1       5       5       0       0
+MscI   TGGCCA  6       2       1       3       3       0       0
+MseI   TTAA    4       2       0       1       3       0       0
+MslI   CAYNNNNRTG      10      2       1       5       5       0       0
+MspI   CCGG    4       2       0       1       3       0       0
+Msp20I TGGCCA  6       2       1       3       3       0       0
+MspA1I CMGCKG  6       2       1       3       3       0       0
+MspCI  CTTAAG  6       2       0       1       5       0       0
+MspR9I CCNGG   5       2       0       2       3       0       0
+MssI   GTTTAAAC        8       2       1       4       4       0       0
+MstI   tgcgca  6       2       1       3       3       0       0
+MunI   CAATTG  6       2       0       1       5       0       0
+MvaI   CCWGG   5       2       0       2       3       0       0
+Mva1269I       GAATGC  6       2       0       7       5       0       0
+MvnI   CGCG    4       2       1       2       2       0       0
+MvrI   CGATCG  6       2       0       4       2       0       0
+MwoI   GCNNNNNNNGC     11      2       0       7       4       0       0
+NaeI   GCCGGC  6       2       1       3       3       0       0
+NarI   GGCGCC  6       2       0       2       4       0       0
+NciI   CCSGG   5       2       0       2       3       0       0
+NcoI   CCATGG  6       2       0       1       5       0       0
+NdeI   CATATG  6       2       0       2       4       0       0
+NdeII  GATC    4       2       0       -1      4       0       0
+NgoAVIII       gacnnnnntga     11      4       0       -13     -15     24      22
+NgoMIV GCCGGC  6       2       0       1       5       0       0
+NhaXI  caagrag 7       0       0       0       0       0       0
+NheI   GCTAGC  6       2       0       1       5       0       0
+NlaIII CATG    4       2       0       4       -1      0       0
+NlaIV  GGNNCC  6       2       1       3       3       0       0
+NlaCI  catcac  6       2       0       25      23      0       0
+Nli3877I       cycgrg  6       2       0       5       1       0       0
+NmeAIII        GCCGAG  6       2       0       27      25      0       0
+NmeDI  rccggy  6       4       0       -13     -8      13      18
+NmuCI  GTSAC   5       2       0       -1      5       0       0
+NotI   GCGGCCGC        8       2       0       2       6       0       0
+NruI   TCGCGA  6       2       1       3       3       0       0
+NsbI   TGCGCA  6       2       1       3       3       0       0
+NsiI   ATGCAT  6       2       0       5       1       0       0
+NspI   RCATGY  6       2       0       5       1       0       0
+NspV   TTCGAA  6       2       0       2       4       0       0
+NspBII cmgckg  6       2       1       3       3       0       0
+OliI   CACNNNNGTG      10      2       1       5       5       0       0
+PabI   gtac    4       2       0       3       1       0       0
+PacI   TTAATTAA        8       2       0       5       3       0       0
+PaeI   GCATGC  6       2       0       5       1       0       0
+PaeR7I CTCGAG  6       2       0       1       5       0       0
+PagI   TCATGA  6       2       0       1       5       0       0
+PalAI  GGCGCGCC        8       2       0       2       6       0       0
+PasI   CCCWGGG 7       2       0       2       5       0       0
+PauI   GCGCGC  6       2       0       1       5       0       0
+PceI   AGGCCT  6       2       1       3       3       0       0
+PciI   ACATGT  6       2       0       1       5       0       0
+PciSI  GCTCTTC 7       2       0       8       11      0       0
+PctI   GAATGC  6       2       0       7       5       0       0
+PdiI   GCCGGC  6       2       1       3       3       0       0
+PdmI   GAANNNNTTC      10      2       1       5       5       0       0
+PfeI   GAWTC   5       2       0       1       4       0       0
+Pfl23II        CGTACG  6       2       0       1       5       0       0
+Pfl1108I       tcgtag  6       0       0       0       0       0       0
+PflFI  GACNNNGTC       9       2       0       4       5       0       0
+PflMI  CCANNNNNTGG     11      2       0       7       4       0       0
+PfoI   TCCNGGA 7       2       0       1       6       0       0
+PhoI   GGCC    4       2       1       2       2       0       0
+PinAI  ACCGGT  6       2       0       1       5       0       0
+PlaDI  catcag  6       2       0       27      25      0       0
+PleI   GAGTC   5       2       0       9       10      0       0
+Ple19I CGATCG  6       2       0       4       2       0       0
+PmaCI  CACGTG  6       2       1       3       3       0       0
+PmeI   GTTTAAAC        8       2       1       4       4       0       0
+PmlI   CACGTG  6       2       1       3       3       0       0
+PpiI   GAACNNNNNCTC    12      4       0       -8      -13     25      20
+PpsI   GAGTC   5       2       0       9       10      0       0
+Ppu10I atgcat  6       2       0       1       5       0       0
+Ppu21I YACGTR  6       2       1       3       3       0       0
+PpuMI  RGGWCCY 7       2       0       2       5       0       0
+PscI   ACATGT  6       2       0       1       5       0       0
+PshAI  GACNNNNGTC      10      2       1       5       5       0       0
+PshBI  ATTAAT  6       2       0       2       4       0       0
+PsiI   TTATAA  6       2       1       3       3       0       0
+Psp03I ggwcc   5       2       0       4       1       0       0
+Psp5II RGGWCCY 7       2       0       2       5       0       0
+Psp6I  CCWGG   5       2       0       -1      5       0       0
+Psp1406I       AACGTT  6       2       0       2       4       0       0
+Psp124BI       GAGCTC  6       2       0       5       1       0       0
+PspCI  CACGTG  6       2       1       3       3       0       0
+PspEI  GGTNACC 7       2       0       1       6       0       0
+PspGI  CCWGG   5       2       0       -1      5       0       0
+PspLI  CGTACG  6       2       0       1       5       0       0
+PspN4I GGNNCC  6       2       1       3       3       0       0
+PspOMI GGGCCC  6       2       0       1       5       0       0
+PspOMII        cgcccar 7       2       0       27      25      0       0
+PspPI  GGNCC   5       2       0       1       4       0       0
+PspPPI RGGWCCY 7       2       0       2       5       0       0
+PspPRI ccycag  6       2       0       21      19      0       0
+PspXI  VCTCGAGB        8       2       0       2       6       0       0
+PsrI   GAACNNNNNNTAC   13      4       0       -8      -13     25      20
+PssI   rggnccy 7       2       0       5       2       0       0
+PstI   CTGCAG  6       2       0       5       1       0       0
+PsuI   RGATCY  6       2       0       1       5       0       0
+PsyI   GACNNNGTC       9       2       0       4       5       0       0
+PteI   GCGCGC  6       2       0       1       5       0       0
+PvuI   CGATCG  6       2       0       4       2       0       0
+PvuII  CAGCTG  6       2       1       3       3       0       0
+RcaI   TCATGA  6       2       0       1       5       0       0
+RceI   catcgac 7       2       0       27      25      0       0
+RgaI   GCGATCGC        8       2       0       5       3       0       0
+RigI   GGCCGGCC        8       2       0       6       2       0       0
+RleAI  cccaca  6       2       0       18      15      0       0
+RpaB5I cgrggac 7       2       0       27      25      0       0
+RruI   TCGCGA  6       2       1       3       3       0       0
+RsaI   GTAC    4       2       1       2       2       0       0
+RsaNI  GTAC    4       2       0       1       3       0       0
+RseI   CAYNNNNRTG      10      2       1       5       5       0       0
+RsrII  CGGWCCG 7       2       0       2       5       0       0
+Rsr2I  CGGWCCG 7       2       0       2       5       0       0
+SacI   GAGCTC  6       2       0       5       1       0       0
+SacII  CCGCGG  6       2       0       4       2       0       0
+SalI   GTCGAC  6       2       0       1       5       0       0
+SanDI  GGGWCCC 7       2       0       2       5       0       0
+SapI   GCTCTTC 7       2       0       8       11      0       0
+SaqAI  TTAA    4       2       0       1       3       0       0
+SatI   GCNGC   5       2       0       2       3       0       0
+SauI   cctnagg 7       2       0       2       5       0       0
+Sau96I GGNCC   5       2       0       1       4       0       0
+Sau3AI GATC    4       2       0       -1      4       0       0
+SbfI   CCTGCAGG        8       2       0       6       2       0       0
+ScaI   AGTACT  6       2       1       3       3       0       0
+SchI   GAGTC   5       2       1       10      10      0       0
+SciI   ctcgag  6       2       1       3       3       0       0
+ScrFI  CCNGG   5       2       0       2       3       0       0
+SdaI   CCTGCAGG        8       2       0       6       2       0       0
+SdeAI  cagrag  6       2       0       27      25      0       0
+SdeOSI gacnnnnrtga     11      4       0       -12     -14     23      21
+SduI   GDGCHC  6       2       0       5       1       0       0
+SecI   ccnngg  6       2       0       1       5       0       0
+SelI   cgcg    4       2       0       -1      4       0       0
+SetI   ASST    4       2       0       4       -1      0       0
+SexAI  ACCWGGT 7       2       0       1       6       0       0
+SfaAI  GCGATCGC        8       2       0       5       3       0       0
+SfaNI  GCATC   5       2       0       10      14      0       0
+SfcI   CTRYAG  6       2       0       1       5       0       0
+SfeI   ctryag  6       2       0       1       5       0       0
+SfiI   GGCCNNNNNGGCC   13      2       0       8       5       0       0
+SfoI   GGCGCC  6       2       1       3       3       0       0
+Sfr274I        CTCGAG  6       2       0       1       5       0       0
+Sfr303I        CCGCGG  6       2       0       4       2       0       0
+SfuI   TTCGAA  6       2       0       2       4       0       0
+SgfI   GCGATCGC        8       2       0       5       3       0       0
+SgrAI  CRCCGGYG        8       2       0       2       6       0       0
+SgrBI  CCGCGG  6       2       0       4       2       0       0
+SgrDI  CGTCGACG        8       2       0       2       6       0       0
+SgsI   GGCGCGCC        8       2       0       2       6       0       0
+SimI   gggtc   5       2       0       2       5       0       0
+SinI   GGWCC   5       2       0       1       4       0       0
+SlaI   CTCGAG  6       2       0       1       5       0       0
+SmaI   CCCGGG  6       2       1       3       3       0       0
+SmiI   ATTTAAAT        8       2       1       4       4       0       0
+SmiMI  CAYNNNNRTG      10      2       1       5       5       0       0
+SmlI   CTYRAG  6       2       0       1       5       0       0
+SmoI   CTYRAG  6       2       0       1       5       0       0
+SmuI   CCCGC   5       2       0       9       11      0       0
+SnaI   gtatac  6       0       0       0       0       0       0
+SnaBI  TACGTA  6       2       1       3       3       0       0
+SpeI   ACTAGT  6       2       0       1       5       0       0
+SphI   GCATGC  6       2       0       5       1       0       0
+SplI   cgtacg  6       2       0       1       5       0       0
+SpoDI  gcggrag 7       0       0       0       0       0       0
+SrfI   GCCCGGGC        8       2       1       4       4       0       0
+Sse9I  AATT    4       2       0       -1      4       0       0
+Sse232I        cgccggcg        8       2       0       2       6       0       0
+Sse8387I       CCTGCAGG        8       2       0       6       2       0       0
+Sse8647I       aggwcct 7       2       0       2       5       0       0
+SseBI  AGGCCT  6       2       1       3       3       0       0
+SsiI   CCGC    4       2       0       1       3       0       0
+SspI   AATATT  6       2       1       3       3       0       0
+SspDI  GGCGCC  6       2       0       1       5       0       0
+SspD5I ggtga   5       2       1       13      13      0       0
+SstI   GAGCTC  6       2       0       5       1       0       0
+SstII  CCGCGG  6       2       0       4       2       0       0
+Sth132I        cccg    4       2       0       8       12      0       0
+Sth302II       ccgg    4       2       1       2       2       0       0
+StrI   CTCGAG  6       2       0       1       5       0       0
+StsI   ggatg   5       2       0       15      19      0       0
+StuI   AGGCCT  6       2       1       3       3       0       0
+StyI   CCWWGG  6       2       0       1       5       0       0
+StyD4I CCNGG   5       2       0       -1      5       0       0
+SwaI   ATTTAAAT        8       2       1       4       4       0       0
+TaaI   ACNGT   5       2       0       3       2       0       0
+TaiI   ACGT    4       2       0       4       -1      0       0
+TaqI   TCGA    4       2       0       1       3       0       0
+TaqII  GACCGA  6       2       0       17      15      0       0
+TaqII  CACCCA  6       2       0       17      15      0       0
+TasI   AATT    4       2       0       -1      4       0       0
+TatI   WGTACW  6       2       0       1       5       0       0
+TauI   GCSGC   5       2       0       4       1       0       0
+TfiI   GAWTC   5       2       0       1       4       0       0
+TliI   CTCGAG  6       2       0       1       5       0       0
+Tru1I  TTAA    4       2       0       1       3       0       0
+Tru9I  TTAA    4       2       0       1       3       0       0
+TscAI  CASTG   5       2       0       7       -3      0       0
+TseI   GCWGC   5       2       0       1       4       0       0
+TsoI   TARCCA  6       2       0       17      15      0       0
+Tsp45I GTSAC   5       2       0       -1      5       0       0
+Tsp509I        AATT    4       2       0       -1      4       0       0
+Tsp4CI acngt   5       2       0       3       2       0       0
+TspDTI ATGAA   5       2       0       16      14      0       0
+TspEI  AATT    4       2       0       -1      4       0       0
+TspGWI ACGGA   5       2       0       16      14      0       0
+TspMI  CCCGGG  6       2       0       1       5       0       0
+TspRI  CASTG   5       2       0       7       -3      0       0
+TssI   gagnnnctc       9       0       0       0       0       0       0
+TstI   CACNNNNNNTCC    12      4       0       -9      -14     24      19
+TsuI   gcgac   5       0       0       0       0       0       0
+Tth111I        GACNNNGTC       9       2       0       4       5       0       0
+Tth111II       caarca  6       2       0       17      15      0       0
+UbaF9I tacnnnnnrtgt    12      0       0       0       0       0       0
+UbaF11I        tcgta   5       0       0       0       0       0       0
+UbaF12I        ctacnnngtc      10      0       0       0       0       0       0
+UbaF13I        gagnnnnnnctgg   13      0       0       0       0       0       0
+UbaF14I        ccannnnntcg     11      0       0       0       0       0       0
+UbaPI  cgaacg  6       0       0       0       0       0       0
+UnbI   ggncc   5       2       0       -1      5       0       0
+Van91I CCANNNNNTGG     11      2       0       7       4       0       0
+Vha464I        CTTAAG  6       2       0       1       5       0       0
+VneI   GTGCAC  6       2       0       1       5       0       0
+VpaK11AI       ggwcc   5       2       0       -1      5       0       0
+VpaK11BI       GGWCC   5       2       0       1       4       0       0
+VspI   ATTAAT  6       2       0       2       4       0       0
+XagI   CCTNNNNNAGG     11      2       0       5       6       0       0
+XapI   RAATTY  6       2       0       1       5       0       0
+XbaI   TCTAGA  6       2       0       1       5       0       0
+XceI   RCATGY  6       2       0       5       1       0       0
+XcmI   CCANNNNNNNNNTGG 15      2       0       8       7       0       0
+XhoI   CTCGAG  6       2       0       1       5       0       0
+XhoII  RGATCY  6       2       0       1       5       0       0
+XmaI   CCCGGG  6       2       0       1       5       0       0
+XmaIII cggccg  6       2       0       1       5       0       0
+XmaCI  CCCGGG  6       2       0       1       5       0       0
+XmaJI  CCTAGG  6       2       0       1       5       0       0
+XmiI   GTMKAC  6       2       0       2       4       0       0
+XmnI   GAANNNNTTC      10      2       1       5       5       0       0
+XspI   CTAG    4       2       0       1       3       0       0
+ZraI   GACGTC  6       2       1       3       3       0       0
+ZrmI   AGTACT  6       2       1       3       3       0       0
+Zsp2I  ATGCAT  6       2       0       5       1       0       0