import Bio.SwissProt as SProt
from Bio import File
parser = E.OptionParser( version = "%prog version: $Id: get_sequences_from_www.py 2782 2009-09-10 11:40:29Z andreas $")
if __name__ == "__main__":
parser.add_option("-f", "--field", dest="fields", type="string",
help="field to write to table." , action="append" )
parser.set_defaults(
fields = [] )
(options, args) = E.Start( parser )
s_parser = SProt.RecordParser()
ninput, nfound, nmissed = 0, 0, 0
for line in sys.stdin:
if line[0] == "#": continue
id = line[:-1].split("\t")[0]
ninput += 1
try:
result = ExPASy.get_sprot_raw(id).read()
except IOError:
options.stdlog.write( "# warning: sequence for id %s not found." % id )
nmissed += 1
continue
