search

GATB / gatb-minia-pipeline

GATB Minia assembly pipeline
29 stars 8 forks source link

GFA output #24

Open ekg opened 4 years ago

ekg commented 4 years ago

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

rchikhi commented 4 years ago

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.