cgi-bin/maps/physical/clone_restriction_gel_image.pl

   1 use strict;
   2 use warnings;
   3
   4 use Bio::Graphics::Gel;
   5
   6 use CXGN::Genomic::Clone;
   7
   8 use CXGN::Page;
   9
  10 my $page = CXGN::Page->new;
  11
  12 my ($clone_id,$enzyme,$fingerprint_id) = $page->get_encoded_arguments(qw/id enzyme fp_id/);
  13 $fingerprint_id += 0; $clone_id += 0; #< enforce numericness
  14
  15 my $clone = CXGN::Genomic::Clone->retrieve($clone_id)
  16   or die "could not retrieve clone with id '$clone_id'";
  17
  18 my $is_frags = $clone->in_silico_restriction_fragment_sizes($enzyme);
  19
  20 my ($iv_frags) = grep {$_->[0] == $fingerprint_id && $_->[1] eq $enzyme } $clone->in_vitro_restriction_fragment_sizes($enzyme);
  21
  22 if($iv_frags) {
  23   shift @$iv_frags; shift @$iv_frags; #shift off the fingerprint id and enzyme name
  24 }
  25
  26 $is_frags || $iv_frags
  27   or die "no restriction fragments available at all for this clone.  user should not have been directed here.";
  28
  29 #warn "drawing fragment\n";
  30 my $gel = Bio::Graphics::Gel->new(
  31                                   $iv_frags ? ('i.v.' => $iv_frags) : (),
  32                                   $is_frags ? ('i.s.' => $is_frags) : (),
  33                                    -min_frag => 1000,
  34                                   -lane_length => 200,
  35                                   -bandcolor => [0xff,0xc6,0x00],
  36                                  );
  37
  38 print "Content-Type: image/png\n\n";
  39 print $gel->img->png;
  40