[Bioperl-l] Re: DNA topology
Heikki Lehvaslaiho
heikki@ebi.ac.uk
15 Nov 2002 15:58:48 +0000
Yes, of course. I fixed that myself at some point. My memory is going,
going.. Now, what were we talking about?
-Heikki
On Fri, 2002-11-15 at 15:33, Jason Stajich wrote:
> I think it's easier than that. Don't we have a function
>
> $seq->is_circular
>
> Which will get set to true if it is circular?
>
> >From Bio/SeqIO/genbank.pm:
> my $circ = shift(@tokens);
> if ($circ eq 'circular') {
> $params{'-is_circular'} = 1;
> $params{'-division'} = shift(@tokens);
> } else {
> # 'linear' or 'circular' may actually be omitted altogether
> $params{'-division'} =
> (CORE::length($circ) == 3 ) ? $circ : shift(@tokens);
> }
>
> -jason
> On 15 Nov 2002, Heikki Lehvaslaiho wrote:
>
> > 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
___ _/_/_/_/_/________________________________________________________