The answer to both question is «Yes». The OP is asking what kind of singularity the foliation of the Riemann sphere induced by the differential equation in Alexandre's answer possesses. It can be rewritten as
$$\dot{z} = \overline{R(z)}=\frac{|R|^2}{R}(z)$$ where $R(z)=\sum_j\frac{a_j}{z-z_j}$. Now this real foliation is the same as the one obtained by cancelling the real term $|R(z)|^2$ altogether (by reparametrization of the time), and is therefore given by a rational differential equation
$$ \dot{z}=\frac{1}{R(z)}$$

Setting $x:= \frac{1}{z}$ and taking into account that $\sum_ja_j=0$ yields the differential equation
$$ \dot{x}=-x^2R^{-1}(x^{-1})=x^{-\nu}(c+\ldots)$$ for some $c\neq 0$ and $\nu\in\mathbb{Z_{\leq 0}}$.

The local conformal type of such foliations is well-known, and is classified by $\nu$: the differential equation is analytically conjugate near $x=0$ to $$\dot{x}=x^{-\nu}$$which can be solved explicitly to understand what happens. In the generic case $\nu=0$, that is $\sum_ja_jz_j\neq0$, the foliation is regular and there is exactly one smooth curve passing through $x=0$. Otherwise you obtain $\nu+1$ curves crossing $x=0$, corresponding to directions $\Im(x^{\nu+1})=0$. Notice in particular that in the original coordinates the OP exceptional lines (called *separatrices*) are tangent to diameters of a regular $2\nu+2$-gone.

Global phase portraits of rational vector fields have been studied intensely since the works of Smale [1] on a differential Newton method to obtain root-finding algorithms.

Of particular interest is the separatrix graph: to determine the graph whose edges link infinity to an electric charge. And the associated inverse problem: to determine which graphs are attained that way (see in particular works by Shub-Tischler-Williams, Berzinger, Sverdlove…). This combinatorial data classifies the global topological type of the phase portrait of field lines. It turns out that pretty much any graph with obvious restrictions (related to the dynamical nature of the question: it is planar since trajectories cannot cross outside $\infty$ and the charges, and its cycles must be compatible with the relative orientation given by the force field) does the job.

The polynomial case (not directly of interest to the OP, but it contains already a lot of information on the rational case) has been generically studied by Douady-Estrada-Sentenac (but their work was unfortunately unpublished). The current complete reference on the topic is Branner-Dias [2]. The recent thesis of Tomasini under the guidance of Lei nearly provides a complete description of the separatrix graph of rational foliations of the Riemann sphere.

It is expected that starting from computer presentable data $R$ there exists an algorithm outputting its separatrix graph.

[1] Smale, Steve, "On the efficiency of algorithms of analysis", Bull. Amer. Math. Soc. (N.S.) 13, 2 (1985), pp. 87--121.

[2] Branner, Bodil and Dias, Kealey, "Classification of complex polynomial vector fields in one complex variable", J. Difference Equ. Appl. 16, 5-6 (2010), pp. 463--517.