生物信息学报错：'global name 'start' is not defined'及TMscanner.py开发需求

下面是完全符合需求的实现，包含两个自定义函数，同时加入了边界处理和友好的结果输出逻辑：

# TMscanner.py
def getProteinRegion(protein_seq, start_idx):
    """
    从蛋白质序列中获取10个氨基酸的滑动窗口
    参数:
        protein_seq: str，完整的蛋白质氨基酸序列（单字母缩写）
        start_idx: int，窗口起始索引（0-based）
    返回:
        str，长度为10的氨基酸窗口；如果剩余长度不足10，返回None
    """
    if start_idx + 10 > len(protein_seq):
        return None
    return protein_seq[start_idx:start_idx+10]

def testForTM(window_seq):
    """
    计算窗口中非极性氨基酸的小数占比
    参数:
        window_seq: str，10个氨基酸的窗口序列
    返回:
        float，非极性氨基酸占比（保留4位小数）
    """
    # 用集合存储非极性氨基酸，查询效率更高
    non_polar_aa = {'A', 'V', 'L', 'I', 'P', 'M', 'F', 'W'}
    count = 0
    for aa in window_seq:
        if aa.upper() in non_polar_aa:  # 兼容小写输入
            count +=1
    return round(count / 10, 4)

def main():
    # 支持用户输入序列，也可替换为从文件读取
    protein_seq = input("请输入蛋白质氨基酸序列（单字母缩写）：").strip()
    if not protein_seq:
        print("错误：输入序列不能为空")
        return
    
    print("\n跨膜结构域检测结果：")
    print("起始位置(1-based)\t窗口序列\t非极性氨基酸占比")
    print("-" * 60)
    
    start_idx = 0
    while True:
        window = getProteinRegion(protein_seq, start_idx)
        if not window:
            break
        tm_ratio = testForTM(window)
        # 转换为1-based起始位置输出，符合生物信息学常用标注习惯
        print(f"{start_idx+1}\t\t\t{window}\t{tm_ratio}")
        # 可根据需求调整阈值，这里以0.6为例标记潜在跨膜区域
        if tm_ratio > 0.6:
            print("\t\t\t⚠️  该窗口为潜在跨膜结构域")
        start_idx +=1

if __name__ == "__main__":
    main()

代码说明： #

• getProteinRegion：做了边界检查，避免索引越界；用0-based索引适配Python字符串切片逻辑。
• testForTM：用集合存储非极性氨基酸提升查询效率，兼容小写输入，返回值保留4位小数更直观。
• 主函数main：支持用户交互输入，输出1-based起始位置（符合生物信息学领域习惯），并添加了阈值判断标记潜在跨膜区域。

使用示例： #

运行程序后输入序列AVLIPMFWXXAVLIPMFW，会输出每个窗口的占比，其中前两个窗口的占比会达到0.8，被标记为潜在跨膜区域。

内容的提问来源于stack exchange，提问作者Seva Galitskiy