Bioinformatics/Sequence mutation: Difference between revisions
Content deleted Content added
Line 946:
i <- randomIndex (length xs)
case m of
Swap -> randomDNA >>= \d -> pure ((Swap, i),
Insert -> randomDNA >>= \d -> pure ((Insert, i),
Delete -> pure ((Delete, i),
where
randomIndex max = getStdRandom (randomR (1, max))
randomDNA = head . randoms <$> newStdGen
(m, ms) <- mutateSequence s
uncurry (printf "%6s @ %d\n") m
main :: IO ()
main = do
putStrLn "\nInitial Sequence:" >> showSequence
putStrLn "\nBase Counts:" >> showBaseCounts
printf "Total: %d\n\n" $ length
ms <-
putStrLn "\nMutated Sequence:" >> showSequence ms
putStrLn "\nBase Counts:" >> showBaseCounts ms
|