I am doing an RNA-seq data analysis. I already have several spesific gene as my target. What I want to do is actually comparing the gene expression level with the SNP mutation. So, for now, I want to get the gene expression level with this RNA-seq data. Because of my analysis focus on SNP and gene expression analysis, I don' have any interest in between group differential expression analysis. I already follow DESeq2 tutorial but I only get the fold change difference for 2 groups, normal and disease. So, the question is, I want to check the fold change to 1 of the sample as the base. I imagine like this. I have 14 sample. I will have fold change for difference between 1 and 2, 1 and 3, and so on, until 1 and 14. So, is this method possible to do in DESeq2 and probably you can share your idea how to do this with DESeq2. Thank you very much for your help.