Cancer cell heterogeneity is a serious problem in the control of tumor progression because it can cause chemoresistance and metastasis. Heterogeneity can be generated by various mechanisms, including genetic evolution of cancer cells, cancer stem cells (CSCs), and niche heterogeneity. Because the genetic heterogeneity of CSCs has been poorly characterized, the genetic mutation status of CSCs was examined using Exome-Seq and RNA-Seq data of liver cancer.Here we show that different surface markers for liver cancer stem cells (LCSCs) showed a unique propensity for genetic mutations. Cluster of differentiation 133 (CD133)-positive cells showed frequent mutations in the IRF2, BAP1, and ERBB3 genes. However, leucine-rich repeat-containing G protein-coupled receptor 5-positive cells showed frequent mutations in the CTNNB1, RELN, and ROBO1 genes. In addition, some genetic mutations were frequently observed irrespective of the surface markers for LCSCs. BAP1 mutations was frequently observed in CD133-, CD24-, CD13-, CD90-, epithelial cell adhesion molecule-, or keratin 19-positive LCSCs. ASXL2, ERBB3, IRF2, TLX3, CPS1, and NFATC2 mutations were observed in more than three types of LCSCs, suggesting that common mechanisms for the development of these LCSCs. The present study provides genetic heterogeneity depending on the surface markers for LCSCs. The genetic heterogeneity of LCSCs should be considered in the development of LCSC-targeting therapeutics.