<p>The PATRICBrc.org is a free resource you can annotate your genome and compare it with publicly available genomes.</p>

<p>Refer this <a href="http://docs.patricbrc.org/tutorial/genome_annotation/annotation.html" rel="nofollow noreferrer" title="tutorial">tutorial</a> for more detail.</p>

<p>Once you have annotated your genome in PATRIC, you will have PATRIC Local families (k-mer based protein family), which will help you identify all genes with a specific function.</p>

CC BY-SA 3.0
2017-11-27T02:59:40.593Z
2017-11-28T04:41:22.757Z
2017-11-28T04:41:22.757Z
Harry Yoo
<p>Im not 100% on what you are looking for - you seem to jump from annotation to detection of genes in your question. </p>

<p>But if you have a database of genes you want to detect (presence/ absence) you could try <a href="https://github.com/katholt/srst2" rel="nofollow noreferrer">srst2</a>
which basically takes your raw reads (fastq) and maps them to the database.  it will give you a report on any hits, and any snps/ indels in the hit. The only problem is making you own database is difficult as you must adhere to the <a href="https://github.com/katholt/srst2#generating-srst2-compatible-clustered-database-from-raw-sequences" rel="nofollow noreferrer">fasta header format</a>.</p>

CC BY-SA 3.0
2017-11-27T16:30:28.680Z
2017-11-27T16:30:28.680Z
AudileF
<p>I'm trying to annotate a genome to find all genes with a specific function. I have a FASTA and the read FASTQs - I'd like to assign the functional group of the identified proteins (e.g. Kegg orthology) automatically.</p>

<p>For more context, I have whole genome sequences of bacteria - I'm trying to find all genes involved in membrane transport. I'd either like a tool to retrieve gene ontologies from a list of gene names (input: genA output: Membrane Transport), or a tool that outputs the function as part of the annotation pipeline (e.g. genbank output with GO function defined). </p>

<p>Or posed in a different way, how would you identify all genes with a specific function from prokaryotic whole genome sequence data?</p>

<p>Are there any scripts that can do this?</p>

CC BY-SA 3.0
2017-11-26T12:09:47.690Z
Please [edit] your question and give us more details about what you have and what you need. Do you have gene models for the bacterium? I mean, do you first have to find the genes, then the gene functions, or do you already have the genes identified? You have FASTA and FASTQ of what? The entire genome? Raw reads? Are they assembled? Is there a reference sequence for the bacterium?
terdon
CC BY-SA 3.0
2017-11-26T14:04:46.727Z
Thanks for your reply - edited in a few more details. I have the raw read data, and relative de novo assemblies of whole genome prokaryotes. I have a few species I'm looking at, a few E. coli, Staph, and Yersinia. Ref seqs available for all, but I'll be annotating draft clinical genomes.

I'm looking for a way to find all genes present in an assembly that match a specific function (e.g. membrane transport). In the genbank format for example I see "/function=" which is exactly what I want to assign - but I can't find a tool that writes that data. Prokka for e.g. does not write this by default
MichaelKirst
<bacteria><sequence-annotation>
Gene function annotation - bacterial genome
# Gene function annotation - bacterial genome

(Asked by: username_0 on 2017-11-26)

I'm trying to annotate a genome to find all genes with a specific function. I have a FASTA and the read FASTQs - I'd like to assign the functional group of the identified proteins (e.g. Kegg orthology) automatically.

For more context, I have whole genome sequences of bacteria - I'm trying to find all genes involved in membrane transport. I'd either like a tool to retrieve gene ontologies from a list of gene names (input: genA output: Membrane Transport), or a tool that outputs the function as part of the annotation pipeline (e.g. genbank output with GO function defined).

Or posed in a different way, how would you identify all genes with a specific function from prokaryotic whole genome sequence data?

Are there any scripts that can do this?

**Tags:** bacteria, sequence-annotation
> **Comments:**
> **username_1**: Please [edit] your question and give us more details about what you have and what you need. Do you have gene models for the bacterium? I mean, do you first have to find the genes, then the gene functions, or do you already have the genes identified? You have FASTA and FASTQ of what? The entire genome? Raw reads? Are they assembled? Is there a reference sequence for the bacterium?  
> **username_0**: Thanks for your reply - edited in a few more details. I have the raw read data, and relative de novo assemblies of whole genome prokaryotes. I have a few species I'm looking at, a few E. coli, Staph, and Yersinia. Ref seqs available for all, but I'll be annotating draft clinical genomes.  I'm looking for a way to find all genes present in an assembly that match a specific function (e.g. membrane transport). In the genbank format for example I see "/function=" which is exactly what I want to assign - but I can't find a tool that writes that data. Prokka for e.g. does not write this by default  

---

## 2 answers

### Answer by username_2 on 2017-11-27. Score: 4

The PATRICBrc.org is a free resource you can annotate your genome and compare it with publicly available genomes.

Refer this [tutorial](http://docs.patricbrc.org/tutorial/genome_annotation/annotation.html "tutorial") for more detail.

Once you have annotated your genome in PATRIC, you will have PATRIC Local families (k-mer based protein family), which will help you identify all genes with a specific function.


---

### Answer by username_3 on 2017-11-27. Score: 1

Im not 100% on what you are looking for - you seem to jump from annotation to detection of genes in your question.

But if you have a database of genes you want to detect (presence/ absence) you could try [srst2](https://github.com/katholt/srst2)
which basically takes your raw reads (fastq) and maps them to the database. it will give you a report on any hits, and any snps/ indels in the hit. The only problem is making you own database is difficult as you must adhere to the [fasta header format](https://github.com/katholt/srst2#generating-srst2-compatible-clustered-database-from-raw-sequences).


---

