get_gene works for one gene

still a problem of initializing gene.seq