Mercurial > repos > petr-novak > repeat_annotation_pipeline3

comparison get_gff_part.R @ 0:ea6a3059a6af draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author petr-novak
date Mon, 18 Oct 2021 11:01:20 +0000
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:ea6a3059a6af
1 #!/usr/bin/env Rscript
2 suppressPackageStartupMessages(library(rtracklayer))
3 g = import(commandArgs(T)[1])
4 gd = mcols(g)
5 col_name=commandArgs(T)[2]
6 col_value=commandArgs(T)[3]
7 inc = gd[,col_name] %in% col_value
8 g_part = sort(sortSeqlevels(g[inc]))
9 export(g_part, format = 'gff3', commandArgs(T)[4])