|
2
|
1 <tool id="crispresso2" name="CRISPResso2" version="0.1.1" python_template_version="3.5">
|
|
0
|
2 <requirements>
|
|
|
3 <requirement type="package" version="2.0.45">crispresso2</requirement>
|
|
|
4 </requirements>
|
|
|
5 <command detect_errors="exit_code"><![CDATA[
|
|
|
6 mkdir crispresso_out;
|
|
|
7 #if $singlePaired.fastq_r1.is_of_type('fastq.gz', 'fastqsanger.gz'):
|
|
|
8 #set $r1 = 'seq_name.fastq.gz'
|
|
|
9 #else:
|
|
|
10 #set $r1 = 'seq_name.fastq'
|
|
|
11 #end if
|
|
|
12 ln -s $singlePaired.fastq_r1 $r1;
|
|
|
13 mkdir -p '${html_file.files_path}' &&
|
|
|
14 #if str($singlePaired.sPaired) == 'paired':
|
|
|
15 #if $singlePaired.fastq_r2.is_of_type('fastq.gz', 'fastqsanger.gz'):
|
|
1
|
16 #set $r2 = 'seq_name2.fastq.gz'
|
|
0
|
17 #else:
|
|
1
|
18 #set $r2 = 'seq_name2.fastq'
|
|
0
|
19 #end if
|
|
|
20 ln -s $singlePaired.fastq_r2 $r2;
|
|
|
21 #end if
|
|
|
22 CRISPResso --fastq_r1 $r1
|
|
|
23 #if str($singlePaired.sPaired) == 'paired':
|
|
|
24 --fastq_r2 $r2
|
|
|
25 #end if
|
|
|
26 --amplicon_seq '$amplicon_seq'
|
|
|
27 -an '$amplicon_name'
|
|
1
|
28 -n reads
|
|
|
29 #if $sgrna_parameters.guide_seq:
|
|
|
30 --guide_seq $sgrna_parameters.guide_seq
|
|
|
31 #end if
|
|
0
|
32 #if $sgrna_parameters.guide_name:
|
|
|
33 --guide_name $sgrna_parameters.guide_name
|
|
|
34 #end if
|
|
|
35 #if $sgrna_parameters.flexiguide:
|
|
|
36 -fg $sgrna_parameters.flexiguide
|
|
|
37 #end if
|
|
|
38 #if $sgrna_parameters.flexiguide_homology:
|
|
|
39 --flexiguide_homology $sgrna_parameters.flexiguide_homology
|
|
|
40 #end if
|
|
|
41 #if $sgrna_parameters.flexiguide_name:
|
|
|
42 --flexiguide_name '$sgrna_parameters.flexiguide_name'
|
|
|
43 #end if
|
|
|
44 #if $coding_seq:
|
|
|
45 --coding_seq '$coding_seq'
|
|
|
46 #end if
|
|
|
47 $sgrna_parameters.discard_guide_positions_overhanging_amplicon_edge
|
|
|
48 $filtering_parameters.split_interleaved_input
|
|
|
49 #if $filtering_parameters.min_average_read_quality:
|
|
|
50 --min_average_read_quality $filtering_parameters.min_average_read_quality
|
|
|
51 #end if
|
|
|
52 #if $filtering_parameters.min_single_bp_quality:
|
|
|
53 --min_single_bp_quality $filtering_parameters.min_single_bp_quality
|
|
|
54 #end if
|
|
|
55 #if $filtering_parameters.min_bp_quality_or_N:
|
|
|
56 --min_bp_quality_or_N $filtering_parameters.min_bp_quality_or_N
|
|
|
57 #end if
|
|
|
58 #if $filtering_parameters.min_paired_end_reads_overlap:
|
|
|
59 --min_paired_end_reads_overlap $filtering_parameters.min_paired_end_reads_overlap
|
|
|
60 #end if
|
|
|
61 #if $filtering_parameters.max_paired_end_reads_overlap:
|
|
|
62 --max_paired_end_reads_overlap $filtering_parameters.max_paired_end_reads_overlap
|
|
|
63 #end if
|
|
2
|
64 #if $window_parameters.cleavage_offset:
|
|
|
65 --cleavage_offset '$window_parameters.cleavage_offset'
|
|
|
66 #end if
|
|
|
67 #if $output_parameters.suppress_report:
|
|
|
68 --suppress_report
|
|
|
69 #end if
|
|
0
|
70 $filtering_parameters.stringent_flash_merging
|
|
|
71 -o '${html_file.files_path}' > $output_log
|
|
1
|
72 && cp '${html_file.files_path}'/*\.html crispresso.html && cp '${html_file.files_path}'/CRISPResso_on_reads/CRISPResso_quantification_of_editing_frequency.txt $quant_editing_freq
|
|
0
|
73 ]]></command>
|
|
|
74 <inputs>
|
|
|
75 <conditional name="singlePaired">
|
|
|
76 <param name="sPaired" type="select" label="Single-end or paired-end reads">
|
|
|
77 <option value="single" selected="true">Single-end</option>
|
|
|
78 <option value="paired">Paired-end</option>
|
|
|
79 </param>
|
|
|
80 <when value="single">
|
|
|
81 <param format="fastq,fastq.gz,fastqsanger.gz" name="fastq_r1" type="data" label="FASTQ file"/>
|
|
|
82 </when>
|
|
|
83 <when value="paired">
|
|
|
84 <param format="fastq,fastq.gz,fastqsanger.gz" name="fastq_r1" type="data" label="FASTQ file, forward reads"/>
|
|
|
85 <param format="fastq,fastq.gz,fastqsanger.gz" name="fastq_r2" type="data" label="FASTQ file, reverse reads"/>
|
|
|
86 </when>
|
|
|
87 </conditional>
|
|
|
88 <param name="amplicon_seq" type="text" label="Amplicon sequence" help="If submitting more than 1 amplicon, use commas to separate them"/>
|
|
|
89 <param name="amplicon_name" type="text" label="Amplicon name" help="--amplicon_name"/>
|
|
|
90 <param name="expected_hdr_amplicon_seq" type="text" label="Amplicon sequence expected after HDR" help="--expected_hdr_amplicon_seq"/>
|
|
|
91 <param argument="--coding_seq" type="text" value="" label="Subsequence/s of the amplicon sequence covering one or more coding sequences for frameshift analysis" help="--coding_seq"/>
|
|
|
92 <section name="sgrna_parameters" expanded="false" title="sgRNA parameters">
|
|
|
93 <param argument="--guide_name" type="text" label="sgRNA names" help="if more than one, please separate by commas. (default: sgRNA)"/>
|
|
1
|
94 <param argument="--guide_seq" type="text" value="" label="sgRNA sequence" help="sgRNA sequence, if more than one, please separate by commas"/>
|
|
|
95 <param argument="--flexiguide_name" type="text" value="" optional="True" label="Names for the flexiguides" help="Similar to --guide_name"/>
|
|
0
|
96 <param argument="--flexiguide" type="text" value="" label="sgRNA sequence (flexible)" help="The flexiguide sequence will be aligned to the amplicon sequence(s), as long as the guide sequence has homology as set by --flexiguide_homology"/>
|
|
|
97 <param argument="--flexiguide_homology" type="integer" value="80" label="Flexiguide homology" help="will yield guides in amplicons with at least this homology to the flexiguide sequence (default:80 meaning 80% homology is required)"/>
|
|
|
98 <param name="discard_guide_positions_overhanging_amplicon_edge" truevalue="" falsevalue="--discard_guide_positions_overhanging_amplicon_edge" type="boolean" default="False" label="Discard guide positions overhanging amplicon edge" help="If set, for guides that align to multiple positions, guide positions will be discarded if plotting around those regions would included bp that extend beyond the end of the amplicon"/>
|
|
|
99 </section>
|
|
|
100 <section name="filtering_parameters" expanded="false" title="Read filtering, trimming, and merging parameters">
|
|
|
101 <param name="split_interleaved_input" type="boolean" default="False" truevalue="--split_interleaved_input" falsevalue="" label="Splits a single fastq file containing paired end reads in two files before running CRISPResso"/>
|
|
|
102 <param argument="--min_average_read_quality" value="0" type="integer" label="Minimum average quality score (phred33) to keep a read"/>
|
|
|
103 <param argument="--min_single_bp_quality" type="integer" value="0" label="Minimum single bp score (phred33) to keep a read"/>
|
|
|
104 <param argument="--min_bp_quality_or_N" type="integer" value="0" label="Bases with a quality score (phred33) less than this value will be set to 'N'"/>
|
|
|
105 <param argument="--min_paired_end_reads_overlap" type="integer" value="10" label="Minimum overlap length between reads" help="Parameter for the FLASH read merging step. Minimum required overlap length between two reads to provide a confident overlap. (default: 10)" />
|
|
|
106 <param argument="--max_paired_end_reads_overlap" type="integer" value="100" label="Maximum overlap length between reads" help="Parameter for the FLASH merging step. Maximum overlap length expected in approximately 90% of read pairs. Please see the FLASH manual for more information. (default: 100)"/>
|
|
|
107 <param name="stringent_flash_merging" truevalue="--stringent_flash_merging" falsevalue="" type="boolean" default="False" label="Use stringent parameters for flash merging" help="In the case where flash could merge R1 and R2 reads ambiguously, the expected overlap is calculated as 2*average_read_length - amplicon_length. The flash parameters for --min-overlap and --max-overlap will be set to prefer merged reads with length within 10bp of the expected overlap. These values override the --min_paired_end_reads_overlap or --max_paired_end_reads_overlap CRISPResso parameters"/>
|
|
|
108 </section>
|
|
2
|
109 <section name="window_parameters" expanded="false" title="Quantification window parameters">
|
|
|
110 <param argument="cleavage_offset" type="integer" optional="True" default="-3" label="Center of quantification window to use within respect to the 3' end of the provided sgRNA sequence." help="-wc or --quantification_window_center or --cleavage_offset (default = -3)"/>
|
|
|
111 </section>
|
|
|
112 <section name="output_parameters" expanded="false" title="Output parameters">
|
|
|
113 <param name="suppress_report" type="boolean" default="False" label="Suppress output report, plots output as .pdf only (not .png)" help="--suppress_report"/>
|
|
|
114 </section>
|
|
0
|
115 </inputs>
|
|
|
116 <outputs>
|
|
|
117 <data format="txt" name="output_log" label="${tool.name} on ${on_string}: log" from_work_dir="Log.final.out"/>
|
|
|
118 <data name="html_file" format="html" from_work_dir="crispresso.html" label="WebPage report"/>
|
|
1
|
119 <data format="txt" name="quant_editing_freq" label="CRISPResso_quantification_of_editing_frequency"/>
|
|
0
|
120 </outputs>
|
|
|
121 <help><![CDATA[
|
|
|
122 TODO: Fill in help.
|
|
|
123 ]]></help>
|
|
|
124 </tool>
|
