dict_genes=[]
with open("findings.txt", mode="r") as file:
    for line in file:
        dict_genes.append(line.strip())
gene=None
with open("genes.fasta", mode="w") as fasta:
    with open("all_proteomes.fa", mode="r") as file:
        for line in file:
            line=line.strip()
            if line.startswith(">"):
                gene=None
                for i in dict_genes:
                    if line.startswith(i):
                        gene=i
                        print(gene, file=fasta)
                        break
            elif gene:
                print(line, file=fasta)




