pip install biopython
conda install biopython
Quick Start
from Bio.Seq import Seq
from Bio.SeqRecord import SeqRecord
from Bio import SeqIO
my_seq = Seq("AGTACACTGGT")
# 输入输出
SeqIO.read("ORF.faa","fasta") #sigal
SeqIO.parse("ls_orchid.fasta", "fasta") #multi
SeqIO.write("w_ORF.faa","fasta")
Tips
from Bio import SeqIO
from Bio.SeqRecord import SeqRecord
def extra_match(genome_sequence, match_sequence):
'''提取序列'''
print(match_sequence)
uniq_in_191_store = []
for i in SeqIO.parse(genome_sequence,'fasta'):
for j in match_sequence:
if i.name.upper() == str(j).upper():
print(i.name)
uniq_in_191_store.append(SeqRecord(id=i.id,\
seq=i.seq ,\
description=i.description))
SeqIO.write(uniq_in_191_store,'result.ffn','fasta')
return None
本文介绍Biopython库的安装与使用,演示如何利用Biopython进行DNA序列的操作,包括序列读取、解析及匹配提取等核心功能。通过具体代码示例,帮助读者快速上手生物信息学的编程实践。
1652

被折叠的 条评论
为什么被折叠?



