I have a dataset consisting of CRISPR gRNA read counts coming from two different samples (something very similar to a RNA-seq experiment output).
I have transformed the read counts in log2 values and computed fold change between treated and non-treated sample. The distribution of the data is normal, but the mean is not = 0. I am plotting these data in a volcano plot, and the plot doesn't look right as I am plotting depleted vs enriched gRNAs but they do not correspond to negative vs positive values. So I thought to transform the values in z-scores. I wonder if this is correct. I have seen it is common to do that for microarray data, but I am not completely sure this applies to my data.
Many thanks for your help!
See below the plot: