A scenario is suggested in which the three light quark flavors are sequentially deconfined under increasing pressure in cold asymmetric nuclear matter as found, for example, in neutron stars. The basis for this analysis is a chiral quark matter model of Nambu-Jona-Lasinio (NJL) type with diquark pairing in the spin-1 single-flavor, spin-0 two-flavor, and three-flavor channels. Nucleon dissociation sets in at about the saturation density, n 0 , when the down-quark Fermi sea is populated (d-quark drip line) because of the flavor asymmetry induced by β equilibrium and charge neutrality. At about 3n 0 , u-quarks appear and a two-flavor color superconducting (2SC) phase is formed. The s-quark Fermi sea is populated only at still higher baryon density, when the quark chemical potential is of the order of the dynamically generated strange quark mass. Two different hybrid equations of state (EOSs) are constructed using the Dirac-Brueckner Hartree-Fock (DBHF) approach and the EOS of Shen et al. [H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Nucl. Phys. A637, 435 (1998)] in the nuclear matter sector. The corresponding hybrid star sequences have maximum masses of 2.1 and 2.0 M , respectively. Twoand three-flavor quark-matter phases exist only in gravitationally unstable hybrid star solutions in the DBHF case, whereas the Shen-based EOSs produce stable configurations with a 2SC phase component in the core of massive stars. Nucleon dissociation via d-quark drip could act as a deep crustal heating process, which apparently is required to explain superbursts and cooling of x-ray transients.