INNER CODE UNIT · C++

col_id

refresh-bio/FAMSA · src/msa_refinement.cpp:53

	size_t col_id = tmp[rnd_rfn() % tmp.size()].first;

	int first_prof_id = 0;
	int second_prof_id = 1;

	if (profile_to_refine->data[0]->GetSymbol(col_id) == GAP)
		swap(first_prof_id, second_prof_id);

	for (size_t j = 0; j < card; ++j)
		if (profile_to_refine->data[j]->GetSymbol(col_id) == GAP)
			dest_prof_id.emplace_back(first_prof_id);
		else
			dest_prof_id.emplace_back(second_prof_id);
}

// *******************************************************************
// Refine alignment
#ifdef DEBUG_MODE

View source record →