Skip to content
Snippets Groups Projects
Commit 317f7fe8 authored by Nathan Johnson's avatar Nathan Johnson
Browse files

corrected 0 start on slice

parent d620e795
No related branches found
No related tags found
No related merge requests found
......@@ -142,10 +142,21 @@ ok( $db , 'Test database instantiated');
my $sa = $db->get_SliceAdaptor;
my $ga = $db->get_GeneAdaptor;
my $slice = $sa->fetch_by_region('chromosome', '1', 0, 10000000 );
#Need to find Slice which genes are on
#Do this via direct mysql or hardcode?
#sql is unlikely to change
my $sql = 'SELECT sr.name, g.seq_region_start from gene g, seq_region sr where g.seq_region_id = sr.seq_region_id order by g.seq_region_start limit 1';
my ($sr_name, $sr_start) = @{$db->dbc->db_handle->selectrow_arrayref($sql)};
my $slice = $sa->fetch_by_region('chromosome', $sr_name, 1, ($sr_start + 1000000));
my @genes = @{$ga->fetch_all_by_Slice($slice)};
#Are there any genes in this test DB?
my $num_genes = scalar(@genes);
ok($num_genes, 'Failed to find genes on test Slice, please Slice redefine in test');
ok($num_genes, "Found $num_genes Gene(s) on test Slice");
SKIP:{
......@@ -175,11 +186,12 @@ SKIP:{
my $gene;
while($gene = $gi->next){
$gene_cnt++;
print $gene->stable_id."\n";
$gene_cnt++;
}
}
ok(($gene_cnt == $num_genes), 'fetch_Iterator_by_Slice returned correct number of features');
ok(($gene_cnt == $num_genes), "fetch_Iterator_by_Slice returned correct number of features($gene_cnt == $num_genes)");
}
}
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment