合并kraken结果

setwd("/mnt/mzy/dairycow/24samples/kraken")
DNAsample<-c("127","134","133","143","140","122","125","117","149","132","130","142","145","152","118","147","121","151","135","136","150","146","141","137")
system("awk -F '\t' '{print $1}' *DNA|sort|uniq >uniq.names")
library(tidyverse)
read.table("uniq.names",header=F,sep="\t",quote="")->a
unique(a)->a
for (i in DNAsample){read.table(paste(i,".DNA",sep=""),header=F,sep="\t",quote="",fill=T)->b; colnames(b)[2]<-i;left_join(a,b,by="V1")->a;a[!duplicated(a[,1]),]->a;print(i);print(date())}
a[is.na(a)]<-0
write.table(a,file="krakenDNA.csv",row.names=F,col.names=T,quote=F,sep=",")