在搜索兩個(gè)正則表達(dá)式的 DNA 序列后，我試圖創(chuàng)建一個(gè) .bed 文件。理想情況下，我想生成一個(gè)制表符分隔的文件，其中包含序列描述、第一個(gè)正則表達(dá)式的起始位置和第二個(gè)正則表達(dá)式的結(jié)束位置。我知道正則表達(dá)式部分有效，它只是創(chuàng)建了我正在努力處理的 \t 分隔文件。我希望我可以打開(kāi)/創(chuàng)建一個(gè)文件并簡(jiǎn)單地為for loop包含此信息的每次迭代打印一個(gè)新行，如下所示：with open("Mimp_hits.bed", "a+") as file_object:    for line in file_object:        print(f'{sequence.description}\t{h.start()}\t{h_rc.end()}')file_object.close()但這似乎不起作用（創(chuàng)建空文件）。我也嘗試過(guò)使用file_object.write，但這同樣會(huì)創(chuàng)建一個(gè)空文件。這是我所有的代碼，包括搜索正則表達(dá)式：import re, sysfrom Bio import SeqIOfrom Bio.SeqRecord import SeqRecordinfile = sys.argv[1]for sequence in SeqIO.parse(infile, "fasta"):    hit = re.finditer(r"CAGTGGG..GCAA[TA]AA", str(sequence.seq))    mimp_length = 400    for h in hit:        h_start = h.start()        hit_rc = re.finditer(r"TT[TA]TTGC..CCCACTG", str(sequence.seq))        for h_rc in hit_rc:            h_rc_end = h_rc.end()            length = h_rc_end - h_start            if length > 0:                if length < mimp_length:                    with open("Mimp_hits.bed", "a+") as file_object:                        for line in file_object:                            print(sequence.description, h.start(), h_rc.end())                    file_object.close()這是所需的輸出：Focub_II5_mimp_1__contig_1.16(656599:656809)    2   208Focub_II5_mimp_2__contig_1.47(41315:41540)  2   223Focub_II5_mimp_3__contig_1.65(13656:13882)  2   224Focub_II5_mimp_4__contig_1.70(61591:61809)  2   216這是示例輸入：>Focub_II5_mimp_1__contig_1.16(656599:656809)TACAGTGGGATGCAAAAAGTATTCGCAGGTGTGTAGAGAGATTTGTTGCTCGGAAGCTAGTTAGGTGTAGCTTGTCAGGTTCTCAGTACCCTATATTACACCGAGATCAGCGGGATAATCTAGTCTCGAGTACATAAGCTAAGTTAAGCTACTAACTAGCGCAGCTGACACAACTTACACACCTGCAAATACTTTTTGCATCCCACTGTA>Focub_II5_mimp_2__contig_1.47(41315:41540)TACAGTGGGAGGCAATAAGTATGAATACCGGGCGTGTATTGTTTTCTGCCGCTAGCCCATTTTAACAGCTAGAGTGTGTATATTAACCTCACACATAGCTATCTCTTATACTAATTGGTTAGGGAAAACCTCTAACCAGGATTAGGAGTCAACATAGCTTGTTTTAGGCTAAGAGGTGTGTGTCAGTACACCAAAGGGTATTCATACTTATTGCCCCCCACTGTA有人能幫忙嗎？