working script for DESeq2 analysis