perl/etc/justonce/sperm/fixkmhist.pl

   1 #!/bin/env perl
   2 use strict;
   3 use warnings;
   4
   5 die "Usage: $0 <input>\n" if @ARGV < 1;
   6 my ($inf)=@ARGV;
   7
   8 open IN,'<',$inf or die "Error opening $inf: $!\n";
   9 open OUTF,'>',$inf.'.fixed' or die "Error opening $inf.fixed: $!\n";
  10
  11 my ($HistSum,$Check1,$Check2)=(0,0,0);
  12 while (<IN>) {
  13         if (/^#/) {
  14                 print OUTF $_;
  15                 if (/\bCntOfFoundKmer\t(\d+)\t/) {
  16                         $HistSum = $1;
  17                 }
  18         } else {
  19                 chomp;
  20                 my @DAT = split /\t/;
  21                 push @DAT,$DAT[2];
  22                 $DAT[2] = $DAT[1]/$HistSum;
  23                 $Check1 += $DAT[2];
  24                 $Check2 += $DAT[4];
  25                 $DAT[3] = $Check1;
  26                 print OUTF join("\t",@DAT),"\n";
  27         }
  28 }
  29 die if $HistSum == 0;
  30 print "$HistSum,$Check1,$Check2\n";
  31 close IN;
  32 close OUTF;
  33
  34 __END__
  35 find *.hist | xargs -n1 ./fixkmhist.pl