@@ -4,6 +4,8 @@ nextflow_process {
44 script "../main.nf"
55 process "GLIMPSE_CONCORDANCE"
66
7+ config "./nextflow.config"
8+
79 tag "modules"
810 tag "modules_nfcore"
911 tag "glimpse"
@@ -12,64 +14,21 @@ nextflow_process {
1214 tag "bcftools/index"
1315
1416 test("test_glimpse_concordance") {
15- setup {
16- run("GLIMPSE_PHASE") {
17- script "../../phase/main.nf"
18- process {
19- """
20- ch_sample = Channel.of('NA12878 2').collectFile(name: 'sampleinfos.txt')
21- region = Channel.fromList([
22- ["chr21:16600000-16750000","chr21:16650000-16700000"]
23- ])
24- input_vcf = Channel.of([
25- [ id:'input'], // meta map
26- file(params.modules_testdata_base_path + "delete_me/glimpse/NA12878.chr21.s.1x.vcf.gz", checkIfExists: true),
27- file(params.modules_testdata_base_path + "delete_me/glimpse/NA12878.chr21.s.1x.vcf.gz.csi", checkIfExists: true)
28- ])
29- ref_panel = Channel.of([
30- file(params.modules_testdata_base_path + "delete_me/glimpse/1000GP.chr21.noNA12878.s.bcf", checkIfExists: true),
31- file(params.modules_testdata_base_path + "delete_me/glimpse/1000GP.chr21.noNA12878.s.bcf.csi", checkIfExists: true)
32- ])
33- ch_map = Channel.of([
34- file(params.modules_testdata_base_path + "delete_me/glimpse/chr21.b38.gmap.gz", checkIfExists: true),
35- ])
36-
37- input[0] = input_vcf
38- | combine(ch_sample)
39- | combine(region)
40- | combine(ref_panel)
41- | combine(ch_map)
42- """
43- }
44- }
45- run("BCFTOOLS_INDEX") {
46- script "../../../bcftools/index/main.nf"
47- process {
48- """
49- input[0] = GLIMPSE_PHASE.out.phased_variants
50- """
51- }
52- }
53- }
5417 when {
5518 process {
5619 """
57- allele_freq = Channel.fromList([
58- file(params.modules_testdata_base_path + "delete_me/glimpse/1000GP.chr21.noNA12878.s.sites.vcf.gz",checkIfExists:true),
59- file(params.modules_testdata_base_path + "delete_me/glimpse/1000GP.chr21.noNA12878.s.sites.vcf.gz.csi",checkIfExists:true)
60- ]).collect()
61- truth = Channel.fromList([
62- file(params.modules_testdata_base_path + "delete_me/glimpse/NA12878.chr21.s.bcf",checkIfExists:true),
63- file(params.modules_testdata_base_path + "delete_me/glimpse/NA12878.chr21.s.bcf.csi",checkIfExists:true)
64- ]).collect()
65- estimate = GLIMPSE_PHASE.out.phased_variants
66- | join (BCFTOOLS_INDEX.out.csi)
67- input[0] = estimate
68- | combine (allele_freq)
69- | combine (truth)
70- | combine (["chr21"])
71- input[1] = []
72- input[2] = []
20+ input[0] = channel.of([
21+ [ id: "NA12878" ],
22+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/illumina/vcf/NA12878.chr21_22.1X.glimpse2.vcf.gz",checkIfExists:true),
23+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/illumina/vcf/NA12878.chr21_22.1X.glimpse2.vcf.gz.csi",checkIfExists:true),
24+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/popgen/1000GP.chr21.sites.vcf.gz",checkIfExists:true),
25+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/popgen/1000GP.chr21.sites.vcf.gz.csi",checkIfExists:true),
26+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/illumina/vcf/NA12878_GIAB.chr21_22.vcf.gz",checkIfExists:true),
27+ file(params.modules_testdata_base_path + "genomics/homo_sapiens/illumina/vcf/NA12878_GIAB.chr21_22.vcf.gz.csi",checkIfExists:true),
28+ "chr21"
29+ ])
30+ input[1] = 0.7
31+ input[2] = 3
7332 input[3] = []
7433 """
7534 }
@@ -78,12 +37,7 @@ nextflow_process {
7837 then {
7938 assertAll(
8039 { assert process.success },
81- { assert snapshot(
82- process.out.errors_cal,
83- process.out.errors_spl,
84- process.out.rsquare_spl,
85- process.out.versions
86- ).match() }
40+ { assert snapshot(process.out).match() }
8741 )
8842 }
8943
0 commit comments