Mercurial > repos > xuebing > sharplabtool

comparison tools/regVariation/compute_motif_frequencies_for_all_motifs.xml @ 0:9071e359b9a3

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author xuebing
date Fri, 09 Mar 2012 19:37:19 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:9071e359b9a3
1 <tool id="compute_motif_frequencies_for_all_motifs" name="Compute Motif Frequencies For All Motifs" version="1.0.0">
2 <description>motif by motif</description>
3
4 <command interpreter="perl">
5 compute_motif_frequencies_for_all_motifs.pl $inputFile1 $inputFile2 $inputWindowSize3 $outputFile1
6 </command>
7
8 <inputs>
9 <param format="tabular" name="inputFile1" type="data" label="Select the motifs file"/>
10 <param format="tabular" name="inputFile2" type="data" label="Select the indel flanking sequences windows file"/>
11 <param type="integer" name="inputWindowSize3" size="6" value="0" label="What is the number of 10bp windows in which the motif frequencies will be computed?" help="'0' = one window only"/>
12 </inputs>
13
14 <outputs>
15 <data format="tabular" name="outputFile1"/>
16 </outputs>
17
18 <tests>
19 <test>
20 <param name="inputFile1" value="motifs2.tabular" />
21 <param name="inputFile2" value="flankingSequencesWindows10_2.tabular" />
22 <param name="inputWindowSize3" value="0" />
23 <output name="outputFile1" file="motifFrequencies_every_indels0.tabular" />
24 </test>
25
26 <test>
27 <param name="inputFile1" value="motifs2.tabular" />
28 <param name="inputFile2" value="flankingSequencesWindows10_2.tabular" />
29 <param name="inputWindowSize3" value="4" />
30 <output name="outputFile1" file="motifFrequencies_every_indels4.tabular" />
31 </test>
32 </tests>
33
34 <help>
35
36 .. class:: infomark
37
38 **What it does**
39
40 This program computes the frequencies of each motif at a window size, determined by the user, in both upstream and downstream sequences flanking indels in all chromosomes.
41
42 - The first input is a TABULAR format file containing the motif names and sequences, one line per motif, such that the file consists of two columns:
43
44 - The left column represents the motif names
45 - The right column represents the motif sequence, as follows::
46
47 dnaPolPauseFrameshift1 GAG
48 dnaPolPauseFrameshift2 ACG
49 xSites1 CCG
50
51 - The second input is a TABULAR format file representing the windows of both upstream and downstream flanking sequences. It consists of multiple left columns representing the windows of the upstream flanking sequences, followed by one column representing the indels, then followed by multiple right columns representing the windows of the downstream flanking sequences, as follows::
52
53 cgaggtcagg agatcgagac catcctggct aacatggtga aatcccgtct ctactaaaaa indel aaatttatat ttataaacaa ttttaataca cctatgttta ttatacattt
54 GCCAGTTTAT GGTCTAACAA GGAGAGAAAC AGGGGGCTGA AGGGGTTTCT TAACCTCCAG indel TTCCGGGCTC TGTCCCTAAC CCCCAGCTAG GTAAGTGGCA AAGCACTTCT
55 CAGTGGGACC AAGCACTGAA CCACTTTGGG GAGAATCTCA CACTGGGGCC CTCTGACACC indel tatatatttt tttttttttt tttttttttt tttttttttg agatggtgtc
56 AGAGCAGCAG CACCCACTTT TGCAGTGTGT GACGTTGGTG GAGCCATCGA AGTCTGTGCT indel GAGCCCTCCC CAGTGCTCCG AGGAGCTGCT GTTCCCCCTG GAGCTCAGAA
57
58 - The third input is an integer number representing the number of windows to be considered starting from the indel and leftward for the upstream flanking sequence and, starting from the indel and rightward for the downstream flanking sequence.
59
60 - The output is a TABULAR format file consisting of three columns:
61
62 - The left column represents the motif name
63 - The middle column represents the motif frequency in the specified windows of the upstream sequence flanking an indel
64 - The right column represents the motif frequency in the specified windows of the downstream sequence flanking an indel
65
66 There is line per indel in the output file, such that the total number of lines in the output file = number of motifs x number of indels.
67
68 Note: The number of windows entered by the user must be a positive integer >= 1. if negative integer or 0 is entered by the user, the program will consider it as 1.
69
70 </help>
71
72 </tool>