Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: GFA to Fasta file
2
gravatar for marcela.uliano
3.2 years ago by
marcela.uliano • 50
European Union
marcela.uliano • 50 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 • 3.3k views
ADD COMMENT • link •
modified 3.2 years ago • written 3.2 years ago by marcela.uliano • 50
1

Awesome!!

 

Thank you so much!

 

 

ADD REPLY • link written 3.2 years ago by marcela.uliano • 50

Hi,

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

ADD REPLY • link written 3.2 years ago by Atu • 0

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 REPLY • link written 6 months ago by ropolocan • 500
7
gravatar for lh3
3.2 years ago by
lh3 ♦ 31k
United States
lh3 ♦ 31k wrote:
awk '/^S/{print ">"$2"\n"$3}' in.gfa | fold > out.fa
ADD COMMENT • link written 3.2 years ago by lh3 ♦ 31k
Please log in to add an answer.

Similar posts • Search »

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