R sushi

library(Sushi)

g1 <- read.delim("genes1.bed")
g1$strand <- as.character(g1$strand)

p1 = plotGenes(g1, chrom = "region_[1_16742]", chromstart = 900, chromend = 16000,
               plotgenetype="arrow",
               bentline=F, arrowlength = 0.025, packrow = T,
               bheight = 0.04, labeloffset = 0.08,
               wigglefactor = 0.001)
labelgenome(chrom = "region_[1_16742]", chromstart = 900, chromend = 16000, n=10, scale="bp", pos = 0.9)