import subprocess as sb
import os

species = ['BRUSU', 'ECOLI', 'NEIMA', 'PARDP', 'POLAQ', 'SACD2', 'SERP5', 'THIDA']

for elem in species:
    dirname = '{}_proteome'.format(elem)
    proteomename = '{}.fasta'.format(elem)
    dbname = '{}_proteome_db'.format(elem)
    logname = '{}_db.log'.format(elem)
    blastoutname = '{}_blastout.tsv'.format(elem)
    replacename = '../{}'.format(blastoutname)
    os.chdir(dirname)
    sb.call(['makeblastdb',
        '-in', proteomename,
        '-dbtype', 'prot',
        '-out', dbname,
        '-logfile', logname])
    sb.call(['blastp',
        '-task', 'blastp',
        '-query', '../CLPX_ECOLI.fasta',
        '-db', dbname,
        '-evalue', '0.001',
        '-word_size', '2',
        '-outfmt', '7',
        '-out', blastoutname])
    os.replace(blastoutname, replacename)
    os.chdir('..')