This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extracting a protein sequence from fasta file

Hi all,

I am currently working on metagenomics and I would like to know if there's any way that I can extract a protein sequence from a fasta file? Thanks.

Cheers

Alan

gene sequence
• 4,371 views • link
updated by Biostar Biostar • written by h.l.wong h.l.wong

What kind of fasta file (single, multi-fasta, DNA. protein)?

• 0 views • link
written by GenoMax GenoMax

Thanks, I have a metagenomics data, do I extract the sequences from the assembled contigs file? Or do I need to extract the sequence out from other files?

And if I have an annotated protein (in KEGG), how can I get the nucleic acid sequence?

Thanks

• 0 views • link
written by h.l.wong h.l.wong

prodigal or genmarks works fine and fast.

• 0 views • link
written by Buffo Buffo

Thanks, should I use prodigal on the assemble contigs file to extract the sequences?

• 0 views • link
written by h.l.wong h.l.wong

yes, read the manual, but yes, it is possible (nucleotide, protein or both), take care about translation table that you use.

http://prodigal.ornl.gov/
• 0 views • link
written by Buffo Buffo

In case you know the sequence name, use: https://github.com/lh3/bioawk

SEQNAME="<insert sequence name>"; bioawk -v x=$SEQNAME -c fastx '{if ($name==x) {print ">"$name"\n"$seq}}' main.fa > selected.fa

This will select the sequence you want from the main.fa file and print it to the selected.fa file.

• 0 views • link
written by Matteo Schiavinato Matteo Schiavinato

It helps if you elucidate your question. Do you mean extract sequences based on header like this?

• 0 views • link
written by Rohit Rohit

0 answers

No answers yet.

Log in to answer this question.

More posts like this