absolute abundance = (number of bases mapped to genome / total number of bases in sample) * 100

def absolute_abundance(coverage, total_bases):
	"""
	absolute abundance = (number of bases mapped to genome / total number of bases in sample) * 100
	"""
	absolute = {}
	for genome in coverage:
		absolute[genome] = []
		index = 0
		for calc in coverage[genome]:
			bases = calc[0]
			total = total_bases[index]
			absolute[genome].append((bases / total) * float(100))
			index += 1
	total_assembled = [0 for i in absolute[genome]]
	for genome in absolute:
		index = 0
		for cov in absolute[genome]:
			total_assembled[index] += cov
			index += 1
	absolute['Unassembled'] = [(100 - i) for i in total_assembled]
	return absolute