python - 串联多个.fasta文件

我试图将数百个.fasta文件连接到一个包含所有序列的大型fasta文件中。我还没有在论坛上找到一个具体的方法来实现这一点。我确实在http://zientzilaria.heroku.com/blog/2007/10/29/merging-single-or-multiple-sequence-fasta-files中遇到了这个代码，我对它进行了一些修改。
Fasta.py包含以下代码：

class fasta:
    def __init__(self, name, sequence):
        self.name = name
        self.sequence = sequence

def read_fasta(file):
    items = []
    index = 0
    for line in file:
        if line.startswith(">"):
           if index >= 1:
               items.append(aninstance)
           index+=1
           name = line[:-1]
           seq = ''
           aninstance = fasta(name, seq)
        else:
           seq += line[:-1]
           aninstance = fasta(name, seq)

    items.append(aninstance)
    return items

下面是连接.fasta文件的改编脚本：

import sys
import glob
import fasta

#obtain directory containing single fasta files for query
filepattern = input('Filename pattern to match: ')

#obtain output directory
outfile = input('Filename of output file: ')

#create new output file
output = open(outfile, 'w')

#initialize lists
names = []
seqs = []

#glob.glob returns a list of files that match the pattern
for file in glob.glob(filepattern):

    print ("file: " + file)

    #we read the contents and an instance of the class is returned
    contents = fasta.read_fasta(open(file).readlines())

    #a file can contain more than one sequence so we read them in a loop
    for item in contents:
        names.append(item.name)
        seqs.append(item.sequence)

#we print the output
for i in range(len(names)):
    output.write(names[i] + '\n' + seqs[i] + '\n\n')

output.close()
print("done")

它可以读取fasta文件，但新创建的输出文件不包含序列。我收到的错误是由于fasta.py，这超出了我处理的能力：

Traceback (most recent call last):
  File "C:\Python32\myfiles\test\3\Fasta_Concatenate.py", line 28, in <module>
    contents = fasta.read_fasta(open(file).readlines())
  File "C:\Python32\lib\fasta.py", line 18, in read_fasta
    seq += line[:-1]
UnboundLocalError: local variable 'seq' referenced before assignment

有什么建议吗？谢谢！

最佳答案

我认为用python来做这项工作是过分了。在命令行中，将单个/多个fasta文件与.fasta或.fa扩展名连接的快速方法是：

cat *.fa* > newfile.txt

关于python - 串联多个.fasta文件，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/11726340/