perl/etc/WoodyMiaoLin/PbeBefore2015/1.make_input.pl

   1 #!/usr/bin/perl
   2 use strict;
   3 use warnings;
   4
   5 open I, "-|", "gzip -dc ../4.G-PhoCS/input.gz";
   6 open O, "|-", "gzip -9c >GPhoCS_PbePvi_input.gz";
   7
   8 while (<I>) {
   9         if (/^chr/) {
  10                 my @a = split /\t/;
  11                 print O "$a[0]\t2\t$a[2]";
  12         } elsif (/^FCA|PVI/) {
  13                 next;
  14         } else {
  15                 print O $_;
  16         }
  17 }
  18 close I;
  19 close O;