diff --git a/scripts/plot.R b/scripts/plot.R
new file mode 100644
index 0000000..3f502a4
--- /dev/null
+++ b/scripts/plot.R
@@ -0,0 +1,15 @@
+library(ggplot2)
+
+gene <- read.csv2(snakemake@input[[1]], dec = ",")
+colnames(gene) <- "id"
+gene$id <- as.numeric(gene$id)
+
+print(paste0("I'm plotting the file ", snakemake@wildcards[["file"]]))
+
+pdf(file = snakemake@output[[1]])
+
+ggplot(gene, aes(x = id)) +
+    geom_density(alpha = .2, fill = "#FF6666") +
+    ggtitle(paste0("Plot "), snakemake@wildcards[["file"]])
+
+dev.off()
\ No newline at end of file