Render blast3

From BioPerl
Jump to: navigation, search
#!/usr/bin/perl
 
# file: render_blast3.pl
# This is code example 3 in the Graphics-HOWTO
# Author: Lincoln Stein
 
use strict;
use lib "$ENV{HOME}/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  => 1,
			      -bgcolor => 'blue',
			      -min_score => 0,
			      -max_score => 1000,
			      -font2color     => 'red',
			      -sort_order     => 'high_score',
			      -description => sub {
				my $feature = shift;
				my $score   = $feature->score;
				return "score=$score";
			       });
 
while (<>) { # read blast file
  chomp;
  next if /^\#/;  # ignore comments
  my($name,$score,$start,$end) = split /\t+/;
  my $feature = Bio::SeqFeature::Generic->new(-score=>$score,
					      -display_name=>$name,
					      -start=>$start,-end=>$end);
  $track->add_feature($feature);
}
 
print $panel->png;
Personal tools
Namespaces
Variants
Actions
Main Links
documentation
community
development
Toolbox