https://github.com/mishashe/HGT
Tip revision: b32b6ebd11b49349893ec69fc4788cf7ede26003 authored by mishashe on 07 June 2021, 16:47:27 UTC
Update README.md
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);
}