File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -125,57 +125,6 @@ process qc_se {
125125}
126126
127127
128- // Quality filtering for pair-end reads
129- process qc_pe {
130-
131- label " main_container"
132-
133- // cpus 10
134-
135- input:
136- path input_R1
137- path input_R2
138-
139- output:
140- path " QC_R1.fq.gz" , emit: filtered_R1
141- path " QC_R2.fq.gz" , emit: filtered_R2
142-
143- script:
144- filter_avgphred = params. qc_avgphred ? " --average_qual ${ params.qc_avgphred} " : " --average_qual 0"
145- filter_phredmin = params. qc_phredmin ? " --qualified_quality_phred ${ params.qc_phredmin} " : " "
146- filter_phredperc = params. qc_phredperc ? " --unqualified_percent_limit ${ params.qc_phredperc} " : " "
147- filter_polyglen = params. qc_polyglen ? " --trim_poly_g --poly_g_min_len ${ params.qc_polyglen} " : " "
148- """
149- echo -e "QC\\ n"
150- echo -e "Input R1: " ${ input_R1}
151- echo -e "Input R2: " ${ input_R2}
152-
153- ## If `filter_phredmin` && `filter_phredperc` are specified,
154- # Filtering based on percentage of unqualified bases
155- # how many percents of bases are allowed to be unqualified (Q < 24)
156-
157- fastp \
158- --in1 ${ input_R1} \
159- --in2 ${ input_R2} \
160- --disable_adapter_trimming \
161- --n_base_limit ${ params.qc_maxn} \
162- ${ filter_avgphred} \
163- ${ filter_phredmin} \
164- ${ filter_phredperc} \
165- ${ filter_polyglen} \
166- --length_required 100 \
167- --thread ${ task.cpus} \
168- --html qc.html \
169- --json qc.json \
170- --out1 QC_R1.fq.gz \
171- --out2 QC_R2.fq.gz
172-
173- echo -e "\\ nQC finished"
174- """
175- }
176-
177-
178-
179128// Validate tags for demultiplexing
180129process tag_validation {
181130
You can’t perform that action at this time.
0 commit comments