annotate get_gff_part.R @ 5:eca9213923b1
draft
"planemo upload commit d13521cdda334d0ce44b9a2bff047168a8bbaae0"
author |
petr-novak |
date |
Fri, 22 Apr 2022 10:26:11 +0000 |
parents |
ea6a3059a6af |
children |
|
rev |
line source |
0
|
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])
|