Immune Atlas data parser

parse_immune_atlas_data.py

"""
Downloaded from: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE118189
"""
with open('GSE118189_ATAC_counts.txt', 'r') as f:
	cell_files = None
	for line in f.readlines():
		line = line.rstrip()
		if not cell_files:
			cell_types = line.split('\t')
			cell_files = [open('./ImmuneAtlasBed/'  + cell + '.bed' , 'w') for cell in cell_types]
		else:
			values = line.split('\t')
			interval = values[0].split('_')
			counts = [int(x) for x in values[1:]]
			for i, file in enumerate(cell_files):
				file.write('\t'.join(interval + [str(counts[i])]) + '\n')
	for file in cell_files:
		file.close()