Mapping reads to genome

To check how the simulation data compares with the input data, we need to map reads to the genome. For this, we can use the STAR aligner (https://github.com/alexdobin/STAR).

Input:

1. Directory with files of reads, one file per cell
2. File with the genome sequence

Output: Bam-format files with alignments of reads to genome