Mercurial > repos > iuc > barrnap

comparison barrnap.xml @ 0:8fd427007e93 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 6d089b07b7a1d50b619e2ac6f33cd2bfe2a793c1
author iuc
date Mon, 28 Aug 2017 14:14:41 -0400
parents
children 7355bbf5e05f
comparison
equal deleted inserted replaced
-1:000000000000 0:8fd427007e93
1 <tool id="barrnap" name="barrnap" version="1.2.0">
2 <description>Locate ribosomal RNA's in a fasta file. (GFF output)</description>
3 <requirements>
4 <requirement type="package" version="0.8">barrnap</requirement>
5 </requirements>
6 <version_command>barrnap --version</version_command>
7 <command detect_errors="exit_code"><![CDATA[
8 barrnap
9 --quiet
10 --threads \${GALAXY_SLOTS:-1}
11 --reject $reject
12 --lencutoff $lencutoff
13 --evalue $evalue
14 #if $incseq
15 --incseq
16 #end if
17 --kingdom $kingdom
18 '$fasta_file'
19 > '$gff'
20 ]]>
21 </command>
22
23 <inputs>
24 <param name="fasta_file" type="data" format="fasta" label="Fasta file" help="Fasta file to look for rRNA's in" />
25 <param name="kingdom" type="select" label="Kingdom">
26 <option value="bac" selected="true">Bacteria</option>
27 <option value="euk">Eukaryote</option>
28 <option value="mito">Mitochondria</option>
29 <option value="arc">Archaea</option>
30 </param>
31 <param name="reject" type="float" value="0.5" label="Reject length threshold" help="Proportional length threshold to reject prediction" />
32 <param name="lencutoff" type="float" value="0.8" label="Length cutoff" help="Proportional length threshold to tag a find as pseudo" />
33 <param name="evalue" type="float" value="1e-06" label="e-value" help="Similarity e-value cutoff" />
34 <param name="incseq" type="boolean" label="Include Sequences in GFF" help="Include the fasta sequences in the GFF output file" />
35 </inputs>
36
37 <outputs>
38 <data name="gff" format="gff3" label="${tool.name} on ${on_string}: rRNA GFF"/>
39 </outputs>
40
41 <tests>
42 <test><!-- Test 1 -->
43 <param name="fasta_file" value="bacteria.fna" ftype="fasta" />
44 <param name="kingdom" value="bac" />
45 <output name="gff" value="test_1.gff" ftype="gff3" />
46 </test>
47 <test> <!-- Test 2 -->
48 <param name="fasta_file" value="small.fna" ftype="fasta" />
49 <param name="kingdom" value="bac" />
50 <param name="incseq" value="true" />
51 <output name="gff" value="test_2.gff" ftype="gff3" />
52 </test>
53 <test> <!-- Test 3 -->
54 <param name="fasta_file" value="small.fna" ftype="fasta" />
55 <param name="kingdom" value="bac" />
56 <param name="evalue" value="0.00000001" />
57 <output name="gff" value="test_3.gff" ftype="gff3" />
58 </test>
59 </tests>
60
61 <help>
62
63 **barrnap**
64
65 barrnap_ predicts the location of 5S, 16S and 23S ribosomal RNA genes in Bacterial genome sequences. Barrnap now supports Archaea, Eukaryota and Mitochondria. 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.
66 This tool is designed to be a substitute for RNAmmer. 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.
67
68 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.
69
70 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.
71
72 .. _barrnap: https://github.com/tseemann/barrnap
73
74 **Inputs/Parameters**
75
76 *Fasta file*: The fasta file whose contents you want to search for ribosomal RNA's.
77
78 *Kingdom*: Select the Kingdom of the organism: Bacteria, Eukarote, Mitochondria or Archaea.
79
80 *Length cutoff*: Proportional length threshold to tag a possible rRNA as pseudo.
81
82 *Reject length threshold*: Proportional length threshold to reject prediction
83
84 *E-value*: Similarity cutoff e-value.
85
86 *Include Sequences in GFF*: This includes the original fasta sequence in the GFF file below the #FASTA tag
87
88 </help>
89 <citations>
90 <citation type="bibtex">@UNPUBLISHED{Seemann2013,
91 author = "Torsten Seemann",
92 title = "barrnap 0.8 : rapid ribosomal RNA prediction",
93 year = "2013",
94 note = "https://github.com/tseemann/barrnap"}
95 </citation>
96 </citations>
97 </tool>