library_construction_protocol: hashtag-stained single cell suspension was pooled and approximately 28678 cells were loaded onto the BD Rhapsody system (BD Biosciences). After loading unique cell barcode-immobilized magnetic beads, cells were lysed, and reverse transcription of polyA RNA and hashtag oligos derived from each single cell were performed according to the manufacture's instructions. Exonuclease I treatment and cDNA/hashtag amplification were performed by using TAS-Seq protocol. Sequencing libraries were prepared by using NEBNext Ultra II FS library prep kit for Illumina (New England Biolabs) and quantified by KAPA library quantification kit (KAPA Biosystems) and a MultiNA system (Shimazu).