def genomes2include(path):
	"""
	Figure out which genomes I want to include in my dataset based on those in selected_kegg.txt
	
	Arguments:
		path -- path to 'selected_kegg.txt'
	
	Returns: 
		tla_to_tnum (dict) -- converts three-letter abbreviations (e.g.: eco) to t-numbers (e.g.: T04989) 
		tnum_to_tla (dict) -- maps tnum to tla for each genome
		keepers (list) -- list of t-numbers to keep (e.g.: T04989)
	"""
	
	# Use phylogenetically thinned list in selected_genomes.txt AND filter out non-bacteria
	path = path+"selected_kegg.txt"
	file = open(path).readlines()
	file = list(map(str.strip, file))
	
	# The keepers list ends up having a few genomes that do not get used
	# This is bc they have no full modules
	keepers = [] 
	# Create dict that converts tla (e.g.: Pea) -> t_num (e.g.: T321890)
	tla_to_tnum = {}
	
	line_counter = 0

	for s in file:
#		if line_counter < 4:
#			line_counter += 1
#			continue
		tla = s.split()[1]
		t_num = s.split()[2]
		tax = s.split()[3]
		   
		if "k__Bacteria" in tax:
			keepers.append(t_num)
			tla_to_tnum[tla] = t_num
	
	tnum_to_tla = {v:k for k,v in tla_to_tnum.items()}
	
	return tla_to_tnum, tnum_to_tla, keepers