generate_bed_file.py
View Options

			# simply generate a bed file with two region on two chromosomes

			import os

			file_bed = os.path.join("data", "toy_data", "peaks.bed")

			# write the cell barcodes
			# add an extra barcode associated with no reads -> empty cell
			with open(file_bed, "wt") as f_out:
			f_out.write("chr1\t990\t1010\n")
			f_out.write("chr2\t990\t1010\n")

scATAC-seq/scripts/toy_data/generate_bed_file.pyfaf6f1d50a68master