fig a

library("ggbio")
library("GenomicRanges")

data <- read.csv("IRN1000003263_TAIR10_12Sep2019_primary_targets.bed", header = FALSE,sep=" ",encoding="utf-8")
gr <- GRanges(seqnames = data$V1,ranges = IRanges(start =data$V2,end = data$V3))
seqlengths(gr)=c(30427671,19698289,23459830,18585056,26975502)
p<-autoplot(gr, layout = "karyogram",lwd=1)+theme(text = element_text(size=25))
	
jpeg(file="fig_a.jpg", width=10, height=4, units="in", res=800)
p
dev.off()