fixed QA_Illumina_report.rb SEQ_LEN > 20

[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 f5b426f0a5e0987b6e73f860f176de2bdfef4212..74d720a2c0a2e0334c8a94e799ec31d28e65af1e 100755 (executable)
--- a/bp_scripts/QA_Illumina_report.rb
+++ b/bp_scripts/QA_Illumina_report.rb
@@ -70,9 +70,9 @@ system(
    progress_meter |
    analyze_vals -k SEQ -o #{analyze_vals_file} |
    trim_seq -l 3 -m 25 |
-   grab -e 'SEQ_LEN > 0' |
+   grab -e 'SEQ_LEN > 20' |
    analyze_vals -k SEQ -o #{analyze_vals_trim_file} |
-   find_adaptor -c 12 -l 6 -L 6 -f ACACGACGCTCTTCCGATCT -r AGATCGGAAGAGCACACGTC |
+   find_adaptor -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} |