画差异基因的染色体分布

颜色：http://www.5tu.cn/colors/yansezhongwenming.html

网站：http://mg2c.iask.in/mg2c_v2.0/

R中取得两个文件：

setwd("D:/biodata/chrplots")

getwd()

allgene <- read.table("gene.good.bed",header=F,sep="\t",row.names = 1)

gene <- read.table("gene.txt",header=F,sep="\t")

result <- allgene[t(gene),]

write.table(result,file="result.txt",quote = F,sep = "\t",col.names=F)

allgene：

gene.txt:

mm10.chrom.sizes:

result.txt: 

截图