Simple Python wrapper for NCBI Blast.
pip install git+https://github.com/MrTomRod/ncbi_blast.git
# load class
from ncbi_blast.Blast import Blast
blast = Blast(blast_path='/local/bin', outfmt=6)
queryfile = open('tests/data/query_single_nucl.fasta').read()
nucl_fasta1 = 'tests/data/ffn/prot_nucl_1.ffn'
nucl_fasta2 = 'tests/data/ffn/prot_nucl_2.ffn'
# create blast database for nucl_fastas
blast.mkblastdb(file=nucl_fasta1, dbtype='nucl')
blast.mkblastdb(file=nucl_fasta2, dbtype='nucl')
# blast database
output = blast.blastn(fasta_string=queryfile, db=nucl_fasta2)
# example output (outfmt=6): 'QUERY_ID\MATCH_ID\t100.000\t1677\t0\t0\t1\t1677\t1\t1677\t0.0\t3097\n'
# blast multiple databases
output = blast.blastn(fasta_string=queryfile, db=[nucl_fasta1, nucl_fasta2])