The DELPHI detector located at LEP accelerator has been used to measure multi-muon bundles originated from cosmic ray interactions. Two subdetectors - Hadron Calorimeter and Time Projection Chamber, are used for this purpose. The 1999 and 2000 data are analyzed over wide range of muon multiplicities. The muon multiplicity distribution is compared with prediction of Monte-Carlo simulation based on CORSIKA/QGSJET. The Monte-Carlo does not describe the large multiplicity part of data. Even the extreme assumption on the cosmic ray composition (pure iron nuclei) hardly predicts number of high multiplicity events comparable with the data. The impact of QGSJET internal parameters onto the result is also studied
