Workflow: 04-peakcall-se.cwl

Fetched 2022-08-09 18:56:17 GMT

ATAC-seq 04 quantification - SE

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


ID Type Title Doc
nthreads Integer
input_bam_files File[]
as_narrowPeak_file File

Definition narrowPeak file in AutoSql format (used in bedToBigBed)

input_genome_sizes File

Two column tab-delimited file with chromosome size information

genome_effective_size String

Effective genome size used by MACS2. It can be numeric or a shortcuts:'hs' for human (2.7e9), 'mm' for mouse (1.87e9), 'ce' for C. elegans (9e7) and 'dm' for fruitfly (1.2e8), Default:hs


ID Runs Label Doc
../spp/spp.cwl (CommandLineTool)
../utils/count-with-output-suffix.cwl (CommandLineTool)

Counts lines in a file and returns a suffixed file with that number

../peak_calling/macs2-callpeak.cwl (CommandLineTool)
../utils/trunk-peak-score.cwl (CommandLineTool)

Trunk scores in ENCODE bed6+4 files

../quant/bedToBigBed.cwl (CommandLineTool)

\"bedToBigBed v. 2.7 - Convert bed file to bigBed. (BigBed version: 4) usage: bedToBigBed in.bed chrom.sizes Where in.bed is in one of the ascii bed formats, but not including track lines and chrom.sizes is two column: <chromosome name> <size in bases> and is the output indexed big bed file. Use the script: fetchChromSizes to obtain the actual chrom.sizes information from UCSC, please do not make up a chrom sizes from your own information. The in.bed file must be sorted by chromosome,start, to sort a bed file, use the unix sort command: sort -k1,1 -k2,2n unsorted.bed > sorted.bed\"

../peak_calling/count-reads-after-filtering.cwl (CommandLineTool)

Count number of dedup-ed reads used in peak calling

../peak_calling/samtools-filter-in-bedfile.cwl (CommandLineTool)

Filter BAM file to only include reads overlapping with a BED file

../spp/extract-best-frag-length.cwl (CommandLineTool)

Extracts best fragment length from SPP output text file

../peak_calling/samtools-extract-number-mapped-reads.cwl (CommandLineTool)

Extract mapped reads from BAM file using Samtools flagstat command


ID Type Label Doc
output_peak_file File[]

peakshift/phantomPeak results file

output_peak_xls_file File[]

Peak calling report file (*_peaks.xls file produced by MACS2)

output_peak_bigbed_file File[]

Peaks in bigBed format

output_spp_x_cross_corr File[]

peakshift/phantomPeak results file

output_peak_summits_file File[]

File containing peak summits

output_extended_peak_file File[]

peakshift/phantomPeak extended fragment results file

output_spp_cross_corr_plot File[]

peakshift/phantomPeak results file

output_filtered_read_count_file File[]

Filtered read count reported by MACS2

output_peak_count_within_replicate File[]

Peak counts within replicate

output_read_in_peak_count_within_replicate File[]

Reads peak counts within replicate