Question

当我尝试修改染色体名称不遵循我的输入文件染色体名称的代码时，我遇到了这个错误。基本上，下面的代码是读取输入文件并通知较短序列的位置，并根据文件中给出的信息输出位置序列。例如，在chr4：154742507-154742714中，值151表示第一个碱基的位置和'CCCAGGCTGG＆＃39;位置是173 - 182.因此，使用下面的代码应该能够通过添加173到154742507返回我的确切位置并获得下面的输出。任何人都可以帮助我吗？

下面是带有示例输入文本文件的代码。

INPUT.TXT

chr4:154742507-154742714


                           CCCAGGCTGG
151  AGTCTTGCTTTTTTTGTCGTTGCCCAGGCTGGAGTGCAGTGGCACCATCTCGGCTCAC


chr9:47303792-47303999

      CCAGCCTGGG
1    TCCAGCCTGGGTGACAGCGTGAGGCTCTTGTCTCAAATAGAAAAAAAACAAAGAACAAAAAACAAAAAACCACCA

输出

chr1    154742680   154742690
chr1    47303794    47303804

预期输出

chr4    154742680   154742690
chr9    47303794    47303804

代码

import re  # regular expressions, not needed (alternatives: the `split` method) but convenient

result = []
output_file=open('output.bed','w')
with open('Input.txt') as f:
    for line in f:
        if line.startswith('chr'):
            label = line.strip()
        elif line[0] == ' ':
            # short sequence
            length = len(line.strip())
            # find the index of the beginning of the short sequence
            for i, c in enumerate(line):
                if c.isalpha():
                    short_index = i
                    break
        elif line[0].isdigit():
            # long sequence
            n = line.split(' ')[0]
            # find the index of the beginning of the long sequence
            for i, c in enumerate(line):
                if c.isalpha():
                    long_index = i
                    break
            start = int(n) + short_index - long_index
            start -= 1
            end = start + length
            result.append('{} {} {}'.format(label, start, end))
            offset, n, start, length = 0, 0, 0, 0
output_line= "\n".join(result)
output_file.write(output_line)
output_file.close()

output_file=open('last_output.bed','w')
with open('output.bed') as fin:
    for line in fin:
        start, _, offset_start, offset_end = re.search(r'[^:]*:(\d+)\D+(\d+)\D+(\d+)\D+(\d+)', line).groups()
        output_line=('chr1\t{}\t{}\n'.format(int(start) + int(offset_start) + 1,int(start) + int(offset_end) + 1))
        output_file.write(output_line)
output_file.close()

Answer 1

如果我正确理解了这个问题，那么您所遇到的问题仅与输出错误的染色体编号（chr##）有关。

这似乎有点明显。在代码的最后，你很难编码：

output_line=('chr1\t{}\t{}\n'.format(stuff))

如果您不希望输出始终显示chr1，则需要更改输出。

前一行的正则表达式似乎与文件中的染色体编号匹配，您只是不在以后可以使用的组中捕获它。尝试：

chromosome, start, _, offset_start, offset_end = re.search(r'([^:]*):(\d+)\D+(\d+)\D+(\d+)\D+(\d+)', line).groups()
output_line=('{}\t{}\t{}\n'.format(chromosome, int(start) + int(offset_start) + 1,int(start) + int(offset_end) + 1))

这仍然相当丑陋，但应该有效。请注意，如果从初始循环中获取正确的输出，而不是写出中间格式然后需要重新分析它，那将会更容易。

Python：如何在输出中输出正确的染色体名称？

1 个答案: