使用python实现两个fasta文件的交集

from Bio import SeqIO

def read_fasta(file_path):
    sequences = set()
    for record in SeqIO.parse(file_path, "fasta"):
        sequences.add(str(record.seq))
    return sequences

def get_intersection(file1, file2):
    sequences1 = read_fasta(file1)
    sequences2 = read_fasta(file2)
    intersection = sequences1.intersection(sequences2)
    return intersection

file1 = "file1.fasta"
file2 = "file2.fasta"
intersection = get_intersection(file1, file2)
for sequence in intersection:
    print(sequence)