[ensembl-dev] is_circular returns 0 for mitochondrial DNA

[Rainer Johannes]

Dear all,

I was puzzled that is_circular on the mitochondrial DNA returns 0. Below is the code I used - eventually there is something wrong with it. I'm using the Perl API version 94 and perl 5, version 18, subversion 0 (v5.18.0).


my $registry = 'Bio::EnsEMBL::Registry';
$registry->load_registry_from_db(-host => $host, -user => $user,
 -pass => $pass, -port => $port);
my $slice_adaptor = $registry->get_adaptor("human", "core", "slice");

my $chr_slice = $slice_adaptor->fetch_by_region("chromosome", "MT");
print $chr_slice->name.": circular? ".$chr_slice->is_circular()."\n";


The output:

chromosome:GRCh38:MT:1:16569:1: circular? 0


thanks for any hints,

jo