Workflow: wf_clipseqcore_pe_1barcode.cwl

Fetched 2024-04-20 09:45:31 GMT

Workflow for handling reads containing one barcode. Returns the bam file containing read2 only. Notes: runs the following steps: - demultiplex - trimfirst_file2string - trimagain_file2string - b1_trim_and_map - view_r2 - index_r2_bam - make_bigwigs

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

Inputs

ID Type Title Doc
read
dataset String
chrom_sizes File
barcodesfasta File
randomer_length String
speciesGenomeDir Directory
repeatElementGenomeDir Directory

Steps

ID Runs Label Doc
view_r2
samtools-viewr2.cwl (CommandLineTool)

samtools-view.cwl is developed for CWL consortium Usage: samtools view [options] <in.bam>|<in.sam>|<in.cram> [region ...]

Options: -b output BAM -C output CRAM (requires -T) -1 use fast BAM compression (implies -b) -u uncompressed BAM output (implies -b) -h include header in SAM output -H print SAM header only (no alignments) -c print only the count of matching records -o FILE output file name [stdout] -U FILE output reads not selected by filters to FILE [null] -t FILE FILE listing reference names and lengths (see long help) [null] -T FILE reference sequence FASTA FILE [null] -L FILE only include reads overlapping this BED FILE [null] -r STR only include reads in read group STR [null] -R FILE only include reads with read group listed in FILE [null] -q INT only include reads with mapping quality >= INT [0] -l STR only include reads in library STR [null] -m INT only include reads with number of CIGAR operations consuming query sequence >= INT [0] -f INT only include reads with all bits set in INT set in FLAG [0] -F INT only include reads with none of the bits set in INT set in FLAG [0] -x STR read tag to strip (repeatable) [null] -B collapse the backward CIGAR operation -s FLOAT integer part sets seed of random number generator [0]; rest sets fraction of templates to subsample [no subsampling] -@ INT number of BAM compression threads [0]

demultiplex
index_r2_bam
samtools-index.cwl (CommandLineTool)

samtools-index.cwl is developed for CWL consortium

make_bigwigs
makebigwigfiles_PE.cwl (CommandLineTool)

Creates strand-specific bigwig files from a BAM file. See original script here: https://github.com/YeoLab/gscripts/blob/master/gscripts/general/make_bigwig_files_pe.py Usage: makebigwigfiles --bam BAM --genome GENOME --dont_flip --bw_pos --bw_neg

b1_trim_and_map

This workflow takes in appropriate trimming params and demultiplexed reads, and performs the following steps in order: trimx1, trimx2, fastq-sort, filter repeat elements, fastq-sort, genomic mapping, sort alignment, index alignment, namesort, PCR dedup, sort alignment, index alignment

trimagain_file2string
file2string.cwl (ExpressionTool)

Returns string expression based on file contents.

trimfirst_file2string
file2string.cwl (ExpressionTool)

Returns string expression based on file contents.

Outputs

ID Type Label Doc
output_neg_bw File
output_pos_bw File
output_r2_bam File
b1_trimx1_fastq File[]
b1_trimx2_fastq File[]
b1_trimx1_metrics File
b1_trimx2_metrics File
b1_mapgenome_stats File
b1_demuxed_fastq_r1 File Barcode1 read1 demultiplexed fastq
b1_demuxed_fastq_r2 File
b1_maprepeats_stats File
b1_output_sorted_bam File
b1_sorted_unmapped_fastq File[]
b1_trimx1_fastqc_stats_R1 File
b1_trimx1_fastqc_stats_R2 File
b1_trimx2_fastqc_stats_R1 File
b1_trimx2_fastqc_stats_R2 File
b1_mapgenome_star_settings File
b1_trimx1_fastqc_report_R1 File
b1_trimx1_fastqc_report_R2 File
b1_trimx2_fastqc_report_R1 File
b1_trimx2_fastqc_report_R2 File
b1_maprepeats_star_settings File
b1_mapgenome_mapped_to_genome File
b1_output_prermdup_sorted_bam File
b1_maprepeats_mapped_to_genome File
b1_output_barcodecollapsepe_bam File
b1_output_barcodecollapsepe_metrics File
Permalink: https://w3id.org/cwl/view/git/c0fffc4979a92371dc0667a03e3d957bf7f77600/cwl/wf_clipseqcore_pe_1barcode.cwl