[Bioperl-l] Sub Seq Feature help
michael watson (IAH-C)
michael.watson at bbsrc.ac.uk
Fri Feb 6 05:32:05 EST 2004
Hello
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
Thanks
Mick
More information about the Bioperl-l
mailing list