def bowtiedb(fa, keepDB):
btdir = '<STR_LIT>' % (os.getcwd())<EOL>if not os.path.exists(btdir):<EOL><INDENT>os.mkdir(btdir)<EOL><DEDENT>btdb = '<STR_LIT>' % (btdir, fa.rsplit('<STR_LIT:/>', <NUM_LIT:1>)[-<NUM_LIT:1>])<EOL>if keepDB is True:<EOL><INDENT>if os.path.exists('<STR_LIT>' % (btdb)):<EOL><INDENT>return btdb<EOL><DEDENT><DEDENT>p = subprocess.Popen('<STR_LIT>'% (fa, btdb), shell = True)<EOL>p.communicate()<EOL>return btdb<EOL>
make bowtie db