颜色: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:
截图
网友评论