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