bioperl读写gz压缩格式的fasta/fastq文件方法如下：
输入文件句柄：
fastq文件
open my $FQ ,"zcat infile.fq.gz|" or die "$!";
my$fq=Bio::SeqIO->new(-fh=>$FQ,-format=>'fastq');
fasta文件
open my $FA ,"zcat infile.fa.gz|" or die "$!";
my$fa=Bio::SeqIO->new(-fh=>$FA,-format=>'fasta');
输出文件句柄：
fastq文件
open my $GZ ,"| gzip >outfile.fq.gz" or die $!;
my$out = Bio::SeqIO->new(-fh => $GZ , -format => 'fastq');
fasta文件
open my $GZ ,"| gzip >outfile.fa.gz" or die $!;
my$out = Bio::SeqIO->new(-fh => $GZ , -format => 'fasta');
更多perl语言知识可观看 Perl语言高级编程 学习！
