There is no universally appropriate basis for delimiting species in protists, including parasites. Many molecular markers used for species delimitation are part of the ribosomal RNA (rRNA) array, with different regions of the array being used for different parasitic protist taxa. However, little is known about sequence variability across the rRNA in most organisms, and there is no standard threshold at which divergence in the sequence of a particular gene can be used as a basis for species delimitation. Here we demonstrate a method to generate the full rRNA array of parasitic protists by amplification of the array in two long, overlapping fragments followed by Illumina and Nanopore sequencing to produce high quality assemblies, to determine variations in sequence variability across the array. We apply this approach to two pairs of closely related ascetosporean parasites of crustaceans and molluscs (respectively Paramarteilia canceri/P. orchestiae and Marteilia cochillia/M. cocosarum (Rhizaria Endomyxa Ascetosporea)) and Bonamia ostreae and demonstrate how full-length rRNA sequences can be used to determine regions of the rRNA array that are most discriminatory, and robustly differentiate between species in combination with other lines of evidence.