Update gene_name_mapping.R

gagneurlab · Jul 27, 2020 · 8679afe · 8679afe
1 parent d254aca
commit 8679afe
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/drop/modules/mae-pipeline/Scripts/MAE/gene_name_mapping.R b/drop/modules/mae-pipeline/Scripts/MAE/gene_name_mapping.R
@@ -25,7 +25,7 @@ gtf_dt <- import(snakemake@input$gtf) %>% as.data.table
 if (!"gene_name" %in% colnames(gtf_dt)) {
   gtf_dt[gene_name := gene_id]
 }
-gtf_dt <- gtf_dt[type == "gene", .(seqnames, start, end, strand, gene_id, gene_name, gene_type, gene_status)]
+gtf_dt <- gtf_dt[type == "gene", .(seqnames, start, end, strand, gene_id, gene_name, gene_type)]
 
 # make gene_names unique
 gtf_dt[, N := 1:.N, by = gene_name] # warning message