I'm trying to construct the variation graph as stated in "Working with a whole genome variation graph" in the github wiki of variation graph toolkit and then convert it to gfa for experimental purposes. I manage to create the .vg graph but when i try to convert it to .gfa format with vg view command, an empty gfa file is created and the process just hangs in the terminal with no result.
The exact command i use is: vg view human.vg >human.gfa
The command used for creating .vg is: vg construct -r ref.fa -v vcf_file.vcf.gz >human.vg
Any suggestions on what the problem might be?
•
written
by
Diomfeas •
1 answer
The .vg format is very inefficient. Loading a whole-genome .vg file is both slow and very memory-intensive. We have replacement formats that are much better and are working to (transparently) integrate them throughout the vg toolchain, including GFA conversion. In the meantime, your safest bet would be to construct a graph for each chromosome as described here: https://github.com/vgteam/vg/wiki/Working-with-a-whole-genome-variation-graph. You can convert them to GFA individually and concatenate the results.
Log in to answer this question.
-
Extracting individual chromosomes from a .gfa graph
written by samuel.michel •
-
Annotate GFA with coverage depth
written by elcortegano
-
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 •
-
vg rna index
written by Yumin •
-
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 •
-
cannot see variations on graph image
written by philgev2 •
-
Issues calling variants with vg when using other tools upstream
written by brian.nadon •
