Safe Haskell | None |
---|---|
Language | Haskell2010 |
Common data types used everywhere. This module is a collection of very basic "bioinformatics" data types that are simple, but don't make sense to define over and over.
- newtype Nucleotide = N {}
- newtype Nucleotides = Ns {}
- newtype Qual = Q {}
- toQual :: (Floating a, RealFrac a) => a -> Qual
- fromQual :: Qual -> Double
- fromQualRaised :: Double -> Qual -> Double
- probToQual :: (Floating a, RealFrac a) => Prob' a -> Qual
- newtype Prob' a = Pr {
- unPr :: a
- type Prob = Prob' Double
- toProb :: Floating a => a -> Prob' a
- fromProb :: Floating a => Prob' a -> a
- qualToProb :: Floating a => Qual -> Prob' a
- pow :: Num a => Prob' a -> a -> Prob' a
- data Word8 :: *
- nucA :: Nucleotide
- nucC :: Nucleotide
- nucG :: Nucleotide
- nucT :: Nucleotide
- nucsA :: Nucleotides
- nucsC :: Nucleotides
- nucsG :: Nucleotides
- nucsT :: Nucleotides
- nucsN :: Nucleotides
- gap :: Nucleotides
- toNucleotide :: Char -> Nucleotide
- toNucleotides :: Char -> Nucleotides
- nucToNucs :: Nucleotide -> Nucleotides
- showNucleotide :: Nucleotide -> Char
- showNucleotides :: Nucleotides -> Char
- isGap :: Nucleotides -> Bool
- isBase :: Nucleotides -> Bool
- isProperBase :: Nucleotides -> Bool
- properBases :: [Nucleotides]
- compl :: Nucleotide -> Nucleotide
- compls :: Nucleotides -> Nucleotides
- type Seqid = ByteString
- data Position = Pos {}
- shiftPosition :: Int -> Position -> Position
- p_is_reverse :: Position -> Bool
- data Range = Range {}
- shiftRange :: Int -> Range -> Range
- reverseRange :: Range -> Range
- extendRange :: Int -> Range -> Range
- insideRange :: Range -> Range -> Range
- wrapRange :: Int -> Range -> Range
- w2c :: Word8 -> Char
- c2w :: Char -> Word8
- findAuxFile :: FilePath -> IO FilePath
Documentation
newtype Nucleotide Source #
A nucleotide base. We only represent A,C,G,T. The contained
Word8
ist guaranteed to be 0..3.
newtype Nucleotides Source #
A nucleotide base in an alignment. Experience says we're dealing with Ns and gaps all the type, so purity be damned, they are included as if they were real bases.
To allow Nucleotides
s to be unpacked and incorporated into
containers, we choose to represent them the same way as the BAM file
format: as a 4 bit wide field. Gaps are encoded as 0 where they
make sense, N is 15. The contained Word8
is guaranteed to be
0..15.
Qualities are stored in deciban, also known as the Phred scale. To
represent a value p
, we store -10 * log_10 p
. Operations work
directly on the "Phred" value, as the name suggests. The same goes
for the Ord
instance: greater quality means higher "Phred"
score, meand lower error probability.
A positive floating point value stored in log domain. We store the
natural logarithm (makes computation easier), but allow conversions
to the familiar "Phred" scale used for Qual
values.
Unbox a => Vector Vector (Prob' a) Source # | |
Unbox a => MVector MVector (Prob' a) Source # | |
Eq a => Eq (Prob' a) Source # | |
(Floating a, Fractional a, Ord a) => Fractional (Prob' a) Source # | |
(Floating a, Ord a) => Num (Prob' a) Source # | |
Ord a => Ord (Prob' a) Source # | |
RealFloat a => Show (Prob' a) Source # | |
Storable a => Storable (Prob' a) Source # | |
Unbox a => Unbox (Prob' a) Source # | |
data MVector s (Prob' a) Source # | |
data Vector (Prob' a) Source # | |
8-bit unsigned integer type
Bounded Word8 | Since: 2.1 |
Enum Word8 | Since: 2.1 |
Eq Word8 | Since: 2.1 |
Integral Word8 | Since: 2.1 |
Data Word8 | Since: 4.0.0.0 |
Num Word8 | Since: 2.1 |
Ord Word8 | Since: 2.1 |
Read Word8 | Since: 2.1 |
Real Word8 | Since: 2.1 |
Show Word8 | Since: 2.1 |
Ix Word8 | Since: 2.1 |
Lift Word8 | |
Storable Word8 | Since: 2.1 |
Bits Word8 | Since: 2.1 |
FiniteBits Word8 | Since: 4.6.0.0 |
Hashable Word8 | |
Prim Word8 | |
Unbox Word8 | |
Vector Vector Word8 | |
MVector MVector Word8 | |
data Vector Word8 | |
data MVector s Word8 | |
nucA :: Nucleotide Source #
nucC :: Nucleotide Source #
nucG :: Nucleotide Source #
nucT :: Nucleotide Source #
nucsA :: Nucleotides Source #
nucsC :: Nucleotides Source #
nucsG :: Nucleotides Source #
nucsT :: Nucleotides Source #
nucsN :: Nucleotides Source #
gap :: Nucleotides Source #
toNucleotide :: Char -> Nucleotide Source #
Converts a character into a Nucleotides
.
The usual codes for A,C,G,T and U are understood, -
and .
become
gaps and everything else is an N.
toNucleotides :: Char -> Nucleotides Source #
Converts a character into a Nucleotides
.
The usual codes for A,C,G,T and U are understood, -
and .
become
gaps and everything else is an N.
nucToNucs :: Nucleotide -> Nucleotides Source #
showNucleotide :: Nucleotide -> Char Source #
showNucleotides :: Nucleotides -> Char Source #
isGap :: Nucleotides -> Bool Source #
Tests if a Nucleotides
is a gap.
Returns true only for the gap.
isBase :: Nucleotides -> Bool Source #
Tests if a Nucleotides
is a base.
Returns True
for everything but gaps.
isProperBase :: Nucleotides -> Bool Source #
Tests if a Nucleotides
is a proper base.
Returns True
for A,C,G,T only.
properBases :: [Nucleotides] Source #
compl :: Nucleotide -> Nucleotide Source #
Complements a Nucleotides.
compls :: Nucleotides -> Nucleotides Source #
Complements a Nucleotides.
type Seqid = ByteString Source #
Sequence identifiers are ASCII strings Since we tend to store them for a while, we use strict byte strings.
Coordinates in a genome. The position is zero-based, no questions about it. Think of the position as pointing to the crack between two bases: looking forward you see the next base to the right, looking in the reverse direction you see the complement of the first base to the left.
To encode the strand, we (virtually) reverse-complement any sequence and prepend it to the normal one. That way, reversed coordinates have a negative sign and automatically make sense. Position 0 could either be the beginning of the sequence or the end on the reverse strand... that ambiguity shouldn't really matter.
shiftPosition :: Int -> Position -> Position Source #
Moves a Position
. The position is moved forward according to the
strand, negative indexes move backward accordingly.
p_is_reverse :: Position -> Bool Source #
Ranges in genomes
We combine a position with a length. In 'Range pos len', pos
is
always the start of a stretch of length len
. Positions therefore
move in the opposite direction on the reverse strand. To get the
same stretch on the reverse strand, shift r_pos by r_length, then
reverse direction (or call reverseRange).
reverseRange :: Range -> Range Source #
Reverses a Range
to give the same Range
on the opposite strand.
extendRange :: Int -> Range -> Range Source #
Extends a range. The length of the range is simply increased.
insideRange :: Range -> Range -> Range Source #
Expands a subrange.
(range1
interprets insideRange
range2)range1
as a subrange of
range2
and computes its absolute coordinates. The sequence name of
range1
is ignored.
wrapRange :: Int -> Range -> Range Source #
Wraps a range to a region. This simply normalizes the start
position to be in the interval '[0,n)', which only makes sense if the
Range
is to be mapped onto a circular genome. This works on both
strands and the strand information is retained.