@prefix LoadListingRequirement: .
@prefix CommandLineBinding: .
@prefix xsd: .
@prefix DockerRequirement: .
@prefix Workflow: .
@prefix sld: .
@prefix rdfs: .
@prefix CommandOutputBinding: .
@prefix ns1: .
@prefix NetworkAccess: .
@prefix rdf: .
@prefix ScatterMethod: .
@prefix cwl: .
rdfs:comment "\tCalculate coverage of 5\" positions (instead of entire interval).\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-5"
] ;
sld:type xsd:boolean , sld:null .
rdfs:comment "--skipNonCoveredRegions, --skipNAs\nThis parameter determines if non-covered regions\n(regions without overlapping reads) in a BAM file\nshould be skipped. The default is to treat those\nregions as having a value of zero. The decision to\nskip non-covered regions depends on the interpretation\nof the data. Non-covered regions may represent, for\nexample, repetitive regions that should be skipped.\n(default: False)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--skipNonCoveredRegions"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--skipNonCoveredRegions"
] ;
sld:type xsd:boolean , sld:null .
a cwl:CommandLineTool ;
rdfs:comment "usage: An example usage is:$ bamCoverage -b reads.bam -o coverage.bw\n\nThis tool takes an alignment of reads or fragments as input (BAM file) and\ngenerates a coverage track (bigWig or bedGraph) as output. The coverage is\ncalculated as the number of reads per bin, where bins are short consecutive\ncounting windows of a defined size. It is possible to extended the length of\nthe reads to better reflect the actual fragment length. *bamCoverage* offers\nnormalization by scaling factor, Reads Per Kilobase per Million mapped reads\n(RPKM), and 1x depth (reads per genome coverage, RPGC).\nRequired arguments:\n --bam BAM file, -b BAM file\n BAM file to process (default: None)\nOutput:\n --outFileName FILENAME, -o FILENAME\n Output file name. (default: None)\n --outFileFormat {bigwig,bedgraph}, -of {bigwig,bedgraph}\n Output file type. Either \"bigwig\" or \"bedgraph\".\n (default: bigwig)\nOptional arguments:\n --help, -h show this help message and exit\n --scaleFactor SCALEFACTOR\n The smooth length defines a window, larger than the\n binSize, to average the number of reads. For example,\n if the –binSize is set to 20 and the –smoothLength is\n set to 60, then, for each bin, the average of the bin\n and its left and right neighbors is considered.\n Any value smaller than –binSize will be ignored and\n no smoothing will be applied. (default: 1.0)\n --MNase Determine nucleosome positions from MNase-seq data.\n Only 3 nucleotides at the center of each fragment are\n counted. The fragment ends are defined by the two mate\n reads. Only fragment lengthsbetween 130 - 200 bp are\n considered to avoid dinucleosomes or other\n artifacts.*NOTE*: Requires paired-end data. A bin size\n of 1 is recommended. (default: False)\n --filterRNAstrand {forward,reverse}\n Selects RNA-seq reads (single-end or paired-end) in\n the given strand. (default: None)\n --version show program's version number and exit\n --binSize INT bp, -bs INT bp\n Size of the bins, in bases, for the output of the\n bigwig/bedgraph file. (default: 50)\n --region CHR:START:END, -r CHR:START:END\n Region of the genome to limit the operation to - this\n is useful when testing parameters to reduce the\n computing time. The format is chr:start:end, for\n example --region chr10 or --region\n chr10:456700:891000. (default: None)\n --blackListFileName BED file, -bl BED file\n A BED file containing regions that should be excluded\n from all analyses. Currently this works by rejecting\n genomic chunks that happen to overlap an entry.\n Consequently, for BAM files, if a read partially\n overlaps a blacklisted region or a fragment spans over\n it, then the read/fragment might still be considered.\n (default: None)\n --numberOfProcessors INT, -p INT\n Number of processors to use. Type \"max/2\" to use half\n the maximum number of processors or \"max\" to use all\n available processors. (default: max/2)\n --verbose, -v Set to see processing messages. (default: False)\nRead coverage normalization options:\n --normalizeTo1x EFFECTIVE GENOME SIZE LENGTH\n Report read coverage normalized to 1x sequencing depth\n (also known as Reads Per Genomic Content (RPGC)).\n Sequencing depth is defined as: (total number of\n mapped reads * fragment length) / effective genome\n size. The scaling factor used is the inverse of the\n sequencing depth computed for the sample to match the\n 1x coverage. To use this option, the effective genome\n size has to be indicated after the option. The\n effective genome size is the portion of the genome\n that is mappable. Large fractions of the genome are\n stretches of NNNN that should be discarded. Also, if\n repetitive regions were not included in the mapping of\n reads, the effective genome size needs to be adjusted\n accordingly. Common values are: mm9: 2,150,570,000;\n hg19:2,451,960,000; dm3:121,400,000 and\n ce10:93,260,000. See Table 2 of http://www.plosone.org\n /article/info:doi/10.1371/journal.pone.0030377 or http\n ://www.nature.com/nbt/journal/v27/n1/fig_tab/nbt.1518_\n T1.html for several effective genome sizes. (default:\n None)\n--ignoreForNormalization IGNOREFORNORMALIZATION [IGNOREFORNORMALIZATION ...]\n A list of space-delimited chromosome names containing\n those chromosomes that should be excluded for\n computing the normalization. This is useful when\n considering samples with unequal coverage across\n chromosomes, like male samples. An usage examples is\n --ignoreForNormalization chrX chrM. (default: None)\n --skipNonCoveredRegions, --skipNAs\n This parameter determines if non-covered regions\n (regions without overlapping reads) in a BAM file\n should be skipped. The default is to treat those\n regions as having a value of zero. The decision to\n skip non-covered regions depends on the interpretation\n of the data. Non-covered regions may represent, for\n example, repetitive regions that should be skipped.\n (default: False)\n --smoothLength INT bp\n The smooth length defines a window, larger than the\n binSize, to average the number of reads. For example,\n if the --binSize is set to 20 and the --smoothLength\n is set to 60, then, for each bin, the average of the\n bin and its left and right neighbors is considered.\n Any value smaller than --binSize will be ignored and\n no smoothing will be applied. (default: None)\nRead processing options:\n --extendReads [INT bp], -e [INT bp]\n This parameter allows the extension of reads to\n fragment size. If set, each read is extended, without\n exception. *NOTE*: This feature is generally NOT\n recommended for spliced-read data, such as RNA-seq, as\n it would extend reads over skipped regions. *Single-\n end*: Requires a user specified value for the final\n fragment length. Reads that already exceed this\n fragment length will not be extended. *Paired-end*:\n Reads with mates are always extended to match the\n fragment size defined by the two read mates. Unmated\n reads, mate reads that map too far apart (>4x fragment\n length) or even map to different chromosomes are\n treated like single-end reads. The input of a fragment\n length value is optional. If no value is specified, it\n is estimated from the data (mean of the fragment size\n of all mate reads). (default: False)\n --ignoreDuplicates If set, reads that have the same orientation and start\n position will be considered only once. If reads are\n paired, the mate's position also has to coincide to\n ignore a read. (default: False)\n --minMappingQuality INT\n If set, only reads that have a mapping quality score\n of at least this are considered. (default: None)\n --centerReads By adding this option, reads are centered with respect\n to the fragment length. For paired-end data, the read\n is centered at the fragment length defined by the two\n ends of the fragment. For single-end data, the given\n fragment length is used. This option is useful to get\n a sharper signal around enriched regions. (default:\n False)\n --samFlagInclude INT Include reads based on the SAM flag. For example, to\n get only reads that are the first mate, use a flag of\n 64. This is useful to count properly paired reads only\n once, as otherwise the second mate will be also\n considered for the coverage. (default: None)\n --samFlagExclude INT Exclude reads based on the SAM flag. For example, to\n get only reads that map to the forward strand, use\n --samFlagExclude 16, where 16 is the SAM flag for\n reads that map to the reverse strand. (default: None)\n" ;
ns1:original_cwlVersion "v1.2" ;
cwl:arguments ( [ CommandLineBinding:position 3 ;
CommandLineBinding:prefix "--outFileName" ;
cwl:valueFrom "${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bdg\"; return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bw\"; }"
]
) ;
cwl:arguments ( [ CommandLineBinding:position 3 ;
CommandLineBinding:prefix "--outFileName" ;
cwl:valueFrom "${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bdg\"; return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bw\"; }"
]
) ;
cwl:baseCommand ( "bamCoverage" ) ;
cwl:baseCommand ( "bamCoverage" ) ;
cwl:cwlVersion ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "reddylab/deeptools:3.0.1"
] ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "reddylab/deeptools:3.0.1"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs , , , , , , , , , , , , , , , , , , , , , , ;
cwl:outputs ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] .
cwl:valueFrom ".raw.bw" .
rdfs:comment "Bigwig file, -b1 Bigwig file\nBigwig file 1. Usually the file for the treatment.\n(default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--bigwig1"
] ;
sld:type cwl:File .
rdfs:comment "Suffix used to replace bigwig1 extension (please, include extension)\n" ;
sld:type xsd:string , sld:null .
rdfs:comment "Suffix used for output file (input BAM filename + suffix)" ;
sld:type xsd:string , sld:null .
rdfs:comment "" ;
cwl:inputBinding [ CommandLineBinding:position 3 ;
CommandLineBinding:prefix "-g"
] ;
sld:type cwl:File .
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
rdfs:comment "{bigwig,bedgraph}, -of {bigwig,bedgraph}\nOutput file type. Either \"bigwig\" or \"bedgraph\".\n(default: bigwig)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--outFileFormat"
] ;
sld:type xsd:string , sld:null .
rdfs:comment "INT bp\nThis parameter allows the extension of reads to\nfragment size. If set, each read is extended, without\nexception. *NOTE*: This feature is generally NOT\nrecommended for spliced-read data, such as RNA-seq, as\nit would extend reads over skipped regions. *Single-\nend*: Requires a user specified value for the final\nfragment length. Reads that already exceed this\nfragment length will not be extended. *Paired-end*:\nReads with mates are always extended to match the\nfragment size defined by the two read mates. Unmated\nreads, mate reads that map too far apart (>4x fragment\nlength) or even map to different chromosomes are\ntreated like single-end reads. The input of a fragment\nlength value is optional. If no value is specified, it\nis estimated from the data (mean of the fragment size\nof all mate reads). (default: False)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--extendReads"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--extendReads"
] ;
sld:type xsd:int , sld:null .
rdfs:comment "{bigwig,bedgraph}, -of {bigwig,bedgraph}\nOutput file type. Either \"bigwig\" or \"bedgraph\".\n(default: bigwig)\n" ;
cwl:default "bigwig" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--outFileFormat"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--outFileFormat"
] ;
sld:type xsd:string .
cwl:valueFrom "${return 200}" .
cwl:default 1 ;
sld:type xsd:int .
rdfs:comment "INT\nNumber of processors to use. Type \"max/2\" to use half\nthe maximum number of processors or \"max\" to use all\navailable processors. (default: max/2)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--numberOfProcessors"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--numberOfProcessors"
] ;
sld:type xsd:int , sld:null .
cwl:outputBinding [ CommandOutputBinding:glob
"${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bdg\"; return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bw\"; }" ] ;
cwl:outputBinding [ CommandOutputBinding:glob
"${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bdg\"; return inputs.bam.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".bw\"; }" ] ;
sld:type cwl:File .
rdfs:comment "show program's version number and exit" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--version"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--version"
] ;
sld:type xsd:boolean , sld:null .
cwl:in , , ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) .
cwl:valueFrom "RPKM" .
cwl:valueFrom ".ctrl_subtracted.bw" .
rdfs:comment "PSEUDOCOUNT\nsmall number to avoid x/0. Only useful when ratio =\nlog2 or ratio (default: 1)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--pseudocount"
] ;
sld:type xsd:int , sld:null .
rdfs:comment "SCALEFACTORS\nSet this parameter to multipy the bigwig values by a\nconstant. The format is scaleFactor1:scaleFactor2. For\nexample 0.7:1 to scale the first bigwig file by 0.7\nwhile not scaling the second bigwig file (default:\nNone)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--scaleFactors"
] ;
sld:type xsd:float , sld:null .
rdfs:comment "If set, reads that have the same orientation and start\nposition will be considered only once. If reads are\npaired, the mate's position also has to coincide to\nignore a read. (default: False)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--ignoreDuplicates"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--ignoreDuplicates"
] ;
sld:type xsd:boolean , sld:null .
cwl:source .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.ibam.path.replace(/^.*[\\\\\\/]/, '').replace(/\\.[^/.]+$/, '') + '.bdg')" ] ;
sld:type cwl:File .
rdfs:comment "\tCalculate coverage of 3\" positions (instead of entire interval).\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-3"
] ;
sld:type xsd:boolean , sld:null .
cwl:source .
rdfs:comment "INT \nExclude reads based on the SAM flag. For example, to\nget only reads that map to the forward strand, use\n--samFlagExclude 16, where 16 is the SAM flag for\nreads that map to the reverse strand. (default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--samFlagExclude"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--samFlagExclude"
] ;
sld:type xsd:int , sld:null .
rdfs:comment "\tScale the coverage by a constant factor.\nEach coverage value is multiplied by this factor before being reported.\nUseful for normalizing coverage by, e.g., reads per million (RPM).\n- Default is 1.0; i.e., unscaled.\n- (FLOAT)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-scale"
] ;
sld:type xsd:float , sld:null .
cwl:valueFrom "${return 1}" .
rdfs:comment "\tReport the depth at each genome position (with zero-based coordinates).\nReports only non-zero positions.\nDefault behavior is to report a histogram.\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-dz"
] ;
sld:type xsd:boolean , sld:null .
cwl:valueFrom "${return 1}" .
cwl:valueFrom "bigwig" .
cwl:in ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) .
rdfs:comment "By adding this option, reads are centered with respect\nto the fragment length. For paired-end data, the read\nis centered at the fragment length defined by the two\nends of the fragment. For single-end data, the given\nfragment length is used. This option is useful to get\na sharper signal around enriched regions. (default:\nFalse)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--centerReads"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--centerReads"
] ;
sld:type xsd:boolean , sld:null .
rdfs:comment "Fragment extended RPKM bigWig (signal) files" ;
cwl:outputSource ;
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
rdfs:comment "INT bp\nThe smooth length defines a window, larger than the\nbinSize, to average the number of reads. For example,\nif the --binSize is set to 20 and the --smoothLength\nis set to 60, then, for each bin, the average of the\nbin and its left and right neighbors is considered.\nAny value smaller than --binSize will be ignored and\nno smoothing will be applied. (default: None)\nRead processing options:\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--smoothLength"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--smoothLength"
] ;
sld:type xsd:int , sld:null .
cwl:source .
rdfs:comment "\tbed_graph is a four column file in the format: \n" ;
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type cwl:File .
rdfs:comment "FILENAME\nOutput file name. (default: input BAM filename with bigwig [*.bw] or bedgraph [*.bdg] extension.)\n" ;
sld:type xsd:string , sld:null .
rdfs:comment "CHR:START:END, -r CHR:START:END\nRegion of the genome to limit the operation to - this\nis useful when testing parameters to reduce the\ncomputing time. The format is chr:start:end, for\nexample --region chr10 or --region\nchr10:456700:891000. (default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--region"
] ;
sld:type xsd:string , sld:null .
rdfs:comment "Bed or bedGraph file to be sorted" ;
cwl:inputBinding [ CommandLineBinding:position
1 ] ;
sld:type cwl:File .
rdfs:comment "Bigwig file, -b2 Bigwig file\nBigwig file 2. Usually the file for the control.\n(default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--bigwig2"
] ;
sld:type cwl:File .
cwl:valueFrom "${return 1}" .
cwl:source .
rdfs:comment "\tReport depth in BedGraph format. For details, see:\ngenome.ucsc.edu/goldenPath/help/bedgraph.html\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-bg"
] ;
sld:type xsd:boolean , sld:null .
rdfs:comment "BED file\nA BED file containing regions that should be excluded\nfrom all analyses. Currently this works by rejecting\ngenomic chunks that happen to overlap an entry.\nConsequently, for BAM files, if a read partially\noverlaps a blacklisted region or a fragment spans over\nit, then the read/fragment might still be considered.\n(default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--blackListFileName"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--blackListFileName"
] ;
sld:type cwl:File , sld:null .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.bed_file.path.replace(/^.*[\\\\\\/]/, '') + \"_sorted\")" ] ;
sld:type cwl:File .
cwl:outputBinding [ CommandOutputBinding:glob
"$(inputs.bed_graph.path.replace(/^.*[\\\\\\/]/, '').replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix)" ] ;
sld:type cwl:File .
cwl:source .
cwl:default ".bw" ;
sld:type xsd:string .
rdfs:comment "INT, -p INT\nNumber of processors to use. Type \"max/2\" to use half\nthe maximum number of processors or \"max\" to use all\navailable processors. (default: max/2)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--numberOfProcessors"
] ;
sld:type xsd:int , sld:null .
cwl:valueFrom "${return true}" .
cwl:valueFrom "RPKM" .
rdfs:comment "INT bp\nSize of the bins, in bases, for the output of the\nbigwig/bedgraph file. (default: 50)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--binSize"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--binSize"
] ;
sld:type xsd:int , sld:null .
cwl:valueFrom "subtract" .
rdfs:comment "BAM file to process " ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--bam"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--bam"
] ;
cwl:secondaryFiles [] ;
cwl:secondaryFiles [] ;
sld:type cwl:File .
cwl:source .
sld:type cwl:File .
rdfs:comment "--verbose \nSet to see processing messages. (default: False)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--verbose"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--verbose"
] ;
sld:type xsd:boolean , sld:null .
rdfs:comment "Adds a UCSC/Genome-Browser track line definition in the first line of the output.\n- See here for more details about track line definition:\nhttp://genome.ucsc.edu/goldenPath/help/bedgraph.html\n- NOTE: When adding a trackline definition, the output BedGraph can be easily\nuploaded to the Genome Browser as a custom track,\nBUT CAN NOT be converted into a BigWig file (w/o removing the first line).\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-trackline"
] ;
sld:type xsd:boolean , sld:null .
a cwl:CommandLineTool ;
rdfs:comment "Tool: bedtools genomecov (aka genomeCoverageBed)\nVersion: v2.25.0\nSummary: Compute the coverage of a feature file among a genome.\n\nUsage: bedtools genomecov [OPTIONS] -i -g \n\nOptions: \n\t-ibam\t\tThe input file is in BAM format.\n\t\t\tNote: BAM _must_ be sorted by position\n\n\t-d\t\tReport the depth at each genome position (with one-based coordinates).\n\t\t\tDefault behavior is to report a histogram.\n\n\t-dz\t\tReport the depth at each genome position (with zero-based coordinates).\n\t\t\tReports only non-zero positions.\n\t\t\tDefault behavior is to report a histogram.\n\n\t-bg\t\tReport depth in BedGraph format. For details, see:\n\t\t\tgenome.ucsc.edu/goldenPath/help/bedgraph.html\n\n\t-bga\t\tReport depth in BedGraph format, as above (-bg).\n\t\t\tHowever with this option, regions with zero \n\t\t\tcoverage are also reported. This allows one to\n\t\t\tquickly extract all regions of a genome with 0 \n\t\t\tcoverage by applying: \"grep -w 0$\" to the output.\n\n\t-split\t\tTreat \"split\" BAM or BED12 entries as distinct BED intervals.\n\t\t\twhen computing coverage.\n\t\t\tFor BAM files, this uses the CIGAR \"N\" and \"D\" operations \n\t\t\tto infer the blocks for computing coverage.\n\t\t\tFor BED12 files, this uses the BlockCount, BlockStarts, and BlockEnds\n\t\t\tfields (i.e., columns 10,11,12).\n\n\t-strand\t\tCalculate coverage of intervals from a specific strand.\n\t\t\tWith BED files, requires at least 6 columns (strand is column 6). \n\t\t\t- (STRING): can be + or -\n\n\t-5\t\tCalculate coverage of 5\" positions (instead of entire interval).\n\n\t-3\t\tCalculate coverage of 3\" positions (instead of entire interval).\n\n\t-max\t\tCombine all positions with a depth >= max into\n\t\t\ta single bin in the histogram. Irrelevant\n\t\t\tfor -d and -bedGraph\n\t\t\t- (INTEGER)\n\n\t-scale\t\tScale the coverage by a constant factor.\n\t\t\tEach coverage value is multiplied by this factor before being reported.\n\t\t\tUseful for normalizing coverage by, e.g., reads per million (RPM).\n\t\t\t- Default is 1.0; i.e., unscaled.\n\t\t\t- (FLOAT)\n\n\t-trackline\tAdds a UCSC/Genome-Browser track line definition in the first line of the output.\n\t\t\t- See here for more details about track line definition:\n\t\t\t http://genome.ucsc.edu/goldenPath/help/bedgraph.html\n\t\t\t- NOTE: When adding a trackline definition, the output BedGraph can be easily\n\t\t\t uploaded to the Genome Browser as a custom track,\n\t\t\t BUT CAN NOT be converted into a BigWig file (w/o removing the first line).\n\n\t-trackopts\tWrites additional track line definition parameters in the first line.\n\t\t\t- Example:\n\t\t\t -trackopts 'name=\"My Track\" visibility=2 color=255,30,30'\n\t\t\t Note the use of single-quotes if you have spaces in your parameters.\n\t\t\t- (TEXT)\n\nNotes: \n\t(1) The genome file should tab delimited and structured as follows:\n\t \n\n\tFor example, Human (hg19):\n\tchr1\t249250621\n\tchr2\t243199373\n\t...\n\tchr18_gl000207_random\t4262\n\n\t(2) The input BED (-i) file must be grouped by chromosome.\n\t A simple \"sort -k 1,1 > .sorted\" will suffice.\n\n\t(3) The input BAM (-ibam) file must be sorted by position.\n\t A \"samtools sort \" should suffice.\n\nTips: \n\tOne can use the UCSC Genome Browser's MySQL database to extract\n\tchromosome sizes. For example, H. sapiens:\n\n\tmysql --user=genome --host=genome-mysql.cse.ucsc.edu -A -e \\\n\t\"select chrom, size from hg19.chromInfo\" > hg19.genome" ;
ns1:original_cwlVersion "v1.0" ;
cwl:baseCommand ( "bedtools" "genomecov" ) ;
cwl:cwlVersion ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "dukegcb/bedtools"
] ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs , , , , , , , , , , , , , ;
cwl:outputs ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] ;
cwl:requirements [ a cwl:ShellCommandRequirement ] ;
cwl:stdout "$(inputs.ibam.path.replace(/^.*[\\\\\\/]/, '').replace(/\\.[^/.]+$/, '') + '.bdg')" .
rdfs:comment " show program's version number and exit\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--version"
] ;
sld:type xsd:boolean , sld:null .
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
rdfs:comment "A list of space-delimited chromosome names containing those chromosomes\nthat should be excluded for computing the normalization. This is useful\nwhen considering samples with unequal coverage across chromosomes, like\nmale samples. An usage examples is –ignoreForNormalization chrX chrM.\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--ignoreForNormalization"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--ignoreForNormalization"
] ;
sld:type sld:null ;
sld:type [ sld:items xsd:string ;
sld:type sld:array
] ;
sld:type [ sld:items xsd:string ;
sld:type sld:array
] .
rdfs:comment "The smooth length defines a window, larger than the\nbinSize, to average the number of reads. For example,\nif the –binSize is set to 20 and the –smoothLength is\nset to 60, then, for each bin, the average of the bin\nand its left and right neighbors is considered.\nAny value smaller than –binSize will be ignored and\nno smoothing will be applied. (default: 1.0)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--scaleFactor"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--scaleFactor"
] ;
sld:type xsd:float , sld:null .
cwl:valueFrom ".rpkm.bw" .
rdfs:comment "INT \nInclude reads based on the SAM flag. For example, to\nget only reads that are the first mate, use a flag of\n64. This is useful to count properly paired reads only\nonce, as otherwise the second mate will be also\nconsidered for the coverage. (default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--samFlagInclude"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--samFlagInclude"
] ;
sld:type xsd:int , sld:null .
cwl:valueFrom "${return 200}" .
rdfs:comment "\tReport depth in BedGraph format, as above (-bg).\nHowever with this option, regions with zero\ncoverage are also reported. This allows one to\nquickly extract all regions of a genome with 0\ncoverage by applying: \"grep -w 0$\" to the output.\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-bga"
] ;
sld:type xsd:boolean , sld:null .
cwl:valueFrom ".rpkm.bw" .
cwl:outputBinding [ CommandOutputBinding:glob
"${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".compared.bdg\"; return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".compared.bw\"; }" ] ;
sld:type cwl:File .
cwl:valueFrom "bigwig" .
a cwl:CommandLineTool ;
rdfs:comment "bedSort - Sort a .bed file by chrom,chromStart\nusage:\n bedSort in.bed out.bed\nin.bed and out.bed may be the same.\n" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( [ CommandLineBinding:position 2 ;
cwl:valueFrom "$(inputs.bed_file.path.replace(/^.*[\\\\\\/]/, '') + \"_sorted\")"
]
) ;
cwl:baseCommand ( "bedSort" ) ;
cwl:cwlVersion ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "dleehr/docker-hubutils"
] ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:inputs ;
cwl:outputs ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] .
cwl:in , , ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) .
rdfs:comment "\tTreat \"split\" BAM or BED12 entries as distinct BED intervals.\nwhen computing coverage.\nFor BAM files, this uses the CIGAR \"N\" and \"D\" operations\nto infer the blocks for computing coverage.\nFor BED12 files, this uses the BlockCount, BlockStarts, and BlockEnds\nfields (i.e., columns 10,11,12).\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-split"
] ;
sld:type xsd:boolean , sld:null .
rdfs:comment "\tCalculate coverage of intervals from a specific strand.\nWith BED files, requires at least 6 columns (strand is column 6).\n- (STRING): can be + or -\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-strand"
] ;
sld:type xsd:string , sld:null .
cwl:source .
rdfs:comment "\tThe input file is in BAM format.\nNote: BAM _must_ be sorted by position\n" ;
cwl:inputBinding [ CommandLineBinding:position 2 ;
CommandLineBinding:prefix "-ibam"
] ;
sld:type cwl:File .
cwl:valueFrom "bigwig" .
a cwl:Workflow ;
rdfs:comment "ChIP-seq - Quantification - samples: treatment and control" ;
ns1:original_cwlVersion "v1.0" ;
Workflow:steps , , , , , ;
cwl:cwlVersion ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs , , , ;
cwl:outputs , , , ;
cwl:requirements [ a cwl:ScatterFeatureRequirement ] ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] ;
cwl:requirements [ a cwl:StepInputExpressionRequirement ] .
rdfs:comment "CHR:START:END\nRegion of the genome to limit the operation to - this\nis useful when testing parameters to reduce the\ncomputing time. The format is chr:start:end, for\nexample --region chr10 or --region\nchr10:456700:891000. (default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--region"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--region"
] ;
sld:type xsd:string , sld:null .
a cwl:CommandLineTool ;
rdfs:comment "usage: bigwigCompare [-h] [--version] [--binSize INT bp]\n [--region CHR:START:END] [--numberOfProcessors INT]\n [--verbose] --outFileName FILENAME\n [--outFileFormat {bigwig,bedgraph}] --bigwig1 Bigwig file\n --bigwig2 Bigwig file [--scaleFactors SCALEFACTORS]\n [--pseudocount PSEUDOCOUNT]\n [--ratio {log2,ratio,subtract,add,reciprocal_ratio}]\n [--skipNonCoveredRegions]\nThis tool compares two bigWig files based on the number of mapped reads. To compare the bigWig files, the genome is partitioned into bins of equal size, then the number of reads found in each BAM file are counted per bin and finally a summary value is reported. This value can be the ratio of the number of readsper bin, the log2 of the ratio, the sum or the difference. optional arguments:\n -h, --help show this help message and exit\n --bigwig1 Bigwig file, -b1 Bigwig file\n Bigwig file 1. Usually the file for the treatment.\n (default: None)\n --bigwig2 Bigwig file, -b2 Bigwig file\n Bigwig file 2. Usually the file for the control.\n (default: None)\n --scaleFactors SCALEFACTORS\n Set this parameter to multipy the bigwig values by a\n constant. The format is scaleFactor1:scaleFactor2. For\n example 0.7:1 to scale the first bigwig file by 0.7\n while not scaling the second bigwig file (default:\n None)\n --pseudocount PSEUDOCOUNT\n small number to avoid x/0. Only useful when ratio =\n log2 or ratio (default: 1)\n --ratio {log2,ratio,subtract,add,reciprocal_ratio}\n The default is to compute the log2(ratio) between the\n two samples. The reciprocal ratio returns the the\n negative of the inverse of the ratio if the ratio is\n less than 0. The resulting values are interpreted as\n negative fold changes. Other possible operations are :\n simple ratio, subtraction, sum (default: log2)\n --skipNonCoveredRegions, --skipNAs\n This parameter determines if non-covered regions\n (regions without a score) in the bigWig files should\n be skipped. The default is to treat those regions as\n having a value of zero. The decision to skip non-\n covered regions depends on the interpretation of the\n data. Non-covered regions in a bigWig file may\n represent repetitive regions that should be skipped.\n Alternatively, the interpretation of non-covered\n regions as zeros may be wrong and this option should\n be used (default: False)\nOptional arguments:\n --version show program's version number and exit\n --binSize INT bp, -bs INT bp\n Size of the bins, in bases, for the output of the\n bigwig/bedgraph file. (default: 50)\n --region CHR:START:END, -r CHR:START:END\n Region of the genome to limit the operation to - this\n is useful when testing parameters to reduce the\n computing time. The format is chr:start:end, for\n example --region chr10 or --region\n chr10:456700:891000. (default: None)\n --numberOfProcessors INT, -p INT\n Number of processors to use. Type \"max/2\" to use half\n the maximum number of processors or \"max\" to use all\n available processors. (default: max/2)\n --verbose, -v Set to see processing messages. (default: False)\nOutput:\n --outFileName FILENAME, -o FILENAME\n Output file name. (default: None)\n --outFileFormat {bigwig,bedgraph}, -of {bigwig,bedgraph}\n Output file type. Either \"bigwig\" or \"bedgraph\".\n (default: bigwig)" ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( [ CommandLineBinding:position 3 ;
CommandLineBinding:prefix "--outFileName" ;
cwl:valueFrom "${ if (inputs.outFileName) return inputs.outFileName; if (inputs.output_suffix) return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix; if (inputs.outFileFormat == \"bedgraph\") return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".compared.bdg\"; return inputs.bigwig1.path.replace(/^.*[\\\\\\/]/, \"\").replace(/\\.[^/.]+$/, \"\") + \".compared.bw\"; }"
]
) ;
cwl:baseCommand ( "bigwigCompare" ) ;
cwl:cwlVersion ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "reddylab/deeptools:3.0.1"
] ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:inputs , , , , , , , , , , , , ;
cwl:outputs ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] .
rdfs:comment "{forward,reverse}\nSelects RNA-seq reads (single-end or paired-end) in\nthe given strand. (default: None)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--filterRNAstrand"
] ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--filterRNAstrand"
] ;
sld:type xsd:string , sld:null .
rdfs:comment "Raw reads bigWig (signal) files" ;
cwl:outputSource ;
sld:type [ sld:items cwl:File ;
sld:type sld:array
] .
cwl:source .
cwl:source .
cwl:in , , , , , , ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) ;
cwl:scatterMethod ScatterMethod:dotproduct .
rdfs:comment "\tCombine all positions with a depth >= max into\na single bin in the histogram. Irrelevant\nfor -d and -bedGraph\n- (INTEGER)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-max"
] ;
sld:type xsd:int , sld:null .
cwl:in , , , , , , ;
cwl:out ;
cwl:run ;
cwl:scatter ( ) .
rdfs:comment "The default is to output the log2ratio of the two\nsamples. The reciprocal ratio returns the the negative\nof the inverse of the ratio if the ratio is less than\n0. The resulting values are interpreted as negative\nfold changes. *NOTE*: Only with --operation subtract\ncan --normalizeUsing RPGC or --normalizeUsing RPKM be\nused. Instead of performing a computation using both\nfiles, the scaled signal can alternatively be output\nfor the first or second file using the '--operation\nfirst' or '--operation second' (default: log2)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "--operation"
] ;
sld:type xsd:string , sld:null .
rdfs:comment "Writes additional track line definition parameters in the first line.\n- Example:\n-trackopts 'name=\"My Track\" visibility=2 color=255,30,30'\nNote the use of single-quotes if you have spaces in your parameters.\n- (TEXT)\n" ;
cwl:inputBinding [ CommandLineBinding:position 1 ;
CommandLineBinding:prefix "-trackopts"
] ;
sld:type xsd:string , sld:null .
a cwl:CommandLineTool ;
rdfs:comment "Tool: bedGraphToBigWig v 4 - Convert a bedGraph file to bigWig format." ;
ns1:original_cwlVersion "v1.0" ;
cwl:arguments ( [ CommandLineBinding:position 3 ;
cwl:valueFrom "$(inputs.bed_graph.path.replace(/^.*[\\\\\\/]/, '').replace(/\\.[^/.]+$/, \"\") + inputs.output_suffix)"
]
) ;
cwl:baseCommand ( "bedGraphToBigWig" ) ;
cwl:cwlVersion ;
cwl:hints [ a ;
LoadListingRequirement:loadListing
"deep_listing"
] ;
cwl:hints [ a cwl:DockerRequirement ;
DockerRequirement:dockerPull "dukegcb/bedgraphtobigwig"
] ;
cwl:hints [ a ;
NetworkAccess:networkAccess true
] ;
cwl:inputs , , ;
cwl:outputs ;
cwl:requirements [ a cwl:InlineJavascriptRequirement ] .