颜色: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:


网友评论