这篇文章主要介绍bioperl如何读写fasta/fastq压缩文件,文中介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们一定要看完!
bioperl读写fasta/fastq压缩文件
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');
以上是“bioperl如何读写fasta/fastq压缩文件”这篇文章的所有内容,感谢各位的阅读!希望分享的内容对大家有帮助,更多相关知识,欢迎关注亿速云行业资讯频道!
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。