Goal: quantify and describe the presence of different gene isoforms (variant transcripts of a gene produced by alternative splicing) in scRNA-seq data
- Creates random sampling of FASTQ files for cells within a specific subclass (for input in STAR script)
- Merges BAM files (the outputs of STAR script) for later visualization and quantification at the subclass-level
- Uses Bioconductor functions to set up and plot a figure visualizing the following:
- read coverage
- junctions (i.e., sashimi plot)
- gene isoform models
- Given FASTQ files, processed genome index directory, and genome annotation file (.gtf or .gff), creates scripts which will produce coordinate-sorted BAMs for gene quantification
- More details and documentation here