Bioinformatics/Sequence mutation: Difference between revisions

Swap -> randomDNA >>= \d -> pure ((Swap, i), swapIndexswapElement i d xs)

Insert -> randomDNA >>= \d -> pure ((Insert, i), insertIndexinsertElement i d xs)

Delete -> pure ((Delete, i), dropIndexdropElement i xs)

dropIndexdropElement i xs = take (pred i) xs <> drop i xs

insertIndexinsertElement i e xs = take i xs <> [e] <> drop i xs

swapIndexswapElement i a xs = take (pred i) xs <> [a] <> drop i xs

mutationsmutate :: Int -> DNASequence -> IO DNASequence

mutationsmutate 0 s = pure s

mutationsmutate n s = do

mutationsmutate (pred n) ms

dnaseqds <- newSequence 200

putStrLn "\nInitial Sequence:" >> showSequence dnaseqds

putStrLn "\nBase Counts:" >> showBaseCounts dnaseqds

printf "Total: %d\n\n" $ length dnaseqds

ms <- mutationsmutate 10 dnaseqds