[Bioperl-l] Sub Seq Feature help

michael watson (IAH-C) michael.watson at bbsrc.ac.uk
Fri Feb 6 05:32:05 EST 2004


I want to manipulate the start and end position of a CDS feature that looks like this:

FT	CDS	join(2307..3221,1..1623)

I have tried:

my @features = $seq->get_all_SeqFeatures;
foreach $f (@features) {
	my @subs = $f->sub_SeqFeature;
	foreach $s (@subs) {
		print $s->start, "-", $s->end, "\n";

However, I get nothing out.  The code doesn't descend into the sub seq features as $f->sub_SeqFeature doesn't return anything.  Nor does $f->get_SeqFeatures.

Clearly I am doing something wrong, but what?  I am using Bioperl-1.2.3


More information about the Bioperl-l mailing list