|
0
|
1 <tool id="barrnap" name="barrnap" version="1.0.0">
|
|
|
2 <description>Locate bacterial ribosomal RNA's in a fasta file. (GFF output)</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="0.1">package_barrnap_0_1</requirement>
|
|
|
5 </requirements>
|
|
|
6 <command>
|
|
|
7 barrnap
|
|
|
8 --lencutoff $lencutoff
|
|
|
9 --evalue $evalue
|
|
|
10 $fasta_file
|
|
|
11 > $gff
|
|
|
12 2> $stderr
|
|
|
13 </command>
|
|
|
14
|
|
|
15 <inputs>
|
|
|
16 <param name="fasta_file" type="data" format="fasta" label="Fasta file" help="Fasta file to look for rRNA's in" />
|
|
|
17 <param name="lencutoff" type="float" value="0.8" label="Length cutoff" help="Proportional length threshold to tag a find as pseudo" />
|
|
|
18 <param name="evalue" type="float" value="1e-06" label="e-value" help="Similarity e-value cutoff" />
|
|
|
19 </inputs>
|
|
|
20
|
|
|
21 <outputs>
|
|
|
22 <data name="stderr" format="tabular" label="${tool.name} on ${on_string}: STDERR"/>
|
|
|
23 <data name="gff" format="gff3" lavel="${tool.name} on ${on_string}: rRNA GFF"/>
|
|
|
24 </outputs>
|
|
|
25
|
|
|
26 <help>
|
|
|
27 **barrnap**
|
|
|
28 ::
|
|
|
29 barrnap_ predicts the location of 5S, 16S and 23S ribosomal RNA genes in Bacterial genome sequences. It takes FASTA DNA sequence as input, and write GFF3 as output. It uses the new NHMMER tool that comes with HMMER 3.1-dev for HMM searching in DNA:DNA style. NHMMER binaries for 64-bit Linux and Mac OS X are included and will be auto-detected. Multithreading is supported and one can expect roughly linear speed-ups with more CPUs.
|
|
|
30 This tool is designed to be a substitute for RNAmmer for bacteria. It was motivated by my desire to remove Prokka's dependency on RNAmmer which is encumbered by an free-for-academic sign-up license, and by the needed legacy HMMER 2.x which conflicts with HMMER 3.x that most people are using now.
|
|
|
31
|
|
|
32 RNAmmer is more sophisticated than Barrnap, and more accurate. because it uses HMMER 2.x in glocal alignment mode, whereas HMMER 3.x currently only supports local alignment (Sean Eddy expects glocal to be supported in 2014). In practice, Barrnap will find all the typical 5/16/23S genes in bacteria, but may get the end points out by a few bases and will probably miss wierd rRNAs. The HMM models it uses are derived from RFAM, Silva, and GreenGenes.
|
|
|
33
|
|
|
34 Barrnap is for bacteria (prokaryota) only. If you want to find rRNA genes in Archaea or Eukaryota then you should continue to use RNAmmer.
|
|
|
35
|
|
|
36 The name Barrnap is derived from BActerial Ribosomal RNA Predictor. It was spawned at CodeFest 2013 in Berlin, Germany by Torsten Seemann and Tim Booth.
|
|
|
37
|
|
|
38 .. _barrnap: http://www.vicbioinformatics.com/software.barrnap.shtml
|
|
|
39
|
|
|
40 ------
|
|
|
41
|
|
|
42 **Inputs/Parameters**
|
|
|
43 ::
|
|
|
44 **Fasta file**
|
|
|
45 The fasta file whose contents you want to search for ribosomal RNA's.
|
|
|
46 ::
|
|
|
47 **Length cutoff**
|
|
|
48 ::
|
|
|
49 Proportional length threshold to tag a possible rRNA as pseudo.
|
|
|
50 ::
|
|
|
51 **E-value**
|
|
|
52 Similarity cutoff e-value.
|
|
|
53
|
|
|
54 </help>
|
|
|
55 </tool>
|
