[Bio] / FigWebServices / show_pegs.cgi Repository:
ViewVC logotype

View of /FigWebServices/show_pegs.cgi

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (annotate)
Wed Jun 2 19:45:53 2010 UTC (10 years, 1 month ago) by olson
Branch: MAIN
display pegs in jbrowse

use FIG;
use CGI;
use strict;
use Data::Dumper;

my $fig = new FIG;

my $cgi = new CGI;

my @pegs = $cgi->param("pegs");

my @expt = $cgi->param("experiment");

my @locs = $fig->feature_location_bulk(\@pegs);

my $genome = $fig->genome_of($pegs[0]);

my($contig, $min, $max) = $fig->boundaries_of($locs[0]->[1]);
($min, $max) = ($max, $min) if $min > $max;

for my $loc (map { $_->[1] } @locs[1..$#locs])
{
    my ($contig, $beg, $end) = $fig->boundaries_of($loc);
    ($beg, $end) = ($end, $beg) if $beg > $end;
    $min = $beg if $beg < $min;
    $max = $end if $end > $max;
}

$min -= 100;
$min = 0 if $min < 0;
$max += 100;

my $expts = join(",", map { "probes_$_" } @expt);
my $url = "http://bioseed.mcs.anl.gov/ross/FIG-Tmp/jbrowse/$genome/?loc=$genome:$contig:$min..$max&tracks=DNA,features,$expts";
print "$url\n";

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3