|
3
|
1 <tool id="chip_seq_ratio_1" name="Chip-Seq Mapper" version="0.1.1">
|
|
|
2 <stdio>
|
|
|
3 <exit_code range="1:" level="fatal" description="Error"/>
|
|
|
4 </stdio>
|
|
|
5 <description></description>
|
|
|
6 <requirements>
|
|
|
7 <requirement type="package">r-base64enc</requirement>
|
|
|
8 <requirement type="package">r-r2html</requirement>
|
|
|
9 <requirement type="package">blast</requirement>
|
|
6
|
10 <!-- <requirement type="package">chip_seq_ration</requirement> -->
|
|
3
|
11 </requirements>
|
|
|
12 <command interpreter="python3">
|
|
|
13 ChipSeqRatioAnalysis.py
|
|
|
14 --ChipSeq=${ChipFile}
|
|
|
15 --InputSeq=${InputFile}
|
|
|
16 --Contigs=${ContigFile}
|
|
|
17 --output=${OutputFile}
|
|
|
18 --html=${ReportFile}
|
|
|
19 --max_cl=${MaxCl}
|
|
|
20 --bitscore=$bitscore
|
|
|
21 --nproc=16
|
|
|
22 </command>
|
|
|
23
|
|
|
24 <inputs>
|
|
|
25 <param name="ChipFile" label="Chip Sequences" type="data" format="fasta" help="NGS data in fasta format"/>
|
|
|
26 <param name="InputFile" label="Input Sequences" type="data" format="fasta" help="NGS data in fasta format"/>
|
|
|
27 <param name="ContigFile" label="Reference - Contig Sequences" type="data" format="fasta"
|
|
|
28 help="Contigs obtained from RepeatExplorer clustering pipeline in fasta file"/>
|
|
|
29 <param name="MaxCl" label="Number of clusters to be shown in graph" type="integer" value="200"/>
|
|
6
|
30 <param name="bitscore" label="Minimum bit score threshold" type="integer" value="50" help="All similarity hits with lower bit score will not be considered for ChIP/Input ratio calculation"/>
|
|
3
|
31 </inputs>
|
|
|
32 <outputs>
|
|
7
|
33 <data name="OutputFile" format="tabular"
|
|
|
34 label="csv table from ChIP-Seq-Mapper on datasets ${InputFile.hid} (Input) ${ChipFile.hid} (ChIP) and ${ContigFile.hid} (reference)"/>
|
|
|
35
|
|
|
36 <data name="ReportFile" format="html"
|
|
|
37 label="HTML report from ChIP-Seq-Mapper on datasets ${InputFile.hid} (Input) ${ChipFile.hid} (ChIP) and ${ContigFile.hid} (reference)"/>
|
|
3
|
38 </outputs>
|
|
|
39
|
|
|
40 <help>
|
|
|
41 **What it does**
|
|
|
42
|
|
|
43 Analysis of NGS sequences from Chromatin Imunoprecipitation. ChiP
|
|
|
44 and Input reads are mapped to contigs obtained from graph based
|
|
|
45 repetitive sequence clustering(`Novak et al. 2013`__) to enriched repeats. Reads from input
|
|
|
46 and ChIP should be ideally short illumina reads with uniform length
|
|
|
47 above 80 nt. It is sufficiant to use about 1 milion of reads for both Input and Chip.
|
|
|
48 This method was first used in (`Neumann et al. 2012`__) for
|
|
|
49 identification of repetitive sequences associated with cetromeric
|
|
|
50 region. If you use this method, reference:
|
|
|
51
|
|
|
52
|
|
|
53 `PLoS Genet. Epub 2012 Jun 21. Stretching the rules: monocentric chromosomes with multiple centromere domains. Neumann P, Navrátilová A, Schroeder-Reiter E, Koblížková A, Steinbauerová V, Chocholová E, Novák P, Wanner G, Macas J.`__.
|
|
|
54
|
|
|
55 .. __: http://bioinformatics.oxfordjournals.org/content/29/6/792.full
|
|
|
56
|
|
|
57 .. __: http://journals.plos.org/plosgenetics/article?id=10.1371/journal.pgen.1002777
|
|
|
58 .. __: http://journals.plos.org/plosgenetics/article?id=10.1371/journal.pgen.1002777
|
|
|
59
|
|
|
60 </help>
|
|
|
61
|
|
|
62 </tool>
|
|
|
63
|
|
|
64
|
