- 博客(1)
- 收藏
- 关注
原创 利用Bio.SeqIO提取基因组中特定序列时注意事项
例如: from Bio import SeqIO # 读取基因组文件 fa_seq = SeqIO.read(fasta_t_path+'\\genome.fasta', "fasta") # 从基因组中提取我想要的基因(基因的位置信息从NCBI上找到为4130142:4131557) # 但不能用起始位置的4130142作为起始,因为seq的从基因组中提取序列的方式类似列表的切片,4130142该位置的碱基是不被包含进去的,所以用4130141 b = fa_seq.seq[4130141:
2022-05-16 11:12:50
543
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人