RemoteRestMap

Módulo para buscar sitios de restricción como cliente del servidor RestrictionMapper, muy sencillo de usar.

#!/usr/bin/perl -w
# program that gets the E.coli restriction maps of a series of E.coli sequences in fasta format
# based on test.pl from http://www.restrictionmapper.org/remote.htm


use lib "/home/compu2/Perl/remoterestmap/";
use strict;
use RemoteRestMap;

my %settings = (
	isoschizomers => "yes",
	enzymelist    => ['EciI','M.Eco67Dam','EcoHAI','EcoNI','Eco29kI','EcoprrI']
	);
	
# instantiation of map object
my $map_obj = new RemoteRestMap("http://www.restrictionmapper.org/cgi-local/sitefind3.pl", "atgc"); 

$map_obj->sequence("ATCGATCGATACGCGATACTAGCGCGATCATCGCGCGCTATATACTCAGCTCGCTCTAGATATCTA");
	
my $map = $map_obj->get_map(\%settings);
my @sites;

while (my $cuts = $map->get_next_enz()) 
{ 
	push( @sites , split(/\,/,$cuts->{CUTLIST}->[0]) );	
}

Bruno Contreras-Moreira
http://www.eead.csic.es/compbio