Skip to content
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

GFA output #24

Open
ekg opened this issue Jul 2, 2020 · 1 comment
Open

GFA output #24

ekg opened this issue Jul 2, 2020 · 1 comment

Comments

@ekg
Copy link

ekg commented Jul 2, 2020

Is it possible to obtain GFA output from the assemblies created by this pipeline?

@rchikhi
Copy link
Member

rchikhi commented Jul 2, 2020

I'd say yes at a certain condition: make sure you skip scaffolding. Please use the https://github.com/GATB/minia/blob/master/scripts/convertToGFA.py script on the final assembly fasta file (it should have the L:x:x elements in headers, indicating links between contigs in the FASTA file, which is equivalent info to GFA). When running that script, you'd need to give as input the final k value ran by the pipeline. I never really looked at that pipeline GFA, and would expect that the GFA would become disconnected at coverage gaps due to the pipeline using large k values.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants