I am comparing the gene expression of a bacterium grown in two conditions (~ 10 replicates each). I am not interested in the mean gene expression differences between the two, rather I am interested in variability of gene expression within each condition. I have seen several papers addressing this, and even this recent package: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5737414/ . However, I would like to stick with DESeq2 if possible, because it is what I use for everything else. Is there a way to determine which genes are 'differentially variable' (i.e. have a higher/lower spread) between two conditions?
I found the discussion here:
https://support.bioconductor.org/p/130742/#130765