GFA to Fasta file
1
4
Entering edit mode
5.6 years ago

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 • 7.8k views
ADD COMMENT
1
Entering edit mode

Awesome!!

Thank you so much!

ADD REPLY
0
Entering edit mode

Hi,

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

ADD REPLY
0
Entering edit mode

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
11
Entering edit mode
5.6 years ago
lh3 32k
awk '/^S/{print ">"$2"\n"$3}' in.gfa | fold > out.fa
ADD COMMENT

Login before adding your answer.

Traffic: 1363 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