def phmmer(query, db, type, out, threads = '<STR_LIT:4>', evalue = '<STR_LIT>'):
if os.path.exists(out) is False:<EOL><INDENT>print('<STR_LIT>' % (query, db))<EOL>os.system('<STR_LIT>' % (out, out, evalue, threads, query, db))<EOL><DEDENT>else:<EOL><INDENT>print('<STR_LIT>' % (query, db))<EOL><DEDENT>phmmer2blast('<STR_LIT>' % out, out)<EOL>
run phmmer