-- Calculate various characteristics on sequence quality

module Metrics where

import Bio.Sequence.SFF
import Bio.Sequence.SeqData

-- import Test.QuickCheck

-- | Take the fractional parts of the flows, and sum their squares (the "K²" metric)
quals :: [Flow] -> Flow
quals q = floor $ ((100::Double) - 2*(sqrt $ (/fromIntegral (length q)) $ sum $ map (fromIntegral . (^(2::Integer)) . (flip (-) 50) . (`mod` 100) . (+50)) $ q))

-- | Count number of n's in the sequence
--   The algorithm for generating Ns is a bit opaque, and appears to depend on the magnitude 
--   of the noise flow values.  We chicken out, and just count the called sequence.
n_count :: ReadBlock -> Int
n_count r = length . filter isN . clip . toStr . bases $ r
    where isN x = x=='N' || x == 'n'
          clip = take (right-left+1) . drop left
          right = fromIntegral $ clip_qual_right (read_header r)
          left = fromIntegral $ clip_qual_left (read_header r)