A weakly-singular, symmetric Galerkin boundary element method (SGBEM) is developed for analysis of fractures in three-dimensional, anisotropic linearly elastic media. The method constitutes a generalization of that by Li et al. [S. Li, M.E. Mear, L. Xiao, Symmetric weak-form integral equation method for three-dimensional fracture analysis, Comput. Methods Appl. Mech. Engrg. 151 (1998) 435-459] for isotropic media, and is based upon a pair of weak-form displacement and traction integral equations recently established by Rungamornrat and Mear [J. Rungamornrat, M.E. Mear, Weakly-singular, weakform integral equations for cracks in three-dimensional anisotropic media, Int. J. Solids Struct. 45 (2008Struct. 45 ( ) 1283Struct. 45 ( -1301]]. The formulation involves only weakly-singular kernels and, as a consequence, standard C 0 elements can be employed in the numerical discretization. The kernels possess a relatively simple form (for general anisotropy) which involves an equatorial line integral like that associated with the displacement fundamental solution. Despite their simple form, it is still necessary to evaluate these kernels efficiently in the numerical implementation; to do so, certain symmetry properties associated with the specific material type under consideration are exploited, and a simple yet accurate interpolation strategy is introduced. Another important aspect of the numerical treatment is the use of a special crack-tip element which allows highly accurate mixed-mode stress intensity factor data to be extracted as a function of position along the crack front. This crack-tip element was originally introduced by Li et al. [S. Li, M.E. Mear, L. Xiao, Symmetric weak-form integral equation method for three-dimensional fracture analysis, Comput. Methods Appl. Mech. Engrg. 151 (1998) 435-459] in their treatment of isotropic media, and here it is extended to allow evaluation of the stress intensity factors for generally anisotropic media. The element involves ''extra" degrees of freedom associated with the nodes along the crack front (with these quantities being directly related to the stress intensity factors), and it is capable of representing the asymptotic behavior in the vicinity of the crack front to a sufficiently high order that relatively large crack-tip elements can be utilized. Various examples are treated for cracks in an unbounded domain and for both embedded and surface-breaking cracks in a finite domain, and it is demonstrated that very accurate results can be obtained even with relatively coarse meshes.