1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81
|
# -*-Perl-*- mode (to keep my emacs happy)
# $Id: RestrictionIO.t,v 1.6 2003/12/15 03:45:41 redwards Exp $
# test for Bio::Restriction::Analysis.pm
# written by Rob Edwards
use strict;
my $NUMTESTS;
my $error;
BEGIN {
eval { require Test; };
if( $@ ) {
use lib 't','..';
}
use Test;
$NUMTESTS = 14;
$error = 0;
plan tests => $NUMTESTS;
}
if( $error == 1 ) {
exit(0);
}
require Bio::Restriction::IO;
use Bio::Root::IO;
require Bio::Restriction::EnzymeCollection;
use Data::Dumper;
ok(1);
#
# default enz set
#
ok my $in = Bio::Restriction::IO->new();
ok my $renzs = $in->read;
ok $renzs->each_enzyme, 532;
ok my $e = $renzs->get_enzyme('AccI');
ok $e->name, 'AccI';
ok my $out = Bio::Restriction::IO->new(-format=>'base', -file=> ">/tmp/r");
#$out->write($renzs);
#map {print $_->name, "\t", $_->site, "\t", $_->overhang, "\n"} $renzs->each_enzyme;
#
# withrefm, 31
#
ok $in = Bio::Restriction::IO->new
(-format=> 31, -verbose => 0,
-file => Bio::Root::IO->catfile("t","data","rebase.withrefm"));
ok $renzs = $in->read;
ok $renzs->each_enzyme, 11;
#
# itype2, 8
#
#enzyme name [tab] prototype [tab] recognition sequence with cleavage site
# [tab] methylation site and type [tab] commercial source [tab] references
ok $in = Bio::Restriction::IO->new
(-format=> 8, -verbose => 0,
-file => Bio::Root::IO->catfile("t","data","rebase.itype2"));
ok $renzs = $in->read;
ok $renzs->each_enzyme, 16;
ok $out = Bio::Restriction::IO->new(-format=>'base');
#$out->write($renzs);
|