Explore Workflows

Sequence reads are first cleaned from adapters and transformed into fully bisulfite-converted forward (C->T) and reverse read (G->A conversion of the forward strand) versions, before they are aligned to similarly converted versions of the genome (also C->T and G->A converted). Sequence reads that produce a unique best alignment from the four alignment processes against the bisulfite genomes (which are running in parallel) are then compared to the normal genomic sequence and the methylation state of all cytosine positions in the read is inferred. A read is considered to align uniquely if an alignment has a unique best alignment score (as reported by the AS:i field). If a read produces several alignments with the same number of mismatches or with the same alignment score (AS:i field), a read (or a read-pair) is discarded altogether. On the next step we extract the methylation call for every single C analysed. The position of every single C will be written out to a new output file, depending on its context (CpG, CHG or CHH), whereby methylated Cs will be labelled as forward reads (+), non-methylated Cs as reverse reads (-). The output of the methylation extractor is then transformed into a bedGraph and coverage file. The bedGraph counts output is then used to generate a genome-wide cytosine report which reports the number on every single CpG (optionally every single cytosine) in the genome, irrespective of whether it was covered by any reads or not. As this type of report is informative for cytosines on both strands the output may be fairly large (~46mn CpG positions or >1.2bn total cytosine positions in the human genome).

https://github.com/datirium/workflows.git

Path: workflows/bismark-methylation-pe.cwl

Branch/Commit ID: 93b844a80f4008cc973ea9b5efedaff32a343895

https://github.com/common-workflow-language/cwltool.git

Path: cwltool/schemas/v1.0/v1.0/count-lines9-wf.cwl

Branch/Commit ID: 4700fbee9a5a3271eef8bc9ee595619d0720431b

https://github.com/ncbi/pgap.git

Path: task_types/tt_cluster_and_qdump.cwl

Branch/Commit ID: 424a01693259a75641dc249d553235aa38a6ce23

https://github.com/genome/analysis-workflows.git

Path: definitions/subworkflows/downsampled_alignment.cwl

Branch/Commit ID: a23f42ef49c10a588fd35a3afaad5de03e253533

https://github.com/datirium/workflows.git

Path: subworkflows/allele-process-reference.cwl

Branch/Commit ID: 09c08f858e91c4cfd387067f07445722ac7e18aa

https://github.com/genome/analysis-workflows.git

Path: definitions/subworkflows/molecular_alignment.cwl

Branch/Commit ID: 93656ed6582073e434eab168c610625a835dce37

https://github.com/common-workflow-language/common-workflow-language.git

Path: v1.0/v1.0/count-lines8-wf-noET.cwl

Branch/Commit ID: 9a23706ec061c5d2c02ff60238d218aadf0b5db9

Inspect provided directory and return filenames. Generate a new directory and return it (including content).

https://github.com/common-workflow-language/cwltool.git

Path: tests/wf/directory.cwl

Branch/Commit ID: 63f539ba60e91f0cb3ce7cda2c5da5c65525c375

Slightly changed original [BioWardrobe's](https://biowardrobe.com) [PubMed ID:26248465](https://www.ncbi.nlm.nih.gov/pubmed/26248465) **RNA-Seq** basic analysis for **strand specific pair-end** experiment. An additional steps were added to map data to mitochondrial chromosome only and then merge the output. Experiment files in [FASTQ](http://maq.sourceforge.net/fastq.shtml) format either compressed or not can be used. Current workflow should be used only with the pair-end strand specific RNA-Seq data. It performs the following steps: 1. `STAR` to align reads from input FASTQ file according to the predefined reference indices; generate unsorted BAM file and alignment statistics file 2. `fastx_quality_stats` to analyze input FASTQ file and generate quality statistics file 3. `samtools sort` to generate coordinate sorted BAM(+BAI) file pair from the unsorted BAM file obtained on the step 1 (after running STAR) 5. Generate BigWig file on the base of sorted BAM file 6. Map input FASTQ file to predefined rRNA reference indices using Bowtie to define the level of rRNA contamination; export resulted statistics to file 7. Calculate isoform expression level for the sorted BAM file and GTF/TAB annotation file using `GEEP` reads-counting utility; export results to file

https://github.com/datirium/workflows.git

Path: workflows/rnaseq-pe-dutp-mitochondrial.cwl

Branch/Commit ID: 4106b7dc96e968db291b7a61ecd1641aa3b3dd6d

https://github.com/junyussh/RunningWorkflows-on-the-GoogleCloud.git

Path: CWL-RNAseq/CWL-RNAseq.cwl

Branch/Commit ID: 6ce4a66ad8847f3214e3ec50194c2dc202e23109