igblast.py 文件源码

def makeblastdb(fasta, database_name):
    with SequenceReader(fasta) as fr:
        sequences = list(fr)
    if not sequences:
        raise ValueError("FASTA file {} is empty".format(fasta))
    process_output = subprocess.check_output(
        ['makeblastdb', '-parse_seqids', '-dbtype', 'nucl', '-in', fasta, '-out', database_name],
        stderr=subprocess.STDOUT
    )
    if b'Error: ' in process_output:
        raise subprocess.SubprocessError()