update

						levels = colnames(design)

					)

	x %>%

	dt <-	x %>%

			voomWithQualityWeights(design, plot = F) %>%

			lmFit(design) %>%

			contrasts.fit(contrasts = contr.matrix) %>%

			topTable(n = Inf) %>%

			as.data.table() %>%

			.[,log10P := -log10(adj.P.Val)] %>%

	merge(gtf_info[,.(gene_id,gene_name,type)])

			merge(gtf_info)

	if("logFC" %in% colnames(dt)){

		setnames(dt,"logFC","log2FC")

	}

	dt

}

## the input data for venn_plot should come from rnaseq123_multi()