在多个fasta fi中计数发生适配器

Question

我在一个文件夹中有30个fastq文件，我想知道在其中的哪个文件中可以找到一个特定的适配器（这样我就可以找出它实际上是哪个示例）。你知道吗

我写了一个很小的biopython脚本，但它只允许我一次查看一个文件，我想计算每个文件同时出现的次数。有人能帮我改进剧本吗？你知道吗

 from Bio import SeqIO
 adaptor = (rec for rec in \
                SeqIO.parse("file.fastq", "fastq") \
                if rec.seq.startswith("TGA"))`
    count = SeqIO.write(adaptor, "adaptor.fastq", "fastq")
    print("Saved %i adaptor" % count)

Answer 1

from Bio import SeqIO
fnames = ["file.fastq", "file1.fastq", "file2.fastq"]
for fname in fnames:
    adaptor = (rec for rec in \
                    SeqIO.parse(fname, "fastq") \
                    if rec.seq.startswith("TGA"))
    count = SeqIO.write(adaptor, "adaptor.fastq", "fastq")
    print("Saved %i adaptor in file %s" %(count, fname))