“…Raw sequences reads (Haenel et al, 2019b , 2021 ) were parsed by library (pool or individual) and aligned to the third‐generation stickleback reference genome assembly (Glazer et al, 2015 ) by using novoalign (version 4.0, http://www.novocraft.com/products/novoalign/ ; alignment settings provided in the Supplementary Codes). From the alignments, we derived nucleotide counts (pileups) for all genome‐wide positions by using the pileup function from the Rsamtools r package (Morgan et al, 2017 ; unless specified otherwise, all analyses were implemented with the r language, version 3.6.0; r Development Core Team, 2019 ).…”