Comments (5)
Can you sweep that under the rug and identify only the single exon genes for now? That'd be better than what I have now. Either identifying spliced genes as partial genes or not at all would be okay. I'm using other pipelines in parallel to try and identify the spliced genes.
from barrnap.
See Cycas taitungensis mitochondrion.
from barrnap.
Bump?
from barrnap.
Not sure I can do this, as rrn26
is spliced into 2 exons in some plants?
# 2 exons
rRNA complement(join(146852..148823,149542..150368))
/gene="rnl"
/locus_tag="MapooMr01"
/gene_synonym="rrn26"
/product="large subunit ribosomal RNA"
# 1 exon
rRNA complement(126101..129506)
/gene="rrn26"
/locus_tag="CQR31_mgr02"
/product="26S ribosomal RNA"
/db_xref="GeneID:35198937"
from barrnap.
Hi, Torsten. I've found that barrnap --kingdom bac
identifies plant 18S and 26S rRNA just fine. It does not however detect plant mitochondrial 5S rRNA (at least in my conifer species). RFAM 5S_rRNA (RF00001) with Infernal detects it just fine though!
Would you consider adding this RFAM family to Barrnap?
You can download an alignment of the family here so that you can build an HMM for hmmer
.
http://rfam.xfam.org/family/RF00001#tabview=tab2
You could also download the prebuilt Infernal covariance model (CM) here:
http://rfam.xfam.org/family/RF00001#tabview=tab9
The covariance model includes an HMM within the model. Search for the line HMMER3/f [i1.1 | October 2013]
. I don't know whether the file format is compatible with hmmer
. If it is, you could just extract this HMM from the covariance model.
Alternatively, would you consider running infernal
as an option instead of hmmer
? It scans first with an HMM and then checks the hit against the RNA covariance model (at least as I understand it).
5S_rRNA (RF00001) is a member of clan 5S_rRNA (CL00113) which includes only one other family, mtPerm-5S (RF02547). I don't need the permuted form, but thought I'd mention it.
http://rfam.xfam.org/clan/CL00113
http://rfam.xfam.org/family/RF02547
The permuted mitochondrial 5S rRNAs in brown algae represent the most unconventional case, where the closing helix I, which otherwise brings together the molecule's 5′ and 3′ ends, is replaced by a (closed) hairpin resulting in an open three-way junction.
https://en.wikipedia.org/wiki/5S_ribosomal_RNA
Here are the RFAM bacterial family and clans for rrnS and rrnL. The clans include the families for other kingdoms. The bacterial families identify plant mitochondrion rrnS and rrnL.
SSU_rRNA_bacteria (RF00177) http://rfam.org/family/RF00177
SSU (CL00111) http://rfam.org/clan/CL00111
LSU_rRNA_bacteria (RF02541) http://rfam.org/family/RF02541
LSU (CL00112) http://rfam.org/clan/CL00112
from barrnap.
Related Issues (20)
- can barrnap-0.9 use hmmer/3.3 HOT 2
- [Syntax error] Outseq argument
- What does column 8 stand for?
- Can not find required 'nhmmer' HOT 1
- Could the barrnap be used to predict the rrna of *.fna.gz file? HOT 1
- possible option for "prokaryotic" search (archaeal + procaryotic)?
- Overlapping 5.8S and 28S annotations in Eukaryotic rDNAs HOT 4
- Adding fungal HMM HOT 4
- Cannot specify a file on cmdline?
- GFF output and Fasta-headers give different start-coordinates of rRNA-genes
- nhmmer failed to run - Error: Invalid alphabet type in target for nhmmer. Expect DNA or RNA HOT 6
- [barrnap] ERROR: nhmmer failed to run HOT 3
- issue with rrna.fa writing and bed tools HOT 2
- calling two overlapping sequences separate partial 16S sequences
- please add option to force overwrite .fai index during, or automatically delete fai file after barrnap run
- barrnap-v0.9 debian package has incomplete database HOT 2
- Fail if Bedtools issues a fasta normalization error
- barrnap doesn't create any output
- barrnap
- Building eukaryotic rrna operon database
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from barrnap.