@@ -441,7 +441,7 @@ test("reads_son [bam]") {
441441 script " ../../../samtools/index/main.nf"
442442 process {
443443 " " "
444- input[0] = [ [id: 'NA12882 ' ], file(" ${projectDir} /data/bam/reads_son.bam" ) ]
444+ input[0] = [ [id: 'reads_son ' ], file(" ${projectDir} /data/bam/reads_son.bam" ) ]
445445 " " "
446446 }
447447 }
@@ -526,7 +526,7 @@ Test Process GATK_HAPLOTYPECALLER
526526 "0": [ "0": [
527527 [ [
528528 { {
529- "id": "NA12882 " "id": "NA12882 "
529+ "id": "reads_son " "id": "reads_son "
530530 }, },
531531 "reads_son.bam.g.vcf:md5,f3583cbbe439469bfc166612e1617694", | "reads_son.bam.g.vcf:md5,428f855d616b34d44a4f0a3bcc1a0b14",
532532 "reads_son.bam.g.vcf.idx:md5,16a78feaf6602adb2a131494e0274f9e" | "reads_son.bam.g.vcf.idx:md5,5a8d299625ef3cd3266229507a789dbb"
@@ -572,8 +572,8 @@ _After:_
572572` ` ` console title=" modules/local/gatk/haplotypecaller/tests/main.nf.test" linenums=" 35"
573573then {
574574 assert process.success
575- assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12882 ' )
576- assert path(process.out[0][0][1]).readLines ().contains(' 20 10040001 . T <NON_REF> . . END=10040048 GT:DP:GQ:MIN_DP:PL 0/0:40:99:37 :0,99,1150 ' )
575+ assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT reads_son ' )
576+ assert path(process.out[0][0][1]).readLines ().contains(' 20_10037292_10066351 3277 . G <NON_REF> . . END=3282 GT:DP:GQ:MIN_DP:PL 0/0:25:72:24 :0,72,719 ' )
577577}
578578` ` `
579579
@@ -619,7 +619,7 @@ test("reads_mother [bam]") {
619619 script " ../../../samtools/index/main.nf"
620620 process {
621621 " " "
622- input[0] = [ [id: 'NA12882 ' ], file(" ${projectDir} /data/bam/reads_mother.bam" ) ]
622+ input[0] = [ [id: 'reads_mother ' ], file(" ${projectDir} /data/bam/reads_mother.bam" ) ]
623623 " " "
624624 }
625625 }
@@ -635,15 +635,15 @@ test("reads_mother [bam]") {
635635 input[1] = file(" ${projectDir} /data/ref/ref.fasta" )
636636 input[2] = file(" ${projectDir} /data/ref/ref.fasta.fai" )
637637 input[3] = file(" ${projectDir} /data/ref/ref.dict" )
638- input[4] = file(" ${projectDir} /data/intervals.list " )
638+ input[4] = file(" ${projectDir} /data/ref/ intervals.bed " )
639639 " " "
640640 }
641641 }
642642
643643 then {
644644 assert process.success
645- assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12878 ' )
646- assert path(process.out[0][0][1]).readLines ().contains(' 20 10040001 . T <NON_REF> . . END=10040013 GT:DP:GQ:MIN_DP:PL 0/0:28:81:27 :0,81,829 ' )
645+ assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT reads_mother ' )
646+ assert path(process.out[0][0][1]).readLines ().contains(' 20_10037292_10066351 3277 . G <NON_REF> . . END=3278 GT:DP:GQ:MIN_DP:PL 0/0:38:99:37 :0,102,1530 ' )
647647 }
648648}
649649` ` `
@@ -658,7 +658,7 @@ test("reads_father [bam]") {
658658 script " ../../../samtools/index/main.nf"
659659 process {
660660 " " "
661- input[0] = [ [id: 'NA12882 ' ], file(" ${projectDir} /data/bam/reads_father.bam" ) ]
661+ input[0] = [ [id: 'reads_father ' ], file(" ${projectDir} /data/bam/reads_father.bam" ) ]
662662 " " "
663663 }
664664 }
@@ -674,15 +674,15 @@ test("reads_father [bam]") {
674674 input[1] = file(" ${projectDir} /data/ref/ref.fasta" )
675675 input[2] = file(" ${projectDir} /data/ref/ref.fasta.fai" )
676676 input[3] = file(" ${projectDir} /data/ref/ref.dict" )
677- input[4] = file(" ${projectDir} /data/intervals.list " )
677+ input[4] = file(" ${projectDir} /data/ref/ intervals.bed " )
678678 " " "
679679 }
680680 }
681681
682682 then {
683683 assert process.success
684- assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12877 ' )
685- assert path(process.out[0][0][1]).readLines ().contains(' 20 10040001 . T <NON_REF> . . END=10040011 GT:DP:GQ:MIN_DP:PL 0/0:30:81:29 :0,81,1025 ' )
684+ assert path(process.out[0][0][1]).readLines ().contains(' #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT reads_father ' )
685+ assert path(process.out[0][0][1]).readLines ().contains(' 20_10037292_10066351 3277 . G <NON_REF> . . END=3281 GT:DP:GQ:MIN_DP:PL 0/0:44:99:42 :0,120,1800 ' )
686686 }
687687}
688688` ` `
@@ -824,7 +824,7 @@ test("family_trio [vcf] [idx]") {
824824 input[2] = file(" ${projectDir} /data/ref/ref.fasta" )
825825 input[3] = file(" ${projectDir} /data/ref/ref.fasta.fai" )
826826 input[4] = file(" ${projectDir} /data/ref/ref.dict" )
827- input[5] = file(" ${projectDir} /data/intervals.list " )
827+ input[5] = file(" ${projectDir} /data/ref/ intervals.bed " )
828828 " " "
829829 }
830830 }
@@ -837,8 +837,8 @@ The output of the joint genotyping step is another VCF file, so we're going to u
837837```groovy title="modules/local/gatk/jointgenotyping/tests/main.nf.test" linenums="25"
838838then {
839839 assert process.success
840- assert path(process.out[0][0]).readLines().contains(' # CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA12877 NA12878 NA12882 ')
841- assert path(process.out[0][0]).readLines().contains('20 10040772 . C CT 1568 .89 . AC=5;AF=0.833;AN=6;BaseQRankSum=0.399 ;DP=82 ;ExcessHet=0.0000;FS=4.291 ;MLEAC=5;MLEAF=0.833;MQ=60.00;MQRankSum=0.00;QD=21.79 ;ReadPosRankSum=-9.150e-01 ;SOR=0.510 GT:AD:DP:GQ:PL 0/1:14 ,16:30 :99:370 ,0,348 1/1:0,17:17:51:487,51 ,0 1/1:0,25:25:75:726,75 ,0')
840+ assert path(process.out[0][0]).readLines().contains(' # CHROM POS ID REF ALT QUAL FILTER INFO FORMAT reads_father reads_mother reads_son ')
841+ assert path(process.out[0][0]).readLines().contains('20_10037292_10066351 3480 . C CT 1625 .89 . AC=5;AF=0.833;AN=6;BaseQRankSum=0.220 ;DP=85 ;ExcessHet=0.0000;FS=2.476 ;MLEAC=5;MLEAF=0.833;MQ=60.00;MQRankSum=0.00;QD=21.68 ;ReadPosRankSum=-1.147e+00 ;SOR=0.487 GT:AD:DP:GQ:PL 0/1:15 ,16:31 :99:367 ,0,375 1/1:0,18:18:54:517,54 ,0 1/1:0,26:26:78:756,78 ,0')
842842}
843843` ` `
844844
0 commit comments