This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Genome Annotation Output Files Naming

I ran this annotation command: for k in *.fasta; do prokka $k --outdir "$k".prokka.output; echo $k; done on hundreds of assembled sequences. Each sequence annotation output produced a folder named "sequence_name.prokka.output" but the actual files in the folders were named according to the date e.g. "PROKKA_03062020.gff" for all the folders.

Now I want the actual output be named with the actual FASTA file ID. Example: sequence ID of D02502 should be like this: "PROKKA_D02502.gff"

What do I add (or remove) from the code above?

annotation genome
• 2,981 views • link
updated by Mensur Dlakic Mensur Dlakic • written by Fid_o • Fid_o

1 answer 

for k in *.fasta; do prokka $k --outdir "$k".prokka.output --prefix PROKKA_$k; echo $k; done
• 0 views • link
written by Mensur Dlakic Mensur Dlakic

Thank you sooooo much Mensur Dlakic.

I am learning Bioinformatics and your help is great.

• 0 views • link
written by Fid_o • Fid_o

Is there a way to put all prokka output from all fasta files in one folder rather than the traditional output (one folder for the output of each fasta)?

• 0 views • link
written by Ahmed Ezat El Zowalaty • Ahmed Ezat El Zowalaty

Log in to answer this question.

More posts like this