Recent Releases of smallvariants
smallvariants - v1.11.0 - Generous Ghent
New features
- Added a new option to the samplesheet (
msi). This option takes a boolean to indicate whether or not the MSI status of the sample should be checked. A baseline should be provided using the--msi_baselineparameter to run MSI calling. - The Multiqc module now generates branded reports.
Fixes
- Fixed PED file writing on cloud platforms.
- Fixed publishing of the MultiQC data directory throwing a
NullPointerException. - QUAL scores of 256 are no longer converted to
.(due to a VCFanno bug).
- Nextflow
Published by nvnieuwk 6 months ago
smallvariants - v1.10.1 - Mighty Mechelen
Fixes
- Fix a bug where the BED processing failed because of unknown chromosomes.
- Nextflow
Published by nvnieuwk 9 months ago
smallvariants - v1.10.0 - Melodic Mechelen
Changes
- Merged the following processes to improve efficiency of the pipeline:
- VCF index creation modules on output VCFs have been merged into the processes that created these VCFs
- The filter modules for
--filterhave been merged - BED filtering and intersecting with Regions Of Interest have been merged
- Added copgt_joint config profile
- Migrate to the third preview of the workflow output definitions
- Bumped minimal Nextflow version to 25.04.0
- Removed
--squash-ploidyfrom the RTG vcfeval process - Small stability fixes for runnning on the nomad cluster
- Added the merged CRAM files to the output of the pipeline. This can be skipped using the
--skip_merged_cram_outputparameter. - Added the filter parameter to the copgt config profiles
- Added mosdepth QC to the multiQC report
- Nextflow
Published by nvnieuwk 9 months ago
smallvariants - v1.9.3 - Nifty Nieuwkerke
Changes
- Fix db postprocess in vcf2db module
- Nextflow
Published by nvnieuwk about 1 year ago
smallvariants - v1.9.2 - New Nieuwkerke
Changes
- The pipeline now also outputs the per-base coverage BED file generated by mosdepth per sample
- Add the
--disable_hc_dict_validationparameter that disables sequence dictionary validation in HaplotypeCaller
- Nextflow
Published by nvnieuwk about 1 year ago
smallvariants - v1.9.1 - Nice Nieuwkerke
Changes
- Updated the ROI file in the WES profile from v5 to v6
- Nextflow
Published by nvnieuwk about 1 year ago
smallvariants - v1.9.0 - Neighborly Nieuwkerke
New features
- Added the
--min_callable_coverageparameter to state what the lowest coverage should for a region to be classified as callable. - Added the
elprepcaller as an alternative to the haplotypecaller. - Added full unit tests for all parts that were missing tests.
Changes
- Added the
--squash-ploidyargument to the RTG vcfeval process. - Update to nf-core v3.0.1
- Completely reworked the output directory structure to a more sensible structure. The pipeline can now be run on the same output directory every time and will incrementally add files to the correct family folder. See the output documentation for more info.
- Migrated to the new workflow output definitions.
- Bumped the minimal Nextflow version to 24.10.0.
- Added the somalier reports to the multiQC report.
- Removed the
--output_suffixparameter - Added some missing required parameters to the
WESandseqplorerprofiles
Fixes
- Validation of all samples now uses an intersect of the golden truth BED files with the BED file used to call the variants. This should fix the WES validation which was broken until this point.
- A couple of small fixes to the vardict flow.
- Only use the standard chromosomes for UPDio analysis.
- Reduced the resources given to some GATK4 modules
- VCF2DB now uses a seqera container to fix some issues when running it in nomad
- Dots in sample and family names are now converted to an underscore automatically.
- Nextflow
Published by nvnieuwk about 1 year ago
smallvariants - v1.8.2 - Outstanding Oostkamp
Fixes
- Fixed some issues where indices were not created
- Updated the docs
- Nextflow
Published by nvnieuwk over 1 year ago
smallvariants - v1.8.1 - Open Oostkamp
Fixes
- Revert VEP version to v105
- Nextflow
Published by nvnieuwk over 1 year ago
smallvariants - v1.8.0 - Optimistic Oostkamp
New features
- Added
watchpathfunctionality to the pipeline. Add thewatch:prefix to a file basename in the samplesheet and the pipeline will automatically wait for the file to be created in the--watchdirdirectory (the lookup happens recursively)
Changes
- Bumped the minimal support nextflow version to
24.04.0 - Bumped all modules to the newest versions
- The pipeline now also outputs
csiindices - Rename the
masterbranch tomain - Low coverage regions (regions with less than 5 reads) are no longer considered for variant calling
Refactors
- Updated the pipeline to the new linting guidelines
- Removed
check_maxin favor ofresourceLimits
- Nextflow
Published by nvnieuwk over 1 year ago
smallvariants - v1.7.0 - Tolerant Tongeren
Fixes
- Automap analysis should now give the correct output files for individuals.
Changes
- Haplotypecaller will not perform phasing by default now. This can still be turned back on using the
--hc_phasingparameter. - Removed the WES and WGS profiles.
- Nextflow
Published by nvnieuwk over 1 year ago
smallvariants - v1.6.0 - Offbeat Ostend
New features
- Added UPDio for Uniparental Disomy detection in family samples. This introduces the
--updioparameter to turn on this detection and--updio_common_cnvsto supply a common CNVs file to UPDio. The family needs to contain at least one child with its mother and father. - Added docs built with MkDocs. See the documentation site here.
- Added AutoMap to find regions of homozygosity from human samples. This introduces the
--automapparameter to turn on this feature and the--automap_repeats,--automap_paneland--automap_panel_nameparameters to configure AutoMap (see the parameters docs for more information)
Changes
- Updated all tests to use snapshots instead.
- Made the pipeline pluggable to enable the use of it in a meta pipeline.
- Nextflow
Published by nvnieuwk almost 2 years ago
smallvariants - v1.5.1 - Great Geraardsbergen
v1.5.1 - Great Geraardsbergen - [March 7 2024]
Fixes
- Fixed an issue with
igenomespaths not being casted correctly to their corresponding parameter
- Nextflow
Published by nvnieuwk almost 2 years ago
smallvariants - v1.5.0 - Amazing Antwerp
Improvements
- Updated to the nf-core template v2.13.0
- Updated all GATK modules to 4.5.0.0
- Moved the pipeline from https://github.com/CenterForMedicalGeneticsGhent/nf-cmgg-germline to https://github.com/nf-cmgg/germline
Changes
- VCF files created with
haplotypecallernow have thehaplotypecallertag in the filename instead ofgatk4-jointto keep naming consistent
- Nextflow
Published by nvnieuwk almost 2 years ago
smallvariants - v1.4.2 - Vibrant Veurne
Fixes
- Set the default ensembl VEP version to 105.0 instead of using dynamic container fetching
- Nextflow
Published by nvnieuwk about 2 years ago
smallvariants - v1.4.1 - Lively Leuven
New Features
- Added the
--output_suffixparameter to add a custom suffix to the basename of the output files. - Implemented files for the alphamissense plugin of VEP.
- Added the
--only_passparameter to only output variants that have thePASSflag in the FILTER column. (This is only applied when--filteris also given) - Added the
--keep_alt_contigsparameter. This will tell the pipeline to not filter out the alternate contigs, which will now be done by default. - Add dbsnp Ids to VCFs coming from vardict. This will be done automatically if a dbsnp VCF is given to the pipeline through the
--dbsnpparameter.
Improvements
- Updated the seqplorer profile so that the output filenames are correct for easy import
- Changed the separator in
--vcfanno_resourcesto;instead of,to allow commas in glob patterns. - Removed the reheader step from the vardict subworkflow and added a simple sed substitution to the vardictjava module
vcf2dbnow uses a python 2 environment to increase it's stability
- Nextflow
Published by nvnieuwk about 2 years ago
smallvariants - v1.4.0 - Kingly Kortrijk
New Features
- Added the
--callersparameter to specify the variant caller to use. Currently onlyhaplotypecallerandvardictare supported. - Added the
vardictvariant caller. - Added the
--vardict_min_afparameter to specify the minimum allele frequency forvardict. This option is also available in the samplesheet asvardict_min_afto set it dynamically per sample. - Added the
--output_genomicsdboption to specify whether a GenomicsDB should be outputted or not. This will betruewhen usingonly_merge. - Added
--normalizeoptions for decomposing and normalizing of variants after calling and genotyping. - Added
WGS,WES,SeqCap,HyperCapandseqplorerprofiles that can be used to set the default parameters for these types of runs.
Improvements
- Refactored the pipeline to accomodate future additions of variant callers and genotypers
- Removed a lot of unnecessary bloat
- Improved GenomicsDBImport (can now be multithreaded and runs a lot faster). This will make very big runs more possible.
- Changed
coverage_fasttomosdepth_slow, reversing the effect of the parameter. By default mosdepth will now be run with--fast-mode. This can be disabled using the newmosdepth_slowparameter. - Automatically merge the regions that are within 150 bps of eachother for the variant calling. This way it's ensured that indel calling happens correctly.
Fixes
- Fixed an issue with the outputting of the validation PNG files, now all three types of PNGs are outputted.
- Fixed a small issue where VCFs without a sample created by the callers could not be used by
bcftools concat, these files will now be filtered from the input of the command. - Removed the
--maxentscanparameter because this file is automatically present in the container
- Nextflow
Published by nvnieuwk about 2 years ago
smallvariants - v1.3.0 - Happy Hasselt
New Features
- Added the
--only_callparameter. Specifying this parameter tells the pipeline to only do variant calling and skip all post-processing. This will only output the GVCFs and files created to help variant calling. - The samplesheet is now also in the output folder.
- Added an option
--only_mergeto tell the pipeline to create genomicsdbs and stop running there - Get regions from the GVCF instead of CRAM for joint genotyping. This removes the need to supply a CRAM file when a GVCF file has been used as input.
Improvements
- Updated
nf-validationto v0.2.1. - Updated the samtools/merge tool to the nf-core version. This increases the efficiency and disk space usage of the tool.
Fixes
- Fixed an error where the truth VCFs caused a join error when the same sample was given multiple times
- Updated some outdated error messages
- Nextflow
Published by nvnieuwk over 2 years ago
smallvariants - v1.2.2 - Benign Brussels
Improvements
Changed the output directory structure to be more bcbio like
- Nextflow
Published by nvnieuwk over 2 years ago
smallvariants - v1.2.1 - Balanced Brussels
New Features
- Added support for the
nf-validationplugin. - Haplotypecaller dragen mode will be automatically disabled when not using a dragstr model.
Bug fixes
- Removed bedtools/jaccard
- Fixed some patterns in the parameter JSON schema (since they are actually used now)
- Fixed a breaking bug where mosdepth didn't output the callable regions (this makes v1.2.0 deprecated, please use v1.2.1 instead)
Improvements
- Genomicsdbs aren't scattered now, this increases the precision of the analyis by almost 3% at the cost of a bit longer runtimes
- Actually do the validation on the output VCFs now instead the freshly called GVCFs
- Improved the efficiency of the VEP run by scattering more efficiently on the amount of variants instead of the chromosomes
- Nextflow
Published by nvnieuwk over 2 years ago
smallvariants - v1.2.0 - Brave Brussels
WARNING: This version contains a fatal bug with Mosdepth which causes the pipeline to only call on low coverage regions. Please use the dev branch in the meantime.
New Features
- Added a
--coverage_fast <true/false>flag which can be used to run mosdepth in fast mode. This flag will also make sure that only the quantized bed from mosdepth is present in the output directory for each WGS individual, otherwise it will output everything - Added the possibility to give GVCF files as inputs and immediately go to the joint-genotyping. This is especially useful for the cases where several samples should be combined. This way the variant calling doesn't need to be re-run. Beware though that a CRAM file should still be given to generate the BED files used for the scatter/gathering. The new header names are
gvcfandtbiwheregvcfis used to give the GVCF andtbiis used to give its index. - Added
bedtools jaccardto the validation. - Added a Dockerfile which creates an image that is able to run a full pipeline run inside of it.
- Added better documentation
Improvements
- Updated the scattering again: it now follows this workflow:
- Sort and merge overlapping intervals of given ROI BED files (WES only)
- Create a BED file with callable regions using mosdepth
- Intersect the callable regions BED with the ROI BED (WES only)
- Split the resulting BED file (or the callable regions BED for WGS) into evenly sized BED files (amount is specified with
--scatter_count) - Run HaplotypeCaller in parallel using these regions
- Merge and sort the BED files of all individuals in a family
- Split the merged BED file into evenly sized BED files (amount is specified with
--scatter_counttimes the family size) - Run GenomicsDBImport and GenotypeGVCFs in parallel using these regions
- Updated the resource requirements of GenomicsDBImport and GenotypeGVCFs to be more efficient (and more cluster friendly)
- Removed ReblockGVCFs (this wasn't worth it and we save the raw GVCFs)
- Added
--merge_distance <integer>to decrease the amount of intervals passed to genomicsdbimport. Increase this value if GenomicsDBImport is running slow. - Renamed
--use_dragstr_modelto--dragstr.
Bug fixes
- Fixed a warning showing up when running with
--dragstr false - Add
--inferflag tosomalier relatewhen no PED file is given
- Nextflow
Published by nvnieuwk almost 3 years ago
smallvariants - v1.1.2 - Groovy Ghent
New features
- Added a parameter for setting the splitting depth threshold
--split_threshold FLOAT
Changes
- Change the default splitting threshold to 0.2 instead of 0.3
- Nextflow
Published by nvnieuwk almost 3 years ago
smallvariants - v1.1.1 - Golden Ghent
Changes
- Set the default of
--validatetofalse
Bug fixes
- Fixed a bug with ensembl VEP. Filenames of the alt contigs should now have a
_altsuffix instead of all alt contigs. - Added file-exist check to the
sdffile - Fixed the scattering when using alt contigs
- Nextflow
Published by nvnieuwk almost 3 years ago
smallvariants - v1.1.0 - Glorious Ghent
! WES support has been temporarily disabled, but is planned to be re-added in the next major release !
New Features
- Added support for samples that aren't part of a family. Just leave the
pedandfamily_idinput fields in the samplesheet empty for a sample to be treated like this. This sample will go through exactly the same workflow but will be emitted as a single-sample VCF. - Added
dumpfunctionality to lots of channels. - Added the
dbsnpoption toGATK HaplotypeCaller. use--dbsnpand--dbsnp_tbito supply these VCFs. - Added the
vcf_extract_somaliersubworkflow to the pipeline. This also creates PED files inferred from the input multi-sample VCF. - Added a validation subworkflow. All files that have a VCF in the
truth_vcfcolumn of the input samplesheet will be validated against this VCF. This can be turned off by supplying the--validate falseflag to the pipeline run.
Improvements
- Improved the scatter/gather logic. This is now done with
goleft indexsplitto define chunks of even coverage. The genotyping scattering now happens withbedtools makewindows. This creates chunks of even regions from the merged BED files for the family. By passing a padding of about 20 bps to the genotype tools, we make sure all variants on the edges of these regions are also genotyped. Duplicates are removed later when runningbcftools concat - Refactored a lot of the code to maintain the same style over the whole pipeline.
- Updated the minimum Nextlow version to
22.10.5to make sure S3 staging works perfectly. - The
post_processingsubworklow has been renamed to the better suitingjoint_genotypingsubworkflow.reblockgvcfhas been moved togermline_variant_callingand thefilterandreheaderinghas been moved to the main workflow. - Merging VCFs of the same family now happens with
GATK GenomicsDBImportinstead ofGATK MergeGVCFsorbcftools merge. This gives more reliable results. - Improved the handling of
vcfanno - The PED headers can now be added to all the output VCFs that are part of a family instead of only those that were given a PED file as input. The PED file used is created using
somalier relate. This feature can be turned on using the--add_ped trueargument. This doesn't happen by default.
Bug fixes
- Fixed some issues when both the
pedandfamily_idwere given for a sample. - Fixed the PED input for
rtgtools_pedfilter(-9isn't recognized as unknown by the tool. Now these will be automatically converted to0before this tool) - Fixed issues with the DBsnp index not being created correctly
- Fixed wrongly formed joins and added checks for mismatches and duplicates
- Nextflow
Published by nvnieuwk almost 3 years ago
smallvariants - v1.0.1 - Happy Hollebeke
Changes
- Upgraded to
nf-corev2.6 template
Fixes
- Fixed the
ensemblvepversion (was 104.1 before and is now 105.0) - Updated the label of
gatk4/calibratedragstrmodeltoprocess_highto match the requirements for bigger inputs
- Nextflow
Published by nvnieuwk over 3 years ago
smallvariants - v1.0.0 - Beautiful Bruges
The first full release of the pipeline!
- Nextflow
Published by nvnieuwk over 3 years ago