• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

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

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

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中的管道, 这样可以方便的读取压缩文件中的内容


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
perl中的循环控制发布时间:2022-07-22
下一篇:
perl读取json格式的文件发布时间:2022-07-22
热门推荐
热门话题
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap