Mercurial > repos > petr-novak > repeat_annotation_pipeline2

diff get_gff_part.R @ 0:cf3cea0a3039 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author petr-novak
date Thu, 07 Oct 2021 06:07:34 +0000
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/get_gff_part.R	Thu Oct 07 06:07:34 2021 +0000
@@ -0,0 +1,9 @@
+#!/usr/bin/env Rscript
+suppressPackageStartupMessages(library(rtracklayer))
+g = import(commandArgs(T)[1])
+gd = mcols(g)
+col_name=commandArgs(T)[2]
+col_value=commandArgs(T)[3]
+inc = gd[,col_name] %in% col_value
+g_part = sort(sortSeqlevels(g[inc]))
+export(g_part, format = 'gff3', commandArgs(T)[4])