I have 2 scRNA-seq data sets which come from different treatment conditions (Drug vs. Control)
I have been using Seurat and running the data through their standard integration pipeline, and at the end can get a result of DEGs between each condition for each cluster of cells.
What I am exploring is if its possible to obtain a fold-change for a given gene on a per-cell basis. For example, I would like to take the expression of gene "X" in cell #1 of cluster 0 for the Drug group, and compare it to all the cells in cluster 0 of the Control group and get a quantitative measure of how the expression in this one cell of the drug group is vs. the entire population of control cells.
Is something like this possible? I feel like it must be since the algorithm is essentially computing this same thing only in batch for all the cells in the cluster right? But Im afraid I can't think of a way to approach this?
Thanks in advance