setwd("/mnt/mzy/dairycow/24samples/mappingrate/binsuper")
fileNames <- read.table("list")
read.table("117.DNA.abd")->first
first[,-(3:4)]->first
colnames(first)<-c("genome","genomelength")
for (i in 1:dim(fileNames)[1]){
read.table(paste(fileNames[i,1]))->tmp
tmp[,3]->first[paste(fileNames[i,1])]
}
write.csv(first,file="mappingrate.csv",row.names = F)
合并samtools idxstats结果
最新推荐文章于 2024-03-11 15:36:22 发布