This is long. I'm fairly sure it is correctI originally thought this would work but there is always the chance for a subtle errorit does not separate all extensions of $f$ and $g$ (see Nik's comment).
Consider the unital embedding $\mathcal A = M_2 \subset M_4$ given by $$A = \left[\begin{array}{cc}a&b\\c&d\end{array}\right] \mapsto A \otimes I_2 = \left[\begin{array}{cccc} a &&b\\&a&&b \\ c&&d\\&c&&d\end{array}\right]$$ and for $\gamma_n = -\sqrt \frac{n-1}{n} + \frac{i}{\sqrt{n}}\in\mathbb T$ let $$\mathcal B_n \equiv \left\{ \left[\begin{array}{cccc}a &&b \\ &a&&\gamma_nb\\c&&d\\&\bar{\gamma}_nc&&d\end{array}\right] : a,b,c,d\in \mathbb C\right\}$$ Then $\mathcal B$ is embedded in $\mathcal B_n$ which is a C$^*$-subalgebra of $M_4$. Finally, define $\mathcal B' \equiv \oplus_{n=1}^\infty \mathcal B_n \subset \oplus_{n=1}^\infty M_4 \equiv \mathcal A'$ with $A\in \mathcal A$ embedded as $\oplus_{n=1}^\infty (A\otimes I_2) \subset \mathcal A'$.
Claim 1: The distance from $P$ to $\mathcal B'$ is 1.
Let $\gamma_\infty = -1 = \lim_{n\rightarrow \infty} \gamma_n$ and define $\mathcal B_\infty$ in the same manner as above. What is the distance between $\mathcal B_\infty$ and $P\otimes I_2$? This is the same as calculating the minimum norm of $\left[\begin{array}{cc} 1+b&\\&1-b\end{array}\right], b\in\mathbb C$ which is 1 since $2 = |1 + b + 1 - b| \leq |1 + b| + |1 -b|$. Because of the definition of these matrix algebras it is straightforward to see that $\lim_{n\rightarrow \infty} d(\mathcal B_n, P\otimes I_2) = d(\mathcal B_\infty, P\otimes I_2) = 1$. Therefore, $$d(\mathcal B', \oplus_{n=1}^\infty (P\otimes I_2)) = 1$$
Claim 2: $\mathcal B'$ separates any two states on $\mathcal A'$ extending $f$ and $g$.
Let $v' = (v_1,v_2,v_3,v_4)\in \mathbb C^4$ such that $\langle (A\otimes I_2) v',v'\rangle = \langle Av,v\rangle$ for all $A\in M_2$. An easy calculation gives that $\|(v_1,v_2)\| = \frac{1}{\sqrt{2}}$ and $(v_1,v_2) = (v_3,v_4)$. The $w'$ case follows similarly except with $(w_1,w_2) = -(w_3,w_4)$.
Now for the general case, $\mathcal A'$ is naturally represented on the Hilbert space $\mathcal H = \oplus_{n=1}^\infty \mathbb C^4$. Let $v', w'\in \mathcal H$ be unit vectors such that $f'(C) = \langle Cv',v'\rangle$ and $g'(C) = \langle Cw',w'\rangle$ are states extending $f$ and $g$ respectively. From the previous paragraph they must have a particular form, namely $v' = (r_1,s_1,r_1,s_1, r_2,s_2,r_2,s_2,\cdots)$ and $w' = (t_1,u_1,-t_1,-u_1, t_2,u_2, -t_2,-u_2,\cdots)$.
Define $B_n = \left[\begin{array}{cccc}0&&1\\&0&&\gamma_n \\ \gamma_n&&0\\ &1&&0\end{array}\right] \in \mathcal B_n$ and let $B' = \oplus_{n=1}^\infty B_n \in \mathcal B'$. Now, $$Im\left\langle B_n\left(\begin{array}{c} r_n\\s_n\\r_n\\s_n\end{array}\right),\left(\begin{array}{c} r_n\\s_n\\r_n\\s_n\end{array}\right)\right\rangle = Im(1+\gamma_n)(\|r_n\|^2 + \|s_n\|^2) = \frac{1}{\sqrt{n}}(\|r_n\|^2 + \|s_n\|^2)$$ whereas $$Im\left\langle B_n\left(\begin{array}{c} t_n\\u_n\\-t_n\\-u_n\end{array}\right),\left(\begin{array}{c} t_n\\u_n\\-t_n\\-u_n\end{array}\right)\right\rangle = Im(1+\gamma_n)(-\|t_n\|^2 - \|u_n\|^2) = -\frac{1}{\sqrt{n}}(\|t_n\|^2 + \|u_n\|^2).$$ Hence, $$Im\langle B'v',v'\rangle = \sum_{n=1}^\infty \frac{1}{\sqrt{n}}(\|r_n\|^2+\|s_n\|^2) > 0$$ and $$Im\langle B'w',w'\rangle = -\sum_{n=1}^\infty \frac{1}{\sqrt{n}} (\|t_n\|^2 + \|u_n\|^2) < 0.$$ Therefore, $f'(B') \neq g'(B')$.