Tumors often show intra-tumor heterogeneity because of genotypic differences between all the cells that compose it and that derive from it. Recent studies have shown significant aspects of neuroblastoma heterogeneity that may affect the diagnostic-therapeutic strategy. Therefore, we developed a laboratory protocol, based on the combination of the advanced dielectrophoresis-based array technology and next-generation sequencing to identify and sort single cells individually and carry out their copy number variants analysis. The aim was to evaluate the cellular heterogeneity, avoiding overestimation or underestimation errors, due to a bulk analysis of the sample. We tested the above-mentioned protocol on two neuroblastoma cell lines, SK-N-BE(2)-C and IMR-32. The presence of several gain or loss chromosomal regions, in both cell lines, shows a high heterogeneity of the copy number variants status of the single tumor cells, even if they belong to an immortalized cell line. This finding confirms that each cell can potentially accumulate different alterations that can modulate its behavior. The laboratory protocol proposed herein provides a tool able to identify prevalent behaviors, and at the same time highlights the presence of particular clusters that deviate from them. Finally, it could be applicable to many other types of cancer.