-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathscrnaseq.nf
36 lines (28 loc) · 995 Bytes
/
scrnaseq.nf
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
// Read and process CTG samplesheet (must be plain .csv - not directly from excel)
samplesheet = file(params.samplesheet)
// Input parsing
include { SPLITSHEET } from "../modules/split_sheet/main"
// Analysis Modules
include { COUNT } from "../modules/cellranger/count-rna/main"
// QC Modules
include { FASTQC } from "../modules/fastqc/main"
// Delivery modules
include { FINISH_PROJECTS } from "../subworkflows/finish.nf"
workflow SCRNASEQ {
main:
// Parse samplesheet
sheet_ch = SPLITSHEET(samplesheet, 'scrna-10x')
// all samplesheet info
sample_info_ch = sheet_ch.data
.splitCsv(header:true)
.map { row -> tuple( row.Sample_ID, row.Sample_Species, row.force, row.agg, row.Sample_Project) }
sample_fastqc_ch = sheet_ch.data
.splitCsv(header:true)
.map { row -> tuple( row.Sample_ID, row.Sample_Project) }
FASTQC(sample_fastqc_ch)
count_ch = COUNT(sample_info_ch)
FINISH_PROJECTS(
count_ch.project_id.collect().flatten().unique(),
'scrna-10x'
)
}