from Bio import Entrez
def get_species_from_ncbi(sequence_id):
Entrez.email = "XX@email.com" # 设置您的电子邮件地址
handle = Entrez.efetch(db="nucleotide", id=sequence_id, retmode="xml")
records = Entrez.read(handle)
handle.close()
species = records[0]['GBSeq_organism']
return species
# 在此处输入您的ID列表文件名
input_file = "XXXid_list.txt"
with open(input_file, 'r') as file:
id_list = file.readlines()
species_list = []
for sequence_id in id_list:
sequence_id = sequence_id.strip()
species_name = get_species_from_ncbi(sequence_id)
species_list.append(species_name)
for species_name in species_list:
print(species_name)
根据NCBI序列ID对应到物种名称
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
相关阅读更多精彩内容
- 最近成为了硕导,克隆基因的基本工作就是获取对应的基因组序列,启动子,内含子等信息。为了降低学生学习获取核酸序列的难...
- 通常我们会使用比对好的fasta文件构建进化树,fasta文件中大于号后的内容就是最终进化树上的文字标签。如果拿到...
- 第一版本:修改名字文件,不对序列处理比如序列名字如下,仅仅需要物种名字 利用python脚本加正则表达处理;思路:...
- 使用场景 假设有一个 fasta 格式的序列文件SRR9620252.faa,我们想要提取其中的一些序列到一个新的...