始める無料で始める

6フレーム内で保存領域を見つける

ここまで到達した皆さん、おめでとうございます!

ここからは、ジカウイルス配列の中から NS5 タンパク質配列を探していきます。NS5 は非常に保存性の高いウイルスタンパク質です。今回は Uniprot から取得し、すでに読み込んであります。

ジカウイルスの DNA 配列は RNAStringSet に転写され、rnaframesZikaSet という名前で用意されています。セットには6つのリーディングフレーム(各配列1フレーム)が含まれており、これらをアミノ酸に翻訳します。検索時には、vcountPattern() の呼び出しで max.mismatch 引数を設定して、検索にある程度の柔軟性を持たせます。

この演習はコースの一部です

Rで始めるBioconductor入門

コースを見る

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Print rnaframesZikaSet 
___

# Translate rnaframesZikaSet 
AAzika6F <- ___(___)
AAzika6F
コードを編集して実行