A three-dimensional nonlinear finite element for thin beams is proposed within the absolute nodal coordinate formulation (ANCF). The deformation of the element is described by means of displacement vector, axial slope and axial rotation parameter per node. The element is based on the Bernoulli-Euler theory and can undergo coupled axial extension, bending and torsion in the large deformation case. Singularities-which are typically caused by such parameterizations-are overcome by a director per element node. Once the directors are properly defined, a cross sectional frame is defined at any point of the beam axis. Since the director is updated during computation, no singularities occur. The proposed element is a three-dimensional ANCF Bernoulli-Euler beam element free of singularities and without transverse slope vectors. Detailed convergence analysis by means of various numerical static and dynamic examples and comparison to analytical solutions shows the performance and accuracy of the element.