您好,登錄后才能下訂單哦!
這篇文章主要介紹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進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。