# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
require 'maasha/biopieces'
+require 'maasha/seq'
require 'gnuplot'
require 'narray'
-require 'pp'
terminals = "dumb,x11,aqua,post,pdf,png,svg"
title = "Mean Quality Scores"
options = Biopieces.options_parse(ARGV, casts)
-SCORES_MAX = 10_000
+SCORES_MAX = 100_000
scores_vec = NArray.int(SCORES_MAX)
count_vec = NArray.int(SCORES_MAX)
plot.title options[:title]
plot.xlabel options[:xlabel]
plot.ylabel options[:ylabel]
- plot.output options[:data_out] if options[:data_out]
+ plot.output options[:data_out] || "/dev/stderr"
plot.xrange "[#{x.min - 1}:#{x.max + 1}]"
plot.yrange "[#{Seq::SCORE_MIN}:#{Seq::SCORE_MAX}]"
plot.style "fill solid 0.5 border"