-2
如何打印第一次出現的蛋白質序列?對於這個查詢,我得到了四個結果,我只想要第一個結果。如何打印每個迭代的搜索序列的第一個結果?
use Bio::DB::GenBank;
use Bio::DB::Query::GenBank;
$query = "LEGK";
$query_obj = Bio::DB::Query::GenBank->new(
-db => 'protein',
-query => $query
);
$gb_obj = Bio::DB::GenBank->new;
$stream_obj = $gb_obj->get_Stream_by_query($query_obj);
while ($seq_obj = $stream_obj->next_seq) {
# do something with the sequence object
print
">$query", ' ',
$seq_obj->display_id, ' ',
$seq_obj->desc, "\n",
$seq_obj->seq[, '\n';
這while
循環應該是這樣的
while ($seq_obj = $stream_obj->next_seq) {
# do something with the sequence object
print $seq_obj->display_id, "\t", $seq_obj->length, "\n";
}
這不會編譯。請發佈足夠的代碼,以便它可以在不修復語法錯誤的情況下運行,除非這些是問題的一部分。這爲任何試圖幫助你的人創造了額外的工作!另外,請在文件的開始部分使用'strict strict;'和'use warnings''。 – bytepusher
嗯,我已經發布了一個答案。如果沒有用,請嘗試更詳細地解釋您要查找的內容。這是一個perl論壇,我不是一個遺傳學家;) – bytepusher
@MTG:請不要編輯你的問題,以便對任何已發佈的答案或評論做出廢話。我已恢復您不正確的Perl代碼並將您的新代碼添加爲更新。 – Borodin