This is a quite simple question: quarks do mix (through the CKM matrix), neutrinos do mix (through the PMNS matrix).
Then why do charged leptons not mix?
This is a quite simple question: quarks do mix (through the CKM matrix), neutrinos do mix (through the PMNS matrix).
Then why do charged leptons not mix?
Short answer: only the down-type quarks mix by the CKM matrix, by convention and without loss of generality. The lepton sector is completely analogous to the quark sector where the up-type quarks play the role of the neutrinos and the down-type quarks play the role of the charged leptons.
Long answer: Let's recall some facts about the standard model. You have three quark weak-doublets $Q_i = (u_L, d_L), (c_L, s_L), (t_L, b_L)$ for $i=1,2,3$, and six quark singlets $u_{Ri}=(u_R, c_R, t_R)$ and similarly for the down-type quarks $d_{Ri}$. The Yukawa couplings for these quarks, which are what end up giving them mass after Higgsing, look like this (I have made no attempt to match this to some standard notation, this should just serve to give the idea):
$$ \mathcal{L} \supset y^U_{ij} Q_i \Phi u_{Rj}+ y^D_{ij} Q_j \tilde{\Phi} d_{Rj} + h.c. $$
If it weren't for these couplings the standard model would have an $(SU(3))^3$ global symmetry under three independent unitaries acting on the generation indices of the doublet, up- and down-type singlets respectively. In other words we could make field redefinitions
$$Q_{i} \to U_{ij}^{Q}Q_{j}\\ u_{Ri} \to U_{ij}^{U}u_{Rj}\\ d_{Ri} \to U_{ij}^{D}d_{Rj} $$
with three independent unitaries $U^{Q,U,D}$ without changing anything else in the Lagrangian apart from the Yukawa terms. That means we could choose $U^Q$ and $U^U$ to diagonalise the up-type quark term (A general fact of matrix algebra: you can diagonalise a matrix if you are free to multiply on the right and left by independent unitaries). So you can diagonalise $y^U$ at the expense of choosing $U^Q$ and $U^U$. Now we try to diagonalise $y^D$, but there's a hitch: we're no longer free to choose $U^Q$. We can only choose $U^D$, so we can't completely diagonalise $y^D$. After symmetry breaking there will inevitably be off-diagonal terms (unless there is some really incredible unexpected flavour symmetry!) which will cause oscillations between the down-type quarks (d,s,b). The CKM matrix, which is related in some way to $y^D, U^Q$ and $U^D$, measures this.
So how are leptons different? You just replace $u\to e, d\to\nu_e,$ etc. In the pure standard model there are no right-handed neutrino singlets analogous to $d_{Ri}$. So you can diagonalise the lepton Yukawas and be done with it. But if you bring in an extension of the standard model to explain neutrino oscillations (which requires at least two right handed neutrinos - I'll assume for simplicity there are three, one for each flavour), then you will have a Yukawa matrix for the neutrinos like there was for the down-type quarks, either appearing directly in your theory or as an effective operator after integrating out some more complicated hidden sector. Now the same logic applies as before. You can always diagonalise the charged lepton Yukawas, no worries, but you don't have enough freedom to completely diagonalize the neutrino Yukawas. There is an analog of the CKM matrix called the PMNS matrix in this context.
This is not to say that neutrino oscillations have no (potentially) observable consequences on the charged leptons. Consider this diagram:
This is clearly forbidden as a real process by energy conservation, but if you radiate a photon then you can get the real lepton flavour violating decay $\mu^- \to e^- \gamma$ which people have looked for. (Check out the Particle Data Group for the current limits on the rate.) Neuneck points out that the standard model rate for this process is very small, so, while it is technically nonzero in SM, any observation of this process will be strong evidence for beyond standard model physics.
It is as @dmckee says, the weak eigenstates are the same as the mass eigenstates for the leptons, therefore the corresponding Cabibo angle is zero.
Without going into the mathematics this plot says all:
The Cabibbo angle represents the rotation of the mass eigenstate vector space formed by the mass eigenstates . θC = 13.04°.
The CKM matrix is an experimental fact:
Note, however, that the specific values of the angles are not a prediction of the standard model: they are open, unfixed parameters. At this time, there is no generally accepted theory that explains why the measured values are what they are.
Thus that the mass eigenstates and the weak eigenstates for leptons have zero angle is an experimental observation too.
Generally physics theories answer how one state implies/transforms into another. "Why" questions can be answered by "how" answers until one reaches the "that is what the experiment says". This is one of those "why" questions.
Somewhat reluctantly, I decided to add to the otherwise excellent technical answers above, since none confronted the fundamental false premise of your question, "why do charged leptons not mix?". Of course they do. Charged weak currents merely link generations.
Let me review its antecedents as you seem to be aware of the phenomenon, when you have all quarks, ups and downs, mix, and not just the downs, as popular convention has it. Consider the early picture of quarks with only 2 generations, some time in the 70s, but without K-M. Schematically, skipping chiral structures and Lorentz indices, the charged current vertex in the weak effective hamiltonian is $$ \overline{\begin{bmatrix} u \\ c \end{bmatrix} }\cdot \begin{bmatrix} d^\prime \\ s^\prime \end{bmatrix} ~W^+= \overline{\begin{bmatrix} u \\ c \end{bmatrix} }\cdot \begin{bmatrix} \cos{\theta_\mathrm{c}} & \sin{\theta_\mathrm{c}} \\ -\sin{\theta_\mathrm{c}} & \cos{\theta_\mathrm{c}}\\ \end{bmatrix} \begin{bmatrix} d \\ s \end{bmatrix} W^+, $$ plus the hermitean conjugate, $$ \overline{\begin{bmatrix} d' \\ s' \end{bmatrix} }\cdot \begin{bmatrix} u \\ c \end{bmatrix} ~W^-= \overline{\begin{bmatrix} d \\ s \end{bmatrix} } \begin{bmatrix} \cos{\theta_\mathrm{c}} & -\sin{\theta_\mathrm{c}} \\ \sin{\theta_\mathrm{c}} & \cos{\theta_\mathrm{c}}\\ \end{bmatrix}\cdot \begin{bmatrix} u \\ c \end{bmatrix} W^- . $$ You did not ask about any up vs. down asymmetries, because there aren't any: if you chose, you could have left the down quarks alone and "mixed" the ups instead, $$ \begin{bmatrix} u' \\ c' \end{bmatrix}= \begin{bmatrix} \cos{\theta_\mathrm{c}} & -\sin{\theta_\mathrm{c}} \\ \sin{\theta_\mathrm{c}} & \cos{\theta_\mathrm{c}}\\ \end{bmatrix} \begin{bmatrix} u \\ c \end{bmatrix}. $$ Since the weak vertices are specified by the same expression, it never occurs to you to ask who mixes: the up or the down quarks, because they all do.
Historically, in 1960, when Gell-Mann and Levy introduced this mixing, only the u, d, s quarks were known, so of course, they were looking only at the upper row of the Cabbibo matrix and naturally they had d and s mix into d' and s'. But, in full fairness, an s weak-couples to a "mixture" of u and c.
So, barring Majorana masses, and making the (monstrous?) assumption the PMNS matrix is as parameter-poor as the CKM, the analog vertex for just 2 generations (think MNS, 1962) is $$ \overline{\begin{bmatrix} e \\ \mu \end{bmatrix} }\cdot \begin{bmatrix} \nu_e \\ \nu_\mu \end{bmatrix} ~W^-= \overline{\begin{bmatrix} e \\ \mu \end{bmatrix} } \begin{bmatrix} \cos{\theta_\mathrm{p}} & -\sin{\theta_\mathrm{p}} \\ \sin{\theta_\mathrm{p}} & \cos{\theta_\mathrm{p}}\\ \end{bmatrix} \begin{bmatrix} \nu_1 \\ \nu_2 \end{bmatrix} W^- . $$ This is a hypothetical world, since the actual PMNS mixing is huge compared to the quark mixing, and no realistic split like this simplistic example could be made; so $\nu_1,\nu_2$ are not quite the actual mass eigenstates $\nu_1,\nu_2$s of today's neutrino experiments: this is just to illustrate the point.
But you see the point: the r.h.side expression, written completely in terms of mass eigenstates of leptons, does not much care who is charged and who is neutral. Since we know so little about the $\nu_i$s it would be perverse to pretend the charged leptons mix into $l_1, ~l_2$, and some of the answers justify the logic of sticking to the present convention.
But, deep down, the formal structures have never failed to be the same. I urge my colleagues to give the $\nu_i$s memorable names (Huey, Dewey, and Ratatouille, anything!) so they are thought of as the real, physical MCoys, would have distinctive flavors if they were quarks!, unlike the odd formal entities of convenience in production and detection, like $\nu_e$, $\nu_\mu$, $\nu_\tau$, with which we are stuck today, an effective confusion machine; but they are reluctant to go there before firmly pinning down the hierarchy of masses...
NB: Looks like, finally!, SM school charts are beginning to clean up their act, listing neutrino L,M,H mass eigenstates in their listings of generations, aligning with quark practice.