Im following this guide: http://www.bioperl.org/wiki/HOWTO:Graphics and what Im tryin to do is create a png file for each hit in the blastoutput file. Also each png file should be named after the hit names. So forexample if there is a hit called 200134, then the resulting png file would be 200134.png etc
sub parseBlastOutput{ use Bio::Graphics; use Bio::SearchIO; use Bio::SeqFeature::Generic; my $searchio = Bio::SearchIO->new( -file => 'CAZyOutputM7.xml', -format => 'blastxml' ) or die +'parse failed'; #my $result = $searchio->next_result() or die "no result"; while (my $result = $searchio->next_result) { my $panel = Bio::Graphics::Panel->new( -length => $result-> +query_length, -width => 800, -pad_left => 10, -pad_right => 10 ); my $full_length = Bio::SeqFeature::Generic->new( -start + => 1, -end + => $result->query_length, -display_name + => $result->query_name ); $panel->add_track( $full_length, -glyph => 'arrow', -tick => 2, -fgcolor => 'black', -double => 1, -label => 1 ); my $track = $panel->add_track( -glyph => 'graded_segment +s', -label => 1, -connector => 'dashed', -bgcolor => 'blue', -font2color => 'red', -sort_order => 'high_score', -description => sub { my $featu +re = shift; return un +less $feature->has_tag('description'); my ($desc +ription) = $feature->each_tag_value('description'); my $score + = $feature->score; "$descrip +tion, score=$score"; } ); while (my $hit = $result->next_hit) { my $hitname = $hit->name; my $feature = Bio::SeqFeature::Generic->new( -score + => $hit->raw_score, -display_name + => $hit->name, -tag + => { description => $hit->description } ); while (my $hsp = $hit->next_hsp) { $feature->add_sub_SeqFeature($hsp,'EXPAND'); } $track->add_feature($feature); open my $blast, '>', $hitname or die "Failed to create $hi +tname\n"; print $blast $panel->png; close $blast or die "Failed to close $blast\n"; } } }
The way i have the script setup now, it doesnt do anything. I would appreciate some help on this.
In reply to BioPerl Graphics problem by Tjuh
| For: | Use: | ||
| & | & | ||
| < | < | ||
| > | > | ||
| [ | [ | ||
| ] | ] |