在Python中搜索genome.txt文件以查找特定实例

我目前正在尝试搜索包含基因组的.txt文件，查找“G”和“C”的特定实例，我可以通过将代码复制粘贴到程序中来让程序正常工作，但我不知道如何从文件中获取它。我猜问题在于代码跨越多行，但不知道如何解决它

以下是我在将代码复制到python中时使用的代码：

sequence='**strings of code**'
found=0

size(400, 100)
background(0)
fill(255)

for letter in sequence:
    if letter=='C' or letter=='G':
        background(0)
        found+=1
        text('The GC content of the sequence is: %', 50, 50)
        text((float(found)/(len(sequence)/100.0)),250,50)

这是我现在无法使用的代码：

f = open('sequence.txt', 'r').read()
sequence=f.split('\n')

found=0


size(400, 100)
background(0)
fill(255)


for letter in sequence:
    if letter=='C' or letter=='G':
    background(0)
    found+=1
    text('The GC content of the sequence is: %', 50, 50)
    text((float(found)/(len(sequence)/100.0)),250,50)

简短答复:

删除此行

sequence=f.split('\n')

和使用：

sequence = open('sequence.txt', 'r').read()

长答覆：

我假设第二段中的错误缩进不是这里的问题，真正的问题是它没有给你正确的答案

我可以在您使用的第一个代码片段中看到这一点

sequence='**strings of code**'

这使序列的值成为字符串，因此当您对其进行迭代时，字母变量将保留单个字符，但是当您尝试从文件中读取它时，您使用了拆分结果字符串，该字符串将序列从字符串转换为字符串列表，因此，这里的字母变成了一个完整字符串，表示“\n”之前的每一行，而不是像第一个代码段中那样的单个字母

---

所以，删除该行将使序列变成一个字符串，像以前一样迭代。

什么是“不能[sic]开始工作”的意思？如果你运行它会发生什么？另外，与此毫不相关的是：那些

大小

、

背景

、

填充

、和

文本

语句与您的实际问题有什么关系？对于这样一个简单的程序，使用

print

要容易得多。从一个基本的调试技术开始：在

For

循环中打印出

字母。这是你期望的吗？