# HG changeset patch # User bgruening # Date 1442524207 14400 # Node ID 42d20ff5da1b455797a739627abdb37d46de5342 planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna/dotknot commit 1973f3035c10db80883d80847ea254289f5cce2a-dirty diff -r 000000000000 -r 42d20ff5da1b dotknot.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/dotknot.xml Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,140 @@ + + pseudoknot prediction in a given RNA sequence + + vienna_rna + dotknot + + + + + + + + dotknot version 1.3.1 + ./output.txt + #if $g: + && + mv ./*.ct ./globalstructure.ct + #end if + ]]> + + + + + + + + + + + + + + + g == True + + + + + + + + + + + + + + + + + + + + followed by the sequence. It can also +contain several consecutive sequences and DotKnot will be executed for each sequence in the file. + +Please ensure that only bases A,C,G,U,T (a,c,g,u,t) are used. + +*Example*:: + + > Arc-Ful-SRP short + GGGGGGUUCGGCGUCCCCUGUAACCCGAAACCGCCGAUACGCGGG + > MMTV + AAAAAACUUGUAAAGGGGCAGUCCCCUAGCCCCGCUCAAAAGGGGGAUG + +----- + +**OPTIONAL ARGUMENTS** + +*Kissing hairpins*:: + Kissing hairpins are complex and biologically relevant types of pseudoknots. + Inclusion of kissing hairpins will lead to increased run time, yet produce more meaningful + results. + +*local pseudoknots*:: + Shows top five near-optimal local pseudoknots in terms of two criteria: estimated free energy + to length ratio and lowest estimated free energy. This can help to identify promising pseudoknot + foldings and may compensate for the limitations of the energy parameters. +*global structure*:: + Shows predicted global structure in addition to predicted pseudoknots + +----- + +**OUTPUT** + +Each pseudoknot is displayed in dot-bracket notation. Unpaired bases are indicated by dots. Base +pairs are written as bracket pairs. The first stem of a pseudoknot is indicated by round brackets , i.e. ( and ). The second stem of a pseudoknot is indicated by square brackets, i.e. [ and ]. + +We also give the start and end positions of the pseudoknot with respect to the input sequence. If the global structure option is chosen (-g:), a CT file will be created which contains the global +structure. The file name will start with the identifier given in the FASTA file. + +----- + +**dotknot** is a Free and Open Source Software, see more details and background information on the dotknot dotknot_ Website. + +.. _dotknot: http://dotknot.csse.uwa.edu.au/ + + ]]> + + 10.1093/nar/gkq021 + 10.1261/rna.2394511 + + diff -r 000000000000 -r 42d20ff5da1b test-data/HCoV229E_short_KISSINGHAIRPINS.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/HCoV229E_short_KISSINGHAIRPINS.fasta Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,2 @@ +> PKB171 human coronavirus 229E, coronavirus; ribosomal frameshifting; polymerase +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG \ No newline at end of file diff -r 000000000000 -r 42d20ff5da1b test-data/TMV.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/TMV.fasta Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,2 @@ +>TMV +UGCAACUUGAGGUAGUCAAGAUGCAUAAUAAAUAACGGAUUGUGUCCGUAAUCACACGUGGUGCGUACGAUAACGCAUAGUGUUUUUCCCUCCACUUAAAUCGAAGGGUUGUGUCUUGGAUCGCGCGGGUCAAAUGUAUAUGGUUCAUAUACAUCCGCAGGCACGUAAUAAAGCGAGGGGUUCGAAUCCCCCCGUUACCCCCGGUAGGGGCCCA diff -r 000000000000 -r 42d20ff5da1b test-data/test1_output.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test1_output.txt Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,30 @@ + +>TMV +UGCAACUUGAGGUAGUCAAGAUGCAUAAUAAAUAACGGAUUGUGUCCGUAAUCACACGUGGUGCGUACGAUAACGCAUAGUGUUUUUCCCUCCACUUAAAUCGAAGGGUUGUGUCUUGGAUCGCGCGGGUCAAAUGUAUAUGGUUCAUAUACAUCCGCAGGCACGUAAUAAAGCGAGGGGUUCGAAUCCCCCCGUUACCCCCGGUAGGGGCCCA +Sequence length: 214 +Using TMV as file identifier + +DotKnot is running... +Predicting pseudoknots... + +Detected pseudoknots: +35 56 -10.97 +ACGGAUUGUGUCCGUAAUCACA +(((((.[[[[)))))...]]]] + +57 78 -12.92 +CGUGGUGCGUACGAUAACGCAU +(((.[[[[[[)))...]]]]]] + +79 108 -8.16 +AGUGUUUUUCCCUCCACUUAAAUCGAAGGG +((((.....[[[[.))))........]]]] + +111 176 -26.94 +GUGUCUUGGAUCGCGCGGGUCAAAUGUAUAUGGUUCAUAUACAUCCGCAGGCACGUAAUAAAGCGA +((((((....[[[[((((((.....(((((((...)))))))))))))))))))........]]]] + +192 210 -13.83 +CCGUUACCCCCGGUAGGGG +(((...[[[[)))..]]]] + diff -r 000000000000 -r 42d20ff5da1b test-data/test2_output.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test2_output.txt Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,56 @@ +Include kissing hairpins. +Show best local pseudoknots. +Show predicted global structure. + +> PKB171 human coronavirus 229E, coronavirus; ribosomal frameshifting; polymerase +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +Sequence length: 54 +Using PKB171 as file identifier + +DotKnot is running... +Predicting pseudoknots... +Predicting kissing hairpins... + +Detected pseudoknots and kissing hairpins: +1 54 -41.2 +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +((((((((((((...[[[[[))))))))))))(((((((..]]]]].))))))) + + +Best 5 pseudoknots and kissing hairpins in terms of energy to length ratio: +1 54 -41.2 +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +((((((((((((...[[[[[))))))))))))(((((((..]]]]].))))))) +1 46 -26.23 +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGC +((((((((((((...[[[[[)))))))))))).........]]]]] +11 54 -18.12 +GUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +(((..((((.....))))....[[[[[[[))).....]]]]]]] +16 37 -7.79 +GCUCGACUAGAGCCCUGUAUGG +((((..[[[))))......]]] +11 48 -13.12 +GUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAG +(((..((((.....)))).[[[[......)))..]]]] + +Best 5 pseudoknots and kissing hairpins in terms of free energy: +1 54 -41.2 +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +((((((((((((...[[[[[))))))))))))(((((((..]]]]].))))))) +1 46 -26.23 +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGC +((((((((((((...[[[[[)))))))))))).........]]]]] +11 54 -18.12 +GUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +(((..((((.....))))....[[[[[[[))).....]]]]]]] +3 54 -16.6 +GGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +(((......[[[.((((.....))))))).(((((((.....]]]))))))) +2 54 -16.2 +GGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +(((...........[[[[[.......)))..(((((((..]]]]].))))))) + +Predicted global structure +CGGGGCUCUAGUGCCGCUCGACUAGAGCCCUGUAUGGACCACGAGCAGUCCAUG +((((((((((((...[[[[[))))))))))))(((((((..]]]]].))))))) diff -r 000000000000 -r 42d20ff5da1b test-data/test2_output2.ct --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test2_output2.ct Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,55 @@ +54 +1 C 0 2 32 1 +2 G 1 3 31 2 +3 G 2 4 30 3 +4 G 3 5 29 4 +5 G 4 6 28 5 +6 C 5 7 27 6 +7 U 6 8 26 7 +8 C 7 9 25 8 +9 U 8 10 24 9 +10 A 9 11 23 10 +11 G 10 12 22 11 +12 U 11 13 21 12 +13 G 12 14 0 13 +14 C 13 15 0 14 +15 C 14 16 0 15 +16 G 15 17 46 16 +17 C 16 18 45 17 +18 U 17 19 44 18 +19 C 18 20 43 19 +20 G 19 21 42 20 +21 A 20 22 12 21 +22 C 21 23 11 22 +23 U 22 24 10 23 +24 A 23 25 9 24 +25 G 24 26 8 25 +26 A 25 27 7 26 +27 G 26 28 6 27 +28 C 27 29 5 28 +29 C 28 30 4 29 +30 C 29 31 3 30 +31 U 30 32 2 31 +32 G 31 33 1 32 +33 U 32 34 54 33 +34 A 33 35 53 34 +35 U 34 36 52 35 +36 G 35 37 51 36 +37 G 36 38 50 37 +38 A 37 39 49 38 +39 C 38 40 48 39 +40 C 39 41 0 40 +41 A 40 42 0 41 +42 C 41 43 20 42 +43 G 42 44 19 43 +44 A 43 45 18 44 +45 G 44 46 17 45 +46 C 45 47 16 46 +47 A 46 48 0 47 +48 G 47 49 39 48 +49 U 48 50 38 49 +50 C 49 51 37 50 +51 C 50 52 36 51 +52 A 51 53 35 52 +53 U 52 54 34 53 +54 G 53 55 33 54 diff -r 000000000000 -r 42d20ff5da1b test-data/tmRNA_alignment.fasta --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/tmRNA_alignment.fasta Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,11 @@ +>Ecoli-tmRNA +GGGGCUGAUUCUGGAUUCGACGGGAUUUGCGAAACCCAAGGUGCAUGCCGAGGGGCGGUUGGCCUCGUAAAAAGCCGCAAAAAAUAGUCGCAAACGACGAAAACUACGCUUUAGCAGCUUAAUAACCUGCUUAGAGCCCUCUCUCCCUAGCCUCCGCUCUUAGGACGGGGAUCAAGAGAGGUCAAACCCAAAAGAGAUCGCGUGGAAGCCCUGCCUGGGGUUGAAGCGUUAAAACUUAAUCAGGCUAGUUUGUUAGUGGCGUGUCCGUCCGCAGCUGGCAAGCGAAUGUAAAGACUGACUAAGCAUGUAGUACCGAGGAUGUAGGAAUUUCGGACGCGGGUUCAACUCCCGCCAGCUCCACCA +>Mesoplasma florum tmRNA +GGGGAUGUCCUGGUUUCGACAGGAUACUCUGGUCAUAUAAUAGCAGUGGUGUGGUAGACCAUAAUACUUCUAGGUUUGAUAAAAUGCAAACAAAAACGAAGAAAACACAAACGAAGUUCCAACAUUUAUGUUAAAUGCUGGACAAGCAAAUUACGCUUUCGCUUAAUUUGACUUAAAUACAAAUAUACGUUUGGCAAUUUUUAAUAGUCUUGGAAUCUGUUAAUUAUUGUACAUACCAAUUCCAAUAGAAAUUAUGAUCGUAAUAAAUAAUUUUGAAACUCAAUUACGGAUGUUAUAUUAGUCUUUUGAUAAUUUUAUUUAAUAUAGCAAAUUAGAAAAUUAUCUAAACUGUAGAAUUUAUAUGAUUCGGUAUUUUGGACGUGGGUUCGAUUCCCACCAUCUCCACCA +>Thermus thermophilus tmRNA, version 1 {HB8}: +GGGGGUGAAACGGUCUCGACGGGGGUCGCCGAGGGCGUGGCUGCGCGCCGAGGUGCGGGUGGCCUCGUAAAAACCCGCAACGGCAUAACUGCCAACACCAACUACGCUCUCGCGGCUUAAUGACCGCGACCUCGCCCGGUAGCCCUGCCGGGGGCUCACCGGAAGCGGGGACACAAACCCGGCUAGCCCGGGGCCACGCCCUCUAACCCCGGGCGAAGCUUGAAGGGGGCUCGCUCCUGGCCGCCCGUCCGCGGGCCAAGCCAGGAGGACACGCGAAACGCGGACUACGCGCGUAGAGGCCCGCCGUAGAGACCUUCGGACGGGGGUUCGACUCCCCCCACCUCCACCA +>Actinobacillus actinomycetemcomitans tmRNA +GGGGCUGAUUCUGGAUUCGACGGGAUUGGCGAAGCCCAAGGUGCACGUCGAGGUGCGGUAGGCCUCGUAAAUAAACCGCAAAAAAAUAGUCGCAAACGACGAACAAUACGCUUUAGCAGCUUAAUAACCUGCAUUUAGCCUUCGCUCCCCAGCUUCCGCUCGUAAGACGGGGAUCAAGCGGAGUCAAACCUAAACGAGAUCGUGUGGACGCUUCGGCUUGGAGAUCGAAACACUAAAUUGAAUCAAGCUAGUUUAUUCUUCGCGUGUCUGUCCGCAGGGGAUAAAUGAAGUUAAAGACCAGACUAAACGUGUAGUGCUGAAGGUAGAGUAAUUUCGGACGCGGGUUCAACUCCCGCCAGCUCCACCA +>marine_actinobacteri.1 AAOB01000009.1/82671-83044 +GGGGAUGAUCGGUUUCGACAUUGCCUGCAAAACUGUGAGAAGCGGGUCGAGAAUGCAGCCUUAUCUCGUUAACGAUGACUGCAAACUAUAAGUGCCAAUUCAAAGCGCACUGACUUCGCCCUCGCUGCCUAAGCGAGCGCACUAAAGAAGUCCGUCAGACCGGGAAUGCUCUCUACCCGGAUCCUGGCGCAAUUUAGAGAGAUUGCUGCGUAGUUACGCCUGAGGGCUACGCGGGACUUGCACUCUGGCUGGGCUUGUUGAUCUAGGUGCUUGUGGCAAAGAUCAGAGCCGAGUAGAACGUCUGCACAAGCUACACCCGUAGAAGGCACAGAAUUACAGCAGUGGACGGGGGUUCAAUUCCCCCCAUCUCCACC + diff -r 000000000000 -r 42d20ff5da1b tool_dependencies.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool_dependencies.xml Thu Sep 17 17:10:07 2015 -0400 @@ -0,0 +1,27 @@ + + + + + + + + + https://github.com/bgruening/download_store/raw/master/dotknot/dotknot_1.3.1.tar.bz2 + + . + $INSTALL_DIR + + + $INSTALL_DIR + $INSTALL_DIR + + + + + + DotKnot:a heuristic method for pseudoknot prediction in a given RNA sequence. + See http://dotknot.csse.uwa.edu.au/ + + + +