search

maickrau / GraphAligner

MIT License
250 stars 27 forks source link

GAM surject incompatible with VG #82

Open jiadong324 opened 1 year ago

jiadong324 commented 1 year ago

Dear developer,

I want to convert the GAM to BAM format via VG surject command, but it failed because of the different attributes between GraphAligner and vg aligner. The differences are shown below:

VG aligner:

image

GraphAligner:

image

It looks like vg surject uses node_id for conversion, and this node_id is the one in the GFA graph. Here, GraphAligner assigns a different value to node_id, as a result, the error shows Attempted to get handle for node 2797 not present in graph.

I am wondering if it convenient for you to change the output format, so that it is compatible with VG's GAM.

Best, Jiadong

maickrau commented 1 year ago

Hi, could you share the graph and the read?