|
0
|
1 <tool id="alignCustomAmplicon" name="Align Custom Amplicon" version="0.0.1">
|
|
|
2 <description>align amplicon to reference with primers</description>
|
|
1
|
3 <requirements>
|
|
|
4 <requirement type="package" version="1.56.0">picard</requirement>
|
|
|
5 <requirement type="package" version="0.1.18">samtools</requirement>
|
|
|
6 <requirement type="package" version="1.3.12">gzip</requirement>
|
|
|
7 <requirement type="perl-module" version="0.42">Inline-CPP</requirement>
|
|
|
8 </requirements>
|
|
0
|
9 <command interpreter="perl">
|
|
|
10 alignCustomAmplicon.pl -s -r -p 4 -o $output $refFile.fields.path $read1 $read2 $primers
|
|
|
11 </command>
|
|
|
12
|
|
|
13 <inputs>
|
|
|
14 <param name="refFile" type="select" label="Select a reference genome">
|
|
|
15 <options from_data_table="all_fasta">
|
|
|
16 <filter type="sort_by" column="2" />
|
|
|
17 <validator type="no_options" message="No indexes are available" />
|
|
|
18 </options>
|
|
|
19 </param>
|
|
|
20 <param name="read1" type="data" format="fastqsanger,fastqillumina, fastq" label="FASTQ read 1 " help="FASTQ with either Sanger-scaled quality values (fastqsanger) or Illumina-scaled quality values (fastqillumina)" />
|
|
|
21 <param name="read2" type="data" format="fastqsanger,fastqillumina, fastq" label="FASTQ read 2" help="FASTQ with either Sanger-scaled quality values (fastqsanger) or Illumina-scaled quality values (fastqillumina)" />
|
|
|
22 <param name="primers" type="data" format="tabular" label="Primers" help="Primers location and length"/>
|
|
|
23 </inputs>
|
|
|
24
|
|
|
25 <outputs>
|
|
|
26 <data type="data" format="bam" name="output"/>
|
|
|
27 </outputs>
|
|
|
28
|
|
|
29 <help>
|
|
|
30
|
|
|
31 .. class:: infomark
|
|
|
32
|
|
|
33 **What it does**
|
|
|
34
|
|
|
35 It is an amplicon aligner that uses primers for higher accuracy.
|
|
|
36
|
|
|
37 Reads with primers are aligned to the reference, then primers are discarded.
|
|
|
38
|
|
|
39 If both reads are long enough, they are aligned with the reference and a consensus alignment is generated.
|
|
|
40
|
|
|
41 Otherwise, each read is aligned separately.
|
|
|
42
|
|
|
43 Sequences with bad quality reads are discarded.
|
|
|
44
|
|
|
45
|
|
|
46
|
|
|
47 **Input**
|
|
|
48
|
|
|
49 ref:
|
|
|
50
|
|
|
51 Fasta file of ref gnome
|
|
|
52
|
|
|
53 read1:
|
|
|
54
|
|
|
55 Fastq file of left to right read
|
|
|
56 (Can also be compressed [fastq.gz])
|
|
|
57
|
|
|
58 read2:
|
|
|
59
|
|
|
60 Fastq file of right to left read
|
|
|
61 (Can also be compressed [fastq.gz])
|
|
|
62
|
|
|
63 primers:
|
|
|
64
|
|
|
65 Text file with primers name and length (see example)
|
|
|
66
|
|
|
67 Example primers format::
|
|
|
68
|
|
|
69 #Name_of_amplicon length_left length_right
|
|
|
70 1:115256345-115256520 23 23
|
|
|
71 1:115256436-115256606 25 22
|
|
|
72 1:115256530-115256724 23 23
|
|
|
73 1:115256532-115256723 23 23
|
|
|
74 4:55151914-55152086 21 23
|
|
|
75 4:55151935-55152132 20 23
|
|
|
76 4:55151991-55152182 23 24
|
|
|
77 4:55591944-55592136 23 24
|
|
|
78 4:55592065-55592263 20 23
|
|
|
79 4:55593504-55593674 24 25
|
|
|
80 ...
|
|
|
81
|
|
|
82
|
|
|
83 </help>
|
|
|
84 </tool>
|
