您好,登錄后才能下訂單哦!
小編給大家分享一下bioperl如何直接讀取輸出gz壓縮格式的fasta序列,相信大部分人都還不怎么了解,因此分享這篇文章給大家參考一下,希望大家閱讀完這篇文章后大有收獲,下面讓我們一起去了解一下吧!
bioperl 直接讀取輸出gz壓縮格式的fasta序列,代碼示例:
die "perl $0 <id><fa><OUT>" unless(@ARGV==3); use Math::BigFloat; use Bio::SeqIO; use Bio::Seq; use Data::Dumper; use PerlIO::gzip; open $FI, "<:gzip", "$ARGV[1]" or die "$!"; open $FO, ">:gzip", "$ARGV[2]" or die "$!"; $in = Bio::SeqIO->new(-fh => $FI , -format => 'Fasta'); $out = Bio::SeqIO->new(-fh => $FO , -format => 'Fasta');
另一種寫法也是可以的:
use Bio::SeqIO; use Bio::Seq; use Data::Dumper; use PerlIO::gzip; open FI, "<:gzip", "$ARGV[1]" or die "$!"; open FO, ">:gzip", "$ARGV[2]" or die "$!"; $in = Bio::SeqIO->new(-fh => \*FI , -format => 'Fasta'); $out = Bio::SeqIO->new(-fh => \*FO , -format => 'Fasta');
以上是“bioperl如何直接讀取輸出gz壓縮格式的fasta序列”這篇文章的所有內容,感謝各位的閱讀!相信大家都有了一定的了解,希望分享的內容對大家有所幫助,如果還想學習更多知識,歡迎關注億速云行業資訊頻道!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。