reconstitute.py 文件源码

def go(args):
    if args.trimmed.endswith('.gz'):
        fh = gzip.open(args.trimmed)
    else:
        fh = open(args.trimmed)

    ids = [rec.id for rec in SeqIO.parse(fh, "fasta")]
    SeqIO.write((seq for seq in SeqIO.parse(args.fasta, "fasta") if seq.id in ids), sys.stdout, "fasta")