Function generateDistFile

ParseAligns/ParseAligns.cpp:261–288

}

static void
generateDistFile()
{
	ofstream distFile(opt::distPath.c_str());
	assert(distFile.is_open());
	for (EstimateMap::iterator mapIt = estMap.begin(); mapIt != estMap.end(); ++mapIt) {
		// Skip empty iterators
		assert(!mapIt->second.estimates[0].empty() || !mapIt->second.estimates[1].empty());
		distFile << mapIt->first;
		for (int refIsRC = 0; refIsRC <= 1; refIsRC++) {
			if (refIsRC)
				distFile << " ;";

			for (Estimates::iterator vecIt = mapIt->second.estimates[refIsRC].begin();
			     vecIt != mapIt->second.estimates[refIsRC].end();
			     ++vecIt) {
				vecIt->second.distance =
				    (int)round((double)vecIt->second.distance / (double)vecIt->second.numPairs);
				if (vecIt->second.numPairs >= opt::c && vecIt->second.numPairs != 0
				    /*&& vecIt->distance > 1 - opt::k*/)
					distFile << ' ' << get(g_contigNames, vecIt->first) << ',' << vecIt->second;
			}
		}
		distFile << '\n';
		assert_good(distFile, opt::distPath);
	}
	distFile.close();
}

static bool
isSingleEnd(const string& id);