"read_sff -m -i #{@sff_file} |
progress_meter |
plot_distribution -k SEQ_LEN -T 'Length Distribution - unclipped' -t png -o #{@plot1} |
- plot_scores -T 'Mean Quality Scores - unclipped' -t png -o #{@plot2} |
+ plot_scores -c -T 'Mean Quality Scores - unclipped' -t png -o #{@plot2} |
clip_seq |
plot_distribution -k SEQ_LEN -T 'Length Distribution - clipped' -t png -o #{@plot3} |
- plot_scores -T 'Mean Quality Scores - clipped' -t png -o #{@plot4} |
+ plot_scores -c -T 'Mean Quality Scores - clipped' -t png -o #{@plot4} |
mean_scores |
bin_vals -k SCORES_MEAN -b 5 |
plot_histogram -s num -k SCORES_MEAN_BIN -T 'Mean score bins' -X 'Bins (size 5)' -Y 'Count' -t png -o #{@plot5} -x"