[Bio] / FigKernelScripts / svr_determine_sets_of_related_contigs.pl Repository:
ViewVC logotype

View of /FigKernelScripts/svr_determine_sets_of_related_contigs.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.4 - (download) (as text) (annotate)
Tue May 18 19:38:33 2010 UTC (9 years, 5 months ago) by parrello
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.3: +8 -0 lines
Documentation fixes.

#!/usr/bin/perl

#
#	This is a SAS Component.
#

=head1 svr_determine_sets_of_related_contigs

    svr_determine_sets_of_related_contigs username password <contig_ids.tbl >output.tbl

This takes as input a list of contig IDs. What is the output?

=cut

use Data::Dumper;
use RASTserver;
use strict;

my $username = shift;
my $password = shift;

my $rast = new RASTserver($username, $password);

my @input_ids = <STDIN>;
chomp @input_ids;

my %seen;
for my $id (@input_ids)
{
    next if $seen{$id};

    my $res = $rast->get_contig_ids_in_project_from_entrez({ -contig_id => $id } );
    # print Dumper($res);
    my $project_ids = $res->{ids};
    my $redundancies = $res->{redundancy_report};

#    my($project_ids, $redundancies) = $rast->get_contig_ids_in_project_from_entrez({ -contig_id => $id } );
    if (@$redundancies)
    {
	for my $redundancy (@$redundancies)
	{
	    print STDERR join("\t", @$redundancy), "\n";
	}
    }
    else
    {
	print join(",", @$project_ids), "\n";
	map { $seen{$_} = 1 } @$project_ids;
    }
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3