https://github.com/mishashe/HGT
Raw File
Tip revision: b32b6ebd11b49349893ec69fc4788cf7ede26003 authored by mishashe on 07 June 2021, 16:47:27 UTC
Update README.md
Tip revision: b32b6eb
MakeHistFromMummer.c
//This program takes as an input mummer file Res.mumm and generates frequency table ResHist.mumm

int main MakeHistFromMumm()
{
	string line,foo;
	long long int L;
	ifstream MummerFile("Res.mumm");
	getline(MummerFile,line);
	fstream OutFile("ResHist.mumm",ios::out);
	OutFile<<line<<endl;
	map<unsigned ,unsigned long long int>  MLD;
	while (!(MummerFile.eof()))
	{
		getline(MummerFile,line);
		if (!(line.find('>')<line.size()))
		{
			istringstream iss(line);
			iss>>foo;iss>>foo;iss>>L;
			MLD[L]++;
		}
	}
	//output result
	for (map<unsigned ,unsigned long long int>::iterator  itMLD=MLD.begin();itMLD!=MLD.end();itMLD++)
	{
		OutFile<<(*itMLD).first<<" "<<(*itMLD).second<<endl;
	}
	MummerFile.close();
	OutFile.close();
	return(0);
}

back to top