[Bioperl-l] Re: DNA topology

[Heikki Lehvaslaiho]

Jing,

RefSeq loosely follows GenBank format. In GenBank, the LOCUS line can
contain the word 'circular'. That word should be part of the string
returned by the Bio::Seq::RichSeq::division method. I am not sure if
there are any circular molecules in RefSeq.

Try:

if ($seq->division =~ /circular/ ) {
   
} else { #linear

}


Yours,
	-Heikki


On Fri, 2002-11-15 at 13:56, jzhao wrote:
> Mr. Lehvaslaiho,
> 
> I am wrapping some BioPerl modules to retrieve genome features from the 
> RefSeq database. Is there any handy methods that can tell whether the DNA 
> is linear or circular?
> 
> Thanks,
> Jing
>    
-- 
______ _/      _/_____________________________________________________
      _/      _/                      http://www.ebi.ac.uk/mutations/
     _/  _/  _/  Heikki Lehvaslaiho          heikki@ebi.ac.uk
    _/_/_/_/_/  EMBL Outstation, European Bioinformatics Institute
   _/  _/  _/  Wellcome Trust Genome Campus, Hinxton
  _/  _/  _/  Cambs. CB10 1SD, United Kingdom
     _/      Phone: +44 (0)1223 494 644   FAX: +44 (0)1223 494 468
___ _/_/_/_/_/________________________________________________________