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

View of /FigWebServices/show_pegs.cgi

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.2 - (download) (annotate)
Wed Jun 2 19:52:14 2010 UTC (9 years, 5 months ago) by olson
Branch: MAIN
CVS Tags: mgrast_dev_08112011, mgrast_dev_08022011, rast_rel_2014_0912, myrast_rel40, mgrast_dev_05262011, mgrast_dev_04082011, rast_rel_2010_0928, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, rast_rel_2014_0729, mgrast_dev_02212011, rast_rel_2010_1206, mgrast_release_3_0, mgrast_dev_03252011, rast_rel_2011_0119, mgrast_release_3_0_4, mgrast_release_3_0_2, mgrast_release_3_0_3, mgrast_release_3_0_1, mgrast_dev_03312011, mgrast_release_3_1_2, mgrast_release_3_1_1, mgrast_release_3_1_0, mgrast_dev_04132011, mgrast_dev_04012011, rast_rel_2010_0827, myrast_33, rast_rel_2011_0928, mgrast_dev_04052011, mgrast_dev_02222011, mgrast_dev_10262011, HEAD
Changes since 1.1: +2 -2 lines
fixe

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

my $fig = new FIG;

my $cgi = new CGI;

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

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 $cgi->redirect($url);

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3