Content Comparison

library(tidyverse)
AmfReads <- read.csv("/Users/gregg/Downloads/AMF1_Test_filtermergestats.csv", header = FALSE)
Map <- read.csv("/Users/gregg/Downloads/Bo_Stevens_Sample_Positions.csv", header=TRUE)
Map <- Map[,c(1,2,4)]
Map$order <- c(1:96, 1:92)
AMFmeta <- data.frame(do.call('rbind', strsplit(as.character(AmfReads$V1),'.',fixed=TRUE)))
AmfReads <- cbind(AMFmeta$X2, AmfReads[,2:4])
names(AmfReads) <- c("Sample.ID", "Reads", "Reads2", "Reads3")
AmfReads <- inner_join(AmfReads, Map, by = "Sample.ID")
LowAmf <- AmfReads[ which(AmfReads$Reads < 500),]
write.csv(LowAmf, "/Users/gregg/Downloads/LowAmf.csv", quote = FALSE, row.names = FALSE)