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
s_iterator = SProt.Iterator(File.StringHandle(result), s_parser)
try:
cur_record = s_iterator.next()
except SyntaxError:
print "# Sequence not found: %s" % id
continue
columns = [ id, ]
for f in options.fields:
if f == "description":
columns.append( cur_record.description )
elif f == "references":
for ref in cur_record.references:
columns.append( ref.authors, ref.title )
elif f == "organism_classification":
