New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
GAF output #3
Comments
Is this feasible to do? I'd love to test astarix, and a standard output format like this is critical to do that! |
Is this still on the agenda? |
Dear @ekg and @ChriKub, please excuse my delays. I was overoptimistic on implementing the GAF output in time but given my upcoming PhD defense, it may be delayed further. In case someone wants to implement this sooner, there are several subtleties to be taken care of:
|
Could you provide GAF output over the nodes of the graph? Would it be possible to do? I could also submit a patch if you can describe how you'd go about setting this up.
It's defined in the minigraph documentation, and GraphAligner and vg both produce it.
It's like PAF but the target is expressed as a walk through nodes in the graph (your GFA input).
The walk is expressed like
>1>2>4>6>7>9>10>11>12>13>14>16>17>19
If you have
<19<17<16<14<13<12<11<10<9<7<6<4<2<1<
then it might represent the reverse complement of the above walk.You'll need to create a CIGAR or md tag to express the base-level alignment, for use downstream in
vg call
for instance.The text was updated successfully, but these errors were encountered: