def deviation(clusters, queryLength)
      hits = clusters.map { |c| c.lengths.map { |x| Array.new(x[1], x[0]) }.flatten }.flatten
      raw_hits = clusters.map { |c| c.lengths.map { |x| Array.new(x[1], x[0]) }.flatten }.flatten.to_s.delete('[').delete(']')
      R.eval("sd = sd(c(#{raw_hits}))")
      sd = R.pull('sd')
      sd = standard_deviation(hits)
      (queryLength - mean).abs / sd
    end
def deviation(clusters, queryLength)
      hits = clusters.map { |c| c.lengths.map { |x| Array.new(x[1], x[0]) }.flatten }.flatten
      raw_hits = clusters.map { |c| c.lengths.map { |x| Array.new(x[1], x[0]) }.flatten }.flatten.to_s.delete('[').delete(']')
      R.eval("sd = sd(c(#{raw_hits}))")
      sd = R.pull('sd')
      sd = standard_deviation(hits)
      (queryLength - mean).abs / sd
    end
Returns the deviation of a value from the values in all clusters
 Params:
 +clusters+: a list of Cluster objects
 +queryLength+: a reference Sequence object
 Output:
 Real number