The recent data suggest that the areas of enhancer action are restricted by partitioning of the genome into topologically-associating domains, TADs. However, most of the observations on the 3D genome organization have been made by conventional methods in cell population where only average characteristics of the so-called typical cell can be identified. On the other hand, FISH-based studies demonstrated that distances between various genomic regions vary significantly in individual cells. To get further insights into the functional significance of 3D genome organization it is necessary to estimate the plasticity of this organization. The high throughput chromosome conformation capture protocol (Hi-C) has been modified recently to allow construction of chromatin contact frequency maps for individual cells. Using this modified protocol we have constructed Hi-C maps for 20 Drosophila cells (line Dm-BG3c2). In the best cell we have captured ~12% of the theoretically available contacts. To analyze these sparse contact matrices, we have developed program tools allowing us to take into account the noise by comparing maps from individual cells with artificially generated random matrices. The results of our analysis demonstrate that contact chromatin domains and chromatin compartments do exist in individual Drosopila cells. We also observed the hierarchical organization of contact chromatin domains in individual cells. Consequently this hierarchical organization does not represent a population average but is an intrinsic feature of chromatin folding. The population TAD borders are well reproduced in individual cells and correlate with certain epigenetic signatures. Presence of TADs in individual cells argues that TADs do not represent a population average as postulated by chromatin loop extrusion model but rather originate due to condensation of chromatin domains possessing specific epigenetic signatures.