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

View of /FigKernelScripts/filter_on_known.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (download) (as text) (annotate)
Sun May 16 00:37:03 2010 UTC (9 years, 6 months ago) by overbeek
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
filter pegs with no raw data

use strict;
use Data::Dumper;

my $usage = "usage: filter_on_known RawData < clusters > filtered";
my($raw_dataF);
(
 ($raw_dataF  = shift @ARGV) && (-s $raw_dataF)
)
    || die $usage;

my %known_pegs = map { ($_ =~ /^(fig\|\d+\.\d+\.[^\.]+\.\d+)/) ? ($1 => 1) : () } `cut -f1 $raw_dataF`;

while (defined($_ = <STDIN>))
{
    my @pegs;
    if (($_ =~ /^(\S+)/) && (@pegs = grep { $known_pegs{$_} } split(/,/,$1)) && (@pegs > 1))
    {
	print join(",",@pegs),"\n";
    }
}


MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3