From a552d97849b94e04ae55674af08dc80d3368efef Mon Sep 17 00:00:00 2001 From: Kevin Menden Date: Tue, 9 Mar 2021 10:04:08 +0100 Subject: [PATCH] Exchange VCF data by sarscov2 data (#261) * added sarscov2 vcf data; adjusted bcftools tests * updated README.md * deleted old VCF files; fixed mergevcfs * fixed tabix --- software/bcftools/consensus/main.nf | 12 ++--- tests/data/README.md | 5 +++ tests/data/dummy/dummy_file2.txt | 0 tests/data/generic/vcf/test.consensus.fa | 20 --------- tests/data/generic/vcf/test.vcf | 24 ---------- tests/data/generic/vcf/test.vcf.gz | Bin 498 -> 0 bytes tests/data/generic/vcf/test.vcf.gz.tbi | Bin 123 -> 0 bytes tests/data/generic/vcf/test2.vcf.gz | Bin 498 -> 0 bytes tests/data/generic/vcf/test2.vcf.gz.tbi | Bin 123 -> 0 bytes tests/data/generic/vcf/test3.vcf.gz | Bin 498 -> 0 bytes tests/data/generic/vcf/test3.vcf.gz.tbi | Bin 123 -> 0 bytes tests/data/generic/vcf/test_a.merge.vcf.gz | Bin 763 -> 0 bytes .../data/generic/vcf/test_a.merge.vcf.gz.tbi | Bin 141 -> 0 bytes tests/data/generic/vcf/test_b.merge.vcf.gz | Bin 769 -> 0 bytes .../data/generic/vcf/test_b.merge.vcf.gz.tbi | Bin 141 -> 0 bytes .../sarscov2/bam/test_single_end.sorted.bam | Bin 0 -> 10732 bytes tests/data/genomics/sarscov2/vcf/test.vcf | 41 +++++++++++++++++ tests/data/genomics/sarscov2/vcf/test.vcf.gz | Bin 0 -> 1312 bytes .../genomics/sarscov2/vcf/test.vcf.gz.tbi | Bin 0 -> 112 bytes tests/data/genomics/sarscov2/vcf/test2.vcf | 41 +++++++++++++++++ tests/data/genomics/sarscov2/vcf/test2.vcf.gz | Bin 0 -> 1312 bytes .../genomics/sarscov2/vcf/test2.vcf.gz.tbi | Bin 0 -> 112 bytes tests/data/genomics/sarscov2/vcf/test3.vcf | 42 ++++++++++++++++++ tests/data/genomics/sarscov2/vcf/test3.vcf.gz | Bin 0 -> 1287 bytes .../genomics/sarscov2/vcf/test3.vcf.gz.tbi | Bin 0 -> 112 bytes tests/software/bcftools/consensus/main.nf | 7 ++- tests/software/bcftools/consensus/test.yml | 2 +- tests/software/bcftools/filter/main.nf | 3 +- tests/software/bcftools/filter/test.yml | 2 +- tests/software/bcftools/isec/main.nf | 13 +++--- tests/software/bcftools/isec/test.yml | 16 +++---- tests/software/bcftools/merge/main.nf | 9 ++-- tests/software/bcftools/merge/test.yml | 2 +- tests/software/bcftools/mpileup/main.nf | 6 +-- tests/software/bcftools/mpileup/test.yml | 7 ++- tests/software/bcftools/stats/main.nf | 3 +- tests/software/bcftools/stats/test.yml | 2 +- tests/software/gatk4/mergevcfs/main.nf | 16 +++---- tests/software/gatk4/mergevcfs/test.yml | 4 +- tests/software/tabix/bgzip/main.nf | 3 +- tests/software/tabix/bgzip/test.yml | 2 +- tests/software/tabix/tabix/main.nf | 2 +- tests/software/tabix/tabix/test.yml | 2 +- 43 files changed, 176 insertions(+), 110 deletions(-) delete mode 100644 tests/data/dummy/dummy_file2.txt delete mode 100644 tests/data/generic/vcf/test.consensus.fa delete mode 100644 tests/data/generic/vcf/test.vcf delete mode 100644 tests/data/generic/vcf/test.vcf.gz delete mode 100644 tests/data/generic/vcf/test.vcf.gz.tbi delete mode 100644 tests/data/generic/vcf/test2.vcf.gz delete mode 100644 tests/data/generic/vcf/test2.vcf.gz.tbi delete mode 100644 tests/data/generic/vcf/test3.vcf.gz delete mode 100644 tests/data/generic/vcf/test3.vcf.gz.tbi delete mode 100644 tests/data/generic/vcf/test_a.merge.vcf.gz delete mode 100644 tests/data/generic/vcf/test_a.merge.vcf.gz.tbi delete mode 100644 tests/data/generic/vcf/test_b.merge.vcf.gz delete mode 100644 tests/data/generic/vcf/test_b.merge.vcf.gz.tbi create mode 100644 tests/data/genomics/sarscov2/bam/test_single_end.sorted.bam create mode 100644 tests/data/genomics/sarscov2/vcf/test.vcf create mode 100644 tests/data/genomics/sarscov2/vcf/test.vcf.gz create mode 100644 tests/data/genomics/sarscov2/vcf/test.vcf.gz.tbi create mode 100644 tests/data/genomics/sarscov2/vcf/test2.vcf create mode 100644 tests/data/genomics/sarscov2/vcf/test2.vcf.gz create mode 100644 tests/data/genomics/sarscov2/vcf/test2.vcf.gz.tbi create mode 100644 tests/data/genomics/sarscov2/vcf/test3.vcf create mode 100644 tests/data/genomics/sarscov2/vcf/test3.vcf.gz create mode 100644 tests/data/genomics/sarscov2/vcf/test3.vcf.gz.tbi diff --git a/software/bcftools/consensus/main.nf b/software/bcftools/consensus/main.nf index 794b854e..2a15136c 100644 --- a/software/bcftools/consensus/main.nf +++ b/software/bcftools/consensus/main.nf @@ -11,19 +11,19 @@ process BCFTOOLS_CONSENSUS { mode: params.publish_dir_mode, saveAs: { filename -> saveFiles(filename:filename, options:params.options, publish_dir:getSoftwareName(task.process), publish_id:meta.id) } - conda (params.enable_conda ? "bioconda::bcftools=1.11" : null) + conda (params.enable_conda ? 'bioconda::bcftools=1.11' : null) if (workflow.containerEngine == 'singularity' && !params.singularity_pull_docker_container) { - container "https://depot.galaxyproject.org/singularity/bcftools:1.11--h7c999a4_0" + container 'https://depot.galaxyproject.org/singularity/bcftools:1.11--h7c999a4_0' } else { - container "quay.io/biocontainers/bcftools:1.11--h7c999a4_0" + container 'quay.io/biocontainers/bcftools:1.11--h7c999a4_0' } input: tuple val(meta), path(vcf), path(tbi), path(fasta) output: - tuple val(meta), path("*.fa"), emit: fasta - path "*.version.txt" , emit: version + tuple val(meta), path('*.fa'), emit: fasta + path '*.version.txt' , emit: version script: def software = getSoftwareName(task.process) @@ -31,7 +31,7 @@ process BCFTOOLS_CONSENSUS { """ cat $fasta | bcftools consensus $vcf $options.args > ${prefix}.fa header=\$(head -n 1 ${prefix}.fa | sed 's/>//g') - sed -i "s/\${header}/${meta.id}/g" ${prefix}.fa + sed -i 's/\${header}/${meta.id}/g' ${prefix}.fa echo \$(bcftools --version 2>&1) | sed 's/^.*bcftools //; s/ .*\$//' > ${software}.version.txt """ diff --git a/tests/data/README.md b/tests/data/README.md index 6e94c988..8d02d790 100644 --- a/tests/data/README.md +++ b/tests/data/README.md @@ -34,6 +34,11 @@ When adding a new module, please check carefully whether the data necessary for * 'test_genomic.gff3.gz': bgzipped-version * paf * 'test_cds_from_genomic.paf': PAF file for MT192765.1 genome + * vcf + * 'test.vcf', 'test2.vcf': generated from 'test_paired_end.sorted.bam' using bcftools mpileup, call and filter + * 'test3.vcf': generated from 'test_single_end.sorted.bam' using bcftools mpileup, call and filter + * '*.gz': generated from VCF files using bgzip + * '.tbi': generated from '.vcf.gz' files using `tabix -p vcf -f ` ### generic diff --git a/tests/data/dummy/dummy_file2.txt b/tests/data/dummy/dummy_file2.txt deleted file mode 100644 index e69de29b..00000000 diff --git a/tests/data/generic/vcf/test.consensus.fa b/tests/data/generic/vcf/test.consensus.fa deleted file mode 100644 index d7688b04..00000000 --- a/tests/data/generic/vcf/test.consensus.fa +++ /dev/null @@ -1,20 +0,0 @@ ->1:2-501 -TACcAtATgTgACAtATAAaAAAGAACATAACCTACGTATCAACTAAAGTGGTTGTTTG -cAGAAAAGGAAGACTTAAAAAGAGTCAGTACTAACCTACATAATATATACAATGTTCATT -AAATAATAAAATGAGCTCATCATACTTAGGTCATCATAAATATATCTGAAATTCACAAAT -ATTGATCAAATGGTAAAATAGACAAGTAGATTTTAATAGGTTAAACAATTACTGATTCTC -TTGAAAGAATAAATTTAATATGAGACCTATTTCATTATAATGAACTCACAAATTAGAAAC -TTCACACTGGGGGCTGGAGAGATGGCTCAGTAGTTAAGAACACTGACTGCTCTTCTGAAG -GTCCTGAGTTCAAATCCCAGCAACCACATGGTGACTTACAACCATCTGTAATGACATCTG -ATGCCCTCTGGTGTGTCTGAAGACAGCTACAGTGTACTTACATAAAATAATAAATAAATC -TTTAAAAACAAAAAAAAAGAA ->2 -gaagatcttttccttattaaggatctgaagctctgtagatttgtattctattaaacatgg -AgagattagtgattttccatattctttaagtcattttagagtaatgtgttcttaagatAa -atcagaaaaacaaaaacttgtgctttcctgtttgaaaaacaaacagctgtggggaatgGt -gtcgggacagcctttttatAaaatttttctaaataatgttgaggctttgatacgtcaaag -ttatatttcaaatggaatcacttagacctcgtttctgagtgtcaatggccatattggggA -tttgctgctgccaatgacaGcacaccctgggaatgccccaactacttactacaaagcagt -gttacatggagaagatcttcaagagtctttttgctagatctttccttggcttttgatgtg -actcctctcaataaaatccacagtaatatagtgagtggtctcctgctccaaaccagtatt -Tcagacacagttaatccagac diff --git a/tests/data/generic/vcf/test.vcf b/tests/data/generic/vcf/test.vcf deleted file mode 100644 index 88895fc7..00000000 --- a/tests/data/generic/vcf/test.vcf +++ /dev/null @@ -1,24 +0,0 @@ -##fileformat=VCFv4.2 -##FORMAT= -##reference=file://some/path/human_g1k_v37.fasta -##INFO= -##INFO= -##ALT= -##contig= -##contig= -#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT NA001 -1 5 . C A . PASS . GT 0/1 -1 5 . C T . PASS . GT 0/1 -1 7 . T A . PASS . GT . -1 10 . G A . PASS . GT 0/1 -1 12 . GACA GA . PASS . GT 0/1 -1 16 . T TAAA . PASS . GT 1/1 -1 19 . A C . PASS . GT 0/1 -1 61 . C A . PASS . GT 0/1 -2 61 . agag aa . PASS . GT 0/1 -2 119 . aaa t . PASS . GT 0/1 -2 179 . g gacgtacgt . PASS . GT 0/1 -2 200 . a . PASS END=210 GT 1/0 -2 300 . a . . PASS END=310;MinDP=10 GT 0/1 -2 320 . a <*> . PASS END=330;MinDP=20 GT 0/1 -2 481 . t c,a . PASS . GT 0/2 diff --git a/tests/data/generic/vcf/test.vcf.gz b/tests/data/generic/vcf/test.vcf.gz deleted file mode 100644 index da39e1d1b6c3ff3b1cc2d06fe8819c91f05a4534..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 498 zcmVG`WmCG#Cz()n-95 z{Fp{~mbuGwwPO-PH@lt~-YUsO=|-%{B|j-%3w$?YqE=eGd(z306-+XbGihIw4&qD_bjGMxA1*P^@I~wpuLq&Z%y=iJZD{)PXq{ zDmfL(*;{bO_6`?X;Ef|$Udz4h>`$Qu7t#77rJ{3HU2QIzxBC&r^A5s}%&aVx*Fy2q zg3zT+!>Zl!ig)3}Ju*?M14g3g_ zPaa##jF4Luw(5h9h;HWJ#4-F&|EQH{NlQR!j|{CV)tc&4Cq^wnN%K;fr-wm-nU{h` zed15Ay1t

5Z_>2bv0BDDN#1!r-gDVz@!&ex+fk&UsdjWaSYsrpHfni5)v2^-86WaP6fT? H8;Ae^YwRQC diff --git a/tests/data/generic/vcf/test2.vcf.gz b/tests/data/generic/vcf/test2.vcf.gz deleted file mode 100644 index da39e1d1b6c3ff3b1cc2d06fe8819c91f05a4534..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 498 zcmVG`WmCG#Cz()n-95 z{Fp{~mbuGwwPO-PH@lt~-YUsO=|-%{B|j-%3w$?YqE=eGd(z306-+XbGihIw4&qD_bjGMxA1*P^@I~wpuLq&Z%y=iJZD{)PXq{ zDmfL(*;{bO_6`?X;Ef|$Udz4h>`$Qu7t#77rJ{3HU2QIzxBC&r^A5s}%&aVx*Fy2q zg3zT+!>Zl!ig)3}Ju*?M14g3g_ zPaa##jF4Luw(5h9h;HWJ#4-F&|EQH{NlQR!j|{CV)tc&4Cq^wnN%K;fr-wm-nU{h` zed15Ay1t

5Z_>2bv0BDDN#1!r-gDVz@!&ex+fk&UsdjWaSYsrpHfni5)v2^-86WaP6fT? H8;Ae^YwRQC diff --git a/tests/data/generic/vcf/test3.vcf.gz b/tests/data/generic/vcf/test3.vcf.gz deleted file mode 100644 index da39e1d1b6c3ff3b1cc2d06fe8819c91f05a4534..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 498 zcmVG`WmCG#Cz()n-95 z{Fp{~mbuGwwPO-PH@lt~-YUsO=|-%{B|j-%3w$?YqE=eGd(z306-+XbGihIw4&qD_bjGMxA1*P^@I~wpuLq&Z%y=iJZD{)PXq{ zDmfL(*;{bO_6`?X;Ef|$Udz4h>`$Qu7t#77rJ{3HU2QIzxBC&r^A5s}%&aVx*Fy2q zg3zT+!>Zl!ig)3}Ju*?M14g3g_ zPaa##jF4Luw(5h9h;HWJ#4-F&|EQH{NlQR!j|{CV)tc&4Cq^wnN%K;fr-wm-nU{h` zed15Ay1t

5Z_>2bv0BDDN#1!r-gDVz@!&ex+fk&UsdjWaSYsrpHfni5)v2^-86WaP6fT? H8;Ae^YwRQC diff --git a/tests/data/generic/vcf/test_a.merge.vcf.gz b/tests/data/generic/vcf/test_a.merge.vcf.gz deleted file mode 100644 index 7ab659d78f2160ed965fd0ef136a501a1d1047b8..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 763 zcmVNEnvctRr>;HdFK00xu!og0*h zBfz~3puwkv$1u#Z3^R!78ZH|`inoL<8>(OSsyCC;RYJvpGA@_P#^WyMLM~fF1$11I zH?^yH>@)E+VVp$qm&YAW6aEYfviTx*=?NEWxWbGuO$`-fPYOmfh*N%?ml4YTVY1ZM zcV$@L)&(!b`cH!x5}I={`NllERJ%@1Nilq4yL`K~p4wk>4JR+AlDU?3UDNiD$M;c~~ zWUKVqSy|UwN@&Dm$Fy(Ew$U+BO>O>@+Im-wZoV3(Xg16>bjop7%zw}mWz29BcrWvi~J7>InyKbpbLYZ;e_5FM)i1?VR|p7?)TPpYl9ib zii_77aa`42`nRgRt`E0&Nx8adN>U;UZ$~*yD2xv0%ub8V9ba^LjW1q|GUpHfm%5)u-ak|cPUP6bFEF;7WgG;CGWii%J)+Lj}5 zD_h7#)MrlPL%SOSd~6{Hn6K8);0=|P4v|Td^ZCayY3AZD8n#~jmll8dup}|nexu1E bzD5Ry_x#Hq9b;o)kVmsnnt>T?KZpPT-=Zl6 diff --git a/tests/data/generic/vcf/test_b.merge.vcf.gz b/tests/data/generic/vcf/test_b.merge.vcf.gz deleted file mode 100644 index 751f5a66d4cf4df0039d538b18e2dc58fe1fd97a..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 769 zcmV+c1OEIUiwFb&00000{{{d;LjnNg0;QB)bD}U5hOg6Kki5ySadTyzi5fa1PDpI!@?1%TprkUv6(_weVH5IT-cE z_P2p&`}cDnjkf8M<~BzD>xS9`CTK`=#PSf9kxgi#zu*Uwv$#Nrj^bbT#v>4oGgxFno%R7~Y`M{`aQB#Ahzs|eM2zgX() zx6-d)>s;jW{KrA`Da(Xh{9vA5@=Yo?oR){EQ@qI~lo=%}=+TXcK14CB{Qc73OEMzP z#4wDXXc9+Rwkk1$p?^OsFwa;&_UyxjOPXXa5bJe={IKG5^2l5kG%xf)#zY(zdW;C? zG+idI_R_pT31y*(Y{R-W+Iq*pHMQ|qYV)KT-+s{z!!q@n+B&1gxNR9X9i^68MRUpm zYOBKRYMPd8xyY#|k4r=lVZE(w^6WPa1lL0Rq)CI8l9)a2XSHaX680#U9(U$VYfZQy z1(&Q+YCA=H8{ZY}eSJ9YCFSZQ2~DUhye$(jW-#1;XMS94&gi0_MVCcu0mkFBQv)^xULIK z>-6ptQ*A8G@47XHfBpcUB8<1&`yT{Whycj6KLtVhg&1QW398+yt}yx@GXyQI)y2p{ z9b}+3GJzVI$hi6g6%Y)lk_P|)ABzYC000000RIL6LPG)o8vp|U0000000000aJOu- diff --git a/tests/data/generic/vcf/test_b.merge.vcf.gz.tbi b/tests/data/generic/vcf/test_b.merge.vcf.gz.tbi deleted file mode 100644 index eadd6a4d3cbef55b424d2b57c6b8042adb39f394..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 141 zcmb2|=3rp}f&Xj_PR>jW1q|GUpHfm%5)u-ak|cPUP6bFEF;7WgG;CGWii%J)x|Snw zD_h7#)MrlPL(Llkd~6{Hn6K7%@P^7thsdPK`TXOUG;{G64O_4NON+mJSdy4(ztLn7 bUn2v(&A<${A4C8E-e)NC diff --git a/tests/data/genomics/sarscov2/bam/test_single_end.sorted.bam b/tests/data/genomics/sarscov2/bam/test_single_end.sorted.bam new file mode 100644 index 0000000000000000000000000000000000000000..e86eb95593797b19c7ba32c98898030b3b96eda5 GIT binary patch literal 10732 zcmVg71b!6!!%~o z1qo@p)x-o&Vn)Lb4XaTygyR9rrII;6P?*G}b}WiXJ1D{if&IzbKmorXQ=lC;me;K? zc9jXGbGgBO<&ERhAbxek+4U7=cy6`cn@3cslaE+}iBEj_=p@A%1?d`jsF5n$i@L2+?!7PcTmM5FUDzM&+F@OzpysBm4csgMcfrqm!EONn03VA81ONa4009360763o z0O2UCT1%`f$64-cJ0ZTdnp2gUOT{fhb*GgSH|FNN=B3!tt)4bi1WBwMi4gIMLSO+n zL^u);0dln}q-3KIz)qyB?gAthEZ87r5d=sfD+FEy$r2lcB7_uy5K35h0r~3DGd+D~ z=8kVt_nvul_slu<_5b+(@BjZji|gi%EF%~1r#G^b?4x-9|7rPi&wTj9)6>~ii{%o~){(O~v;unAE;YV*i{p-KF%Gw)or`G|#~M{LX&Y*J<+x(pbwm+r6{UuhreRiVQ`MHp($?+E)Km2UG^{>i)FzL;?I_;?)uD_P7J9`12G$$sttj+gW0a*_QjOvV|fN=a1PT4Q8mE3GJ3 z4Ohf9t~F(?sO~HwsxcOa>hP!Lx$xlRqR8MMgCKZYHDXtd+ra{h0y<3pA+oh~mggXLJl4RLJ=tAM>w*%8|+t}DWw zHLjzM@UAqxWYoBlmaKJ(FhXEcyp@i(nh{rXsU@$lXZS%O@=5Rpc++H(r~4PSi4$RM z3?n>>2wz0e0~SFN^%}hdT(lwo-IK`%vSQn`Ch6ImRYLIwET{68Qw=SL07si}{EL_V z<-H!qA5A#E{gI6ae(?Z~i_7WddG;ybIKr*k#*G}pjk!ikd{FU*;7Q+2Acakfd~GEI zqw#2^w{@_qya;cKU}br|b|HT>g6uzE0QmT|xU2{<&fBGTzaoz~3Lm`ca2-wvAN>5? z@^%d*9|Dq3cqE?~BKbIw{6s?X)%K^J@<@I&`QZBpNdC+LBxlR%*+urZ@IitIC3y!Z zTJ2qpkj5zlX=7WXc&9qsnXc5f?r2Hu2CMJg%99l(2mG?ST77x2-5lHf>5&6oCc`0* zZ5;5W14(K+U!0$1fB6<3HzjLY7Eig4e(&y5N}Uz6FPb5666($KySccRZ4JhDhE5k}vN`0dLKLWF0N=5R&JYi`g_FDXr6O zQyH!eL-7hlw<>WYvoFj5A1C9*>zA?Jh3VE+@J(xew)G3goGzTZ zpGAItWRY1y@|lxOO#I~oNX}-<`82!tW>_Q=rRTGyq~IrwX%upRMf%iDNg4c^b5c=9 zt5%n-MydtNr#rW`V-1{na8)4^y*Rvv#O&*^rmv6E;Lq=gZ$C+ovZl+~^l};r1(gb> zW6DNzVzucUb?`|bQ&yIMm3pd5DlSV2Uv$zK9}YXjI}2Nq)K@*4gJg7$9BobeuNlR- z{MwoI3b~_aec|OK{caISNzrqeI#ILQTT?x>rndsg-(M4o*WNpfZx6AfC@wht;&dJ; zt|XR&)k})3${OdCrl#@$sZKKIB)FKST$c@)aOdkHqFcggMTs`BECb+i21!?$NLAO{ zi*AvJG9?757;diPHhTAlJ$tWd4m0lc;^+r};~$PJ@_h-%JChBLf4m>Z=_#trStfv^ zX)SZiM>mFon&1s}8swHYB}X3akQ<27(vliQr6xpLm?C994}%UP8Zy#Y$2E;0h5bzG zV!f}T-_qK?^8QeL+!8xq9$rHPuUn$$EPpgYaFr1Jo5{wfUpau_#rbl+%$@~;C5)() za8^qI8K`Y*9jppM1#Gi5q-2s3VmhNKDO=W2=>V9r9V;s;n~vdw(y-F%lYK6oKZtql zSb506@un2%9(hZT)mij{L^d}UBp7dd7;ldq`m=jreB#!|5D&q4c{)9hFe+VE&S_)4 z(^AAk33jYdsx}-ox-&352@YsVkQ8nzWnPD^ zI=`5o&O%(&^eD{u2s4mEUxpA}b&TNIC!O$Tk6H4AuYC96`?zm!RpP?B^=-P zA0G_`>pw>rvp*Y{;|C7mIJ;OZvR6^wHk{GcN<>4{`JF5k$Fr^sq&Ti!XS}vT*$aWw zZ(>^dlZWpH$JdYETDE=K=hL2od`VrczBJg_O>FGxk%9bRGLUEAw%OS4?%&w#?Ck9P zB761wvJ8||8d9;67zm5Vi_lF$Yh8=%$vDz#V-d;>_y!`DM6HGTf#NXunPk}T_T3P@ zx71>De>WHGyyxOUg4vLxy{@iSuMYq^0YIM~0rXVwL2r)$Jp_P0cK{#=6&H)F#nF4H z-zoG!1)5NCkTKYRb>Q`Fhg9gpysxGF;{Lb^I#h?10BmZRa84u%*tRqIAf6P1kzS z(L(Sj6>NTN^$S6(7Lm*byDECGO}h$aBZg>;`y(W_p5EBeNn$^RcD&O>O9H75#TJrj z>r*Iz6sc-n?_=v}5c0(4Sp?|USM=E-7d`y-4!p(SP^?gruI_JEl#rvtA7E%6#?Opk z{NM=2hYuI8i|K57KF>~2z@ysDKKAefz-eUZYQs7x=Y~6JA&SB(5h2UHjQ7wQ3U%SR zSl^Q#k-C2?LZFQ5lBaq8br7;Z`LV7ambDcJbdC6|t`+lJRYi z=#NCh{LleJ7gKaIv;PLk7#kQACe;3|Hm*gAF}5V!S|eYGLS`Q9AYcCv{enCm5QRxb|r_9pCYp7!=X>YoCE`iW#MeA8?~0}7vshk5sD}^4!t!B=y@VF1sYv-P#Yu1DT0O@S~S6^kbe6@7I+Ha z`ETlK^}nM+-HdFdPX~|u@KNbtdA6L*&q6xj@FPo+O*=@m4*o?;W6?aQh-?{T$*wQx=W}mRgobx$wHy>m%HhcL(^Sj4;`bJ0b1+XbyF@|H@M{5ZBupJOsVqqQN53ce3H^;36@k)Z?a5tMsdVLkC)v?~=_;!1& zg40Z*g5Nhb6V>%y1mCPIVys<#LtU-DF~H?f;PT@Mmv`NJG~n_~!sQ>|jh4s}vlNiI zX3L9M7fa3?*%)nG%OStQmk@wjjt6C+7~0N4*k{sHd`wROrLtqvM<;5+lttOVS`dwG ztcz2^nZ!Ov9(&s&jRn+`6qI|a8Nuw+`4HCri-H6*CcE*dw>>{XQH0tdHhaaSNL~{< zF!h09id{0z+Y*qku7Nzse)|xR=cn`eJi7-A`>t?|e|ppv@*pCMKH#Q~f00;dKQNN7 zv;7*HeDK(>$?#M!;Y{M;Vd?lSb+!7|0G#g!IM0p1Q3)LVuq`mI zJLfte<=agW5%I~_-x&apb!MmzHk7F!02&mNn?BJyy@kDVtYD34=8q*1|N6;|g?;fL zi02R~LZ$;NXeD|dT%uJ1oeI6R5~Dk4n^3gQwS!c_&<4T87wRGujd)o}%s{ykjX4eV zN_GsBK4>xd!l$Rr^ zetK~pq8?gH67z9T>HWla83sA|VtB|D(V8+6Sr`ogk=V1XFgpg;sU;m9$a zB~XT>Z{R_b!(BiewGH&4Yq75o+bV35h{U#siqZ3qQ^&56Ha_(1>uK0?*q8PN$@d%= zB+sU2^NU$52Vh%HVU*F9O$o^jYK6tv0S1Ru2Y@lbfuS=<{j4zVl$-nRdu!c>S!3@BYRp!r)S3nHKf_a=^_L*D1;8(M^&N- zf>dMBHL5^Ctf6p6(xL;5ksysm0XhT*PU+AH(H4PBLAC81=7cJYT_NbvP?@JR5MoSG z56Wv1n`;;xjj<|@W92cLrManaLGJl*|0C)Sa&*{|e=(xWKv*SPE>+RvUH7tzkAiz%$pM>Gc!8(m1&LMOtMc83{eblr%b zU)LD#$B0#p29@?b%?5Loro?pMcUn8UP+798j-O=4a6xEpm(?>BD~!9Etu)bxbGpwt z-A6}qt8Y0|(F6x>9)|caY9ML^M@TdXf!wMz+EF9H81>D7s31XnQSd?WIjBS$r*Ob4ro@UC;dUvt2!C{`kZJ{eQ zk0epXzYc=Y;%Y%yXnL2SfZ60R&?J`-)u4MoZ>G_;q78i9(~!Z$1%a5B=Y3=)J~G`A$svTOCur+nMt3*7e#xjG1zC zvD@98J-8>hK`=L9bs|h&D&ItSl+?ZDvBVR}S|l(W$oHguo=TIHmV1^3P@O?L2_S0L zzowP$x=7~+Yv0(XFCHG44cKfbpvYe7taVSUb+7Aqzc}D{|3206AaOW;+a#C5b_KOr zmeDwc(aWxxLsF&xcVJBmRe?@xxWDgdSEfN~EZBI$hL^Qc4P&8`JLnKzdE(V^$9uM? zXXQ=#1DlDa+^+FGOBEE~fldVNRPaeETmGM93CRkTFo73Oj55H3dAT}fAPte`))jXx za~3*qaRSUqr!Y((8+LB>4Be+{eRH~;?f0`=-aLMdn^NX~J~e>Klnqmkp)!{j>)pnp zGQuNG!{wYX4TD!ovOEKpqWUfAuI%-3PyUydGxPMu1gk%1ADTEjs zQ%GKxu%7dyH2)y1CgN+lJm?(ciUPNQI07noh^UUWo-gT0lp)X`1xgZWu$sB1=aP1i1a2?xz9s)Xwa0sQKmEwY}WlXy%xz49I;@x|#Qk`s1$*pfYb7D>B?PJNObWhjO~T?UaEZ4mS@f-Er?hB+%gN&#UYOD~hOfP@6{8bF&q z3d*T05lvtO2^R6fvvUPzoh@7eS9P7uA3h?RukHuuu+ETiy`l%L??1vu-Z)+a2R%rk zH{{2!r-Ay!nu{B2$S5kbE(%Q#3?3^KRrz)=tQ-p8LQ-!qYGxpzn;!%)1Vvn>36QId zBb7{d(JRY}Lv+-6D+_$gnS;$XM>-O^1`o;hJrLDOTZU$DJRn}eC^dyHz*w?%4a$8Z1c+98WkJKk)r z;hK)3QCtF&649>$oo0NH;@EKG$E2PWCzPYB0Lf5vGu^J zbF{GKUKP6=Z8zWhXk^g#cVo9FnoYHazlOcdhavBA($;(}#i)rGBLfM@89r-B2@OyQhKB!&>F>_`iV9ulBKY)Zt=;W`BZ zs6wy^!BX==hB-4@S4-?l=ziRI1pg%dOEm@G=t}wBQB`Ja{j%15CcoQt8i)9efw1z| z``VNvVP$>3*<7Dr+4j>OCeyOm<6A4nMOOO~S~aaXa2{TZ_8ME`owV`b`A7;>pyYrX zvbGDErDnoz&t6KV9xA&x#qJTz+AO98JF0=Q8Yc?xqmmyai=z;Pn=X%8ga}te70qx8-8rXETkq z*z%^7%b)!|sOo%C7$vtd3|}qh_OWI$!1JO%}6-+gc;wytF<32C8Sm1#70WO@V%W`eSbYxj*MpBpV>|IDpR&ek^>SEB~f%uN7U)%A++Kyvu8@%@F#MeL-P;zOlHs_*SAR0ryZHLK!9t zna`L@@h83boB++hc2E>qN#E{pT$XfQ4rQj1!XY9_wFo>!2wyeVcKb_h3e|bJc|)vq z-)>jm7b{28@mp+_OQ{<9jT^sWZyJp^pBacZ-<3^^H#9|ej_wfb1GsEV6J&ehavIP- zprnbQgW!z9z=8hA5D)C+yj&ssa?JG}5QuIAT|Fsd@|4US0#lTl%%^_x&J3vez8g92 zbXi}V{<~^C-bfDg>5d`a?+p3S8S<$Ui+jy19>aK_y=T+PFN-@ru`5FArLZQSl+i}U zInt(UHgJRLyTQDs+huQkS(v&d?sMB(f9STEV{dx5TRtZSeQt<|aW`a~B2EAjBa z$aM%t9Sb069!JPP3$^fVQL(P{%jJ#OWy!H2?9NoZ%IfcpiB13AvcRM-V|;TF$6j2; zUACKVt?E9B+KBHGBfhIUkL|mkpZ@!QZDl5Z@qN?0aecWtUt5k5X&jWIGAI#%mfvqtDUVnnwB^&ne$;jF9URz26Bo|LW{w& zqESKS2asJNVgv$6yClgH`IBK{vk}arJCANOu7Qk2-M**&uKl=dCSAYpcyji`8Yf~7 zksJQ(g-M{odc8eA-&%45I6=TbyM#wmJ0O7~0aC0(zPzWt0OhB{Ld_8}e70Hv|I1j6 znY+a$MJhwV zG_TxTuD92=5d{~EH5_422sJQ5fCED$EXey%4cIdb6NwCin7cqu24xC5psRAltPCu4 zZ6sBAO>~hCPcnX{)Dte?CYr{NHc@qGES19AI5U!1m6Zq|%OTAe@^DJdGe zA`>y7r?fZ$|Ml{af}*mnS-+8}7fBN|3y#ZUCgreb*MyP0MxRAmSJDnFm9BY~s(w%h zxx^K%uns56mtv}cin;ejoO9}RdZ!!$75?YNX`sU9`r>-Gv9SMJU_prg0R#5tzuOh| zHR9~`-)bW@6QX9@v!2K|WSqElTQg#r49JUs7~rO2S}g3pZT>Kz&M55v>jMt?`}s6b z;qu~YcVQh;1fkdx>5PC-gSP_2F!`760b0o{*)LIL5V)24qD+9i6Ej#rn<|4v%y2Zl z`}7zv>w`l+!O|XlCqD6NCnVJo$Ny>R0wbWpqGw9S8}g4P8gjF{x^jk0i8K*G1bV*= zf=6MxJcK%th8{R|isFZGMCeK(JV>8D5WXuWhVzOp(jf)55uz9qJNeP(k&D&?4qg@X zYO?c%JVx8Xk9grYG7|6BLs@5;+x{Og@Rl z=te=?p?z>LbP2jHkjTP{+}UG9)%SBFw)?a)cAKcugCB-VZ5;dEPqVg?AlLA>6P}}x z!cdF(UnYCx?q=sSY45gtwZI2dI*33R5k)KsJP((VZreHL!MHe@+#mgd{_ z2jPCIn`6}?YA`DsveuTN0 z9E5KtTnNU^QzYbaNQ(I~Fy%R7LHQHxmDTm6zvnIoAS&zK+8MbiO&f|q!|qHG??;8-N_7~558bdxUhJ3ndNKqDg$r)GwNF-f+VU!>Qi&Dj0`UWbteR%U_xVxdU0!Z4O=TS{Be6ivj3_w_h<6|Y-^|#q^qwvskdEXw z2|8`kO%4S64aCUQJSF%gGrK?>z-fo&edUM_1`Bs7w128tdBS4rc^zZSRVNO$umQcc zXEQyV)i2z(Ya6Xo{v?sJL0QOdLFk7Z|G#0!pBT5}Zo4r_YAp*KjPEnapZNoUGHw7!8(TA&j~}0zzH71yac1~dt`tUwtS)`NMXXf3V@7^ z1pv7Mt?wK-Rp)E#>Ln+8gI3g! zTVbIqz5F(3(r16Pj;$XVGwE)7b#-l&ycu_rW{C!%F*z_=U@(8-#q5!AWrmwPgd_%L zjnQIY*pcD9?1{qYosnpEBU8=5Nl8(A2&7zKK08F$i!2zyl(HuI zww~dX!#Lqt)x5X!8x0btd%-C-G!k70WjhyjJJIbU9)g}}9FO0JKR&G$LSDALxwL~# zGRq{xoF%p~Z}*Q|-!cs?7`aZ|-3*S?6!pC+&ptG*6>@WBl&r**Fw$slnCW72E_%iq4-_m) z7`<}@Q4)|Ts8-O-0Zu4kX(cg;7nA^4B3`OB30x8&R-L~2G*|V%*8g|9k0CZVcWjbe zH(6z)XzZ%5aZ`3R!FJA+fBd1^D?c-4%I#Waa=vWlJW)G2>}50?&2#oLBtA5!{a`Zc zRiwh?j1iO%_lSTb(uoHd$mx$F)DK_S_%nK)j~7Oc$;*Q&fBKj0Q`jaxu-|WDJg|3qf=Z0_Zoqul%zCQPdKV`PUr=0EH z^WoYJe`eTroAV1!yV;%n3Gsk3e!hgqmvSD)o*7w;A^RdDOh%2G7M#dw9+|-lf$@+t zMR?frP#5~nE3kcy$5;D?%QRj~A3c1nV}@4-O3sgdcv8s;>0rHPV)NWfZ}&eeW(h=n z4R1L@)(}~J(KA}1{3MZ6uS%Md^o4+y_t3`gx)#-hSE_4}G^yG@)jumeQgZeUe)v4~ z|9;XM5h=*CSH_HZwMAj^_fvOD^APc=$!U=ZlR87lFDLIuD+B@A>I^1N!91gI_?ngW~;PqRX8~k7g;V z?C>b?4<-%LV1Jb9D>nEH=?KlaA$bU9A(BocTF?fkT; zd({I~<&3#&=HUjd+8efV+!ywH;l7|Mz0+9z<+ zb~4BylO;m9RBo`_J*3rtow=v>PxP_$(Zi25o=A=XUArgx-M@MF_cpV=_Kva&{#>3@Ct1w&y4>3UP>PmD^BW zpp{N!UKTjLwvKD(5Ms*{RkPUV681GG-V;{}W!4EC6Lscr6sJS4ETZDa9h~Zq#po8w z|1{HezETNFncMvCAsvc$h}*y8Z5>nk!@Pd^nSXk|Y04jdbyC8Nb8Ig!?Wt%?uCkSi zG1y3GL8N*v@t7SogJda!>CCJrWO{IQvFj=%dq!<=zCS=AS{9<`NEH`IyYjxOx7pU! z3Y)DBpmTSvuv@on+P_=vEdC%n?Z!=9$4s^U4f|<~4zd6005<%|S0`b^8-X>?*9IHr zcm(5Xy5Vc{%GZR^jQA1900`}6g^-&d>>T4qIt{Unyg4PJajXzfvC~zyhZ)vK?h0Gg z2&vjrXN~TY12|_hXD<~y*0Vgnojg+hy_Vrs%lvN79)_me%6=4x6@wv!{rr7@&~nIc z{^+D)c(VZ+zOrOCxVm(BXnZq&O|e0}2Ebk>2NpQZ2U-;g^K){JZ<5OCRkoFD4q~wDbdH?_) eiwFb&00000{{{d;LjnLB00RI3000000000NZpUE& literal 0 HcmV?d00001 diff --git a/tests/data/genomics/sarscov2/vcf/test.vcf b/tests/data/genomics/sarscov2/vcf/test.vcf new file mode 100644 index 00000000..02a8a1f6 --- /dev/null +++ b/tests/data/genomics/sarscov2/vcf/test.vcf @@ -0,0 +1,41 @@ +##fileformat=VCFv4.2 +##FILTER= +##bcftoolsVersion=1.11+htslib-1.11 +##bcftoolsCommand=mpileup -Ou -f data/genomics/sarscov2/fasta/test_genome.fasta data/genomics/sarscov2/bam/test_paired_end.sorted.bam +##reference=file://data/genomics/sarscov2/fasta/test_genome.fasta +##contig= +##ALT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##FORMAT= +##FORMAT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##bcftools_callVersion=1.11+htslib-1.11 +##bcftools_callCommand=call -Ou -mv; Date=Mon Mar 8 10:55:13 2021 +##FILTER=0"> +##bcftools_filterVersion=1.11+htslib-1.11 +##bcftools_filterCommand=filter -s LowQual -e '%QUAL<1 || DP>0'; Date=Mon Mar 8 10:55:13 2021 +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT data/genomics/sarscov2/bam/test_paired_end.sorted.bam +MT192765.1 197 . G T 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 4788 . C T 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 8236 . C A 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 10506 . TTATGACTGTGTCTCTTTTTGTTACATGCACCATATG TTATG 18.4617 LowQual INDEL;IDV=1;IMF=0.333333;DP=3;VDB=0.5;SGB=-0.379885;MQSB=1;MQ0F=0;ICB=1;HOB=0.5;AC=1;AN=2;DP4=1,1,1,0;MQ=60 GT:PL 0/1:51,0,54 +MT192765.1 11037 . T C 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 15009 . G A 30.4183 LowQual DP=2;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:60,3,0 +MT192765.1 18807 . T C 136 LowQual DP=6;VDB=0.618918;SGB=-0.590765;MQSB=1;MQ0F=0;AC=2;AN=2;DP4=0,0,1,4;MQ=60 GT:PL 1/1:166,15,0 +MT192765.1 23813 . T C 4.38466 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:32,3,0 +MT192765.1 24103 . A G 30.4183 LowQual DP=2;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:60,3,0 diff --git a/tests/data/genomics/sarscov2/vcf/test.vcf.gz b/tests/data/genomics/sarscov2/vcf/test.vcf.gz new file mode 100644 index 0000000000000000000000000000000000000000..9625700417c45a46633ad5bdd03c9b8d4fe6a81b GIT binary patch literal 1312 zcmV+*1>gD~iwFb&00000{{{d;LjnK;1-({lZ`w!{{p|e}qeQE=yTKQ{pwZbY*s~5dGRtCd2jglT zHZcBsCq>NGE%p3piJv72Nh3GO0eqjM)_jjz8x#>qF19qy5*CUJLAVIBgLSbX0>&j3 z@@+9mjbh|W%hx0+7w3d=8oi}yWQdGQ8W}JJEI8dzPScRO%J-{_i|@k)3t^T@wk_nD z2DoQkb~^^vW14Q|&b4}u)l-6alb~4g*QtCDg`2OOk?_7yg-@UQI6z}q$pbp z$`8P{P<1>T&I>E!+2CfP&-Tfha@W*@k2!URG1;CPAEyxu32;`nB>YfxAe@m@A_jOV zD$LYP5AKeJacOv*O1h=|)b^BoV97o~X|WV#8^pO2egccaZ@bC2B$UGteh+)7u}1o02?NBHZpt%;d*rfY(0(cufTT zKD6oW>eu;q5NIVi6trg8t5E~!72R^WB}Ks2;9i6hDme>{x*ukPyc#~rC)iQt({qui z0`!FFaK4;+W$9YnSDevTGHSv~T^2j2w=EHp80e~!DPaqW*}pVqyDW=D-TWx1aR(Ve z*i-AH%22JI_zgrHr)y4Wq+;b~bhBY0gSrb*Ql*~6WkA0buGjMhNJZcwy;Z?|sE@Db zQ*nI#lG^t;o_~b#zM{2PTp1|L_Nt_(_cE6Sh>~!J;NGcj7X3&`kc6U)M*t5zV~AIE zX+NOPn$*r9lOWrJgC$Jid8nQr!Y=BixqYJn}8yP)neS0BTY6tOoV$o1XZum-DH%n6I?)KwI7nH8`T0GNP3( zYAQZjSs1h*YwI5RFb(%EHA5R|fp%%wri0sBMP3^$T-=B9>$Xh8zU(_dAYxX_g zwff!+e&D1+du8faH(}E4n%XG1T1+&2fv@bYZtLb_toEhj0GR*3fH{`kEx^2=0t2`P zVA$-K3M~k{VC4D12);gi)n^2QJ`DSw4}b7buF|k$w7dB7o>kFo_MyE5a?k*}z^(WI z71!=VvjroaC$v}BE-dRqCiv2ULA}#HaR!^VIurqj_xt07 zJEqwyj-aR6rqRZZ{TR*q8MLnX1kG{Gqi+}v$Rng~rEwQKJ?tD|J3SM2)`b2Xdix1F z?sj$DImNeZ2ip~X+pwK>xBKg~wd&qlZ3rcxd7$oZ()K_5lSwHQ4*&okiwFb&00000 W{{{d;LjnLB00RI3000000000;MS$r5 literal 0 HcmV?d00001 diff --git a/tests/data/genomics/sarscov2/vcf/test.vcf.gz.tbi b/tests/data/genomics/sarscov2/vcf/test.vcf.gz.tbi new file mode 100644 index 0000000000000000000000000000000000000000..f6e862edbfb216e8d168096f47b282f554caff4e GIT binary patch literal 112 zcmb2|=3rp}f&Xj_PR>jW!3^AmpHfm%5)u-ak|cPUP6f=8?x=O<(d|8R*5@nF!AtX; w1eU7rRPZYM;FO`i!lLK=fdhw5F_;R)>N5P@`}t5g&?tE{Go=}r!B&6>0LW+{tpET3 literal 0 HcmV?d00001 diff --git a/tests/data/genomics/sarscov2/vcf/test2.vcf b/tests/data/genomics/sarscov2/vcf/test2.vcf new file mode 100644 index 00000000..02a8a1f6 --- /dev/null +++ b/tests/data/genomics/sarscov2/vcf/test2.vcf @@ -0,0 +1,41 @@ +##fileformat=VCFv4.2 +##FILTER= +##bcftoolsVersion=1.11+htslib-1.11 +##bcftoolsCommand=mpileup -Ou -f data/genomics/sarscov2/fasta/test_genome.fasta data/genomics/sarscov2/bam/test_paired_end.sorted.bam +##reference=file://data/genomics/sarscov2/fasta/test_genome.fasta +##contig= +##ALT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##FORMAT= +##FORMAT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##bcftools_callVersion=1.11+htslib-1.11 +##bcftools_callCommand=call -Ou -mv; Date=Mon Mar 8 10:55:13 2021 +##FILTER=0"> +##bcftools_filterVersion=1.11+htslib-1.11 +##bcftools_filterCommand=filter -s LowQual -e '%QUAL<1 || DP>0'; Date=Mon Mar 8 10:55:13 2021 +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT data/genomics/sarscov2/bam/test_paired_end.sorted.bam +MT192765.1 197 . G T 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 4788 . C T 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 8236 . C A 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 10506 . TTATGACTGTGTCTCTTTTTGTTACATGCACCATATG TTATG 18.4617 LowQual INDEL;IDV=1;IMF=0.333333;DP=3;VDB=0.5;SGB=-0.379885;MQSB=1;MQ0F=0;ICB=1;HOB=0.5;AC=1;AN=2;DP4=1,1,1,0;MQ=60 GT:PL 0/1:51,0,54 +MT192765.1 11037 . T C 7.30814 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 15009 . G A 30.4183 LowQual DP=2;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:60,3,0 +MT192765.1 18807 . T C 136 LowQual DP=6;VDB=0.618918;SGB=-0.590765;MQSB=1;MQ0F=0;AC=2;AN=2;DP4=0,0,1,4;MQ=60 GT:PL 1/1:166,15,0 +MT192765.1 23813 . T C 4.38466 LowQual DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:32,3,0 +MT192765.1 24103 . A G 30.4183 LowQual DP=2;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:60,3,0 diff --git a/tests/data/genomics/sarscov2/vcf/test2.vcf.gz b/tests/data/genomics/sarscov2/vcf/test2.vcf.gz new file mode 100644 index 0000000000000000000000000000000000000000..9625700417c45a46633ad5bdd03c9b8d4fe6a81b GIT binary patch literal 1312 zcmV+*1>gD~iwFb&00000{{{d;LjnK;1-({lZ`w!{{p|e}qeQE=yTKQ{pwZbY*s~5dGRtCd2jglT zHZcBsCq>NGE%p3piJv72Nh3GO0eqjM)_jjz8x#>qF19qy5*CUJLAVIBgLSbX0>&j3 z@@+9mjbh|W%hx0+7w3d=8oi}yWQdGQ8W}JJEI8dzPScRO%J-{_i|@k)3t^T@wk_nD z2DoQkb~^^vW14Q|&b4}u)l-6alb~4g*QtCDg`2OOk?_7yg-@UQI6z}q$pbp z$`8P{P<1>T&I>E!+2CfP&-Tfha@W*@k2!URG1;CPAEyxu32;`nB>YfxAe@m@A_jOV zD$LYP5AKeJacOv*O1h=|)b^BoV97o~X|WV#8^pO2egccaZ@bC2B$UGteh+)7u}1o02?NBHZpt%;d*rfY(0(cufTT zKD6oW>eu;q5NIVi6trg8t5E~!72R^WB}Ks2;9i6hDme>{x*ukPyc#~rC)iQt({qui z0`!FFaK4;+W$9YnSDevTGHSv~T^2j2w=EHp80e~!DPaqW*}pVqyDW=D-TWx1aR(Ve z*i-AH%22JI_zgrHr)y4Wq+;b~bhBY0gSrb*Ql*~6WkA0buGjMhNJZcwy;Z?|sE@Db zQ*nI#lG^t;o_~b#zM{2PTp1|L_Nt_(_cE6Sh>~!J;NGcj7X3&`kc6U)M*t5zV~AIE zX+NOPn$*r9lOWrJgC$Jid8nQr!Y=BixqYJn}8yP)neS0BTY6tOoV$o1XZum-DH%n6I?)KwI7nH8`T0GNP3( zYAQZjSs1h*YwI5RFb(%EHA5R|fp%%wri0sBMP3^$T-=B9>$Xh8zU(_dAYxX_g zwff!+e&D1+du8faH(}E4n%XG1T1+&2fv@bYZtLb_toEhj0GR*3fH{`kEx^2=0t2`P zVA$-K3M~k{VC4D12);gi)n^2QJ`DSw4}b7buF|k$w7dB7o>kFo_MyE5a?k*}z^(WI z71!=VvjroaC$v}BE-dRqCiv2ULA}#HaR!^VIurqj_xt07 zJEqwyj-aR6rqRZZ{TR*q8MLnX1kG{Gqi+}v$Rng~rEwQKJ?tD|J3SM2)`b2Xdix1F z?sj$DImNeZ2ip~X+pwK>xBKg~wd&qlZ3rcxd7$oZ()K_5lSwHQ4*&okiwFb&00000 W{{{d;LjnLB00RI3000000000;MS$r5 literal 0 HcmV?d00001 diff --git a/tests/data/genomics/sarscov2/vcf/test2.vcf.gz.tbi b/tests/data/genomics/sarscov2/vcf/test2.vcf.gz.tbi new file mode 100644 index 0000000000000000000000000000000000000000..f6e862edbfb216e8d168096f47b282f554caff4e GIT binary patch literal 112 zcmb2|=3rp}f&Xj_PR>jW!3^AmpHfm%5)u-ak|cPUP6f=8?x=O<(d|8R*5@nF!AtX; w1eU7rRPZYM;FO`i!lLK=fdhw5F_;R)>N5P@`}t5g&?tE{Go=}r!B&6>0LW+{tpET3 literal 0 HcmV?d00001 diff --git a/tests/data/genomics/sarscov2/vcf/test3.vcf b/tests/data/genomics/sarscov2/vcf/test3.vcf new file mode 100644 index 00000000..a62ed968 --- /dev/null +++ b/tests/data/genomics/sarscov2/vcf/test3.vcf @@ -0,0 +1,42 @@ +##fileformat=VCFv4.2 +##FILTER= +##bcftoolsVersion=1.11+htslib-1.11 +##bcftoolsCommand=mpileup -Ou -f ../fasta/test_genome.fasta ../bam/test_single_end.sorted.bam +##reference=file://../fasta/test_genome.fasta +##contig= +##ALT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##FORMAT= +##FORMAT= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##INFO= +##bcftools_callVersion=1.11+htslib-1.11 +##bcftools_callCommand=call -Ou -mv; Date=Mon Mar 8 13:13:34 2021 +##FILTER=1"> +##bcftools_filterVersion=1.11+htslib-1.11 +##bcftools_filterCommand=filter -s LowQual -e '%QUAL<1 || DP>1'; Date=Mon Mar 8 13:13:34 2021 +#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT ../bam/test_single_end.sorted.bam +MT192765.1 3744 . A T 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=36 GT:PL 1/1:36,3,0 +MT192765.1 4788 . C T 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 8236 . C A 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 10506 . TTATGACTGTGTCTCTTTTTGTTACATGCACCATATG TTATG 21.4353 LowQual INDEL;IDV=1;IMF=0.5;DP=2;VDB=0.78;SGB=-0.379885;MQSB=1;MQ0F=0;ICB=1;HOB=0.5;AC=1;AN=2;DP4=0,1,1,0;MQ=60 GT:PL 0/1:54,0,54 +MT192765.1 11037 . T C 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 15009 . G A 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 15636 . AG A 8.99861 PASS INDEL;IDV=1;IMF=1;DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,0,1;MQ=60 GT:PL 1/1:38,3,0 +MT192765.1 18807 . T C 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 +MT192765.1 18820 . TCC TC 6.51141 PASS INDEL;IDV=1;IMF=1;DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:35,3,0 +MT192765.1 24103 . A G 7.30814 PASS DP=1;SGB=-0.379885;MQ0F=0;AC=2;AN=2;DP4=0,0,1,0;MQ=60 GT:PL 1/1:36,3,0 diff --git a/tests/data/genomics/sarscov2/vcf/test3.vcf.gz b/tests/data/genomics/sarscov2/vcf/test3.vcf.gz new file mode 100644 index 0000000000000000000000000000000000000000..21e819ded311af93bcd9408e3e21f07711978a5c GIT binary patch literal 1287 zcmV+i1^D_OiwFb&00000{{{d;LjnNm1ie;mbJ{o%{_Oh|nq+R;-hpF*0Tb%ZIHnM0 zcni?n$BwWU(41eWB+_=KKmJy>NrICm^wKtDz_O%$o@ceY+SAh&k6@LGm?`(#@7>sD z?T^#b-e9y6)7Q&nsP!f#6%k)5JaTqKM!w1wxW7lM=#{cerO6tW=st@*3ak zb=0Kv?~RfXUsm<=y(NAc$1Dlmc#Gh>EvZg+q`D%edA4Fwu`>lyeOyD5#$aY|^wg5Y z`Gn-jI)aaogr-b|g3!b%v@T!;0+Ik+t?B&i?ECg`b&w{CuQQ3G1#Q<_&4x)U5hQE1 zacga-*4CoD;UZh|*P#n$uoWO7QIfDIf(TyAHzZXX5JYWQLeixKaf8e<&kV-BNoHj* z?!F&Z#=Ce4!nG=kuUl|?5nCS`A0#0U7;@IO6n+7@VS=+n5sr8X$;^~Zcdz$`X>NFs zC|H9yv^`>9c)W{Ak}W0aiU{_Z2w0Khe7ygeM_HC5%z8gwcyqrXELF2^Q!# zKZwzt9wTE`7$hNV)#lK`b@$?Vme++tyBwndLxx;%CJE+25!src+(tZrs^kg?!#G=^ zAc>sGZ0uDfd?F#k$b7xzYviYoW(f+nzB#GvQ{be2NNqMfg4&2BN%fx%R|$M2ABc`! z?ZcUtJ57;|c4|iV5$FTB(bfNl8?o(HJK<`_BCft32YB%UzzZhf$IwPs^Pi{hMv}P_ zSd@;#p7)P{p2J$enq>#JM0p{JK?xq1WgXTLdDXj@Pq?GVr`Iyp1?U0M-eflN^3pZE ztvLNBGHS?kzf3m*!m12dgmiVuRJesi{9lOpCQUAxT71j&&DL39C70QGTG$9oOaL2&C-bF}-!cy{nIxlabsHzvt~+9M9jw zcu)A+%dQL(q&r>GlUq0E1qfKMA$ac;H%nztNwOG7Dnf)uo)Ml`aq-@vkB_OHi&UZ8 z?gZvAnesqC-yORwmu7o!-+sdW_}211XMK;k0N3mylp%115uh<+tow{KeXbKbGUHyzT&i`pGwT8@|? z#38hPj{kL=)T|mk?4*Y2XI*W|H3CCmMU>dVIr;PIgExFj$+vH$JAFqFx#kzsF>Ldd zd%pRIR3#}0kSdT@C8$@=D)BF8laVo<%#A_Un7!{A7@~$YV&oDHeLhBB7>wuJwY#3n zsJCp}Fb&UG7%j7IIn*|E(`j_4F705Ub*q+HZ?zq#(b1*SwK|^f);iu8zc41TkjW!3^AmpHfm%5)u-ak|cPUP6f=8?x=O<(d|8R*5@nF!AtX; w1eQ*X5l{`UP~0LOB6BP~q2lZz22-I}U51-~_nDOgjgm(*Q<{MpYz2q_08}I%9RL6T literal 0 HcmV?d00001 diff --git a/tests/software/bcftools/consensus/main.nf b/tests/software/bcftools/consensus/main.nf index 8b9ce4fe..325de143 100644 --- a/tests/software/bcftools/consensus/main.nf +++ b/tests/software/bcftools/consensus/main.nf @@ -5,12 +5,11 @@ nextflow.enable.dsl = 2 include { BCFTOOLS_CONSENSUS } from '../../../../software/bcftools/consensus/main.nf' addParams( options: [:] ) workflow test_bcftools_consensus { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf.gz", checkIfExists: true) ], - [ file("${launchDir}/tests/data/vcf/test.vcf.gz.tbi", checkIfExists: true) ], - [ file("${launchDir}/tests/data/vcf/test.consensus.fa", checkIfExists: true) ] ] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf.gz", checkIfExists: true) ], + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf.gz.tbi", checkIfExists: true) ], + [ file("${launchDir}/tests/data/genomics/sarscov2/fasta/test_genome.fasta", checkIfExists: true) ] ] BCFTOOLS_CONSENSUS ( input ) } diff --git a/tests/software/bcftools/consensus/test.yml b/tests/software/bcftools/consensus/test.yml index 26384b3e..9d239bc2 100644 --- a/tests/software/bcftools/consensus/test.yml +++ b/tests/software/bcftools/consensus/test.yml @@ -5,4 +5,4 @@ - bcftools_consensus files: - path: output/bcftools/test.fa - md5sum: c9e7ac4537756a0b33bcf17117f9a065 + md5sum: e57d3e4d113f989bc069c5fd61627091 diff --git a/tests/software/bcftools/filter/main.nf b/tests/software/bcftools/filter/main.nf index bdf21a59..5ca7e47a 100644 --- a/tests/software/bcftools/filter/main.nf +++ b/tests/software/bcftools/filter/main.nf @@ -6,10 +6,9 @@ nextflow.enable.dsl = 2 include { BCFTOOLS_FILTER } from '../../../../software/bcftools/filter/main.nf' addParams( options: ['args': '--no-version'] ) workflow test_bcftools_filter { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf", checkIfExists: true) ]] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf", checkIfExists: true) ]] BCFTOOLS_FILTER ( input ) } diff --git a/tests/software/bcftools/filter/test.yml b/tests/software/bcftools/filter/test.yml index a47d3cd5..b1dfdfba 100644 --- a/tests/software/bcftools/filter/test.yml +++ b/tests/software/bcftools/filter/test.yml @@ -5,4 +5,4 @@ - bcftools_filter files: - path: output/bcftools/test.vcf.gz - md5sum: 16947ce72a127938d881113a1e6e696b + md5sum: 9d491cfa84067450342ba8e66c75e5b8 diff --git a/tests/software/bcftools/isec/main.nf b/tests/software/bcftools/isec/main.nf index 7d22f437..ab214019 100644 --- a/tests/software/bcftools/isec/main.nf +++ b/tests/software/bcftools/isec/main.nf @@ -5,16 +5,13 @@ nextflow.enable.dsl = 2 include { BCFTOOLS_ISEC } from '../../../../software/bcftools/isec/main.nf' addParams( options: ['args': '--nfiles +2 --output-type z --no-version'] ) workflow test_bcftools_isec { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf.gz", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test2.vcf.gz", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test3.vcf.gz", checkIfExists: true)], - [ file("${launchDir}/tests/data/vcf/test.vcf.gz.tbi", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test2.vcf.gz.tbi", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test3.vcf.gz.tbi", checkIfExists: true) ]] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf.gz", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf.gz", checkIfExists: true)], + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf.gz.tbi", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf.gz.tbi", checkIfExists: true)] + ] BCFTOOLS_ISEC ( input ) - } diff --git a/tests/software/bcftools/isec/test.yml b/tests/software/bcftools/isec/test.yml index 8955521e..9fe06207 100644 --- a/tests/software/bcftools/isec/test.yml +++ b/tests/software/bcftools/isec/test.yml @@ -5,18 +5,14 @@ - bcftools_isec files: - path: output/bcftools/test/0000.vcf.gz - md5sum: 4e84e3b6903fa44e8bd1acdeff9b265d + md5sum: fc178eb342a91dc0d1d568601ad8f8e2 - path: output/bcftools/test/0000.vcf.gz.tbi - md5sum: 8484b151ef902e25e54f7713d46ed90e + md5sum: 36e11bf96ed0af4a92caa91a68612d64 - path: output/bcftools/test/0001.vcf.gz - md5sum: 4e84e3b6903fa44e8bd1acdeff9b265d + md5sum: fc178eb342a91dc0d1d568601ad8f8e2 - path: output/bcftools/test/0001.vcf.gz.tbi - md5sum: 8484b151ef902e25e54f7713d46ed90e - - path: output/bcftools/test/0002.vcf.gz - md5sum: 4e84e3b6903fa44e8bd1acdeff9b265d - - path: output/bcftools/test/0002.vcf.gz.tbi - md5sum: 8484b151ef902e25e54f7713d46ed90e + md5sum: 36e11bf96ed0af4a92caa91a68612d64 - path: output/bcftools/test/README.txt - md5sum: 63ef64134d2685bc5d50332ef20389d2 + md5sum: 10fc33b66522645600d44afbd41fb792 - path: output/bcftools/test/sites.txt - md5sum: 01bb949ed7825ecf692bf0640e363647 + md5sum: 1cea3fbde7f6d3c97f3d39036f9690df diff --git a/tests/software/bcftools/merge/main.nf b/tests/software/bcftools/merge/main.nf index d9fe395d..a466e387 100644 --- a/tests/software/bcftools/merge/main.nf +++ b/tests/software/bcftools/merge/main.nf @@ -6,13 +6,12 @@ nextflow.enable.dsl = 2 include { BCFTOOLS_MERGE } from '../../../../software/bcftools/merge/main.nf' addParams( options: ['args': '--force-samples --no-version'] ) workflow test_bcftools_merge { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test_a.merge.vcf.gz", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test_b.merge.vcf.gz", checkIfExists: true) ], - [ file("${launchDir}/tests/data/vcf/test_a.merge.vcf.gz.tbi", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test_b.merge.vcf.gz.tbi", checkIfExists: true) ]] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf.gz", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test3.vcf.gz", checkIfExists: true) ], + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf.gz.tbi", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test3.vcf.gz.tbi", checkIfExists: true) ]] BCFTOOLS_MERGE ( input ) } diff --git a/tests/software/bcftools/merge/test.yml b/tests/software/bcftools/merge/test.yml index 471ea52b..d31ed1b9 100644 --- a/tests/software/bcftools/merge/test.yml +++ b/tests/software/bcftools/merge/test.yml @@ -5,4 +5,4 @@ - bcftools_merge files: - path: output/bcftools/test.vcf.gz - md5sum: 23eee1acfd99581889be3fd2b86b6950 + md5sum: 52165fc4149c61547e63800b876c9661 diff --git a/tests/software/bcftools/mpileup/main.nf b/tests/software/bcftools/mpileup/main.nf index 1c28df55..fed52804 100644 --- a/tests/software/bcftools/mpileup/main.nf +++ b/tests/software/bcftools/mpileup/main.nf @@ -6,12 +6,10 @@ include { BCFTOOLS_MPILEUP } from '../../../../software/bcftools/mpileup/main.nf 'args3': '--no-version' ] ) workflow test_bcftools_mpileup { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/bam/test-sc2-artic-v3.bam", checkIfExists: true) ]] - fasta = [ file("${launchDir}/tests/data/fasta/sarscov2/MN908947.3.fa", checkIfExists: true) ] + [ file("${launchDir}/tests/data/genomics/sarscov2/bam/test_paired_end.sorted.bam", checkIfExists: true) ]] + fasta = [ file("${launchDir}/tests/data/genomics/sarscov2/fasta/test_genome.fasta", checkIfExists: true) ] BCFTOOLS_MPILEUP ( input, fasta ) - } diff --git a/tests/software/bcftools/mpileup/test.yml b/tests/software/bcftools/mpileup/test.yml index 0c55f614..3703286a 100644 --- a/tests/software/bcftools/mpileup/test.yml +++ b/tests/software/bcftools/mpileup/test.yml @@ -5,9 +5,8 @@ - bcftools_mpileup files: - path: output/bcftools/test.vcf.gz - md5sum: f25141161cb935ef7719cfb41563ffb0 + md5sum: 518f770a9d04004c4493fc2167280653 - path: output/bcftools/test.vcf.gz.tbi - md5sum: 4423be9f13fa07e5eadb13b1cbe8baf7 + md5sum: f4be4d707434b3b960fde100b7ec3bf3 - path: output/bcftools/test.bcftools_stats.txt - md5sum: b8d6b6927f7bf024acf11ba438b54fa0 - + md5sum: 2d506e32837a53a01fea0fc90402632a diff --git a/tests/software/bcftools/stats/main.nf b/tests/software/bcftools/stats/main.nf index ee7833da..829e3ca8 100644 --- a/tests/software/bcftools/stats/main.nf +++ b/tests/software/bcftools/stats/main.nf @@ -5,10 +5,9 @@ nextflow.enable.dsl = 2 include { BCFTOOLS_STATS } from '../../../../software/bcftools/stats/main.nf' addParams( options: [:] ) workflow test_bcftools_stats { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf", checkIfExists: true) ]] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf", checkIfExists: true) ]] BCFTOOLS_STATS ( input ) } diff --git a/tests/software/bcftools/stats/test.yml b/tests/software/bcftools/stats/test.yml index fdab4904..fa097efc 100644 --- a/tests/software/bcftools/stats/test.yml +++ b/tests/software/bcftools/stats/test.yml @@ -5,4 +5,4 @@ - bcftools_stats files: - path: output/bcftools/test.bcftools_stats.txt - md5sum: abfc6a90f84e24b2cc7e92cbce06200a + md5sum: c4c5938add12a20050eec3782c8ad623 diff --git a/tests/software/gatk4/mergevcfs/main.nf b/tests/software/gatk4/mergevcfs/main.nf index 48b9b15e..49763ae5 100644 --- a/tests/software/gatk4/mergevcfs/main.nf +++ b/tests/software/gatk4/mergevcfs/main.nf @@ -5,27 +5,23 @@ nextflow.enable.dsl = 2 include { GATK4_MERGEVCFS } from '../../../../software/gatk4/mergevcfs/main.nf' addParams( options: [:] ) workflow test_gatk4_mergevcfs { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test2.vcf.gz", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test3.vcf.gz", checkIfExists: true) ] ] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf", checkIfExists: true) ] ] - ref_dict = file("tests/data/fasta/test.consensus.for_vcf.dict", checkIfExists: true) + ref_dict = file('tests/data/genomics/sarscov2/fasta/test_genome.dict', checkIfExists: true) GATK4_MERGEVCFS ( input, ref_dict, false ) } workflow test_gatk4_mergevcfs_refdict { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/vcf/test.vcf", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test2.vcf.gz", checkIfExists: true), - file("${launchDir}/tests/data/vcf/test3.vcf.gz", checkIfExists: true) ] ] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf", checkIfExists: true), + file("${launchDir}/tests/data/genomics/sarscov2/vcf/test2.vcf", checkIfExists: true) ] ] - ref_dict = file("tests/data/fasta/test.consensus.for_vcf.dict", checkIfExists: true) + ref_dict = file('tests/data/genomics/sarscov2/fasta/test_genome.dict', checkIfExists: true) GATK4_MERGEVCFS ( input, ref_dict, true ) } diff --git a/tests/software/gatk4/mergevcfs/test.yml b/tests/software/gatk4/mergevcfs/test.yml index 4f053cc1..2bd2e50c 100644 --- a/tests/software/gatk4/mergevcfs/test.yml +++ b/tests/software/gatk4/mergevcfs/test.yml @@ -5,7 +5,7 @@ - gatk4_mergevcfs files: - path: output/gatk4/test.merged.vcf.gz - md5sum: f25850b7bd4d362b5ea67d4453e9df55 + md5sum: ff48f175e26db2d4b2957762f6d1c715 - name: gatk4 mergevcfs refdict command: nextflow run ./tests/software/gatk4/mergevcfs -entry test_gatk4_mergevcfs_refdict -c tests/config/nextflow.config @@ -14,4 +14,4 @@ - gatk4_mergevcfs files: - path: output/gatk4/test.merged.vcf.gz - md5sum: f25850b7bd4d362b5ea67d4453e9df55 + md5sum: ff48f175e26db2d4b2957762f6d1c715 diff --git a/tests/software/tabix/bgzip/main.nf b/tests/software/tabix/bgzip/main.nf index 9e32e24a..7042a283 100644 --- a/tests/software/tabix/bgzip/main.nf +++ b/tests/software/tabix/bgzip/main.nf @@ -5,10 +5,9 @@ nextflow.enable.dsl = 2 include { TABIX_BGZIP } from '../../../../software/tabix/bgzip/main.nf' addParams( options: [:] ) workflow test_tabix_bgzip { - def input = [] input = [ [ id:'test' ], // meta map - [ file("${launchDir}/tests/data/generic/vcf/test.vcf", checkIfExists: true) ]] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf", checkIfExists: true) ]] TABIX_BGZIP ( input ) } diff --git a/tests/software/tabix/bgzip/test.yml b/tests/software/tabix/bgzip/test.yml index d75918e9..7e284124 100644 --- a/tests/software/tabix/bgzip/test.yml +++ b/tests/software/tabix/bgzip/test.yml @@ -5,4 +5,4 @@ - tabix_bgzip files: - path: ./output/tabix/test.vcf.gz - md5sum: 40419fb7562475d1c8ec4ab725796de2 + md5sum: 0f1c94af3aa3e7e203d9e034ef6f8f4d diff --git a/tests/software/tabix/tabix/main.nf b/tests/software/tabix/tabix/main.nf index d89f6aee..45aa9b4d 100644 --- a/tests/software/tabix/tabix/main.nf +++ b/tests/software/tabix/tabix/main.nf @@ -25,7 +25,7 @@ workflow test_tabix_tabix_gff { workflow test_tabix_tabix_vcf { def input = [] input = [ [ id:'test.vcf' ], // meta map - [ file("${launchDir}/tests/data/generic/vcf/test.vcf.gz", checkIfExists: true) ] ] + [ file("${launchDir}/tests/data/genomics/sarscov2/vcf/test.vcf.gz", checkIfExists: true) ] ] TABIX_VCF ( input ) } diff --git a/tests/software/tabix/tabix/test.yml b/tests/software/tabix/tabix/test.yml index 288ef78c..6ee3bbb8 100644 --- a/tests/software/tabix/tabix/test.yml +++ b/tests/software/tabix/tabix/test.yml @@ -24,4 +24,4 @@ - vcf files: - path: output/tabix/test.vcf.gz.tbi - md5sum: a03f56d3e968f32256ffb9f6b7d01812 + md5sum: bbec39fd53cf2834909d52094980d094