I am new to Deseq analysis. Recently i conducted metatranscriptomics analysis of bacterial community in two different conditions. My aim to explore the overall significent transcript expression and also some selective transcripts (related to specific function). I used raw featurecounts output to run Deseq and created heatmap on VSD transformed data on to 20 features.
But, due to vast number of transcripts, this does not cover my transcripts of interest. Is it ok to extract raw counts of all transcripts related to a particular function and run Deseq?
Also, the heatmap generated via vsd transformed data shows scale as -1, -0.5, 0.5 and 1.0. Are this suppossed to be the Log2F change?