Strand-specific paired end reads were filtered for rRNA sequences. Filtered reads were then aligned using TopHat (v1.4.1) against the D. melanogaster genome (FlyBase r5.44) and a maximum of six mismatches. Introns between 20 and 150,000 bp were allowed. snRNA, rRNA, tRNA, snoRNA and pseudogenes were masked.