Do vg filter out some variants during the construction of the pan-genome, and if so, what are the criteria for filtering? The number of variations in the pan-genome is reduced compared to the variations in the input VCF file. The command I used is "vg autoindex --workflow giraffe -v sv.vcf -r ref.genome -p sv -t 64" The number of variants decreased by almost half after undergoing vg deconstruct compared to the number of variants in the VCF used to construct the pan-genome graph. What could be the reason for this?
•
written
by
Wenke •
1 answer
If variants overlap in the genome, vg deconstruct will combine them into one locus with multiple alleles. If your input VCF has a lot of deletions, that's probably what the cause is.
Log in to answer this question.
-
How to parallelize running vg autoindex to create graph references from VCFs with vg?
written by Ricky •
-
The number of SVs called by `vg call` is much smaller than the number of SVs in the VCF used to con…
written by Maxine •
-
Comparing SV calls produced by VG call on a Minigraph-Cactus produced Pangenome Graph vs. a linear …
written by peterdfields •
-
How to identify which genome the variations in the pan-genome graph originate from
written by Wenke •
-
Indexing the human pangenome draft
written by lisle.mose •
-
unable to find most of SVs in constructed graph.vg
written by Maxine •
-
Construct a pantranscriptome reference with two haplotypes from a single sample.
written by Juhyun •
-
vg: Mapping fewer reads to a genome graph than a linear reference
written by John Lloyd •
-
Using VG to genotype a GAM using a graph constructed from GFA
written by Michael •
-
Creating a variation graph for Giraffe alignment from assemblies
written by Michael •
