#!/usr/bin/perl
# Generate an enzyme list for RestrictionEnzyme.pm from rebase
# From Ryan Brinkman
my $strider = $ARGV[0]; #commercial_version_rebase_strider_format
open (FILEIN,"$strider") or die "can't open $strider: $!\n";
while (<FILEIN>){
chomp;
if ( /^[A-Z]\S+,\S+/ ){
($enzyme,$cutsite)=split(',');
if ($cutsite =~ m#/#){
$match=$-[0];
}
($seqfixed=$cutsite) =~ s/\///g;
$seqfixed=uc $seqfixed;
print " \'$enzyme\'\t=> \'".$seqfixed." ".$match."\'\,\n";
}
}