How to do annotation after gene prediction?

[garethboy]

Dear all,

I am pretty new in bioinformatics. I want to ask, after I perform gene prediction using PRODIGAL, I got the following:

>final_001_R1_00001 # 2 # 898 # 1 # ID=1_1;partial=10;start_type=Edge;rbs_motif=None;rbs_spacer=None;gc_cont=0.644
ATGCCAGCAGGTCGATTGGCCAAGCCACTGGAGGCAATCGCGGGCCCGCTGCGGCGGTCA
TCGCCGGAGTCCGTTGACCGGATGCTCATCGGCTTGTGCGCGGTCATCTGGTTGGCGTTC
GTCGGGATGCTGGTCGGTGCCATCGTGGTGATCGCCGGAATGGGGGAGGACTCGGGTTCG
TCGGGGTCTTCGTGGGGCATCTACATCGTCATCGGGGTATCGCTGGCCGTCATCGTCGGT
GCCGTGCCGCTGTTGCTACGTGCCCGCAACACCGCTGCCCAGCGCACCTCGAAAACGGCC
CCTAAAACCAAGGGCGCCAATTCCTCCGGTGGGGCCACTGCGCCGTCTGGACCGCCTCCC
GCGGAGGCGTCGACCGAGAAGCTCAAGACCTTCGGTGGTCTCGCCGATCACGTGAGCCAT
GTTCCCAAGGACTACAACGGTCCTGGGTCACCCGCGTACCGGGCAGCCACCGAGGCGGCC
GCGAAGAGGGCGGTGGTCAACCGGCTGTGGTTGCGGGCCACGGTGGGGATCGCCGGCGCG
ATCGGCCTGGCACTGCTGGCGGTGGTGACCGGTACCTATCTGCTGGCAGTGGAAAGCGAC
ACCGCCGCCATCGTCGCGTTCGTCTTCGCCGGGGTGATCACCTGCGCGATGGGGGCCATT
CCATGGGTGACGCTCAAGAAGCTGCGCGAGCTCGGCTCGTCGGGGTCTTAG 

After I perform a local BLAST, how can I get something like the following?

>final_001_R1_00001 hypothetical protein
ATGCCAGCAGGTCGATTGGCCAAGCCACTGGAGGCAATCGCGGGCCCGCTGCGGCGGTCA
TCGCCGGAGTCCGTTGACCGGATGCTCATCGGCTTGTGCGCGGTCATCTGGTTGGCGTTC
GTCGGGATGCTGGTCGGTGCCATCGTGGTGATCGCCGGAATGGGGGAGGACTCGGGTTCG
TCGGGGTCTTCGTGGGGCATCTACATCGTCATCGGGGTATCGCTGGCCGTCATCGTCGGT
GCCGTGCCGCTGTTGCTACGTGCCCGCAACACCGCTGCCCAGCGCACCTCGAAAACGGCC
CCTAAAACCAAGGGCGCCAATTCCTCCGGTGGGGCCACTGCGCCGTCTGGACCGCCTCCC
GCGGAGGCGTCGACCGAGAAGCTCAAGACCTTCGGTGGTCTCGCCGATCACGTGAGCCAT
GTTCCCAAGGACTACAACGGTCCTGGGTCACCCGCGTACCGGGCAGCCACCGAGGCGGCC
GCGAAGAGGGCGGTGGTCAACCGGCTGTGGTTGCGGGCCACGGTGGGGATCGCCGGCGCG
ATCGGCCTGGCACTGCTGGCGGTGGTGACCGGTACCTATCTGCTGGCAGTGGAAAGCGAC
ACCGCCGCCATCGTCGCGTTCGTCTTCGCCGGGGTGATCACCTGCGCGATGGGGGCCATT
CCATGGGTGACGCTCAAGAAGCTGCGCGAGCTCGGCTCGTCGGGGTCTTAG 

In other words, after performing gene prediction, how to get a file which project the sequences as following?

>gene_0001 GENE_NAME_A
SEQUENCE...
>gene_0002 GENE_NAME_B
SEQUENCE...
>gene_0003 GENE_NAME_C
SEQUENCE...
>gene_0004 GENE_NAME_D
SEQUENCE...
>gene_0005 GENE_NAME_E
SEQUENCE...
>gene_0006 GENE_NAME_F
SEQUENCE... 

[Desai]

Dear garethboy

I dont't know how to get sequence like that but yes u can use BLAST2GO tool for mapping and annotation. could I know what is your conclusion?

[boetsie]

Use a tool like Prokka (http://www.vicbioinformatics.com/software.prokka.shtml) to predict and annotate your sequences. After the genbank or other files are generated, write a script to add the annotation to your fastA file.

Boetsie