Comments (11)
Done! nf-core/configs@f84bc30
Are you able to test it out with the latest version of the pipeline and let us know if everything is working as expected please?
from viralrecon.
I test it with the following configuration:
nextflow run nf-core/viralrecon
-profile sbc_sharc
-resume
-with-report
-with-trace
-with-timeline
--input '/fastdata/md1mbdx/Covid_and_ITS2_PA_PZ_260123/SARScov2/viralrecon/Sars_cov_2.csv'
--outdir results
--platform nanopore
--artic_minion_caller medaka
--artic_minion_medaka_model r941_min_hac_g507
--genome 'MN908947.3'
--primer_set_version 5.3.2
--fastq_dir '/fastdata/md1mbdx/Covid_and_ITS2_PA_PZ_260123/SARScov2/input/'
-c viralrecon_resources.conf
Which is the same configuration I always use, the only difference is between 4.1 and 5.3.2 but I get this error:
ERROR: Validation of pipeline parameters failed!
* --primer_set_version: expected type: Number, found: String (5.3.2)
As I mentioned, this used to work with 4.1 and never complained about it being a string. I assume this is because of the two dots instead of one that makes the pipeline interpret it as a string? Perhaps this could be simply fixed by changing the nextflow_schema.json primer_set_version type to string instead of number, but not sure whether this would have any further consequences.
from viralrecon.
Yep, this is a more complex issue with the way parameters are coerced by Nextflow and then mismatch the validation.
Changing the line below to string
does fix error above:
viralrecon/nextflow_schema.json
Line 112 in 3731dd3
However, the pipeline fails when you provide older primer sets that are integers:
ERROR: Validation of pipeline parameters failed!
* --primer_set_version: expected type: String, found: Integer (1)
Any suggestions @ewels ?
For now, I think the best solution is to provide the parameters via a -params-file params.yml
explicitly:
fasta: 'https://github.com/artic-network/artic-ncov2019/raw/master/primer_schemes/nCoV-2019/V5.3.2/SARS-CoV-2.reference.fasta'
gff: 'https://github.com/nf-core/test-datasets/raw/viralrecon/genome/MN908947.3/GCA_009858895.3_ASM985889v3_genomic.200409.gff.gz'
primer_bed: 'https://github.com/artic-network/artic-ncov2019/raw/master/primer_schemes/nCoV-2019/V5.3.2/SARS-CoV-2.scheme.bed'
scheme: 'SARS-CoV-2'
Haven't tested this so please let me know if the pipeline complains elsewhere!
from viralrecon.
In general though, if you plan on using this primer set even more often it's worth downloading the files locally and updating the -params-file
so you don't download the files every time you run the pipeline e.g.
fasta: '/my/local/path/nCoV-2019/V5.3.2/SARS-CoV-2.reference.fasta'
gff: '/my/local/path/MN908947.3/GCA_009858895.3_ASM985889v3_genomic.200409.gff.gz'
primer_bed: '/my/local/path/nCoV-2019/V5.3.2/SARS-CoV-2.scheme.bed'
scheme: 'SARS-CoV-2'
from viralrecon.
When using the params-file do I just need to specify the genome and leave primer set etc as defaults? I actually never used it before but this could be useful for testing new primer sets that arent commercially available yet.
from viralrecon.
Yep, this is a more complex issue with the way parameters are coerced by Nextflow and then mismatch the validation.
Changing the line below to
string
does fix error above:viralrecon/nextflow_schema.json
Line 112 in 3731dd3
However, the pipeline fails when you provide older primer sets that are integers:
ERROR: Validation of pipeline parameters failed! * --primer_set_version: expected type: String, found: Integer (1)
Any suggestions @ewels ?
For now, I think the best solution is to provide the parameters via a
-params-file params.yml
explicitly:fasta: 'https://github.com/artic-network/artic-ncov2019/raw/master/primer_schemes/nCoV-2019/V5.3.2/SARS-CoV-2.reference.fasta' gff: 'https://github.com/nf-core/test-datasets/raw/viralrecon/genome/MN908947.3/GCA_009858895.3_ASM985889v3_genomic.200409.gff.gz' primer_bed: 'https://github.com/artic-network/artic-ncov2019/raw/master/primer_schemes/nCoV-2019/V5.3.2/SARS-CoV-2.scheme.bed' scheme: 'SARS-CoV-2'
Haven't tested this so please let me know if the pipeline complains elsewhere!
I would add that in the ignore_params thingy
from viralrecon.
Ok. Ignore everything I said about using a params file. Can you try with your original command in this comment and add this parameter too please?
--schema_ignore_params 'genomes,primer_set_version'
Thanks @maxulysse !
from viralrecon.
Hello @drpatelh
I'd like to hijack this thread, if I may.
My lab has generated tiled-PCR amplfication assays using PrimalScheme for a couple of other viruses. I was wondering if there is a workaround to allow Viralrecon to take in completely custom schemes? They should be functionally the same as those being generated for SARS-CoV-2.
I wasn't sure if the thoughts in this thread might be applicable to my scenario.
Thanks in advance.
-Cody
from viralrecon.
Hi @Codes1985 . You might be able to use similar parameters to the ones suggested for the SWIFT panel with customisations for your genome reference and primer sets. You might need to check that the pipeline is using the primer sets as expected by checking some of the heatmaps / plots etc:
https://nf-co.re/viralrecon/2.6.0/docs/usage#swift-primer-sets
Might be good to have something like this as generic docs in that section too.
from viralrecon.
Sorry to further the hijacking but did @Codes1985 have any luck with custom schemes? I am looking to do the same kind of analysis with custom primer sets for other viruses as well.
from viralrecon.
Sorry to further the hijacking but did @Codes1985 have any luck with custom schemes? I am looking to do the same kind of analysis with custom primer sets for other viruses as well.
Hello @kneubehl
I'm so sorry - I haven't had a chance to try it out yet! I will reply here as soon as I can. However, if you beat me to it, I 'd love to hear.
Take care!
from viralrecon.
Related Issues (20)
- [bug] Fails to find FASTQ files when stored on S3
- epi2me-labs/wf-cnv failed during analysis HOT 7
- Allow viralrecon to take gtf as input annotation file
- primer_set is not taking param
- Artic v5 mismatched primer names in artic-ncov2019 repo cause certain amplicons to be erroneously filtered/removed
- Add QIAseq DIRECT SARS-CoV-2 Kit amplicons
- Argument input-fasta is missing in NFCORE_VIRALRECON:ILLUMINA:CONSENSUS_BCFTOOLS:CONSENSUS_QC:NEXTCLADE_RUN" HOT 4
- Problem installing viralrecon
- Adding "aggregate" and "plot" methods in the freyja subworkflow HOT 4
- Error running version 2.6.0 with Nanopore data in process NFCORE_VIRALRECON:NANOPORE:ARTIC_MINION
- Unable to download the python script HOT 3
- Split authors in generated RO crate
- Adding "--grouplineages" parameter in the nf-core/viralrecon HOT 3
- Make cutadapt primer's position ext.arg
- Properly deal with multiqc in the config files before the next release HOT 3
- Temp file problem in VARIANTS_IVAR:BCFTOOLS_SORT HOT 1
- Allow skipping `freyja boot` HOT 1
- `MOSDEPTH_AMPLICON` doesn't run in `dev` branch HOT 1
- Non-SCV2 amplicon run returns consensus genomes with no low-coverage masking HOT 5
- nf-core/viralrecon run halted due to R version clash HOT 6
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from viralrecon.