Newer
Older
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
use strict;
use lib 't';
use TestUtils qw(test_getter_setter debug);
BEGIN { $| = 1;
use Test;
plan tests => 8;
}
use MultiTestDB;
use Bio::EnsEMBL::AssemblyExceptionFeature;
our $verbose = 0;
my $multi = MultiTestDB->new;
# get a core DBAdaptor
#
my $dba = $multi->get_DBAdaptor("core");
my $aefa = $dba->get_AssemblyExceptionFeatureAdaptor();
ok($aefa);
#
# 1 create a new AssemblyExceptionFeature
#
my $aef = new Bio::EnsEMBL::AssemblyExceptionFeature;
ok($aef);
#
# test the basic getter and setters
#
# start
ok(test_getter_setter($aef,'start',10));
# end
ok(test_getter_setter($aef,'end',14));
# type
ok(test_getter_setter($aef,'type', 'HAP'));
# check adaptor attaching
$aef->adaptor($aefa);
ok($aef->adaptor->isa('Bio::EnsEMBL::DBSQL::AssemblyExceptionFeatureAdaptor'));
# fetch all
my $chr_slice = $dba->get_SliceAdaptor->fetch_by_region('chromosome',
'20_HAP1');
my @features = @{$aefa->fetch_all_by_Slice($chr_slice)};
ok(@features);
foreach my $f (@features) {
debug( "Feature: " . $f->slice->seq_region_name . " " .
$f->start . " " . $f->end . " " . $f->type);
my $as = $f->alternate_slice();
debug(" Alternate slice: " . $as->seq_region_name . " " .
$as->start . " " . $as->end);
}
my ($f) = @features;
ok($f->display_id eq $f->alternate_slice->seq_region_name);