Workflow: TOPMed_RNA-seq

Fetched 2021-03-04 15:38:37 GMT

TOPMed RNA-seq CWL workflow. Documentation on the workflow can be found [here]( Example input files: [Dockstore.json]( and [rnaseq_pipeline_fastq-example.yml]( Quickstart instructions are [here]( Start). [GitHub Repo]( Pipeline steps: 1. Align RNA-seq reads with [STAR v2.5.3a]( 2. Run [Picard]( [MarkDuplicates]( 2a. Create BAM index for MarkDuplicates BAM with [Samtools 1.6]( index. 3. Transcript quantification with [RSEM 1.3.0]( 4. Gene quantification and quality control with [RNA-SeQC 1.1.9](

children parents
Workflow as SVG
  • Selected
  • Default Values
  • Nested Workflows
  • Tools
  • Inputs/Outputs


ID Type Title Doc
rsem_ref_dir Directory
star_index Directory
fastqs File[]
prefix_str String
genome_fasta File
is_stranded Boolean
max_frag_len Integer
rnaseqc_flags String[]
paired_end Boolean
estimate_rspd Boolean
genes_gtf File


ID Runs Label Doc
indexbam.cwl (CommandLineTool)

A wrapper for running `samtools index <bam>`.

star.cwl (CommandLineTool)

A CWL wrapper for [](

Runs [STAR v2.5.3a](

This CWL Tool was developed as step 1 of the TOPMed RNA-seq workflow.

[GitHub Repo](

indexbam.cwl (CommandLineTool)

A wrapper for running `samtools index <bam>`.

rsem.cwl (CommandLineTool)

A CWL wrapper for [](

Runs [RSEM 1.3.0](

This CWL Tool was developed as step 3 of the TOPMed RNA-seq workflow.

markduplicates.cwl (CommandLineTool)

A CWL wrapper for [](

Runs [Picard]( [MarkDuplicates](

This CWL Tool was developed as step 2 of the TOPMed RNA-seq workflow.

[GitHub Repo](

rna_seqc.cwl (CommandLineTool)

A CWL wrapper for []( duplicated from []( with minor modifications.

Runs [RNA-SeQC 1.1.9](

This CWL Tool was developed as step 4 of the TOPMed RNA-seq workflow.

[GitHub Repo](

samtools-sort.cwl (CommandLineTool)

Sort alignments by leftmost coordinates, or by read name when -n is used. An appropriate @HD-SO sort order header tag will be added or an existing one updated if necessary.

Usage: samtools sort [-l level] [-m maxMem] [-o out.bam] [-O format] [-n] -T out.prefix [-@ threads] [in.bam]

Options: -l INT Set the desired compression level for the final output file, ranging from 0 (uncompressed) or 1 (fastest but minimal compression) to 9 (best compression but slowest to write), similarly to gzip(1)'s compression level setting.

If -l is not used, the default compression level will apply.

-n Sort by read names (i.e., the QNAME field) rather than by chromosomal coordinates.

-o FILE Write the final sorted output to FILE, rather than to standard output.

-O FORMAT Write the final output as sam, bam, or cram.

By default, samtools tries to select a format based on the -o filename extension; if output is to standard output or no format can be deduced, -O must be used.

-T PREFIX Write temporary files to PREFIX.nnnn.bam. This option is required.


ID Type Label Doc
rna-seqc_output_count_metrics File
markduplicates_output_bam File
markduplicates_bam_index File
rsem_output_isoforms_results File
rna-seqc_output_gene_rpkm File
markduplicates_output_metrics File
rsem_output_gene_results File
rna-seqc_output_count_outputs File
star_output_transcriptome_bam File
star_output_junctions_pass1 File
star_output_read_counts File
star_output_junctions File
star_output_chimeric_junctions File
star_output_bam_index File
star_output_logs File[]
star_output_bam File
rna-seqc_output_gene_counts File
rna-seqc_output_exon_counts File