Question: GFA to Fasta file
3
gravatar for marcela.uliano
5.0 years ago by
European Union
marcela.uliano60 wrote:

Hey guys,

Does anyone have a script (or a awk line) to extract fasta reads from a GFA file?

 

I've runned this miniasm in some Pacbio reads I have, and I would like to extract the final unitigs from the final file I have (which is a GFA file) in the fasta format!

 

Thanks a lot guys!

sequence • 6.3k views
ADD COMMENTlink modified 5.0 years ago • written 5.0 years ago by marcela.uliano60
1

Awesome!!

Thank you so much!

ADD REPLYlink modified 18 months ago by h.mon31k • written 5.0 years ago by marcela.uliano60

Hi,

could you describe what does GFA stand for? I don't know about this format.

ADD REPLYlink modified 11 months ago by _r_am31k • written 5.0 years ago by Atu0

GFA is the Graphical Fragment Assembly format. Here is a post by Heng Li (@lh3) on this: https://lh3.github.io/2014/07/19/a-proposal-of-the-grapical-fragment-assembly-format.

ADD REPLYlink written 2.3 years ago by ropolocan680
10
gravatar for lh3
5.0 years ago by
lh332k
United States
lh332k wrote:
awk '/^S/{print ">"$2"\n"$3}' in.gfa | fold > out.fa
ADD COMMENTlink modified 11 months ago by _r_am31k • written 5.0 years ago by lh332k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2357 users visited in the last hour