First of all, thank you very much for your attention. This is my first time participating in a forum of this kind, and any help is greatly appreciated.
I am new to the field of bioinformatics and currently using the Mageck tool (I’ve been working with Mageck count, Mageck test, and Mageck mle) within the Galaxy platform to analyze a single-cell screen experiment I conducted.
My main challenge lies in how to structure the samples (I have extracted the counts from NGS) when performing the corresponding comparisons. Specifically, my experiment has two genotypes (WT and KO) that were infected with an sgRNA library. After a positive selection, I collected samples on day 5. Additionally, I have two replicates for each experiment in both genotypes, resulting in a total of 8 samples (see the diagram).
My primary interest is identifying which guides are overrepresented in the KO condition compared to the standard WT distribution. I can handle this comparison using Mageck test (RRA) by setting the WT genotype as the control sample.
However, I would like to account for the samples from day 0 to control for any potential biases in the initial guide distribution between both genotypes, ensuring that they are not mistakenly interpreted as enriched guides in the KO vs. WT comparison.
How could I achieve this initial standardization? Is it possible to do so using Mageck? Is there any other tool that would be better suited for my experimental design?
Thank you very much for your help!