Tags: biowdl/RNA-seq
Tags
+ Default docker images for various tools have been updated. + Changes in Cromwell 48 made it impossible to use the wide array of inputs in our documentation (such as `RNAseq.sampleJobs.qc.Cutadapt.minimumLength`). Fixes have been made upstream and in the pipeline. From Cromwell 52 onwards these options will be available again. + Bam files are no longer indexed after alignment, saving compute time. + WDL files and imports zip packages are now provided each release to make running pipelines easier. + The output directory was simplified. All files related to a sample are now in the `samples/<sample_id>` directory. + Use the scatter-regions tool to replace biopet-scatterregions. + The pipeline was renamed from "pipeline" to "RNAseq". + Tasks were updated to contain the `time_minutes` runtime attribute and associated `timeMinutes` input, describing the maximum time the task will take to run. + Added a step to generate the STAR index, if neither a STAR nor a Hisat2 index is provided. + Document the use of cromwell's `final_workflow_outputs_dir` feature which makes the RNA-seq pipeline usable on all of Cromwell's supported backends. Users are encouraged to use this feature. `outputDir` references are removed from the documentation. + Make the MultiQC task suitable for use with a `final_workflow_outputs_dir` so it can be used on all of Cromwell's supported backends. + Added a picard markduplicates step after UMI deduplication. + Move common optional inputs to top-level workflow, so nested inputs are not required anymore for the majority of configurations. + The pipeline has been altered so it starts the variant calling jobs in a more efficient way. + Major bug fix: The --dont-use-soft-clipped-reads is now used on HaplotypeCaller in concordance with GATK best practices.
In version 3.0.0 the RNA-seq pipeline was brought up to date with the…
… GATK best
practices pipeline. Several errors in the variant calling part of the pipeline
were fixed.
+ UMI deduplication stats collection is now optional.
+ Add scatterSize option to centrally control the scatter size
+ Multisample VCFs are no longer generated for RNA-seq by default as this is
not GATK best practice. It can optionally be turned on again.
+ Add proper copyright headers to WDL files. So the free software license
is clear to end users who wish to adapt and modify.
+ Added UMI based deduplication as an optional step.
+ Major fix: Joint genotyping is not performed on RNA-seq reads. Instead, each
sample is genotyped seperately by the HaplotypeCaller.
+ Major fix: base recalibration is now applied after splitting the cigar reads.
previously this was done in an erroneous manner.
+ Structs are removed from the input.
+ Fixed a typo in the inputs: detectNovelTranscipts -> detectNovelTranscripts.
+ Removed a number of unused inputs ("dbsnp", "targetIntervals",
"ampliconIntervals" and "variantCalling") from the sample workflow.
+ Added input overview to docs.
+ Added WDL-AID to linting.
+ Update default htseq image version to 0.11.2
version 2.0.0 --------------------------- + Replace the bam-to-gvcf and jointgenotyping pipelines with a gatk-variantcalling pipeline. This reduces complexity, allows for more efficient running of the pipeline and uses less filesystem storage. + Update biowdl-input-converter tot the latest bugfix release (0.2.1) + Update STAR version to 2.7.3a + Simplify the pipeline so it uses much less subworkflows. This reduces the complexity for cromwell and reduces inefficiencies that are caused by waiting for the subworkflows to finish. It also makes configuring memory or cpu requirements for tasks in the workflow a lot easier, as these are not as deeply nested anymore.