C++ (Cpp) Genome::init_genomeの例

コード例 #1

ファイル: gene_tools.cpp プロジェクト: JonasBehr/trsk

vector<Region*> GeneTools::init_regions(const char* gio_fname)
{
	Genome* gio = new Genome(); 
	gio->init_genome((char*) gio_fname);
	const char* strands = "+-";

	vector<Region*> regions; 
	for (int i=0; i<gio->num_contigs; i++)
	{
		for (int s=0; s<2; s++)
		{
			int start = 1; 
			int stop = gio->contig_len(i)-1;
			Region* reg = new Region(start, stop, i, strands[s]);
			reg->set_gio(gio);
			regions.push_back(reg);
		}
	}
	return regions;
}