In Némethi's book "Normal surface singularities", Example 5.1.17, there is a formula to find the Seifert invariants of a Brieskorn complete intersection $\Sigma(a_1,...,a_n)$. I am interested in the case $n=3$. I was trying to find a counterexample of a Brieskorn intersection that is not a rational homology sphere according to Nemethi's Example 5.1.17. The conditions to be a rational homology sphere are the following:
(i) $(a_1,...,a_n) = (kp_1,kp_2,p_3,...,p_n)$ where $k\geq 1$ and the integers $\{p_j\}_{j=1}^n$ are pairwise coprime, and $\gcd(k,p_j) = 1$ for any $j\geq 3$.
(ii) $(a_1,..,a_n) = (2^c p_1,2p_2, 2p_3, p_4,...,p_n)$ where the integers $\{p_j\}_{j=1}^n$ are pairwise odd and coprime, and $c\geq 1$.
I picked $\Sigma(4,4,4)$, which should not satisfy either condition to be a rational homology sphere. However, I computed the Seifert invariants and got $(b_0;q_1(\alpha_1,\omega_1), q_2(\alpha_2,\omega_2), q_3(\alpha_3,\omega_3)) = (0;4(1,-1), 4(1,0), 4(1,0))$ where $q_k(\alpha_k,\omega_k)$ means that Seifert invariant is repeated $q_k$ times. The matrix whose cokernel is isomorphic to $H_1(M,\mathbb{Z})$ is
$$\begin{pmatrix} \alpha_1 & 0 & \cdots & 0 & \omega_1 \\ 0 & \alpha_2 & \ & 0 & \omega_2\\ \vdots & \ & \ddots & \ & \vdots\\ 0 & 0 & \ & \alpha_n & \omega_n\\ 1 & 1 & \cdots & 1 & b_0 \end{pmatrix}$$
Plugging in the Seifert invariants I found, this matrix has finite cokernel. This would mean $H_1(M,\mathbb{Z})$ is finite, hence $M$ is a rational homology sphere, which it shouldn't. What went wrong?
P.S. I tried for $\Sigma(4,8,12)$ as well and got the invariants $(0;4(1,0), 4(2,-1), 4(3,1))$, or the normalized version $(-4;4(1,0), 4(2,1), 4(3,1))$. The cokernel of the matrix above is still finite even though the tuple $(4,8,12)$ doesn't satisfy conditions (i) nor (ii).
Edit: these are the conditions to compute the Seifert invariants listed in Nemethi's book. $$A := \prod_{i=1}^n a_i,\ a:=lcm\{a_i:1\leq i\leq n\},\ \alpha_i = a/lcm\{a_j: j\neq i\}$$ $$\ q_i = \frac{A\alpha_i}{aa_i},\ \ \ \ \omega_ia/a_i \equiv -1 \mod \alpha_i$$ $$ -A/a^2 = -b_0 + \sum_j \omega_j/\alpha_j$$