24
|
1 <tool id="PhageTerm" name="PhageTerm" version="1.0.12">
|
|
2 <description>Determine Bacteriophage Termini and Packaging Mode using randomly fragmented NGS data.</description>
|
|
3 <requirements>
|
|
4 <requirement type="package" version="2.7">python</requirement>
|
|
5 <requirement type="package" version="1.16.5">numpy</requirement>
|
|
6 <requirement type="package" version="2.2.5">matplotlib</requirement>
|
|
7 <requirement type="package" version="0.24.2">pandas</requirement>
|
|
8 <requirement type="package" version="0.20.4">scikit-learn</requirement>
|
|
9 <requirement type="package" version="1.2.1">scipy</requirement>
|
|
10 <requirement type="package" version="0.10.2">statsmodels</requirement>
|
|
11 <requirement type="package" version="3.5.42">reportlab</requirement>
|
|
12 </requirements>
|
|
13 <edam_operations>
|
|
14 <edam_operation>operation_3198</edam_operation>
|
|
15 </edam_operations>
|
|
16 <edam_topics>
|
|
17 <edam_topic>topic_0781</edam_topic>
|
|
18 </edam_topics>
|
|
19 <stdio>
|
|
20 <!-- Assume anything other than zero is an error -->
|
|
21 <exit_code range="1:" />
|
|
22 </stdio>
|
|
23 <command interpreter="python2.7">PhageTerm.py -c \${GALAXY_SLOTS} -f $fastq -r $reference
|
|
24 #if $phagename
|
|
25 -n $phagename
|
|
26 #end if
|
|
27 #if $paired
|
|
28 -p $paired
|
|
29 #end if
|
|
30 #if $seed
|
|
31 -s $seed
|
|
32 #end if
|
|
33 #if $surrounding
|
|
34 -d $surrounding
|
|
35 #end if
|
|
36 #if $host
|
|
37 -g $host
|
|
38 #end if
|
|
39 #if $coverage
|
|
40 -m $coverage
|
|
41 #end if
|
|
42 </command>
|
|
43 <inputs>
|
|
44 <param format="fastq" name="fastq" type="data" label="Phage Reads (FASTQ)" help="Phage sequencing reads file in fastq format from randomly fragmented NGS data (e.g. Illumina TruSeq) but NOT Nextera."/>
|
|
45 <param format="fasta" name="reference" type="data" label="Phage Genome (FASTA)" help="Phage reference genome file in fasta format, Multi-fasta NOT accepted."/>
|
|
46 <param name="phagename" type="text" value="my_phage" size="20" optional="true" label="Prefix" help="Prefix for the output file name."/>
|
|
47 <param format="fastq" name="paired" type="data" optional="true" label="Phage Paired-End Reads (FASTQ)" help="Phage sequencing paired-end reads file in fastq format."/>
|
|
48 <param format="fasta" name="host" type="data" optional="true" label="Bacterial Host Genome (FASTA)" help="Host reference genome file in fasta format, Multi-fasta NOT accepted. Warning : increases process time."/>
|
|
49 <param name="seed" type="integer" optional="true" value="20" label="Seed length" help="Seed length value for alignment of reads."/>
|
|
50 <param name="surrounding" type="integer" optional="true" value="20" label="Peak surrounding region" help="Length of the surrounding region defining close peaks to be merged in the analysis process."/>
|
|
51 <param name="coverage" type="integer" optional="true" value="250" label="Limit Coverage" help="Phage upper limit coverage."/>
|
|
52 </inputs>
|
|
53 <outputs>
|
|
54 <data name="Statistics" format="csv" label="${phagename}_statistics.csv" from_work_dir="*_statistics.csv" />
|
|
55 <data name="CohesiveSequence" format="fasta" label="${phagename}_cohesive-DTR.fasta" from_work_dir="*_cohesive-sequence.fasta"/>
|
|
56 <data name="PhageSequence" format="fasta" label="${phagename}_sequence.fasta" from_work_dir="*_sequence.fasta"/>
|
|
57 <data name="ReportPDF" format="pdf" label="${phagename}_PhageTerm_report.pdf" from_work_dir="*_PhageTerm_report.pdf"/>
|
|
58 </outputs>
|
|
59 <tests>
|
|
60 <test>
|
|
61 <param name="phagename" value="Test-cohesive-5'" />
|
|
62 <param name="fastq" value="COS-5.fastq" />
|
|
63 <param name="reference" value="COS-5.fasta" />
|
|
64 </test>
|
|
65 <test>
|
|
66 <param name="phagename" value="Test-cohesive-3'" />
|
|
67 <param name="fastq" value="COS-3.fastq" />
|
|
68 <param name="reference" value="COS-3.fasta" />
|
|
69 </test>
|
|
70 <test>
|
|
71 <param name="phagename" value="Test-short-direct-terminal-repeats" />
|
|
72 <param name="fastq" value="DTR-short.fastq" />
|
|
73 <param name="reference" value="DTR-short.fasta" />
|
|
74 </test>
|
|
75 <test>
|
|
76 <param name="phagename" value="Test-long-direct-terminal-repeats" />
|
|
77 <param name="fastq" value="DTR-long.fastq" />
|
|
78 <param name="reference" value="DTR-long.fasta" />
|
|
79 </test>
|
|
80 <test>
|
|
81 <param name="phagename" value="Test-Headful" />
|
|
82 <param name="fastq" value="Headful.fastq" />
|
|
83 <param name="reference" value="Headful.fasta" />
|
|
84 <param name="surrounding" value="0" />
|
|
85 </test>
|
|
86 <test>
|
|
87 <param name="phagename" value="Test-Mu-like" />
|
|
88 <param name="fastq" value="Mu-like_R1.fastq" />
|
|
89 <param name="paired" value="Mu-like_R2.fastq" />
|
|
90 <param name="reference" value="Mu-like.fasta" />
|
|
91 <param name="host" value="Mu-like_host.fasta" />
|
|
92 <param name="surrounding" value="0" />
|
|
93 </test>
|
|
94 </tests>
|
|
95 <help>
|
|
96
|
|
97 **What it does**
|
|
98
|
|
99 PhageTerm software is a tool to determine phage termini and packaging mode from high throughput sequences
|
|
100 that rely on the random fragmentation of DNA (e.g. Illumina TruSeq, but NOT Nextera).
|
|
101
|
|
102 -----
|
|
103
|
|
104 **Methods**
|
|
105
|
|
106 Our methods is based on the detection of biases in the number of reads that start at natural DNA termini compared to elsewhere along the phage genome.
|
|
107 Preprint manuscript: http://biorxiv.org/content/early/2017/02/16/108100
|
|
108 Source Code: https://sourceforge.net/projects/phageterm.
|
|
109
|
|
110 -----
|
|
111
|
|
112 **Inputs**
|
|
113
|
|
114 - Raw reads file in fastq format
|
|
115 - Phage genome in fasta format
|
|
116 - Raw paired-ends reads file in fastq format (optional)
|
|
117 - Host genome in fasta format (optional)
|
|
118
|
|
119 -----
|
|
120
|
|
121 **Options**
|
|
122
|
|
123 Name of the phage being analyzed by the user (Default: my_phage)
|
|
124 Length of the seed used for reads in the mapping process (Default: 20)
|
|
125 Length of the surrounding used to merge close peaks in the analysis process (Default: 20)
|
|
126 Define phage target mean coverage (Default: 250)
|
|
127
|
|
128 -----
|
|
129
|
|
130 **Outputs**
|
|
131
|
|
132 - Report (.pdf)
|
|
133 - Statistical table (.csv)
|
|
134 - Phage Sequence (.fasta)
|
|
135 - Cohesive or DTR sequence (.fasta)
|
|
136
|
|
137
|
|
138 </help>
|
|
139 <citations>
|
|
140 <citation type="doi">10.1038/s41598-017-07910-5</citation>
|
|
141 </citations>
|
|
142 </tool>
|