Suppose $U$ is a coordinate neighbourhood of $p$ in $M$, contained in the preimage under $F$ of a coordinate neighbourhood of $F(p)$ in $N$. We have a continuous map
$$ x \in U \subset \mathbb R^n \mapsto (dF)_x \in {\rm M}(n \times m , \mathbb R).$$
(Of course, I used the coordinate charts to identify $U$ with $\mathbb R^n$ and to identify the respective tangent spaces with $\mathbb R^n$ and $\mathbb R^m$.)
Since $F$ is a submersion at $x = p$, we know that the $n \times m$ matrix $(dF)_x$ is of maximal rank at $x = p$. Said another way, at least one of the $m \times m$ submatrices of $(dF)_p$ has non-vanishing determinant.
This suggests that we should consider the map $g : U \to \mathbb R^{\frac{n!}{m!(n-m)!}}$
that sends each point $x \in U$ to an $\frac{n!}{m!(n-m)!}$-component vector containing the determinants of all $m \times m$ submatrices of $(dF)_x$. Clearly $g$ is also continuous.
Finally, $g(p) \neq 0$. By continuity of $g$, there is an open neighbourhood $V \subset U$ of $p$ such that $g(x) \neq 0$ for all $x \in V$. This means that, for all $x \in V$, the matrix $(dF)_x$ has at least one $m \times m$ submatrix with non-vanishing determinant. So the original $F$ is a submersion on $V$.