根据NCBI序列ID对应到物种名称

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)
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容