[Bioperl-l] Problems with Bioperl graphics
Wacki
wackattack at gmail.com
Tue Jul 5 03:23:04 EDT 2005
I followed the tutorial here:
http://bioperl.org/HOWTOs/Graphics-HOWTO/gettingStarted.html
And ran this exact code:
http://biokdd.informatics.indiana.edu/jnowacki/render_blast1.txt
The image produced is shown here:
http://biokdd.informatics.indiana.edu/jnowacki/test.png
It doesn't have the name of the hits. What is wrong? The code is
exactly the same as the tutorial is it not?
code:
#!/usr/bin/perl
# This is code example 2 in the Graphics-HOWTO
use strict;
use lib '/home/lstein/projects/bioperl-live';
use Bio::Graphics;
use Bio::SeqFeature::Generic;
my $panel = Bio::Graphics::Panel->new(-length => 1000,
-width => 800,
-pad_left => 10,
-pad_right => 10,
);
my $full_length = Bio::SeqFeature::Generic->new(-start=>1,-end=>1000);
$panel->add_track($full_length,
-glyph => 'arrow',
-tick => 2,
-fgcolor => 'black',
-double => 1,
);
my $track = $panel->add_track(-glyph => 'graded_segments',
-label => 4,
-bgcolor => 'blue',
-min_score => 0,
-max_score => 1000);
while (<>) { # read blast file
chomp;
next if /^\#/; # ignore comments
my($name,$score,$start,$end) = split /\t+/;
my $feature = Bio::SeqFeature::Generic->new(-display_name=>$name,-score=>$score,
-start=>$start,-end=>$end);
$track->add_feature($feature);
}
binmode(STDOUT);
print $panel->png;
More information about the Bioperl-l
mailing list