Open biozzq opened 2 years ago
vg call
can use any paths in the graph as references via the -p
option. But cycles in the reference path (which PGGB can produce) will be collapsed in the VCF as well, which makes them hard to interpret in vg call
's output (ie wrong).
For example, if there are two variants on path chr1
chr1 10 A T
chr1 20 A G
but positions 10 and 20 are on the same node in the graph (ie a cycle), then vg call
will just collapse them into the first position found
chr1 10 A T,G
I guess in theory if you're calling with a GBWT call -g
it would be possible to port over the unfolding logic from deconstruct
to resolve cycles -- but we haven't done that yet.
One thing you might have more luck with now is odgi untangle
Dear all,
When learning how to detect SVs based on the
pggb
constructed graph genome usingvg giraffe
, I found that the results are a bit strange. And this may be due to the absense of the backbone reference genome when running vg. I wonder that if we can detect SVs on the specified paths. Also, if we could detect SVs on the specified paths, such as the backbone full genome, I think the SV type could be also annotated in the final VCF file. Thank you in advance.Sincerely, Zheng zhuqing