Hello,
I am trying to find why I have to use files like .vg.gbwt and .vg.xg to feed sequenceTubeMap? Wouldn't it be more convenient to use human readable file types like json or txt?
Also, I want to know how I can produce these files when I have something human readable with me.
Thank you for your response!
•
written
by
Mayank •
1 answer
It's mostly a question of efficiency. JSON representations of a eukaryotic pangenome graph are monstrously large. There is a human-readable, text-based format in GFA, which is frequently used as an interchange format between tools. However, naively loading GFA into memory leads to data structures that either are impractically large or don't support the basic queries that you need to compute on a graph. This is especially true for storing a large number of haplotype paths (which are necessary for sequenceTubeMap's visualizations)
You would use the vg toolkit. The vg convert subcommand can convert a GFA to an XG index, and the vg gbwt subcommand can make GBWT indexes.
I was successfully able to read the xg format, but am still having issues with opening up my gbwt file in .gfa or any other readable formats. I suspect I need to edit this file to change the number of tracks in a plot. The menu that pops up when I type vg gbwt isnt very informative on how I can do so. Could you throw some light on this?
GBWT isn't convertible into a GFA on its own. The GBWT index contains only the haplotype paths, so it lacks the node sequences that are required to fully specify the graph. If you want to make a full graph with the GBWT, you should augment it into a GBZ, which is essentially a GBWT with node sequences added. You can do that in vg gbwt as well.
Could you clarify what you mean by the "number of tracks"?
By tracks I mean the thick lines that represent each sequence that go through the various nodes. By being able to change 'number of tracks' I meant increase or decrease these lines while making necessary corresponding changes to the .xg file.
In short: I want to add another sequence (something I think can achieved by editing the gbwt file.)
GBWTs are not particularly easy to edit, but if you have two GBWT files over a graph with the same node IDs, you can combine them with vg gbwt --merge. You can also remove haplotypes with vg gbwt --remove-sample.
This wiki article is probably the best resource. The academic papers on the GBWT are more focused on the underlying algorithmics. https://github.com/vgteam/vg/wiki/VG-GBWT-Subcommand
Log in to answer this question.
-
Making Changes to SequenceTubeMap
written by Mayank •
-
trouble using sequenceTubeMap in chromium installed on WSL
written by ayushm64 •
-
Reverse complement only sequences resulted from BLASTN minus
written by Oscar •
-
Ipyrad de novo for variant calling
written by giulia.trauzzi •
-
bcftools isec: tandem mutations
written by whb •
-
Retrieving GO mappings for GOseq analysis with non-native organism
written by thomas.welch •
-
Intersectbed of Homer and MAcs2 peak caller results
written by bulkybogany •
-
Error with String data base CYTOSCAPE
written by sebas827 •
-
Why annotations are not encoded using standard informatic format ?
written by jsgounot
-
How/when to use cuffmerge and/or cuffcompare?
written by awilli84 •
