2021-02-03 16:02:58 +00:00
|
|
|
#!/usr/bin/env nextflow
|
|
|
|
|
|
|
|
nextflow.enable.dsl = 2
|
|
|
|
|
2021-02-03 16:20:02 +00:00
|
|
|
include { SALMON_INDEX } from '../../../../software/salmon/index/main.nf' addParams( options: [:] )
|
|
|
|
include { SALMON_QUANT } from '../../../../software/salmon/quant/main.nf' addParams( options: [args: '--minAssignedFrags 1'] )
|
2021-02-03 16:02:58 +00:00
|
|
|
|
|
|
|
workflow test_salmon_quant_single_end {
|
|
|
|
|
|
|
|
def genome_fasta = file("${launchDir}/tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.fna", checkIfExists: true)
|
|
|
|
def transcript_fasta = file("${launchDir}/tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_cds_from_genomic.fna", checkIfExists: true)
|
|
|
|
def gtf = file("${launchDir}/tests/data/gff/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.gtf", checkIfExists: true)
|
|
|
|
def input = [ [ id:'test', single_end:true ], // meta map
|
|
|
|
file("${launchDir}/tests/data/fastq/rna/sarscov2/EPI_ISL_486436_1.fastq.gz", checkIfExists: true) ]
|
|
|
|
|
|
|
|
SALMON_INDEX ( genome_fasta, transcript_fasta )
|
|
|
|
SALMON_QUANT( input, SALMON_INDEX.out.index, gtf, transcript_fasta, false)
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
workflow test_salmon_quant_paired_end {
|
|
|
|
|
|
|
|
def genome_fasta = file("${launchDir}/tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.fna", checkIfExists: true)
|
|
|
|
def transcript_fasta = file("${launchDir}/tests/data/fasta/sarscov2/GCA_011545545.1_ASM1154554v1_cds_from_genomic.fna", checkIfExists: true)
|
|
|
|
def gtf = file("${launchDir}/tests/data/gff/sarscov2/GCA_011545545.1_ASM1154554v1_genomic.gtf", checkIfExists: true)
|
|
|
|
def input = [ [ id:'test', single_end:false ], // meta map
|
|
|
|
[ file("${launchDir}/tests/data/fastq/rna/sarscov2/EPI_ISL_486436_1.fastq.gz", checkIfExists: true),
|
|
|
|
file("${launchDir}/tests/data/fastq/rna/sarscov2/EPI_ISL_486436_2.fastq.gz", checkIfExists: true) ] ]
|
|
|
|
|
|
|
|
SALMON_INDEX ( genome_fasta, transcript_fasta )
|
|
|
|
SALMON_QUANT( input, SALMON_INDEX.out.index, gtf, transcript_fasta, false)
|
|
|
|
|
|
|
|
}
|