perl/etc/justonce/fqtrim50.pl

   1 #!/bin/env perl
   2 use strict;
   3 use warnings;
   4
   5 unless (@ARGV){
   6         print "perl $0 <in>\n";
   7         exit;
   8 }
   9
  10 my ($fq1) = @ARGV;
  11 open FQ,'-|',"gzip -dc $fq1";
  12
  13 while (<FQ>) {
  14         my $a=$_;
  15         my $b=<FQ>;
  16         my $c=<FQ>;
  17         my $d=<FQ>;
  18         $b=substr $b,0,50;
  19         $d=substr $d,0,50;
  20         print "$a$b\n$c$d\n";
  21 }