hg.stevelosh.com > rosalind

src/problems/tran.lisp @ 32944636e2d3 

TRAN
author Steve Losh <steve@stevelosh.com>
date Mon, 24 Dec 2018 01:25:44 -0500
parents (none)
children 2735aa6aab79 
(in-package :rosalind)

(defparameter *input-tran* ">Rosalind_0209
GCAACGCACAACGAAAACCCTTAGGGACTGGATTATTTCGTGATCGTTGTAGTTATTGGA
AGTACGGGCATCAACCCAGTT
>Rosalind_2200
TTATCTGACAAAGAAAGCCGTCAACGGCTGGATAATTTCGCGATCGTGCTGGTTACTGGC
GGTACGAGTGTTCCTTTGGGT")

(defparameter *output-tran* "1.21428571429")

(defun transitionp (x y)
  (and (char/= x y)
       (= (rings x) (rings y))))

(defun transversionp (x y)
  (and (char/= x y)
       (/= (rings x) (rings y))))

(define-problem tran (data stream)
    *input-tran*
    *output-tran*
  (destructuring-bind (x y)
      (mapcar #'cdr (read-fasta-into-alist data))
    (format nil "~,11F" (coerce (/ (mapcount #'transitionp x y)
                                   (mapcount #'transversionp x y))
                                'double-float))))

;; (problem-tran "2")
;; (solve tran)