[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