Screening of intracellular parasites
To explore the presence of relevant intracellular parasites that could
explain their speciation process (60), we collapsed unique reads from
the butterfly-free sequences and removed from the dataset low complexity
sequences using Prinseq (61). Afterwards, we applied kraken2 to assign
reads against standard contaminant databases (bacteria, archaea, fungi,
protozoa and viral)(62).