more power to the QA

[biopieces.git] / bp_scripts / QA_Illumina_report.rb

diff --git a/bp_scripts/QA_Illumina_report.rb b/bp_scripts/QA_Illumina_report.rb
index 2655b10a4856076fa41a88254db65e8e7e811459..3139d49127b1e3a796c4e30cc819f3314df94f14 100755 (executable)
--- a/bp_scripts/QA_Illumina_report.rb
+++ b/bp_scripts/QA_Illumina_report.rb
@@ -72,7 +72,7 @@ system(
    trim_seq -l 3 -m 25 |
    grab -e 'SEQ_LEN > 0' |
    analyze_vals -k SEQ -o #{analyze_vals_trim_file} |
-   find_adaptor -l 6 -L 6 -f ACACGACGCTCTTCCGATCT -r AGATCGGAAGAGCACACGTC |
+   find_adaptor -c 12 -l 6 -L 6 -f ACACGACGCTCTTCCGATCT -r AGATCGGAAGAGCACACGTC |
    clip_adaptor |
    grab -e 'SEQ_LEN > 0' |
    analyze_vals -k SEQ -o #{analyze_vals_trim_noadapt_file} |
@@ -112,7 +112,7 @@ template = %{
       </table>
       <p>Sequence trimming was performed by removing from the ends all residues until 3 consecutive</p>
       <p>residues with quality score larger than or equal to 25.</p>
-      <p>All plots are after sequence trimming.</p>
+      <p>All plots are after sequence trimming and adaptor removal.</p>
       <h2>Sequence length distribution</h2>
       <p><img src="<%= png2base64(lendist_file) %>" width="600" /></p>
       <p><img src="<%= png2base64(lendist_bin_file) %>" width="600" /></p>