module Biobase.Primary.AA where
import Control.Arrow ((***),first)
import Data.Hashable
import Data.Ix (Ix(..))
import Data.Map.Strict (Map)
import Data.Primitive.Types
import Data.Tuple (swap)
import Data.Vector.Unboxed.Deriving
import GHC.Base (remInt,quotInt)
import GHC.Generics (Generic)
import qualified Data.ByteString.Char8 as BS
import qualified Data.ByteString.Lazy.Char8 as BSL
import qualified Data.Map.Strict as M
import qualified Data.Text as T
import qualified Data.Vector.Generic as VG
import qualified Data.Vector.Generic.Mutable as VGM
import qualified Data.Vector.Unboxed as VU
import qualified Data.Bijection.Map as B
import Biobase.Primary.Letter
data AA
pattern Stop = Letter 0 :: Letter AA
pattern A = Letter 1 :: Letter AA
pattern B = Letter 2 :: Letter AA
pattern C = Letter 3 :: Letter AA
pattern D = Letter 4 :: Letter AA
pattern E = Letter 5 :: Letter AA
pattern F = Letter 6 :: Letter AA
pattern G = Letter 7 :: Letter AA
pattern H = Letter 8 :: Letter AA
pattern I = Letter 9 :: Letter AA
pattern K = Letter 10 :: Letter AA
pattern L = Letter 11 :: Letter AA
pattern M = Letter 12 :: Letter AA
pattern N = Letter 13 :: Letter AA
pattern P = Letter 14 :: Letter AA
pattern Q = Letter 15 :: Letter AA
pattern R = Letter 16 :: Letter AA
pattern S = Letter 17 :: Letter AA
pattern T = Letter 18 :: Letter AA
pattern V = Letter 19 :: Letter AA
pattern W = Letter 20 :: Letter AA
pattern X = Letter 21 :: Letter AA
pattern Y = Letter 22 :: Letter AA
pattern Z = Letter 23 :: Letter AA
pattern Undef = Letter 24 :: Letter AA
aa :: Int -> Letter AA
aa = Letter
aaRange = [Stop .. pred Undef]
charAA :: Char -> Letter AA
charAA = B.findWithDefaultL Undef charBaa
aaChar :: Letter AA -> Char
aaChar = B.findWithDefaultR '?' charBaa
charBaa :: B.Bimap Char (Letter AA)
charBaa = B.fromList
[ ('*',Stop)
, ('A',A)
, ('B',B)
, ('C',C)
, ('D',D)
, ('E',E)
, ('F',F)
, ('G',G)
, ('H',H)
, ('I',I)
, ('K',K)
, ('L',L)
, ('M',M)
, ('N',N)
, ('P',P)
, ('Q',Q)
, ('R',R)
, ('S',S)
, ('T',T)
, ('V',V)
, ('W',W)
, ('X',X)
, ('Y',Y)
, ('Z',Z)
, ('?',Undef)
]
instance Show (Letter AA) where
show n = [aaChar n]
instance Read (Letter AA) where
readsPrec p [] = []
readsPrec p (x:xs)
| x==' ' = readsPrec p xs
| aa <- charAA x = [(aa,xs)]
| otherwise = []
instance Enum (Letter AA) where
succ Undef = error "succ/Undef:AA"
succ (Letter x) = Letter $ x+1
pred Stop = error "pred/Stop:AA"
pred (Letter x) = Letter $ x1
toEnum k | k>=0 && k<=(getLetter Undef) = Letter k
toEnum k = error $ "toEnum/Letter RNA " ++ show k
fromEnum (Letter k) = k
instance MkPrimary (VU.Vector Char) AA where
primary = VU.map charAA