Question

我有一个作业，其中我必须编写一个python类来表示和操纵生物序列。我几乎完成了该类，但是我需要从输入该类的FASTA文件中导入序列，以便它可以运行它们。为了在命令行中运行该类，我需要创建一个run_me.py文件，当在命令行中执行该文件时，应先请求FASTA文件，然后运行该类。

这是类的构造函数：

class Bioseq:
    """Biological sequence class"""
    def __init__(self, seq):
        self.seq = seq

这是我对run_me.py文件的了解：

def read_fasta():
        with open(file_name, "r") as fh:
            line = fh.readline()
            header = ""
            seq = ""
            while line:
                line = line.rstrip('\n')
                if ">" in line:
                    header = line
                else:
                    seq = seq + line
                line = fh.readline()
        return seq
        print(seq)


from TPC2_MR2 import Bioseq
file_name = input("Please enter sequence file name:")
Bioseq = Bioseq(read_fasta())
Bioseq.main()

TPC2_MR2是类模块的名称。现在的问题是，当我执行run_me.py文件时，要求我提供FASTA文件的名称，但是之后，该序列未导入到类中。我该怎么办？

请注意，我不能使用外部Pyhton模块，例如Biopython。

如何读取FASTA文件并将序列插入另一个调用类的函数中？

0 个答案: