[Bioperl-l] Bio::Seq issue

Paolo Pavan paolo.pavan at gmail.com
Thu May 31 11:34:11 EDT 2012


Hello everyone again,
By the way, I think I have encountered a minor issue in the method
Bio::Seq->is_circular(), that is defined in the file Bio/Seq.pm as a pure
getter as:

sub is_circular { shift->primary_seq->is_circular }

while it's counterpart in Bio/PrimarySeq.pm is defined as a getter/setter
as:

sub is_circular{
    my $self = shift;
    return $self->{'is_circular'} = shift if @_;
    return $self->{'is_circular'};
}

the result is that if you have a Bio::Seq object, for instance read by
Bio::SeqIO you can't change any more the property (well, unless you do
$seq->primary_seq->is_circular($is_circular) ).

Do someone agree that they should have the same behaviour?
It the case, attached patch file applied to Bio/Seq.pm file should do the
job.

Best regards,
Paolo
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Seq.pm.patch
Type: application/octet-stream
Size: 239 bytes
Desc: not available
URL: <http://lists.open-bio.org/pipermail/bioperl-l/attachments/20120531/a7e39139/attachment.obj>


More information about the Bioperl-l mailing list