[ensembl-dev] GERP Conservation Score via Compara API

Rashid Mamunur rm8 at sanger.ac.uk
Wed Aug 20 09:55:47 BST 2014


Hello everyone,
I have been trying to extract conservation score for single base pair using compare API and following code block. 
I was wondering if anyone confirm if MethodLinkSpeciesSet is the right way to extract conservation score. 

Is there any way to double check these scores via Biomart. ? 

Thanks,
Mamun
==================================================================================================================================

my $registry = "Bio::EnsEMBL::Registry";
$registry->load_registry_from_db(
		-host => "ensembldb.ensembl.org",
		-user => "anonymous",
		-port => 5306 );

my $mlss_adaptor = $registry->get_adaptor("Multi", "compara", "MethodLinkSpeciesSet");
my $mlss = $mlss_adaptor->fetch_by_method_link_type_species_set_name("GERP_CONSERVATION_SCORE", "mammals");

#create slice
my $slice_adaptor = $registry->get_adaptor($species, 'core', 'Slice');
my $slice = $slice_adaptor->fetch_by_region('toplevel', $seq_region, $seq_region_start, $seq_region_end);

#get conservation score adaptor
my $cs_adaptor = $registry->get_adaptor("Multi", 'compara', 'ConservationScore');			
my $display_size = $slice->end - $slice->start + 1; 

## Get score array 
my $scores = $cs_adaptor->fetch_all_by_MethodLinkSpeciesSet_Slice($method_linked_species_Set, $slice, $display_size);

==================================================================================================================================
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20140820/a5b3c0e8/attachment.html>


More information about the Dev mailing list