gene = cell.loci.get_or_create(id='gene_{}_{}'.format(
                    i_chr + 1, i_gene + 1), __type=wc_kb.prokaryote.GeneLocus)