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

View of /FigKernelScripts/make_md5_figfams_file.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Fri Nov 14 02:55:03 2008 UTC (11 years, 3 months ago) by arodri7
Branch: MAIN
CVS Tags: mgrast_dev_08112011, rast_rel_2009_05_18, mgrast_dev_08022011, rast_rel_2014_0912, myrast_rel40, mgrast_dev_05262011, rast_rel_2008_12_18, mgrast_dev_04082011, rast_rel_2010_0928, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, rast_rel_2009_0925, rast_rel_2010_0526, rast_rel_2014_0729, mgrast_dev_02212011, rast_rel_2010_1206, mgrast_release_3_0, mgrast_dev_03252011, rast_rel_2010_0118, rast_rel_2009_02_05, 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_2009_07_09, rast_rel_2010_0827, myrast_33, rast_rel_2011_0928, mgrast_dev_04052011, mgrast_dev_02222011, rast_rel_2009_03_26, mgrast_dev_10262011, HEAD
script makes the md5 to figfams file

use FIG;
my $fig = new FIG;

my $usage = "usage: make_md5_figfams_file figfams_dir > md5.figfams";

($dir = shift @ARGV) || die $usage;

my $md5Hash = {};
open (FH, "$dir/families.2c") || die "could not open families.2c file\n";

while (my ($ff, $id) = <FH> =~ m/(.*)\t(.*)\n/){
    $prevFF = $ff;
    my $md5 = $fig->md5_of_peg($id);
    push (@{$md5Hash->{$md5}}, $ff);
#    push (@{$idHash->{$md5}}, $id);
}
close FH;

foreach my $md5 (sort keys %{$md5Hash}){
    my %saw;
    @saw{@{$md5Hash->{$md5}}} = ();
    my @array = sort keys %saw;  # remove sort if undesired

    print join(",", @array) . "\t$md5\n";
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3