reference.py 文件源码

def write_genome_fasta(self, out_fasta_fn):
        if len(self.genomes) > 1:
            with open(out_fasta_fn, 'w') as f:
                for genome_prefix, in_fasta_fn in itertools.izip(self.genome_prefixes, self.in_fasta_fns):
                    with open(in_fasta_fn, 'r') as g:
                        for line in g:
                            line = line.strip()
                            if line.startswith('>'):
                                line = '>' + genome_prefix + '_' + line[1:]
                            f.write(line + '\n')
        else:
            cr_utils.copy(self.in_fasta_fns[0], out_fasta_fn)