Newer
Older
#!/usr/bin/env Rscript
args = commandArgs(trailingOnly=TRUE)
source("https://bioconductor.org/biocLite.R")
biocLite("karyoploteR")
bed = read.table(args[1], sep="\t", header=TRUE)
#Read output plot files?
#chrom as "chr1"
kp <- plotKaryotype(chromosomes=c(chrom))
#kpDataBackground(kp)
kpAxis(kp)
kpRect(kp, chr=chrom, x0=start, x1=end, y0=0.2, y1=0.4)