modernized final annotation filtering

dge_workflow/featcounts_deseq.R

@@ -265,6 +265,8 @@ normCounts <- counts(dds,normalized=TRUE) %>%
     set_names(colnames(countMatrix)) %>% rownames2column("ensembl_gene_id")
 
 
+normCounts %>% write.delim("normalized_counts_by_replicate.txt")
+
 ## .. should be same as input
 #filter(countData, ensembl_gene_id=="ENSDARG00000000001")