Creating assembly graphs .gfa
3
1
Entering edit mode
4.9 years ago

Hello all,

Is it possible to create a .gfa file out of a fasta file?

This is outputted by some assembly programs and is great for assembly visualization in bandage. But I have now scaffolded my assembly and would like to visualize that as well.

Cheers,
Ricardo

assembly gfa graph bandage visualization • 6.7k views
ADD COMMENT
0
Entering edit mode

I have the same question. It would be great to get some tips on how to do this.

ADD REPLY
2
Entering edit mode
4.8 years ago
Asaf 10k

These files are generated in the assembly process using software like e.g. flye. Usually assemblers do the process of reads -> assembly graph -> assembly (contigs or scaffolds). The graphs are intermediates and represent repeats and branches, basically uncertainties in the assembly process. You can't go back from an assembly to a graph.

ADD COMMENT
0
Entering edit mode

What if I try to do an assembly using my scaffolds as reads. Wont it build an assembly graph for them?

ADD REPLY
0
Entering edit mode
4.1 years ago
tothepoint ▴ 940

Yes, by SibeliaZ you can get .gfa file out of a fasta file. They recently updated there old version and now large genome file can also be used for the same.

ADD COMMENT
0
Entering edit mode
10 months ago
Adam Taranto ▴ 40

You might want to convert a FASTA file to GFA format to edit it in a GFA-viz tool like Bandage.

You can use fa2gfa from GraphTagger to generate a basic gfa with segment length tags. You can add additional tags such as depth "DP" using csv2tag.

This graph won't have any links between contigs like you might see in a GFA output by an assembler. If you have long-reads available you can update the basic gfa with links using tools such as ntLINK or SGTK (skip scaffolding steps to just update the links if you want to manually curate the graph before scaffolding). You can also manually add link lines to a GFA using the format from the GFA-Spec.

ADD COMMENT

Login before adding your answer.

Traffic: 1817 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6