|
0
|
1 <tool id="fastqFilter" name="Filter FASTQ files">
|
|
|
2 <description></description>
|
|
|
3 <command>
|
|
|
4 fastq_filter.pl -v
|
|
|
5 #if $sampleIndex.filterBySampleIndex == "yes":
|
|
|
6 -index $sampleIndex.sequence
|
|
|
7 #end if
|
|
|
8 -maxN $maxN -if sanger -f $filterString -of $outputFormat $inputfile $outputfile
|
|
|
9 </command>
|
|
|
10
|
|
|
11 <inputs>
|
|
|
12 <param name="inputfile" format="fastq" type="data" label="Input Sanger FASTQ file (.gz file accepted; see help below for more information)" />
|
|
|
13
|
|
|
14 <conditional name="sampleIndex">
|
|
|
15 <param name="filterBySampleIndex" type="select" label="Filter by sample index (see help below for parameter suggestion)" >
|
|
|
16 <option value="yes">Yes</option>
|
|
|
17 <option value="no" selected="true">No</option>
|
|
|
18 </param>
|
|
|
19 <when value="yes">
|
|
|
20 <param name="sequence" type="text" value="" label="Index position and sequence" />
|
|
|
21 </when>
|
|
|
22 <when value="no">
|
|
|
23 </when>
|
|
|
24 </conditional>
|
|
|
25
|
|
|
26 <param name="filterString" type="text" value="" label="Quality score filter string; format: Method:Start-End:Score (zero-based; see help below for parameter suggestion)" />
|
|
|
27 <param name="maxN" type="integer" value="-1" label="Max number of N in sequence (default off - value less than 0) " />
|
|
|
28 <param name="outputFormat" type="select" label="Output data type">
|
|
|
29 <option value="fastq">FASTQ</option>
|
|
|
30 <option value="fasta">FASTA</option>
|
|
|
31 </param>
|
|
|
32
|
|
|
33 </inputs>
|
|
|
34
|
|
|
35 <outputs>
|
|
|
36 <data name="outputfile" format="fastq" label="Read quality filtering on ${on_string}">
|
|
|
37 <change_format>
|
|
|
38 <when input="OutputFormat" value="fasta" format="fasta" />
|
|
|
39 </change_format>
|
|
|
40 </data>
|
|
|
41 </outputs>
|
|
|
42
|
|
|
43 <help>
|
|
|
44
|
|
|
45 .. class:: infomark
|
|
|
46
|
|
|
47 **What this tool does**
|
|
|
48
|
|
|
49 This tool extracts reads passing quality filters.
|
|
|
50
|
|
|
51 It takes as input Sanger FASTQ files and output FASTQ/A files of filtered reads.
|
|
|
52
|
|
|
53 -----
|
|
|
54
|
|
|
55 **FASTQ format**
|
|
|
56
|
|
|
57 Check quality score in the FASTQ file for the right format.
|
|
|
58
|
|
|
59 Reference https://en.wikipedia.org/wiki/FASTQ_format#Quality :
|
|
|
60
|
|
|
61 * Sanger format can encode a Phred quality score from 0 to 93 using ASCII 33 to 126.
|
|
|
62 * Solexa/Illumina 1.0 format can encode a Solexa/Illumina quality score from -5 to 62 using ASCII 59 to 126.
|
|
|
63
|
|
|
64 See http://www.asciitable.com/ for ASCII table.
|
|
|
65
|
|
|
66 -----
|
|
|
67
|
|
|
68 **Filter by sample index (optional)**
|
|
|
69
|
|
|
70 For users who would like to start from a FASTQ file consisting of multiple libraries.
|
|
|
71
|
|
|
72 For example:
|
|
|
73
|
|
|
74 If you have six samples with indexes GTCA, GCAT, ACTG, AGCT, GCAT, TCGA, you can extract reads for each library with indicated index sequences (e.g. GTCA, etc.) starting from position 0 in the read. For example, you could specify 0:GTCA, etc.
|
|
|
75
|
|
|
76 -----
|
|
|
77
|
|
|
78 **How to set the filter**
|
|
|
79
|
|
|
80 You can apply multiple filtering criteria based on the quality scores for each read. They are separated by commas.
|
|
|
81
|
|
|
82 Each critieron is composed of four components (e.g. method1:start1-end1:score1,method2:start2-end2:score2)
|
|
|
83
|
|
|
84 1. Method: min or mean, which means requirement on minimal or mean score of a region
|
|
|
85 2. Start: the first nucleotide to consider (0-based)
|
|
|
86 3. End: the last nucleotide to consider (0-based)
|
|
|
87 4. score: the threshold required
|
|
|
88
|
|
|
89 **Parameter suggestion**
|
|
|
90
|
|
|
91 For example:
|
|
|
92
|
|
|
93 * For Standard CLIP protocol filtering: mean:0-29:20 (this specifies a mean score of 20 or above in the first 30 bases, which includes 5 positions with sample indexes and the random barcode, followed by 25 positions with the actual CLIP tag).
|
|
|
94 * For iCLIP/BrdU CLIP filtering: mean:0-38:20 (this specifies a mean score of 20 or above in the first 39 bases, which includes 14 positions with sample indexes and the random barcode, followed by 25 positions with the actual CLIP tag).
|
|
|
95
|
|
|
96 The reason to filter as such is because low quality reads can introduce mapping errors and background. They will inflate the number of unique tags after removal of PCR duplicates.
|
|
|
97
|
|
|
98
|
|
|
99
|
|
|
100 </help>
|
|
|
101 </tool>
|
