Workflow: scRNA-seq pipeline using Salmon and Alevin
Fetched 2023-01-14 18:52:26 GMT
Verified with cwltool version 3.1.20221201130942
- Selected
- |
- Default Values
- Nested Workflows
- Tools
- Inputs/Outputs
This workflow is Open Source and may be reused according to the terms of:
GNU General Public License v3.0 only
Note that the tools invoked by the workflow may have separate licenses.
Inputs
| ID | Type | Title | Doc |
|---|---|---|---|
| assay | String | scRNA-seq assay | |
| threads | Integer | Number of threads for Salmon | |
| fastq_dir | Directory[] | Directory containing FASTQ files |
Steps
| ID | Runs | Label | Doc |
|---|---|---|---|
| fastqc |
steps/fastqc.cwl
(CommandLineTool)
|
Runs fastQC on each fastq file in fastq directory | |
| salmon |
steps/salmon.cwl
(CommandLineTool)
|
Run Salmon Alevin tool on FASTQ input | |
| annotate_cells |
steps/annotate-cells.cwl
(CommandLineTool)
|
Assay-specific annotation of cell barcodes after quantification | |
| adjust_barcodes |
steps/adjust-barcodes.cwl
(CommandLineTool)
|
Assay-specific adjustment of cell barcodes | |
| scanpy_analysis |
steps/scanpy-analysis.cwl
(CommandLineTool)
|
Dimensionality reduction and clustering | |
| alevin_to_anndata |
steps/alevin-to-anndata.cwl
(CommandLineTool)
|
Convert Alevin sparse output to anndata.AnnData object, save as h5ad |
Outputs
| ID | Type | Label | Doc |
|---|---|---|---|
| umap_pdf | File | UMAP dimensionality reduction plot | |
| fastqc_dir | Directory[] | Directory of FastQC output files, mirroring input directory structure | |
| qc_results | File | Quality control metrics | |
| count_matrix | File | Unfiltered count matrix from Alevin, converted to H5AD | |
| filtered_data | File | Full data set of filtered results: expression matrix, coordinates in dimensionality-reduced space (PCA and UMAP), cluster assignments via the Leiden algorithm, and marker genes for one cluster vs. rest | |
| salmon_output | Directory | Full output of `salmon alevin` | |
| full_count_matrix | File | Unfiltered count matrix from Alevin, converted to H5AD, with intronic regions | |
| marker_gene_plot_logreg | File | Cluster marker genes, logreg method | |
| marker_gene_plot_t_test | File | Cluster marker genes, t-test |
Permalink:
https://w3id.org/cwl/view/git/278890c1f502467a04f94e9c38f43cf1464e5cf5/pipeline.cwl
