[Bioperl-l] wormbase115/Bio::DB::GFF::Aggregator
wormbase_transcript problem
Philip MacMenamin
pm66 at nyu.edu
Tue Dec 23 11:01:13 EST 2003
Hi,
Previously I ran the following code to draw a curated gene with UTRs
hanging on the ends (attached):
my $aggregator = Bio::DB::GFF::Aggregator->new(-method => 'transcript',
-sub_parts => ['UTR:UTR','exon:curated','CDS:curated']
);
my $db = new Bio::DB::GFF(-adaptor=>'dbi::mysqlopt',
-dsn=>'dbi:mysql:wormbase115Mod;host=localhost',
-user=>'philip',
-pass=> $passwd,
-aggregator =>$aggregator
) or die();
my @all_transcripts = $searchSeg->features('transcript');
if (scalar @all_transcripts )
{
$panel->add_track(wormbase_transcript=>\@all_transcripts,
-bgcolor => 'wheat',
-fgcolor => 'black',
-forwardcolor => 'blue',
-reversecolor => 'blue',
-spacing => 0,
-utr_color => '#D0D0D0',
-font2color => 'blue',
-height => 10,
-description => 1,
-label => 1,
-key => "Curated genes");
}
This worked fine.
With the latest release of wormbase this code does not work.
So I have changed the -sub_parts arg for the Aggregator object to
-sub_parts => ['UTR:UTR','exon:curated','CDS:curated']
which gets rid of some things that I dont want. However, I cannot get it to
draw the UTRs actually hanging on the ends of the gene, it stacks/bumps them
now. They do not over lap in their start /stop co-ords.
???
Thanks.
--
Philip MacMenamin
More information about the Bioperl-l
mailing list