biopythonでswissprotのレコードをパースする

python

意味無くSwissprotの
ftp://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.dat.gz
からデータをダウンロード，解凍しパースしましょう．

from Bio import SwissProt

records = SwissProt.parse(open("uniprot_sprot.dat"))

for rec in records:
    if "subtilis" in rec.organism:
        print rec.sequence

Class Recordのメンバーはここで確認できます．
http://www.biopython.org/DIST/docs/api/Bio.SwissProt.Record-class.html
本当は目的を持ってパースしましょう．
さよなら．