test.fasta内容如下:
>heavy
EVQLVESGGGLVQPGGSLRLSCAASGFTLSGDWIHWVRQAPGKGLEWLGEISAAGGYTDYADSVKGRFTISADTSKNTAYLQMNSLRAEDTAVYYCARESRVSFEAAMDYWGQGT
>light
DIQMTQSPSSLSASVGDRVTITCRASQDLATDVAWYQQKPGKAPKLLIYSASFLYSGVPSRFSGSGSGTDFTLTISSLQP
编写代码
from Bio import SeqIO
#fasta文件解析
for seq_record in SeqIO.parse("test.fasta","fasta"):
print(seq_record.id)
print(repr(seq_record.seq))
print(len(seq_record))
##GenBank文件解析,唯一的区别是改变上述代码的文件名和格式字符串
输出结果
heavy
Seq('EVQLVESGGGLVQPGGSLRLSCAASGFTLSGDWIHWVRQAPGKGLEWLGEISAA...QGT')
115
light
Seq('DIQMTQSPSSLSASVGDRVTITCRASQDLATDVAWYQQKPGKAPKLLIYSASFL...LQP')
80