File tree Expand file tree Collapse file tree 1 file changed +9
-1
lines changed
Expand file tree Collapse file tree 1 file changed +9
-1
lines changed Original file line number Diff line number Diff line change @@ -5,7 +5,7 @@ BASH_TAP_ROOT=../deps/bash-tap
55
66PATH=../bin:$PATH # for vg
77
8- plan tests 45
8+ plan tests 46
99
1010vg msga -f GRCh38_alts/FASTA/HLA/V-352962.fa -t 1 -k 16 | vg mod -U 10 - | vg mod -c - > hla.vg
1111vg index hla.vg -x hla.xg
@@ -255,3 +255,11 @@ vg paths --compute-augref --min-augref-len 0 -x cyclic_ref_nested.pg -Q x > cycl
255255vg deconstruct cyclic_ref_nested.augref.pg -p x -a > cyclic_ref_nested.vcf
256256is $( grep -v ^# cyclic_ref_nested.vcf | wc -l) 1 "cyclic reference with nesting produces single variant"
257257rm -f cyclic_ref_nested.pg cyclic_ref_nested.augref.pg cyclic_ref_nested.vcf
258+
259+ # Test 5: Cyclic reference outputs multiple variants (one per reference traversal)
260+ # Reference x visits the snarl twice (via node 2 then 3), alt a visits twice (via node 3 then 4)
261+ # With -c 0 to disable context-jaccard (which doesn't work well on tiny graphs), we should get 2 variants
262+ vg convert -g nesting/cyclic_ref_multiple_variants.gfa -p > cyclic_ref_multi.pg
263+ vg deconstruct cyclic_ref_multi.pg -p x -a -c 0 > cyclic_ref_multi.vcf
264+ is $( grep -v ^# cyclic_ref_multi.vcf | wc -l) 2 "cyclic reference with -a outputs variant for each reference traversal"
265+ rm -f cyclic_ref_multi.pg cyclic_ref_multi.vcf
You can’t perform that action at this time.
0 commit comments