我是python编程的新手，并且有一个fasta文件，我想解析它以便在特定软件中使用。该文件包含两行：1）在分类学的序列标识符和用空格隔开的分类法，和最后种名也可以包含空格，和2）的DNA序列（参见下面的实施例）： >123876987 Bacteria;test;test;test test test
ATCTGCTGCATGCATGCATCGACTGCATGAC
>239847239 Ba
我需要为以下术语解析FASTA头文件：叶，芽，茎和嫩芽，如果序列包含任何一个术语，则打开一个文件并将其放在那里使用Biopython。 所以我让他们转换为使用SeqIO.to_dict字典： from Bio import SeqIO
records_dict = SeqIO.to_dict(SeqIO.parse("my_example.fasta","fasta"))
但现在我不知道如何
我尝试使用头文件中的增量ID从.csv文件创建多个fasta。但是我的脚本运行不正常。任何建议？ M=open('/home/anna/Scrivania/db_mat/matk_2db_c.csv','r')
M_out=open('/home/anna/Scrivania/db_mat/db_matk_bronx.fas', 'w')
for i in range[1,92]:
我试图从Ensembl FASTA文件中找到蛋白质图案。我已经完成了大部分脚本，比如检索序列ID和序列本身，但是我收到了一些有趣的结果。 #!/usr/bin/perl
use strict;
use warnings;
use autodie;
my $motif1 = qr/(HE(\D)(\D)H(\D{18})E)/x;
my $motif2 = qr/(AMEN)/x;
m