perl/etc/WoodyMiaoLin/PbeBefore2015/missing_sample.pl

   1 #!/usr/bin/perl
   2 use strict;
   3 use warnings;
   4
   5 open I, "<", "CALB-CHRNA-CMA-GNB.fasta";
   6 open O1, ">", "CALB-CHRNA-CMA-GNB_missMore100.fasta";
   7 open O2, ">", "CALB-CHRNA-CMA-GNB_missLess100.fasta";
   8
   9 while (<I>) {
  10         my $seq = <I>;
  11         my $miss = @{[$seq =~ /-/g]};
  12         if ($miss > 100) {
  13                 print O1 $_, $seq;
  14         } else {
  15                 print O2 $_, $seq;
  16         }
  17 }
  18 close I;
  19 close O1;
  20 close O2;