Conversation
Some fixes
…omap is used with PE reads
CHANGELOG.md
Outdated
| The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/) | ||
| and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html). | ||
|
|
||
| ## v2.1.0 - [date] |
There was a problem hiding this comment.
Update date when releasing!
There was a problem hiding this comment.
I also missed in mine.....how do you do this? just direct commit to dev ?
| label 'process_low' | ||
|
|
||
| conda (params.enable_conda ? "conda-forge::sed=4.7" : null) | ||
| conda "conda-forge::sed=4.7" |
There was a problem hiding this comment.
should the conda be reading from a environment.yml file ?
There was a problem hiding this comment.
Yes, you are right but this now will need to refactor all the local modules and I want to release as it is long time we have not done it. If you don't mind I opened an issue to do it for next release, see here
There was a problem hiding this comment.
why have local versions of these nf-core modules ?
There was a problem hiding this comment.
The reason is this: the indexes produced by the nf-core modules will be incompatible with iGenomes.
nextflow.config
Outdated
| version = '2.0.0' | ||
| nextflowVersion = '!>=23.04.0' | ||
| version = '2.1.0' | ||
| doi = '' |
| @@ -100,45 +105,53 @@ | |||
| }, | |||
There was a problem hiding this comment.
For clarity I'd write for both gtf and gff
"Either this parameter or gtf is mandatory if --genome is not specified."
bin/check_samplesheet.py
Outdated
| print_error( | ||
| f"Control identifier and replicate has to match a provided sample identifier and replicate!", | ||
| "Control", | ||
| val[4], |
There was a problem hiding this comment.
| val[4], | |
| val[-1], |
There was a problem hiding this comment.
I believe val[5] should also be fine, but val[4] refers to the antibody
README.md
Outdated
| These scripts were originally written by Chuan Wang ([@chuan-wang](https://github.com/chuan-wang)) and Phil Ewels ([@ewels](https://github.com/ewels)) for use at the [National Genomics Infrastructure](https://portal.scilifelab.se/genomics/) at [SciLifeLab](http://www.scilifelab.se/) in Stockholm, Sweden. The pipeline was re-implemented by Harshil Patel ([@drpatelh](https://github.com/drpatelh)) from [Seqera Labs, Spain](https://seqera.io/) and converted to Nextflow DSL2 by Jose Espinosa-Carrasco ([@JoseEspinosa](https://github.com/JoseEspinosa)) from [The Comparative Bioinformatics Group](https://www.crg.eu/en/cedric_notredame) at [The Centre for Genomic Regulation, Spain](https://www.crg.eu/). | ||
|
|
||
| Many thanks to others who have helped out and contributed along the way too, including (but not limited to): [@apeltzer](https://github.com/apeltzer), [@bc2zb](https://github.com/bc2zb), [@crickbabs](https://github.com/crickbabs), [@drejom](https://github.com/drejom), [@houghtos](https://github.com/houghtos), [@KevinMenden](https://github.com/KevinMenden), [@mashehu](https://github.com/mashehu), [@pditommaso](https://github.com/pditommaso), [@Rotholandus](https://github.com/Rotholandus), [@sofiahaglund](https://github.com/sofiahaglund), [@tiagochst](https://github.com/tiagochst) and [@winni2k](https://github.com/winni2k). | ||
| The pipeline workflow diagram was designe by Sarah Guinchard ([@G-Sarah](https://github.com/G-Sarah)). |
There was a problem hiding this comment.
| The pipeline workflow diagram was designe by Sarah Guinchard ([@G-Sarah](https://github.com/G-Sarah)). | |
| The pipeline workflow diagram was designed by Sarah Guinchard ([@G-Sarah](https://github.com/G-Sarah)). |
CHANGELOG.md
Outdated
| | `multiqc` | 1.13 | 1.23 | | ||
| | `picard` | 2.27.4 | 3.2.0 | | ||
| | `samtools` | 1.15.1 | 1.20 | | ||
| | `samtools` | 1.15.1 | 1.20 | |
There was a problem hiding this comment.
| | `samtools` | 1.15.1 | 1.20 | |
docs/usage.md
Outdated
| ### Multiple replicates | ||
|
|
||
| The `sample` identifiers have to be the same when you have re-sequenced the same sample more than once e.g. to increase sequencing depth. The pipeline will perform the alignments in parallel, and subsequently merge them before further analysis. Below is an example where the samples called `WT_BCATENIN_IP_REP2` and `WT_INPUT_REP2` have been re-sequenced multiple times: | ||
| The `sample` identifier should be identical when you have multiple replicates from the same experimental group, just increment the `replicate` identifier appropriately. The first replicate value for any given experimental group must be 1. |
There was a problem hiding this comment.
| The `sample` identifier should be identical when you have multiple replicates from the same experimental group, just increment the `replicate` identifier appropriately. The first replicate value for any given experimental group must be 1. | |
| The `sample` identifier should be identical when you have multiple replicates from the same experimental group; just increment the `replicate` identifier appropriately. The first replicate value for any given experimental group must be 1. |
docs/usage.md
Outdated
| It is a good idea to specify a pipeline version when running the pipeline on your data. This ensures that a specific version of the pipeline code and software are used when you run your pipeline. If you keep using the same tag, you'll be running the same version of the pipeline, even if there have been changes to the code since. | ||
|
|
||
| First, go to the [nf-core/chipseq releases page](https://github.com/nf-core/chipseq/releases) and find the latest version number - numeric only (eg. `1.2.2`). Then specify this when running the pipeline with `-r` (one hyphen) - eg. `-r 1.2.2`. | ||
| First, go to the [nf-core/chipseq releases page](https://github.com/nf-core/chipseq/releases) and find the latest pipeline version - numeric only (eg. `2.0.0`). Then specify this when running the pipeline with `-r` (one hyphen) - eg. `-r 2.0.0`. Of course, you can switch to another version by changing the number after the `-r` flag. |
There was a problem hiding this comment.
| First, go to the [nf-core/chipseq releases page](https://github.com/nf-core/chipseq/releases) and find the latest pipeline version - numeric only (eg. `2.0.0`). Then specify this when running the pipeline with `-r` (one hyphen) - eg. `-r 2.0.0`. Of course, you can switch to another version by changing the number after the `-r` flag. | |
| First, go to the [nf-core/chipseq releases page](https://github.com/nf-core/chipseq/releases) and find the latest pipeline version - numeric only (eg. `2.1.0`). Then specify this when running the pipeline with `-r` (one hyphen) - eg. `-r 2.1.0`. Of course, you can switch to another version by changing the number after the `-r` flag. |
| ch_peak_count_header = file("$projectDir/assets/multiqc/peak_count_header.txt", checkIfExists: true) | ||
| ch_frip_score_header = file("$projectDir/assets/multiqc/frip_score_header.txt", checkIfExists: true) | ||
| ch_peak_annotation_header = file("$projectDir/assets/multiqc/peak_annotation_header.txt", checkIfExists: true) | ||
| ch_deseq2_pca_header = Channel.value(file("$projectDir/assets/multiqc/deseq2_pca_header.txt", checkIfExists: true)) |
There was a problem hiding this comment.
I can't spot a difference in the usage of this header vs e.g. the ch_peak_annotation_header. What's the reason for defining some as paths and some as value channels?
JoseEspinosa
left a comment
There was a problem hiding this comment.
I also missed in mine.....how do you do this? just direct commit to dev ?
I will open a small PR to dev instead, ping me if you need a review
Apply review suggestions
Fix how replicates are set for downstream analysis
No description provided.