Mercurial > repos > petr-novak > repeat_annotation_pipeline3
view get_gff_part.R @ 6:b53f5a456d01 draft
"planemo upload commit 3aefb0555456837d10fe69e4ad25de08d5972cb2"
author | petr-novak |
---|---|
date | Thu, 19 May 2022 08:23:55 +0000 |
parents | ea6a3059a6af |
children |
line wrap: on
line source
#!/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])