Tumor complexity and intratumor heterogeneity contribute to subclonal diversity. Despite advances in next-generation sequencing (NGS) and bioinformatics, detecting rare mutations in primary tumors and metastases contributing to subclonal diversity is a challenge for precision genomics. Here, in order to identify rare mutations, we adapted a recently described epithelial reprograming assay for short-term propagation of epithelial cells from primary and metastatic tumors. Using this approach, we expanded minor clones and obtained epithelial cell-specific DNA/RNA for quantitative NGS analysis. Comparative Ampliseq Comprehensive Cancer Panel sequence analyses were performed on DNA from unprocessed breast tumor and tumor cells propagated from the same tumor. We identified previously uncharacterized mutations present only in the cultured tumor cells, a subset of which has been reported in brain metastatic but not primary breast tumors. In addition, whole-genome sequencing identified mutations enriched in liver metastases of various cancers, including Notch pathway mutations/chromosomal inversions in 5/5 liver metastases, irrespective of cancer types. Mutations/rearrangements in FHIT, involved in purine metabolism, were detected in 4/5 liver metastases, and the same four liver metastases shared mutations in 32 genes, including mutations of different HLA-DR family members affecting OX40 signaling pathway, which could impact the immune response to metastatic cells. Pathway analyses of all mutated genes in liver metastases showed aberrant tumor necrosis factor and transforming growth factor signaling in metastatic cells. Epigenetic regulators including KMT2C/MLL3 and ARID1B, which are mutated in >50% of hepatocellular carcinomas, were also mutated in liver metastases. Thus, irrespective of cancer types, organ-specific metastases may share common genomic aberrations. Since recent studies show independent evolution of primary tumors and metastases and in most cases mutation burden is higher in metastases than primary tumors, the method described here may allow early detection of subclonal somatic alterations associated with metastatic progression and potentially identify therapeutically actionable, metastasis-specific genomic aberrations.