Mercurial > repos > iuc > meme_meme
changeset 19:7ce0fe11b411 draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/meme commit 93fe7345e2b6a97000975d978bf06513b3a442d9
| author | iuc |
|---|---|
| date | Thu, 29 Aug 2024 10:20:10 +0000 |
| parents | 19237efe98d1 |
| children | |
| files | macros.xml meme.xml test-data/dreme_output_test1.html test-data/dreme_output_test1.txt test-data/dreme_output_test1.xml test-data/dreme_output_test2.html test-data/dreme_output_test2.txt test-data/dreme_output_test2.xml test-data/fimo_output_test1.gff test-data/fimo_output_test1.html test-data/fimo_output_test1.txt test-data/fimo_output_test1.xml test-data/fimo_output_test2.gff test-data/fimo_output_test2.html test-data/fimo_output_test2.txt test-data/fimo_output_test2.xml test-data/fimo_output_test3.html test-data/fimo_output_test3.txt test-data/fimo_output_test3.xml |
| diffstat | 19 files changed, 68 insertions(+), 686 deletions(-) [+] |
line wrap: on
line diff
--- a/macros.xml Sat Apr 09 08:34:02 2022 +0000 +++ b/macros.xml Thu Aug 29 10:20:10 2024 +0000 @@ -1,7 +1,7 @@ -<?xml version='1.0' encoding='UTF-8'?> <macros> <token name="@VERSION_SUFFIX@">0</token> - <token name="@TOOL_VERSION@">5.4.1</token> + <token name="@TOOL_VERSION@">5.5.6</token> + <token name="@PROFILE@">23.0</token> <xml name="requirements"> <requirements> <requirement type="package" version="@TOOL_VERSION@">meme</requirement>
--- a/meme.xml Sat Apr 09 08:34:02 2022 +0000 +++ b/meme.xml Thu Aug 29 10:20:10 2024 +0000 @@ -1,12 +1,12 @@ -<tool id="meme_meme" name="MEME" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@"> +<tool id="meme_meme" name="MEME" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="@PROFILE@"> <description>- Multiple EM for Motif Elicitation</description> + <macros> + <import>macros.xml</import> + </macros> <xrefs> <xref type="bio.tools">meme_meme</xref> </xrefs> - <macros> - <import>macros.xml</import> - </macros> - <expand macro="requirements" /> + <expand macro="requirements"/> <command detect_errors="exit_code"><![CDATA[ @CHECK_NON_COMMERCIAL_USE@ export TMPDIR=\${TMPDIR:-.}; @@ -81,15 +81,15 @@ && mv '${html_outfile.files_path}/meme.xml' '${xml_outfile}' ]]></command> <inputs> - <param format="fasta" name="input1" type="data" label="Sequences"/> + <param name="input1" type="data" format="fasta" label="Sequences"/> <conditional name="options_type"> <param name="options_type_selector" type="select" label="Options Configuration"> <option value="basic" selected="true">Basic</option> <option value="advanced">Advanced</option> </param> - <when value="basic" /> + <when value="basic"/> <when value="advanced"> - <param name="sf" type="text" value="Galaxy FASTA Input" label="Name of sequence set" argument="-sf"/> + <param argument="-sf" type="text" value="Galaxy FASTA Input" label="Name of sequence set"/> <conditional name="alphabet_type"> <param name="alphabet_type_selector" type="select" label="Sequence Alphabet"> <option value="protein">Protein</option> @@ -97,7 +97,7 @@ </param> <when value="protein"> <conditional name="prior_type"> - <param name="prior_type_selector" type="select" label="Choice of prior" argument="-prior"> + <param name="prior_type_selector" argument="-prior" type="select" label="Choice of prior"> <option value="dirichlet">simple Dirichlet prior</option> <option value="dmix" selected="true">mixture of Dirichlets prior</option> <option value="mega">extremely low variance dmix</option> @@ -105,45 +105,45 @@ <option value="addone">add +1 to each observed count</option> </param> <when value="dirichlet"> - <param name="prior_b" type="float" value="0.01" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0.01" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="dmix"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="mega"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="megap"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> - <when value="addone" /> + <when value="addone"/> </conditional> <conditional name="spmap_type"> - <param name="spmap_type_selector" type="select" label="EM starting points" argument="-spmap"> + <param name="spmap_type_selector" argument="-spmap" type="select" label="EM starting points"> <option value="uni">uni</option> <option value="pam" selected="true">pam</option> <option value="cons">Use starting point from string</option> </param> <when value="uni"> - <param name="spfuzz" type="float" value="0.5" label="Fuzziness of the mapping" argument="-spfuzz"/> + <param argument="-spfuzz" type="float" value="0.5" label="Fuzziness of the mapping"/> </when> <when value="pam"> - <param name="spfuzz" type="integer" value="120" label="Fuzziness of the mapping" argument="-spfuzz"/> + <param argument="-spfuzz" type="integer" value="120" label="Fuzziness of the mapping"/> </when> <when value="cons"> - <param name="cons" type="text" value="" label="Starting point from string" argument="-cons"/> + <param argument="-cons" type="text" value="" label="Starting point from string"/> </when> </conditional> </when> <when value="dna"> - <param name="revcomp" label="Check reverse complement" type="boolean" truevalue="-revcomp" falsevalue="" checked="False"/> - <param name="pal" label="Check for palindromes" type="boolean" truevalue="-pal" falsevalue="" checked="False"/> + <param argument="-revcomp" type="boolean" truevalue="-revcomp" falsevalue="" checked="False" label="Check reverse complement"/> + <param argument="-pal" type="boolean" truevalue="-pal" falsevalue="" checked="False" label="Check for palindromes"/> <conditional name="prior_type"> - <param name="prior_type_selector" type="select" label="Sequence Alphabet" argument="-prior"> + <param name="prior_type_selector" argument="-prior" type="select" label="Sequence Alphabet"> <option value="dirichlet" selected="true">simple Dirichlet prior</option> <option value="dmix">mixture of Dirichlets prior</option> <option value="mega">extremely low variance dmix</option> @@ -151,50 +151,50 @@ <option value="addone">add +1 to each observed count</option> </param> <when value="dirichlet"> - <param name="prior_b" type="float" value="0.01" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0.01" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="dmix"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="mega"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> <when value="megap"> - <param name="prior_b" type="float" value="0" label="strength of prior on model parameters" argument="-b"/> - <param name="plib" type="data" format="txt" optional="True" label="Dirichlet prior file" argument="-plib"/> + <param name="prior_b" argument="-b" type="float" value="0" label="strength of prior on model parameters"/> + <param argument="-plib" type="data" format="txt" optional="True" label="Dirichlet prior file"/> </when> - <when value="addone" /> + <when value="addone"/> </conditional> <conditional name="spmap_type"> - <param name="spmap_type_selector" type="select" label="EM starting points" argument="-spmap"> + <param name="spmap_type_selector" argument="-spmap" type="select" label="EM starting points"> <option value="uni" selected="true">uni</option> <option value="pam">pam</option> <option value="cons">Use starting point from string</option> </param> <when value="uni"> - <param name="spfuzz" type="float" value="0.5" label="Fuzziness of the mapping" argument="-spfuzz"/> + <param argument="-spfuzz" type="float" value="0.5" label="Fuzziness of the mapping"/> </when> <when value="pam"> - <param name="spfuzz" type="integer" value="120" label="Fuzziness of the mapping" argument="-spfuzz"/> + <param argument="-spfuzz" type="integer" value="120" label="Fuzziness of the mapping"/> </when> <when value="cons"> - <param name="cons" type="text" value="" label="Starting point from string" argument="-cons"/> + <param argument="-cons" type="text" value="" label="Starting point from string"/> </when> </conditional> </when> </conditional> - <param name="nmotifs" type="integer" value="1" label="Number of different motifs to search" argument="-nmotifs" /> - <param name="evt" type="float" value="inf" label="E-value to stop looking for motifs" argument="-evt"/> + <param argument="-nmotifs" type="integer" value="1" label="Number of different motifs to search"/> + <param argument="-evt" type="float" value="inf" label="E-value to stop looking for motifs"/> <conditional name="mod_type"> - <param name="mod_type_selector" type="select" label="Expected motif distribution" argument="-mod"> + <param name="mod_type_selector" argument="-mod" type="select" label="Expected motif distribution"> <option value="oops">One Occurrence Per Sequence</option> <option value="zoops" selected="true">Zero or One Occurrence Per Sequence</option> <option value="anr">Any Number of Repetitions</option> </param> - <when value="oops" /> + <when value="oops"/> <when value="zoops"> <conditional name="motif_occurrence_type"> <param name="motif_occurrence_type_selector" type="select" label="Number of motif occurrences"> @@ -202,13 +202,13 @@ <option value="nsites">nsites</option> <option value="min_max_sites">min and max sites</option> </param> - <when value="default" /> + <when value="default"/> <when value="nsites"> - <param name="nsites" type="integer" value="1" label="Search nsites number of occurrences" argument="-nsites"/> + <param argument="-nsites" type="integer" value="1" label="Search nsites number of occurrences"/> </when> <when value="min_max_sites"> - <param name="minsites" type="integer" value="1" label="minsites" argument="-minsites"/> - <param name="maxsites" type="integer" value="50" label="maxsites" argument="-maxsites"/> + <param argument="-minsites" type="integer" value="1" label="minsites"/> + <param argument="-maxsites" type="integer" value="50" label="maxsites"/> </when> </conditional> </when> @@ -219,29 +219,29 @@ <option value="nsites">nsites</option> <option value="min_max_sites">min and max sites</option> </param> - <when value="default" /> + <when value="default"/> <when value="nsites"> - <param name="nsites" type="integer" value="1" label="Search nsites number of occurrences" argument="-nsites"/> + <param argument="-nsites" type="integer" value="1" label="Search nsites number of occurrences"/> </when> <when value="min_max_sites"> - <param name="minsites" type="integer" value="1" label="minsites" argument="-minsites"/> - <param name="maxsites" type="integer" value="50" label="maxsites" argument="-maxsites"/> + <param argument="-minsites" type="integer" value="1" label="minsites"/> + <param argument="-maxsites" type="integer" value="50" label="maxsites"/> </when> </conditional> </when> </conditional> - <param name="wnsites" type="float" value="0.8" label="Weight on the prior on nsites" argument="-wnsites"/> + <param argument="-wnsites" type="float" value="0.8" label="Weight on the prior on nsites"/> <conditional name="motif_width_type"> <param name="motif_width_type_selector" type="select" label="Motif width type"> <option value="exact">Exact width</option> <option value="range" selected="true">Specify a range</option> </param> <when value="exact"> - <param name="width" type="integer" value="10" label="Width of motif to search" argument="-w"/> + <param name="width" argument="-w" type="integer" value="10" label="Width of motif to search"/> </when> <when value="range"> - <param name="minw" type="integer" value="8" label="Min width of motif to search" argument="-minw"/> - <param name="maxw" type="integer" value="50" label="Max width of motif to search" argument="-maxw"/> + <param argument="-minw" type="integer" value="8" label="Min width of motif to search"/> + <param argument="-maxw" type="integer" value="50" label="Max width of motif to search"/> </when> </conditional> <conditional name="motif_trim_type"> @@ -249,31 +249,31 @@ <option value="nomatrim">No motif trim</option> <option value="trim" selected="true">Trim motif</option> </param> - <when value="nomatrim" /> + <when value="nomatrim"/> <when value="trim"> - <param name="wg" type="integer" value="11" label="Gap cost" argument="-wg"/> - <param name="ws" type="integer" value="1" label="Space cost" argument="-ws"/> - <param name="noendgaps" label="Do not penalize endgaps" type="boolean" truevalue="-noendgaps" falsevalue="" checked="False"/> + <param argument="-wg" type="integer" value="11" label="Gap cost"/> + <param argument="-ws" type="integer" value="1" label="Space cost"/> + <param argument="-noendgaps" type="boolean" truevalue="-noendgaps" falsevalue="" checked="False" label="Do not penalize endgaps"/> </when> </conditional> - <param name="bfile" type="data" format="txt" optional="True" label="Background Model" argument="-bfile"/> - <param name="pspfile" type="data" format="txt" optional="True" label="Position-Specific Prior" argument="-psp"/> - <param name="maxiter" type="integer" value="50" label="Number of iterations of EM to run" argument="-maxiter"/> - <param name="distance" type="float" value="0.001" label="Convergence criterion" argument="-distance"/> + <param argument="-bfile" type="data" format="txt" optional="True" label="Background Model"/> + <param name="pspfile" argument="-psp" type="data" format="txt" optional="True" label="Position-Specific Prior"/> + <param argument="-maxiter" type="integer" value="50" label="Number of iterations of EM to run"/> + <param argument="-distance" type="float" value="0.001" label="Convergence criterion"/> <conditional name="branching_type"> <param name="branching_type_selector" type="select" label="x-branching type"> <option value="x_branch">Perform x-branching</option> <option value="no_x_branch" selected="true">No x-branching</option> </param> - <when value="no_x_branch" /> + <when value="no_x_branch"/> <when value="x_branch"> - <param name="bfactor" type="integer" value="3" label="Number of iterations of branching" argument="-bfactor"/> - <param name="heapsize" type="integer" value="64" label="Maximum number of heaps to use" argument="-heapsize"/> + <param argument="-bfactor" type="integer" value="3" label="Number of iterations of branching"/> + <param argument="-heapsize" type="integer" value="64" label="Maximum number of heaps to use"/> </when> </conditional> </when> </conditional> - <param name="non_commercial_use" label="I certify that I am not using this tool for commercial purposes." type="boolean" truevalue="NON_COMMERCIAL_USE" falsevalue="COMMERCIAL_USE" checked="False"> + <param name="non_commercial_use" type="boolean" truevalue="NON_COMMERCIAL_USE" falsevalue="COMMERCIAL_USE" checked="False" label="I certify that I am not using this tool for commercial purposes."> <validator type="expression" message="This tool is only available for non-commercial use.">value == True</validator> </param> </inputs>
--- a/test-data/dreme_output_test1.html Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,100 +0,0 @@ - </div> - <div class="box"> - <h4>Sequences</h4> - <table id="seq_info" class="inputs"> - <tr><th>Source <div class="help" data-topic="pop_seq_source"></div></th> - <th>Alphabet <div class="help" data-topic="pop_seq_alph"></div></th> - <th>Sequence Count <div class="help" data-topic="pop_seq_count"></div></th> - </tr> - <tr> - <td id="ins_seq_source"></td> - <td id="ins_seq_alphabet"></td> - <td id="ins_seq_count"></td> - </tr> - </table> - <script> - { - var db = data.sequence_db; - $("ins_seq_source").innerHTML = db.file; - $("ins_seq_alphabet").innerHTML = current_alphabet.get_alphabet_name(); - $("ins_seq_count").innerHTML = db.count; - } - </script> - <h4>Control Sequences</h4> - <table id="seq_info" class="inputs"> - <tr><th>Source <div class="help" data-topic="pop_seq_source"></div></th> - <th>Sequence Count <div class="help" data-topic="pop_seq_count"></div></th> - </tr> - <tr> - <td id="ins_cseq_source"></td> - <td id="ins_cseq_count"></td> - </tr> - </table> - <script> - { - var db = data.control_db; - if (db.from == "shuffled") { - $("ins_cseq_source").innerHTML = "Shuffled Sequences"; - } else { - $("ins_cseq_source").innerHTML = db.file; - } - $("ins_cseq_count").innerHTML = db.count; - } - </script> - <h4>Background</h4> - <span id="alpha_bg"></span> - <script> - { - $("alpha_bg").appendChild(make_alpha_bg_table(current_alphabet, data.control_db.freqs)); - } - </script> - <h4>Other Settings</h4> - <table id="tbl_settings" class="inputs hide_advanced"> - <tr> - <th>Strand Handling</th> - <td id="opt_strand"> - <span class="strand_none">This alphabet only has one strand</span> - <span class="strand_given">Only the given strand is processed</span> - <span class="strand_both">Both the given and reverse complement strands are processed</span> - </td> - </tr> - <tr><th># REs to Generalize</th><td id="opt_ngen"></td></tr> - <tr><th>Shuffle Seed</th><td id="opt_seed"></td></tr> - <tr><th>E-value Threshold</th><td id="opt_stop_evalue"></td></tr> - <tr><th>Max Motif Count</th><td id="opt_stop_count"></td></tr> - <tr><th>Max Run Time</th><td id="opt_stop_time"></td></tr> - </table> - <script> - { - $("opt_strand").className = (current_alphabet.has_complement() ? (data.options.revcomp ? "both" : "given") : "none"); - $("opt_ngen").innerHTML = data.options.ngen; - $("opt_seed").innerHTML = data.options.seed; - $("opt_stop_evalue").innerHTML = data.options.stop.evalue; - $("opt_stop_count").innerHTML = (typeof data.options.stop.count == "number" ? data.options.stop.count : "No maximum motif count."); - $("opt_stop_time").innerHTML = (typeof data.options.stop.time == "number" ? data.options.stop.time + " seconds." : "No maximum running time."); - } - </script> - </div> - <!-- list information on this program --> - <div id="info_sec" class="bar" style="position:relative"> - <div style="position: absolute; right: 0;"><a href="#inputs_sec">Previous</a> <a href="#">Top</a></div> - <div class="subsection"> - <h5 id="version">DREME version</h5> - <span id="ins_version"></span> - (Release date: <span id="ins_release"></span>)<br> - </div> - <script> - $("ins_version").innerHTML = data["version"]; - $("ins_release").innerHTML = data["release"]; - </script> - <div class="subsection" id="reference"> <script>print_citation("reference", "DREME");</script> </div> - <div class="subsection"> - <h5 id="command">Command line</h5> - <textarea id="cmd" rows="3" style="width:100%;" readonly="readonly"> - </textarea> - <script>$("cmd").value = data["cmd"].join(" ");</script> - </div> - </div> - - </body> -</html>
--- a/test-data/dreme_output_test1.txt Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,64 +0,0 @@ -MOTIF UUYUCY DREME-1 - -# Word Pos Neg P-value E-value -# BEST UUYUCY 387 210 2.6e-018 1.2e-013 -# UUUUCC 147 75 1.8e-007 8.1e-003 -# UUUUCU 155 94 2.2e-005 1.0e+000 -# UUCUCU 94 51 1.3e-004 6.1e+000 -# UUCUCC 75 42 1.1e-003 5.0e+001 - -letter-probability matrix: alength= 4 w= 6 nsites= 459 E= 1.2e-013 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.294118 0.000000 0.705882 -0.000000 0.000000 0.000000 1.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.474946 0.000000 0.525054 - - -MOTIF YAGG DREME-2 - -# Word Pos Neg P-value E-value -# BEST YAGG 600 416 1.1e-016 5.1e-012 -# CAGG 441 304 1.5e-010 6.6e-006 -# UAGG 232 165 1.1e-004 4.7e+000 - -letter-probability matrix: alength= 4 w= 4 nsites= 793 E= 5.1e-012 -0.000000 0.692308 0.000000 0.307692 -1.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 - - -MOTIF GAAGAW DREME-3 - -# Word Pos Neg P-value E-value -# BEST GAAGAW 81 22 8.2e-010 3.4e-005 -# GAAGAU 45 7 2.4e-008 9.9e-004 -# GAAGAA 40 16 7.9e-004 3.3e+001 - -letter-probability matrix: alength= 4 w= 6 nsites= 89 E= 3.4e-005 -0.000000 0.000000 1.000000 0.000000 -1.000000 0.000000 0.000000 0.000000 -1.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -1.000000 0.000000 0.000000 0.000000 -0.494382 0.000000 0.000000 0.505618 - - -MOTIF SMUGGA DREME-4 - -# Word Pos Neg P-value E-value -# BEST SMUGGA 110 47 9.1e-008 3.7e-003 -# GAUGGA 22 6 1.7e-003 7.1e+001 -# GCUGGA 33 14 3.6e-003 1.5e+002 -# CCUGGA 32 15 8.6e-003 3.5e+002 -# CAUGGA 29 13 9.1e-003 3.7e+002 - -letter-probability matrix: alength= 4 w= 6 nsites= 119 E= 3.7e-003 -0.000000 0.529412 0.470588 0.000000 -0.428571 0.571429 0.000000 0.000000 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -1.000000 0.000000 0.000000 0.000000
--- a/test-data/dreme_output_test1.xml Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,44 +0,0 @@ - <motifs> - <motif id="m01" alt="DREME-1" seq="UUYUCY" length="6" nsites="459" p="387" n="210" pvalue="2.6e-018" evalue="1.2e-013" unerased_evalue="1.2e-013"> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="0.294118" G="0.000000" U="0.705882"/> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="1.000000" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.474946" G="0.000000" U="0.525054"/> - <match seq="UUUUCC" p="147" n="75" pvalue="1.8e-007" evalue="8.1e-003"/> - <match seq="UUUUCU" p="155" n="94" pvalue="2.2e-005" evalue="1.0e+000"/> - <match seq="UUCUCU" p="94" n="51" pvalue="1.3e-004" evalue="6.1e+000"/> - <match seq="UUCUCC" p="75" n="42" pvalue="1.1e-003" evalue="5.0e+001"/> - </motif> - <motif id="m02" alt="DREME-2" seq="YAGG" length="4" nsites="793" p="600" n="416" pvalue="1.1e-016" evalue="5.1e-012" unerased_evalue="2.4e-012"> - <pos A="0.000000" C="0.692308" G="0.000000" U="0.307692"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <match seq="CAGG" p="441" n="304" pvalue="1.5e-010" evalue="6.6e-006"/> - <match seq="UAGG" p="232" n="165" pvalue="1.1e-004" evalue="4.7e+000"/> - </motif> - <motif id="m03" alt="DREME-3" seq="GAAGAW" length="6" nsites="89" p="81" n="22" pvalue="8.2e-010" evalue="3.4e-005" unerased_evalue="3.5e-004"> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <pos A="0.494382" C="0.000000" G="0.000000" U="0.505618"/> - <match seq="GAAGAU" p="45" n="7" pvalue="2.4e-008" evalue="9.9e-004"/> - <match seq="GAAGAA" p="40" n="16" pvalue="7.9e-004" evalue="3.3e+001"/> - </motif> - <motif id="m04" alt="DREME-4" seq="SMUGGA" length="6" nsites="119" p="110" n="47" pvalue="9.1e-008" evalue="3.7e-003" unerased_evalue="2.6e-005"> - <pos A="0.000000" C="0.529412" G="0.470588" U="0.000000"/> - <pos A="0.428571" C="0.571429" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <match seq="GAUGGA" p="22" n="6" pvalue="1.7e-003" evalue="7.1e+001"/> - <match seq="GCUGGA" p="33" n="14" pvalue="3.6e-003" evalue="1.5e+002"/> - <match seq="CCUGGA" p="32" n="15" pvalue="8.6e-003" evalue="3.5e+002"/> - <match seq="CAUGGA" p="29" n="13" pvalue="9.1e-003" evalue="3.7e+002"/> - </motif> - </motifs> \ No newline at end of file
--- a/test-data/dreme_output_test2.html Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,50 +0,0 @@ - <h4>Other Settings</h4> - <table id="tbl_settings" class="inputs hide_advanced"> - <tr> - <th>Strand Handling</th> - <td id="opt_strand"> - <span class="strand_none">This alphabet only has one strand</span> - <span class="strand_given">Only the given strand is processed</span> - <span class="strand_both">Both the given and reverse complement strands are processed</span> - </td> - </tr> - <tr><th># REs to Generalize</th><td id="opt_ngen"></td></tr> - <tr><th>Shuffle Seed</th><td id="opt_seed"></td></tr> - <tr><th>E-value Threshold</th><td id="opt_stop_evalue"></td></tr> - <tr><th>Max Motif Count</th><td id="opt_stop_count"></td></tr> - <tr><th>Max Run Time</th><td id="opt_stop_time"></td></tr> - </table> - <script> - { - $("opt_strand").className = (current_alphabet.has_complement() ? (data.options.revcomp ? "both" : "given") : "none"); - $("opt_ngen").innerHTML = data.options.ngen; - $("opt_seed").innerHTML = data.options.seed; - $("opt_stop_evalue").innerHTML = data.options.stop.evalue; - $("opt_stop_count").innerHTML = (typeof data.options.stop.count == "number" ? data.options.stop.count : "No maximum motif count."); - $("opt_stop_time").innerHTML = (typeof data.options.stop.time == "number" ? data.options.stop.time + " seconds." : "No maximum running time."); - } - </script> - </div> - <!-- list information on this program --> - <div id="info_sec" class="bar" style="position:relative"> - <div style="position: absolute; right: 0;"><a href="#inputs_sec">Previous</a> <a href="#">Top</a></div> - <div class="subsection"> - <h5 id="version">DREME version</h5> - <span id="ins_version"></span> - (Release date: <span id="ins_release"></span>)<br> - </div> - <script> - $("ins_version").innerHTML = data["version"]; - $("ins_release").innerHTML = data["release"]; - </script> - <div class="subsection" id="reference"> <script>print_citation("reference", "DREME");</script> </div> - <div class="subsection"> - <h5 id="command">Command line</h5> - <textarea id="cmd" rows="3" style="width:100%;" readonly="readonly"> - </textarea> - <script>$("cmd").value = data["cmd"].join(" ");</script> - </div> - </div> - - </body> -</html>
--- a/test-data/dreme_output_test2.txt Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,55 +0,0 @@ -ALPHABET "RNA" RNA-LIKE -A "Adenine" CC0000 -C "Cytosine" 0000CC -G "Guanine" FFB300 -U "Uracil" 008000 -N "Any base" = ACGU -X = ACGU -. = ACGU -V "Not U" = ACG -H "Not G" = ACU -D "Not C" = AGU -B "Not A" = CGU -M "Amino" = AC -R "Purine" = AG -W "Weak" = AU -S "Strong" = CG -Y "Pyrimidine" = CU -K "Keto" = GU -T = U -END ALPHABET - -Background letter frequencies (from dataset): -A 0.221 C 0.245 G 0.221 U 0.312 - - -MOTIF UUYUCY DREME-1 - -# Word Pos Neg P-value E-value -# BEST UUYUCY 387 210 2.6e-018 3.3e-013 -# UUUUCC 147 75 1.8e-007 2.2e-002 -# UUUUCU 155 94 2.2e-005 2.8e+000 -# UUCUCU 94 51 1.3e-004 1.7e+001 -# UUCUCC 75 42 1.1e-003 1.4e+002 - -letter-probability matrix: alength= 4 w= 6 nsites= 459 E= 3.3e-013 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.000000 0.000000 1.000000 -0.000000 0.294118 0.000000 0.705882 -0.000000 0.000000 0.000000 1.000000 -0.000000 1.000000 0.000000 0.000000 -0.000000 0.474946 0.000000 0.525054 - - -MOTIF YAGG DREME-2 - -# Word Pos Neg P-value E-value -# BEST YAGG 600 416 1.1e-016 1.4e-011 -# CAGG 441 304 1.5e-010 1.8e-005 -# UAGG 232 165 1.1e-004 1.3e+001 - -letter-probability matrix: alength= 4 w= 4 nsites= 793 E= 1.4e-011 -0.000000 0.692308 0.000000 0.307692 -1.000000 0.000000 0.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 -0.000000 0.000000 1.000000 0.000000 \ No newline at end of file
--- a/test-data/dreme_output_test2.xml Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,22 +0,0 @@ - <motifs> - <motif id="m01" alt="DREME-1" seq="UUYUCY" length="6" nsites="459" p="387" n="210" pvalue="2.6e-018" evalue="3.3e-013" unerased_evalue="3.3e-013"> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="0.294118" G="0.000000" U="0.705882"/> - <pos A="0.000000" C="0.000000" G="0.000000" U="1.000000"/> - <pos A="0.000000" C="1.000000" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.474946" G="0.000000" U="0.525054"/> - <match seq="UUUUCC" p="147" n="75" pvalue="1.8e-007" evalue="2.2e-002"/> - <match seq="UUUUCU" p="155" n="94" pvalue="2.2e-005" evalue="2.8e+000"/> - <match seq="UUCUCU" p="94" n="51" pvalue="1.3e-004" evalue="1.7e+001"/> - <match seq="UUCUCC" p="75" n="42" pvalue="1.1e-003" evalue="1.4e+002"/> - </motif> - <motif id="m02" alt="DREME-2" seq="YAGG" length="4" nsites="793" p="600" n="416" pvalue="1.1e-016" evalue="1.4e-011" unerased_evalue="6.3e-012"> - <pos A="0.000000" C="0.692308" G="0.000000" U="0.307692"/> - <pos A="1.000000" C="0.000000" G="0.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <pos A="0.000000" C="0.000000" G="1.000000" U="0.000000"/> - <match seq="CAGG" p="441" n="304" pvalue="1.5e-010" evalue="1.8e-005"/> - <match seq="UAGG" p="232" n="165" pvalue="1.1e-004" evalue="1.3e+001"/> - </motif> - </motifs> \ No newline at end of file
--- a/test-data/fimo_output_test1.gff Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,12 +0,0 @@ -##gff-version 3 -chrM fimo nucleotide_motif 6529 6536 46.6 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-1-chrM;pvalue=2.18e-05;qvalue=0.142;sequence=TACTAACA; -chrM fimo nucleotide_motif 7741 7748 46.6 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-2-chrM;pvalue=2.18e-05;qvalue=0.142;sequence=TACTAACA; -chrM fimo nucleotide_motif 13656 13663 46.6 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-3-chrM;pvalue=2.18e-05;qvalue=0.142;sequence=TACTAACA; -chrM fimo nucleotide_motif 13740 13747 46.6 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-4-chrM;pvalue=2.18e-05;qvalue=0.142;sequence=TACTAACA; -chrM fimo nucleotide_motif 2299 2306 46.6 - . Name=TACTAAYM_chrM-;Alias=MEME-1;ID=TACTAAYM-MEME-1-5-chrM;pvalue=2.18e-05;qvalue=0.142;sequence=TACTAACA; -chrM fimo nucleotide_motif 861 868 44 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-6-chrM;pvalue=3.96e-05;qvalue=0.185;sequence=TACTAACC; -chrM fimo nucleotide_motif 9346 9353 44 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-7-chrM;pvalue=3.96e-05;qvalue=0.185;sequence=TACTAACC; -chrM fimo nucleotide_motif 3767 3774 41.8 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-8-chrM;pvalue=6.62e-05;qvalue=0.216;sequence=TACTAATA; -chrM fimo nucleotide_motif 5497 5504 41.8 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-9-chrM;pvalue=6.62e-05;qvalue=0.216;sequence=TACTAATA; -chrM fimo nucleotide_motif 10105 10112 41.8 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-10-chrM;pvalue=6.62e-05;qvalue=0.216;sequence=TACTAATA; -chrM fimo nucleotide_motif 10959 10966 40.6 + . Name=TACTAAYM_chrM+;Alias=MEME-1;ID=TACTAAYM-MEME-1-11-chrM;pvalue=8.79e-05;qvalue=0.261;sequence=TACTAACT;
--- a/test-data/fimo_output_test1.html Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,34 +0,0 @@ - <td style="padding-right: 2em">background file name = --nrdb--</td> - <td style="padding-left: 5em; padding-right: 2em">alphabet = DNA</td> - <td style="padding-left: 5em; padding-right: 2em">max stored scores = 100000</td> - </tr> <tr> - <td style="padding-right: 2em">allow clobber = false</td> - <td style="padding-left: 5em; padding-right: 2em">compute q-values = true</td> - <td style="padding-left: 5em; padding-right: 2em">parse genomic coord. = false</td> - </tr> - <tr> - <td style="padding-right: 2em">text only = false</td> - <td style="padding-left: 5em; padding-right: 2em">scan both strands = true</td> - <td style="padding-left: 5em; padding-right: 2em">max strand = false</td> - </tr> - <tr> - <td style="padding-right: 2em">threshold type = p-value</td> - <td style="padding-left: 5em; padding-right: 2em">output theshold = 0.0001</td> - <td style="padding-left: 5em; padding-right: 2em">pseudocount = 0.1</td> - </tr> - <tr> - <td style="padding-right: 2em">alpha = 1</td> - <td style="padding-left: 5em; padding-right: 2em">verbosity = 2</td> - <td style="padding-left: 5em; padding-right: 2em"></td> - </tr> - -</table> -</pre> -<p> -This information can be useful in the event you wish to report a -problem with the FIMO software. -</p> -<hr> -<span style="background-color: #DDDDFF"><a href="#top_buttons"><b>Go to top</b></a></span> -</body> -</html>
--- a/test-data/fimo_output_test1.txt Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,12 +0,0 @@ -motif_id motif_alt_id sequence_name start stop strand score p-value q-value matched_sequence -TACTAAYM MEME-1 chrM 6529 6536 + 12.9701 2.18e-05 0.142 TACTAACA -TACTAAYM MEME-1 chrM 7741 7748 + 12.9701 2.18e-05 0.142 TACTAACA -TACTAAYM MEME-1 chrM 13656 13663 + 12.9701 2.18e-05 0.142 TACTAACA -TACTAAYM MEME-1 chrM 13740 13747 + 12.9701 2.18e-05 0.142 TACTAACA -TACTAAYM MEME-1 chrM 2299 2306 - 12.9701 2.18e-05 0.142 TACTAACA -TACTAAYM MEME-1 chrM 861 868 + 12.2836 3.96e-05 0.185 TACTAACC -TACTAAYM MEME-1 chrM 9346 9353 + 12.2836 3.96e-05 0.185 TACTAACC -TACTAAYM MEME-1 chrM 3767 3774 + 11.7164 6.62e-05 0.216 TACTAATA -TACTAAYM MEME-1 chrM 5497 5504 + 11.7164 6.62e-05 0.216 TACTAATA -TACTAAYM MEME-1 chrM 10105 10112 + 11.7164 6.62e-05 0.216 TACTAATA -TACTAAYM MEME-1 chrM 10959 10966 + 11.6567 8.79e-05 0.261 TACTAACT
--- a/test-data/fimo_output_test1.xml Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,39 +0,0 @@ -<setting name="background file name">--nrdb--</setting> -<setting name="allow clobber">false</setting> -<setting name="compute q-values">true</setting> -<setting name="parse genomic coord.">false</setting> -<setting name="text only">false</setting> -<setting name="scan both strands">true</setting> -<setting name="output threshold">0.0001</setting> -<setting name="threshold type">p-value</setting> -<setting name="max stored scores">100000</setting> -<setting name="pseudocount">0.1</setting> -<setting name="verbosity">2</setting> -</settings> -<sequence-data num-sequences="1" num-residues="16569" /> -<alphabet name="DNA" like="dna"> -<letter id="A" symbol="A" complement="T" name="Adenine" colour="CC0000"/> -<letter id="C" symbol="C" complement="G" name="Cytosine" colour="0000CC"/> -<letter id="G" symbol="G" complement="C" name="Guanine" colour="FFB300"/> -<letter id="T" symbol="T" aliases="U" complement="A" name="Thymine" colour="008000"/> -<letter id="N" symbol="N" aliases="X." equals="ACGT" name="Any base"/> -<letter id="V" symbol="V" equals="ACG" name="Not T"/> -<letter id="H" symbol="H" equals="ACT" name="Not G"/> -<letter id="D" symbol="D" equals="AGT" name="Not C"/> -<letter id="B" symbol="B" equals="CGT" name="Not A"/> -<letter id="M" symbol="M" equals="AC" name="Amino"/> -<letter id="R" symbol="R" equals="AG" name="Purine"/> -<letter id="W" symbol="W" equals="AT" name="Weak"/> -<letter id="S" symbol="S" equals="CG" name="Strong"/> -<letter id="Y" symbol="Y" equals="CT" name="Pyrimidine"/> -<letter id="K" symbol="K" equals="GT" name="Keto"/> -</alphabet> -<motif name="TACTAAYM" alt="MEME-1" width="8" best-possible-match="TACTAACA"/> -<background source="--nrdb--"> -<value letter="A">0.275</value> -<value letter="C">0.225</value> -<value letter="G">0.225</value> -<value letter="T">0.275</value> -</background> -<cisml-file>cisml.xml</cisml-file> -</fimo>
--- a/test-data/fimo_output_test2.gff Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,12 +0,0 @@ -##gff-version 3 -chrM fimo nucleotide_motif 440 446 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-1-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 2093 2099 40.9 - . Name=ACTAAYH_chrM-;Alias=DREME-1;ID=ACTAAYH-DREME-1-2-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 2299 2305 40.9 - . Name=ACTAAYH_chrM-;Alias=DREME-1;ID=ACTAAYH-DREME-1-3-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 5186 5192 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-4-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 6530 6536 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-5-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 7742 7748 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-6-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 13657 13663 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-7-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 13741 13747 40.9 + . Name=ACTAAYH_chrM+;Alias=DREME-1;ID=ACTAAYH-DREME-1-8-chrM;pvalue=8.2e-05;qvalue=0.332;sequence=ACTAACA; -chrM fimo nucleotide_motif 510 516 43.8 + . Name=CCAGCAY_chrM+;Alias=DREME-5;ID=CCAGCAY-DREME-5-1-chrM;pvalue=4.15e-05;qvalue=0.683;sequence=CCAGCAC; -chrM fimo nucleotide_motif 5137 5143 43.8 + . Name=CCAGCAY_chrM+;Alias=DREME-5;ID=CCAGCAY-DREME-5-2-chrM;pvalue=4.15e-05;qvalue=0.683;sequence=CCAGCAC; -chrM fimo nucleotide_motif 4241 4247 40.3 + . Name=CCAGCAY_chrM+;Alias=DREME-5;ID=CCAGCAY-DREME-5-3-chrM;pvalue=9.37e-05;qvalue=1;sequence=CCAGCAT;
--- a/test-data/fimo_output_test2.html Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,33 +0,0 @@ - <td style="padding-left: 5em; padding-right: 2em">alphabet = DNA</td> - <td style="padding-left: 5em; padding-right: 2em">max stored scores = 100000</td> - </tr> <tr> - <td style="padding-right: 2em">allow clobber = false</td> - <td style="padding-left: 5em; padding-right: 2em">compute q-values = true</td> - <td style="padding-left: 5em; padding-right: 2em">parse genomic coord. = false</td> - </tr> - <tr> - <td style="padding-right: 2em">text only = false</td> - <td style="padding-left: 5em; padding-right: 2em">scan both strands = true</td> - <td style="padding-left: 5em; padding-right: 2em">max strand = false</td> - </tr> - <tr> - <td style="padding-right: 2em">threshold type = p-value</td> - <td style="padding-left: 5em; padding-right: 2em">output theshold = 0.0001</td> - <td style="padding-left: 5em; padding-right: 2em">pseudocount = 0.1</td> - </tr> - <tr> - <td style="padding-right: 2em">alpha = 1</td> - <td style="padding-left: 5em; padding-right: 2em">verbosity = 2</td> - <td style="padding-left: 5em; padding-right: 2em"></td> - </tr> - -</table> -</pre> -<p> -This information can be useful in the event you wish to report a -problem with the FIMO software. -</p> -<hr> -<span style="background-color: #DDDDFF"><a href="#top_buttons"><b>Go to top</b></a></span> -</body> -</html>
--- a/test-data/fimo_output_test2.txt Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,12 +0,0 @@ -motif_id motif_alt_id sequence_name start stop strand score p-value q-value matched_sequence -CCAGCAY DREME-5 chrM 510 516 + 13.5843 4.15e-05 0.683 CCAGCAC -CCAGCAY DREME-5 chrM 5137 5143 + 13.5843 4.15e-05 0.683 CCAGCAC -ACTAAYH DREME-1 chrM 440 446 + 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 2093 2099 - 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 2299 2305 - 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 5186 5192 + 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 6530 6536 + 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 7742 7748 + 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 13657 13663 + 11.7286 8.2e-05 0.332 ACTAACA -ACTAAYH DREME-1 chrM 13741 13747 + 11.7286 8.2e-05 0.332 ACTAACA -CCAGCAY DREME-5 chrM 4241 4247 + 13.1461 9.37e-05 1 CCAGCAT \ No newline at end of file
--- a/test-data/fimo_output_test2.xml Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,34 +0,0 @@ -<setting name="compute q-values">true</setting> -<setting name="parse genomic coord.">false</setting> -<setting name="text only">false</setting> -<setting name="scan both strands">true</setting> -<setting name="output threshold">0.0001</setting> -<setting name="threshold type">p-value</setting> -<setting name="max stored scores">100000</setting> -<setting name="pseudocount">0.1</setting> -<setting name="verbosity">2</setting> -</settings> -<sequence-data num-sequences="1" num-residues="16569" /> -<alphabet name="DNA" like="dna"> -<letter id="A" symbol="A" complement="T" name="Adenine" colour="CC0000"/> -<letter id="C" symbol="C" complement="G" name="Cytosine" colour="0000CC"/> -<letter id="G" symbol="G" complement="C" name="Guanine" colour="FFB300"/> -<letter id="T" symbol="T" aliases="U" complement="A" name="Thymine" colour="008000"/> -<letter id="N" symbol="N" aliases="X." equals="ACGT" name="Any base"/> -<letter id="V" symbol="V" equals="ACG" name="Not T"/> -<letter id="H" symbol="H" equals="ACT" name="Not G"/> -<letter id="D" symbol="D" equals="AGT" name="Not C"/> -<letter id="B" symbol="B" equals="CGT" name="Not A"/> -<letter id="M" symbol="M" equals="AC" name="Amino"/> -<letter id="R" symbol="R" equals="AG" name="Purine"/> -<letter id="W" symbol="W" equals="AT" name="Weak"/> -<letter id="S" symbol="S" equals="CG" name="Strong"/> -<letter id="Y" symbol="Y" equals="CT" name="Pyrimidine"/> -<letter id="K" symbol="K" equals="GT" name="Keto"/> -</alphabet> -<motif name="ACTAAYH" alt="DREME-1" width="7" best-possible-match="ACTAACA"/> -<motif name="YTAACA" alt="DREME-2" width="6" best-possible-match="TTAACA"/> -<motif name="TCTGT" alt="DREME-3" width="5" best-possible-match="TCTGT"/> -<motif name="SCCAGG" alt="DREME-4" width="6" best-possible-match="CCCAGG"/> -<motif name="CCAGCAY" alt="DREME-5" width="7" best-possible-match="CCAGCAC"/> -<motif name="GMATGT" alt="DREME-6" width="6" best-possible-match="GAATGT"/> \ No newline at end of file
--- a/test-data/fimo_output_test3.html Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,23 +0,0 @@ - <td style="padding-right: 2em">background file name = --uniform--</td> - <td style="padding-left: 5em; padding-right: 2em">alphabet = DNA</td> - <td style="padding-left: 5em; padding-right: 2em">max stored scores = 100000</td> - </tr> <tr> - <td style="padding-right: 2em">allow clobber = false</td> - <td style="padding-left: 5em; padding-right: 2em">compute q-values = true</td> - <td style="padding-left: 5em; padding-right: 2em">parse genomic coord. = false</td> - </tr> - <tr> - <td style="padding-right: 2em">text only = false</td> - <td style="padding-left: 5em; padding-right: 2em">scan both strands = true</td> - <td style="padding-left: 5em; padding-right: 2em">max strand = false</td> - </tr> - <tr> - <td style="padding-right: 2em">threshold type = p-value</td> - <td style="padding-left: 5em; padding-right: 2em">output theshold = 0.01</td> - <td style="padding-left: 5em; padding-right: 2em">pseudocount = 0.1</td> - </tr> - <tr> - <td style="padding-right: 2em">alpha = 1</td> - <td style="padding-left: 5em; padding-right: 2em">verbosity = 2</td> - <td style="padding-left: 5em; padding-right: 2em"></td> - </tr> \ No newline at end of file
--- a/test-data/fimo_output_test3.txt Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,40 +0,0 @@ -motif_id motif_alt_id sequence_name start stop strand score p-value q-value matched_sequence -ACTAAYH DREME-1 chrM 440 446 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 2093 2099 - 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 2299 2305 - 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 5186 5192 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 6530 6536 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 7742 7748 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 13657 13663 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 13741 13747 + 12.1831 6.1e-05 0.239 ACTAACA -ACTAAYH DREME-1 chrM 3768 3774 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 5498 5504 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 7736 7742 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 9872 9878 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 10106 10112 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 10313 10319 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 11818 11824 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 15903 15909 + 11.5915 0.000122 0.239 ACTAATA -ACTAAYH DREME-1 chrM 7732 7738 + 11.4507 0.000183 0.294 ACTAACT -ACTAAYH DREME-1 chrM 10960 10966 + 11.4507 0.000183 0.294 ACTAACT -ACTAAYH DREME-1 chrM 862 868 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 1832 1838 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 8679 8685 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 8770 8776 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 9347 9353 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 9359 9365 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 10302 10308 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 14765 14771 + 11.2394 0.000244 0.294 ACTAACC -ACTAAYH DREME-1 chrM 4466 4472 + 10.8592 0.000305 0.302 ACTAATT -ACTAAYH DREME-1 chrM 11248 11254 + 10.8592 0.000305 0.302 ACTAATT -ACTAAYH DREME-1 chrM 475 481 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 7930 7936 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 8649 8655 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 8670 8676 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 10768 10774 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 11053 11059 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 11101 11107 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 12720 12726 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 13149 13155 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 14216 14222 + 10.6479 0.000366 0.302 ACTAATC -ACTAAYH DREME-1 chrM 678 684 - -1.07042 0.000549 0.366 ACTAAGA \ No newline at end of file
--- a/test-data/fimo_output_test3.xml Sat Apr 09 08:34:02 2022 +0000 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,32 +0,0 @@ -<sequence-data num-sequences="1" num-residues="16569" /> -<alphabet name="DNA" like="dna"> -<letter id="A" symbol="A" complement="T" name="Adenine" colour="CC0000"/> -<letter id="C" symbol="C" complement="G" name="Cytosine" colour="0000CC"/> -<letter id="G" symbol="G" complement="C" name="Guanine" colour="FFB300"/> -<letter id="T" symbol="T" aliases="U" complement="A" name="Thymine" colour="008000"/> -<letter id="N" symbol="N" aliases="X." equals="ACGT" name="Any base"/> -<letter id="V" symbol="V" equals="ACG" name="Not T"/> -<letter id="H" symbol="H" equals="ACT" name="Not G"/> -<letter id="D" symbol="D" equals="AGT" name="Not C"/> -<letter id="B" symbol="B" equals="CGT" name="Not A"/> -<letter id="M" symbol="M" equals="AC" name="Amino"/> -<letter id="R" symbol="R" equals="AG" name="Purine"/> -<letter id="W" symbol="W" equals="AT" name="Weak"/> -<letter id="S" symbol="S" equals="CG" name="Strong"/> -<letter id="Y" symbol="Y" equals="CT" name="Pyrimidine"/> -<letter id="K" symbol="K" equals="GT" name="Keto"/> -</alphabet> -<motif name="ACTAAYH" alt="DREME-1" width="7" best-possible-match="ACTAACA"/> -<motif name="YTAACA" alt="DREME-2" width="6" best-possible-match="TTAACA"/> -<motif name="TCTGT" alt="DREME-3" width="5" best-possible-match="TCTGT"/> -<motif name="SCCAGG" alt="DREME-4" width="6" best-possible-match="CCCAGG"/> -<motif name="CCAGCAY" alt="DREME-5" width="7" best-possible-match="CCAGCAC"/> -<motif name="GMATGT" alt="DREME-6" width="6" best-possible-match="GAATGT"/> -<background source="--uniform--"> -<value letter="A">0.250</value> -<value letter="C">0.250</value> -<value letter="G">0.250</value> -<value letter="T">0.250</value> -</background> -<cisml-file>cisml.xml</cisml-file> -</fimo>