Short answer: Given a prime $p>2$, an infinite extra-special $p$ group is pseudofinite, and the quantifier-free formula $xy=yx$ witnesses the independence property (and so witnesses the order property too).
Details: I am basically just quoting from the Appendix in Definable envelopes in groups having a simple theory by Milliet. Let $p>2$ be prime. An extra-special $p$ group is a group $G$ such that $g^p=1$ for all $g\in G$, $Z(G)$ is cyclic of order $p$, and $Z(G)=[G,G]$.
It is well-known that there is a unique countably infinite extra-special $p$-group. In fact, Felgner showed that the theory of infinite extra-special $p$-groups is a well-defined complete $\aleph_0$-categorical first-order theory in the language of groups. The following is a concrete construction of the countable model.
Let $V$ be a vector space over $\mathbb{F}_p$ of dimension $\aleph_0$, and let $\langle\cdot,\cdot\rangle$ be a non-degenerate skew-symmetric bilinear form on $V$. Let $G=V\times\mathbb{F}_p$ and define the group operation $$ (u,a)\ast (v,b)=(u+v,a+b+\langle u,v\rangle). $$ Then $(G,\ast)$ is the countably infinite extra-special $p$-group.
Now, since the bilinear form is skew-symmetric, and $p>2$, it follows that $(u,a)\ast(v,b)=(v,b)\ast (u,a)$ if and only if $\langle u,v\rangle=0$. So we can witness the independence property for the formula $\varphi(x,y)$ given by $x\ast y=y\ast x$. Specifically, we fix $n$ and find $g_1,\ldots,g_n\in G$ such that for any $X\subseteq\{1,\ldots,n\}$ there is some $h_X\in G$ so that $g_i\ast h_X=h_X\ast g_i$ if and only if $i\in X$. Let $u_1,\ldots,u_n$ be linearly independent vectors in $V$ and set $g_i=(u_i,0)$. Since the bi-linear form is non-degenerate, the maps $v\mapsto \langle u_i,v\rangle$ are linearly independent. So we can find $v_X\in V$ such that $\langle u_i,v_X\rangle =0$ if $i\in X$ and $\langle u_i,v_X\rangle =1$ if $i\not\in X$. Now let $h_X=(v_X,0)$ and we have the desired elements of $G$.
Note that we are really just using the fact that $G$ interprets the structure $V$ (with the bilinear form), which has the independence property witnessed by the formula $\langle x,y\rangle =0$. (Actually these theories are bi-interpretable.)
Finally, to see that this theory is pseudofinite, one only needs to know that there are arbitrarily large finite extra-special $p$-groups, which is also well-known. Then any non-principal ultraproduct is a model of the theory. Actually this is a $1$-dimensional asymptotic class in the sense of Macpherson and Steinhorn (see Proposition 3.11 of One-dimensional asymptotic classes of finite structures).