質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

88.06%

Biopython GenBank形式のファイルを読込みその配列を翻訳後、他のファイル翻訳配列と同じか比べたい。

受付中

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 52

score 0

GenBank形式のファイルを読込みその配列を翻訳後、他のファイル翻訳配列と同じか比べたい。

ここに質問の内容を詳しく書いてください。
(例)PHP(CakePHP)で●●なシステムを作っています。
Biopythonで、GenBank形式ファイルをseqIOで読み込んで、翻訳しましたが、その後のステップがわかりません。
しかし、print実行で、翻訳が見られる状態で、実際の配列翻訳したデータだけはどこに保存されているのでしょうか? Seq (MKFRP....)などの状態で表記されます。でも、他のファイルと比べるときは、Seqはいらなくて、MKFRP...の部分だけで比べたいです。また、比べるときは、MKFRPという部分は文字列にして、リスト化する必要があるのでしょうか?その場合は、numpyとか使うのですか?SeqIOだけで翻訳し、MKFRP部分を他の配列比べられるのでしょうか?初心者なのでどうしてよいかわかりません。結局は、GenBank形式から呼びだして、翻訳し、それを他のGenBank形式の配列と同じかどうか調べたい。(その際、文字列化やリスト形式にするのかなと他でも良い方法があれば。)
■■な機能を実装中に以下のエラーメッセージが発生しました。

発生している問題・エラーメッセージ

GenBankファイルをSeqIOで読み込み後の処理がわからない。

該当のソースコード

python, biopython

試したこと

nuc_translate = nuc.seq[351:666].translate()と入力したが、Seq(MKFRP...)と出力される。prinにしたら、このSeqは消えるが。MKFRPを文字列としてどう設定するか、リスト化して、他の配列と同じか判定するやり方がわからない。初心者なのでnumpyやpandasとか使う方法しかわからない。

!pip install biopython
from Bio import SeqIO

nuc = SeqIO.read("sequence1.gb", "genbank")
nuc

print(nuc)
print(nuc.seq)
len(nuc.seq)

print(nuc.seq[351:666].translate())
nuc_translate = nuc.seq[351.666].transalte()

補足情報(FW/ツールのバージョンなど)

JupyterLab,  biopython ver 19.2.3

ここにより詳細な情報を記載してください。

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

0

Biopythonについては全く知らないので一般論です。

nuc_translateからその部分的な情報を取り出したいとのことですね。
部分的な情報は属性になっていることが多いです。
属性名を見てそれらしいものを見てみると見つかる場合が多いですね。

print(dir(nuc_translate))


で、メソッド名を含む属性名の一覧が出てきます。
仮に属性名がabcなら、

print(nuc_translate.abc)


で属性の値をみることができます。アンダーバーで始まる属性名は無関係の場合が多いでしょう。
nuc_translate.abcでメソッドだというような表示が出てきたときは、

print(nuc_translate.abc.__doc__)


を行うとメソッドの使い方が出てきます。
これでだめな場合は、英語の公式ドキュメントを見るしかないでしょう。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2021/05/05 10:38

    print(nuc_translate.encode)と入力したら、<bound method Seq.encode of Seq('MKFRPLHDRVVVKRIDAEEKTAGGIIIPDTVKEKPSQGEVIAVGPGGRDESGKL...AA*')>
    と出ました。これが使えるのですか?この後、どうやったら、リスト化でできますか?

    キャンセル

  • 2021/05/05 13:34

    nuc_list = list(range(0:10)
    としても、invalid syntax としかでません。

    キャンセル

  • 2021/05/05 13:55

    nuc_list = list("MKFRPLHDRVVVKRIDAEEKTAGGIIIPDTVKEKPSQGEVIAVGPGGRDESGKL...")とコピペしたらできますけど、コピペしない方法教えてください!

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 88.06%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

同じタグがついた質問を見る

  • トップ
  • Pythonに関する質問
  • Biopython GenBank形式のファイルを読込みその配列を翻訳後、他のファイル翻訳配列と同じか比べたい。