Dmel R6.28
Reads were aligned to the Drosophila melanogaster reference genome using STAR v2.4.2a (Dobin et al., 2013). Uniquely mapped reads that overlap with genes were counted using HTSeq-count v0.7.1 (Anders et al., 2015), with cells having fewer than 300,000 uniquely mapped reads removed. Principal component analysis was performed using the top 500 most variable genes; typically 7 to 12 principal components were used for further analysis.