diff --git a/dge_workflow/featcounts_deseq.R b/dge_workflow/featcounts_deseq.R
index 089d4fc4363a42db6b9d7f1c2dbcd908bf7c93a3..157c0a7f0a8eb6b8eb42398e3e70538e9e5fcc77 100755
--- a/dge_workflow/featcounts_deseq.R
+++ b/dge_workflow/featcounts_deseq.R
@@ -406,8 +406,11 @@ degs <- deAnnot %>% filter(is_hit)
 #    inner_join(normCounts) %>%
 #    left_join(geneInfo) %>%
 #    left_join(bindCats)
-degs %>% write.delim(file=paste0(resultsBase, "diffex_genes.txt"))
-#'  [degs](`r paste0(resultsBase, "diffbind_genes.txt")`)
+degs %>% write_tsv(paste0(resultsBase, "diffex_genes.txt"))
+#' [Differentially Expressed Genes](`r paste0(resultsBase, "diffbind_genes.txt")`)
+
+#write.delim(degs, file="degs.txt")
+# degs <- read.delim("degs.txt")
 
 
 #unloadNamespace("dplyr"); require(dplyr)
@@ -450,9 +453,6 @@ degs %>% count( sample_1, sample_2, s1_overex) %>% kable()
 
 
 #'  DEGs (differentially expressed genes) are contained in the following table
-write.delim(degs, file="degs.txt")
-# degs <- read.delim("degs.txt")
-#' [Differentially Expressed Genes](degs.txt)
 
 
 #' DEGs can be browsed in Excel using the exported table or via the embedded table browser. To enable the IGV links, you need to set the port in the IGV preferences to 3334.