0

我有读取 VCF 文件的功能

def readVCF(VCF):
    VCF = VCF.readlines()
    header = [] 
    data = []
    for line in VCF:
        if line.startswith('#'):
            header.append(line) 
        else:
            data.append(line)
    VCF_contents =[]  
    for line in data:
        VCF_contents.extend(line.strip().split('\t'))
    for i in range(0, len(VCF_contents), 10): 
        lines.append(VCF_contents[i : i+10])

该功能可以按我的意愿工作，但是当我多次使用它时就会出现问题。

我想在不同的文件上多次使用这个函数，并同时使用所有这些不同文件的输出。我一直在使用这样的功能：

VCF = open("frogs.txt", 'r')
lines = []
readVCF(VCF)
frogs = lines
VCF.close

VCF = open("rabbits.txt.", 'r')
readVCF(VCF)
rabbits = lines
VCF.close

VCF = open("lizards.txt", 'r')
readVCF(VCF)
lizards = lines
VCF.close

问题是，我的所有输出都是相同的，并且都有来自 的输出，即frogs.txt我使用该函数的第一个文件。如何在不获取相同值的情况下重用此函数？

4

1 回答 1

1

可能，您要添加：

lines = []

一开始在你的函数内部，并且：

return lines

在函数的末尾，并像这样调用它：

VCF = open("frogs.txt", 'r')
frogs = readVCF(VCF)
VCF.close()

或者（可能更好）：

with open("frogs.txt", 'r') as VCF:
    frogs = readVCF(VCF)
于 2019-10-23T21:21:10.210 回答