-
Notifications
You must be signed in to change notification settings - Fork 0
/
createImputationPanel.nf
43 lines (36 loc) · 1.03 KB
/
createImputationPanel.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
37
38
39
40
41
42
43
#!/usr/bin/env nextflow
nextflow.enable.dsl = 2
include {
getChromosomes;
getHapmapGeneticMap;
getPlinkGeneticMap;
getThousandGenomesReference;
getVcf;
splitVcfByChrom;
alignGenotypesToReference;
beaglephase;
eaglePhaseWithoutRef;
getEagleHapmapGeneticMap;
getVcfIndex;
} from "${projectDir}/modules/phasing_and_imputation.nf"
include {
getPhasedVcf;
validateVcf;
getValidationExitStatus;
createLegendFile;
getm3vcf;
prepareChrXPanel;
} from "${projectDir}/modules/custom_panel.nf"
workflow {
chromosome = getChromosomes()
vcf = getPhasedVcf()
validateVcf(vcf).map { chr, vcf_file -> tuple(chr.baseName, vcf_file) }.set { chrom_vcf }
getVcfIndex(chrom_vcf).view().set { vcf_fileset }
/*
hapmapGeneticMap = getEagleHapmapGeneticMap()
vcf_fileset.combine(hapmapGeneticMap).set { eagle_no_ref_phase_input }
phased_vcf = eaglePhaseWithoutRef(eagle_no_ref_phase_input)
createLegendFile(phased_vcf).view()
*/
m3vcf = getm3vcf(vcf_fileset).view()
}