写一个简单的python程序完成病毒的基因获取并绘制饼图

好的，以下是一个简单的 Python 程序，用于获取病毒基因序列并绘制饼图：

# 导入所需的库
from Bio import Entrez
import matplotlib.pyplot as plt

# 设置邮箱地址
Entrez.email = "your.email@example.com"

# 搜索病毒基因序列
handle = Entrez.esearch(db="nucleotide", term="HIV-1[Organism] AND pol gene", retmax=100)
record = Entrez.read(handle)

# 获取基因长度
gene_ids = record["IdList"]
gene_len = []
for gene_id in gene_ids:
    handle = Entrez.efetch(db="nucleotide", id=gene_id, rettype="fasta", retmode="text")
    gene_seq = handle.read().strip().split("\n")[1]
    gene_len.append(len(gene_seq))

# 绘制饼图
labels = ['<5000', '5000-10000', '10000-15000', '>15000']
sizes = [0, 0, 0, 0]
for length in gene_len:
    if length < 5000:
        sizes[0] += 1
    elif length < 10000:
        sizes[1] += 1
    elif length < 15000:
        sizes[2] += 1
    else:
        sizes[3] += 1

plt.pie(sizes, labels=labels, autopct='%1.1f%%')
plt.title("HIV-1 pol 基因长度分布")
plt.show()

该程序使用 Biopython 库获取病毒 HIV-1 的 pol 基因序列，首先连接 PubMed 数据库，然后搜索 HIV-1 的 pol 基因序列，并获取序列数据。接着，程序使用 Matplotlib 库绘制饼图，展示 HIV-1 pol 基因序列长度的分布情况。你可以根据自己的需求修改病毒名称和基因名称，以及绘制的图形类型和样式。