perl 利用管道读取压缩文件内容

perl的文件句柄不仅支持普通文件， 还支持管道，今天需要统计一个fastq文件中的序列数和碱基数，而NGS的fastq文件一般都是gzip压缩的，所以

需要读取压缩文件中的内容，代码如下：

my ($fastq) = @ARGV:
my ($reads, $bases) = cal_sequence_info($fastq);
print qq{$reads\t$bases\n};

sub cal_sequence_info {
my $fastq = shift;
my $file_handle =   $fastq =~ /gz$/ ? qq{$fastq} : qq{gzip -dc $fastq |};
open FASTQ, $file_handle or die "Can't open $fastq\n";
my ($reads, $bases);
while (my $readid = <FASTQ>) {
my $quality   = <FASTQ>;
my $comment   = <FASTQ>;
my $sequence  = <FASTQ>;
chomp($sequence);
$reads++;
$bases += length $sequence;
}
close FASTQ;
return ($reads, $bases);
}

直接利用linux中的管道， 这样可以方便的读取压缩文件中的内容