main.py

def filter_by_min_positives(blast_results, minPositivesPercent):
    correct_indices = []
    for index, sequence in enumerate(blast_results.alignments):
        if sequence.hsps[0].positives * 100.0 / sequence.length >= minPositivesPercent:
            correct_indices.append(index)

    filtered_results = Bio.Blast.Record.Blast()
    filtered_results.multiple_alignment = blast_results.multiple_alignment

    for index in correct_indices:
        filtered_results.descriptions.append(blast_results.descriptions[index])
        filtered_results.alignments.append(blast_results.alignments[index])

    return filtered_results