Simon Pearce
07/19/2024, 11:16 AMnextflow run main.nf \
--input <https://raw.githubusercontent.com/nf-core/test-datasets/fastquorum/testdata/samplesheets/samplesheet.tiny.csv> \
--genome GRCh38 \
--outdir test \
-profile docker
It won't align anything.
It seems to have all the required boilerplate in https://github.com/nf-core/fastquorum/blob/da1148fed20864d5e041d3201560fde74705bcb3/main.nf#L33
params.fasta = getGenomeAttribute('fasta')
params.fasta_fai = getGenomeAttribute('fasta_fai')
params.dict = getGenomeAttribute('dict')
params.bwa = getGenomeAttribute('bwa')
// Initialize fasta file with meta map:
fasta = params.fasta ? Channel.fromPath(params.fasta).map{ it -> [ [id:it.baseName], it ] }.collect() : Channel.empty()
and when I try with some println action getGenomeAttribute seems to be working absolutely fine.
But setting params.fasta just doesn't do anything:
params.fasta = getGenomeAttribute('fasta')
println params.fasta
gives null , but
params.fasta2 = getGenomeAttribute('fasta')
println params.fasta2
gives <s3://ngi-igenomes/igenomes//Homo_sapiens/NCBI/GRCh38/Sequence/WholeGenomeFasta/genome.fa>
What is silently stopping params.fasta from being reassigned here? It seems to stop the other fasta_fai, dict and bwa parameters on the next lines too.Maxime Garcia
07/19/2024, 11:20 AMgetGenomeAtribute function coming from a Nextflow script instead of the groovy libMaxime Garcia
07/19/2024, 11:20 AMSimon Pearce
07/19/2024, 11:23 AMSimon Pearce
07/19/2024, 11:23 AMSimon Pearce
07/19/2024, 11:23 AMparams.fasta is defined multiple times -- Assignments following the first are ignored`Simon Pearce
07/19/2024, 11:24 AMNicolas Vannieuwkerke
07/19/2024, 11:41 AMgetGenomeAttribute code to before you include all other subworkflows/functions 😛Nicolas Vannieuwkerke
07/19/2024, 11:41 AMSimon Pearce
07/19/2024, 11:43 AMMaxime Garcia
07/19/2024, 11:43 AMSimon Pearce
07/19/2024, 11:43 AMSimon Pearce
07/19/2024, 11:43 AMNicolas Vannieuwkerke
07/19/2024, 11:44 AMSimon Pearce
07/19/2024, 11:46 AMAccess to undefined parameter thoughNicolas Vannieuwkerke
07/19/2024, 11:49 AMSimon Pearce
07/19/2024, 12:01 PMSimon Pearce
07/19/2024, 12:04 PMSimon Pearce
07/19/2024, 12:04 PMSimon Pearce
07/19/2024, 12:05 PMNicolas Vannieuwkerke
07/19/2024, 12:05 PMSimon Pearce
07/19/2024, 12:05 PMNicolas Vannieuwkerke
07/19/2024, 12:06 PMSimon Pearce
07/19/2024, 12:06 PMMaxime Garcia
07/19/2024, 12:14 PMMaxime Garcia
07/19/2024, 12:14 PMNicolas Vannieuwkerke
07/19/2024, 12:15 PMSimon Pearce
07/22/2024, 2:43 PM// Check input path parameters to see if they exist
def checkPathParamList = [
params.bwa,
params.dict,
params.fasta,
params.fasta_fai
]
Is this bit of code outdated now?
I can't find checkPathParamList used anywhere, and I can't see it in a few other pipelines.Nicolas Vannieuwkerke
07/22/2024, 2:55 PMexists: true keyword in your schema to check those parametersSimon Pearce
07/22/2024, 3:19 PMcheckPathParamList can be removed from sarek too?Nicolas Vannieuwkerke
07/22/2024, 3:20 PMexists to the schema 🙂Simon Pearce
07/22/2024, 3:22 PMexists in the schema checks that the file exists if it is present?Nicolas Vannieuwkerke
07/22/2024, 3:23 PMSimon Pearce
07/22/2024, 3:23 PMNicolas Vannieuwkerke
07/22/2024, 3:24 PMNicolas Vannieuwkerke
07/22/2024, 3:24 PM