[Bio] / FigMetagenomeTools / correctphageparsing.pl Repository:
ViewVC logotype

View of /FigMetagenomeTools/correctphageparsing.pl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1.1.1 - (download) (as text) (annotate) (vendor branch)
Mon Feb 19 17:15:26 2007 UTC (12 years, 9 months ago) by olson
Branch: x, MAIN
CVS Tags: mgrast_dev_08112011, mgrast_dev_08022011, mgrast_dev_05262011, mgrast_dev_04082011, mgrast_version_3_2, mgrast_dev_12152011, mgrast_dev_06072011, mgrast_rel_2008_0806, mgrast_dev_10262011, mgrast_dev_02212011, mgrast_rel_2008_0923, mgrast_release_3_0, mgrast_dev_03252011, mgrast_rel_2008_0924, mgrast_rel_2008_1110_v2, mgrast_rel_2008_0625, 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, mgrast_rel_2008_0919, mgrast_rel_2008_1110, myrast_33, mgrast_rel_2008_0917, mgrast_dev_04052011, mgrast_dev_02222011, y, HEAD
Changes since 1.1: +0 -0 lines
Initial import

#!/usr/bin/perl -w

use strict;
my $genomefile=shift || die "phage database file";
my %genome;
open(IN, $genomefile) || die "CAn't open $genomefile";
while (<IN>) {
 my @a=split /\t/;
 $genome{$a[1]}=$a[0];
}

die "other file(s) to correct?" unless ($ARGV[0]);

foreach my $f (@ARGV) {
 open(IN, $f) || die "Can't opnbe $f";
 open(OUT, ">/home2/rob/temp.txt") || die "Can't open /home2/rob/temp.txt";
 while (<IN>)
 {
  my @a=split /\t/;
  if ($genome{$a[6]}) {$a[7]=$genome{$a[6]}}
  else {print STDERR "No genome for $a[6]\n"}
  print OUT join("\t", @a);
 }
 
 `mv $f $f.old`;
 `cp -i /home2/rob/temp.txt $f`;
}

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3