批量计算测序文件reads长度

输入文件（包含测序材料号的文件，测序材料后缀需要是_1.fq.ft.gz，_2.fq.ft.gz ）

SRR13661968
SRR13661969
SRR13661970
SRR13661971
SRR13661972
SRR13661973
SRR13661974

使用的perl脚本如下：

#!/usrr/bin/perl -w
use strict;

my $in0 = $ARGV[0]; ##- samples 

my $out = "read_length.lst";

open OUT, ">$out";
open IN0, $in0;
while(<IN0>){
  chomp;
  my $samid = $_;

  my $left  = $samid."_1.fq.ft.gz";
  my $right = $samid."_2.fq.ft.gz";
 
  my $leftLen = &countLen($left);
  my $rightLen = &countLen($right);
     print OUT join("\t", $samid, $leftLen, $rightLen), "\n";
}
close IN0;
close OUT;


sub countLen {
 
    my ($fileIn) = @_;
 
    my $count = 0;
    my $length = 0;
 
    open IN1, "gzip -dc $fileIn | ";
    while(<IN1>){
      chomp;
      $count += 1;
      if($count == 2){
         $length = length($_);    
         last;
      }
    }
    close IN1;
    return($length);

}

输出结果(材料对应的左端和右端reads长度)：

SRR13661968     100     100
SRR13661969     100     100
SRR13661970     100     100
SRR13661971     100     100
SRR13661972     100     100
SRR13661973     100     100