Comments (1)
I think your problem is that autoindex imports the GFA, and then you map to that, and then you re-import the GFA again to make a different graph, and try to interpret the mappings against the second graph. The node IDs are not the same between the two graphs; vg has several ways of importing GFAs and they don't always name the nodes in the result the same way.
You should bge able to use the same .xg file you used for mapping as the input to vg augment
. If it doesn't like that, you can use vg convert
to convert the .xg to .vg (which won't change the node IDs) and then use that file with vg augment
.
You can check if vg thinks a file of reads is interpretable against a particular graph with vg validate [graph] -a [gam]
.
The problem should also go away if you make sure vg doesn't need to do any node ID modifications to your GFA graph: give all the S lines numerical names starting at 1
, and make sure that the edges don't specify any overlaps, and that none of the nodes are too long. 32 bp is definitely a safe limit, and I think 1024 bp might be the real upper limit.
from vg.
Related Issues (20)
- Fastq analysis in euka HOT 2
- Protobuf errors mapping HOT 4
- vg augment crashes in combination with giraffe .gam outputs HOT 5
- when i do the code :vg stats -a var1.gam > var1.stats ,show me a bug HOT 2
- vg surject and pack fail on GraphAligner produced .gam HOT 4
- About the usage of vg deconstruct HOT 5
- Right method to annotate genes on pangenome graph HOT 1
- vg call report missing allele HOT 2
- vg call error HOT 2
- How do different methods affect precision? HOT 4
- Giraffe alignment is very slow and produces warning[vg::Watchdog] messages unless rescue is disabled HOT 13
- Giraffe needs to better deal with systems where multiple processes opening the distance index causes slowdowns
- Any option like BWA MEM '-C'? HOT 1
- vg convert failes while converting gbz to xg with signal 6 error HOT 8
- Release vg v1.60.0
- `vg convert --vg-algorithm` loses start coordinates of paths in its W lines HOT 3
- The Amazing Uninjectable SAM HOT 1
- Throw an error "error: trying to access a snarl tree node of the wrong type" in `haplotypes` HOT 11
- Get input file error HOT 5
- vg pack error
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from vg.