File tree 2 files changed +4
-4
lines changed
2 files changed +4
-4
lines changed Original file line number Diff line number Diff line change @@ -14,8 +14,8 @@ vg construct -r small/x.fa -v small/x.vcf.gz > x2.vg
14
14
vg index -x x.xg -G x.gbwt -v small/x.vcf.gz x.vg
15
15
16
16
# List path/thread names from various input formats
17
- is " $( vg paths --list -v x2.vg) " " x" " path listing works from vg"
18
- is " $( vg paths --list -x x.xg) " " x" " path listing works from XG"
17
+ is " $( vg paths --list -v x2.vg | grep -v _alt_ ) " " x" " path listing works from vg"
18
+ is " $( vg paths --list -x x.xg | grep -v _alt_ ) " " x" " path listing works from XG"
19
19
is $( vg paths --list -g x.gbwt | wc -l) 2 " thread listing works from GBWT"
20
20
21
21
# Select threads by name
Original file line number Diff line number Diff line change @@ -42,8 +42,8 @@ is $(vg chunk -x x.xg -r 1 -c 0 | vg view - -j | jq .node | grep id | wc -l) 1 "
42
42
43
43
# Check that traces work on a GBWT
44
44
is $( vg chunk -x x.xg -G x.gbwt -r 1:1 -c 2 -T | vg view - -j | jq .node | grep id | wc -l) 5 " id chunker traces correct chunk size"
45
- is " $( vg chunk -x x.xg -r 1:1 -c 2 -T | vg view - -j | jq -c ' .path[] | select(.name != "x")' | wc -l) " 0 " chunker extracts no threads from an empty gPBWT"
46
- is " $( vg chunk -x x.xg -G x.gbwt -r 1:1 -c 2 -T | vg view - -j | jq -c ' .path[] | select(.name != "x")' | wc -l) " 2 " chunker extracts 2 local threads from a gBWT with 2 locally distinct threads in it"
45
+ is " $( vg chunk -x x.xg -r 1:1 -c 2 -T | vg view - -j | jq -c ' .path[] | select(.name != "x")' | grep -v _alt_ | wc -l) " 0 " chunker extracts no threads from an empty gPBWT"
46
+ is " $( vg chunk -x x.xg -G x.gbwt -r 1:1 -c 2 -T | vg view - -j | jq -c ' .path[] | select(.name != "x")' | grep -v _alt_ | wc -l) " 2 " chunker extracts 2 local threads from a gBWT with 2 locally distinct threads in it"
47
47
is " $( vg chunk -x x.xg -G x.gbwt -r 1:1 -c 2 -T | vg view - -j | jq -r ' .path[] | select(.name == "thread_0") | .mapping | length' ) " 3 " chunker can extract a partial haplotype from a GBWT"
48
48
49
49
# check that n-chunking works
You can’t perform that action at this time.
0 commit comments