# genome coverage

setwd("/Users/christineewers/Dropbox/Chelonibia/1.5_Chelonibia_structure/Data_files/singleCopyGene_coverage")

acc <- scan("nuccore_result.txt", what="integer")

cl <- c("no.nuc", "no.seq", "per.ident.sites", "per.pair.identity", "max.seq.length", "mean.cov", "min.seq.length", "acc.ID", "seq.length")

cov <- 0
for(i in 1:length(acc)) {
  cov[i] <- read.table(acc[i], sep=",", header=F)[1,6]
}

summary(cov)
hist(cov, breaks=10)

genome.size <- 6.2/cov
hist(genome.size)
mean(genome.size)
summary(genome.size)