“…RNA-seq datasets were first quality assessed using FastQC (Andrews, 2010), trimmed with BBDuk (parameters ktrim = r, k = 21, mink = 10, hdist = 1) (Bushnell, 2015), aligned using HISAT2 (option -U -S), and sorted and indexed with Samtools (Li et al, 2009). Splicing events were visualized in IGV (sashimi plot function) to enumerate the intron-spanning reads between avian exons 4, 5, and 6 (Katz et al, 2015). …”