91超碰碰碰碰久久久久久综合_超碰av人澡人澡人澡人澡人掠_国产黄大片在线观看画质优化_txt小说免费全本

溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

bioperl如何直接讀取輸出gz壓縮格式的fasta序列

發布時間:2022-02-24 11:46:16 來源:億速云 閱讀:160 作者:小新 欄目:開發技術

小編給大家分享一下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序列”這篇文章的所有內容,感謝各位的閱讀!相信大家都有了一定的了解,希望分享的內容對大家有所幫助,如果還想學習更多知識,歡迎關注億速云行業資訊頻道!

向AI問一下細節

免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。

AI

都江堰市| 珠海市| 开化县| 黑龙江省| 清水河县| 甘谷县| 德钦县| 云霄县| 双城市| 廉江市| 临安市| 普定县| 高密市| 包头市| 清新县| 马边| 萝北县| 丹巴县| 沽源县| 海口市| 宜良县| 定边县| 织金县| 南木林县| 清水河县| 阳泉市| 逊克县| 天镇县| 新乐市| 孟村| 绿春县| 分宜县| 集安市| 通河县| 嵩明县| 吉水县| 东乌| 安阳市| 正镶白旗| 惠水县| 明星|