1 Introduction

In this work we focus on multi-state systems modeled by means of a particular class of non-homogeneous Markov processes introduced in Vergne (2008), called drifting Markov processes. Most of the estimation methods and reliability results are developed in Barbu and Vergne (2019).

Note that in many mathematical models it is assumed the homogeneity with respect to time, which is inappropriate in most of the applications. But, considering general non-homogeneous processes could be unrealistic from a practical point of view. For this reason, the drifting Markov chains assume that the Markov transition matrix is a linear/polynomial function of two/several Markov transition matrices. Thus we obtain a “smooth” non-homogeneity, with sensibly less parameters than in the general case.

Few R packages have been developed to handle similar types of models of Markov, semi-Markov or hidden semi-Markov type, useful in reliability or DNA analysis. For semi-Markov models we have semiMarkov R package Król and Saint-Pierre (2015) that performs maximum likelihood estimation for parametric continuous-time semi-Markov processes, smm R package Barbu, Bérard, et al. (2018) which performs parametric and non-parametric estimation and simulation for multi-state discrete-time semi-Markov processes. Two R packages are also dedicated to hidden semi-Markov models, implementing estimation and prediction methods: the hsmm R package Bulla et al. (2010) and the mhsmm R package O’Connell and Højsgaard (2011).

Note that there is no R package developed for drifting Markov models (DMMs). Thus the purpose of this paper is to present an R package that we have developed, called drimmR, which performs estimation and simulation for such models, as well as the estimation of associated reliability measures. The aim of this paper is to describe the different possibilities of this package. To summarize, the package drimmR that we present deals with different problems:

We consider one or several sample paths; for several sample paths, two estimation methods are proposed: one is the usual LSE, the other one is the average of LSEs obtained on each sample;
The samples paths are complete or incomplete;
The sample paths come from drifting Markov chains that are of the same length or of different lengths (between the same Markov transition probability matrices);
We derive exact computations of reliability/survival analysis measures (reliability or survival function, availability, maintainability, failure rates);
We compute the probabilities of appearance of a word in a given sequence under a given model.

We would like to mention that a web interface called WebDRIMM has been also developed [cf. Barbu, Brelurut, et al. (2018)] for simulating and estimating drifting Markov models, as well as associated reliability indicators; it is available at http://bioinfo.univ-rouen.fr/WebDRIMM
The paper is organized as follows. Section 2 describes the drifting Markov models used in this package, present associated reliability indicators and corresponding estimation results and techniques. Section 3 illustrates the different functions of the drimmR package and we end the paper by presenting some concluding remarks on this R package in Section 4.

2 Drifting Markov models, estimation and associated reliability measures

Let us consider a random system with finite state space $E=\{1,\ldots, s\}$, $s < \infty.$ We assume that the time evolution of a system is governed by a discrete-time stochastic process with values in $E.$ In the following we will introduce a class of processes according to Vergne (2008) in Section 2.1, we will briefly present associated reliability indicators according to Barbu and Vergne (2019) in Section 2.2 and estimation of the parameters in Section 2.3, according to Vergne (2008) Barbu and Vergne (2019).

2.1 Drifting Markov models

Let $\boldsymbol{\Pi}_0 = (\Pi_0(u, v))_{u,v \in E}$ and $\boldsymbol{\Pi}_1 = (\Pi_1(u, v))_{u,v \in E}$ be two Markov transition matrices of order 1 over the state space $E.$

Definition 1 (linear drifting Markov chain of order $1$ and of length $n$). A sequence $X_0,$ $X_1, \ldots, X_n$ with state space $E=\mathbb \{1, 2, \ldots, s\}$ is said to be a linear drifting Markov chain (of order 1)* of length $n$ between the Markov transition matrices $\boldsymbol{\Pi}_0$ and $\boldsymbol{\Pi}_1$ if the distribution of $X_t,$ $t = 1, \ldots, n,$ is defined by \[\label{eq_DriftLin1} \mathbb P(X_t=v \mid X_{t-1} = u, X_{t-2}, \ldots ) = \Pi_{\frac{t}{n}}(u, v), \; u, v \in E,\] where \[\label{eq_DriftLin2} \Pi_{\frac{t}{n}}(u, v) = \left( 1 - \frac{t}{n} \right) \Pi_0(u, v) + \frac{t}{n} \Pi_1(u, v), \; u, v \in E.\]*

Let us denote by $\boldsymbol{\alpha} = (\alpha(1), \ldots, \alpha(s))$ the initial distribution of the chain, that is the distribution of $X_0,$ $\alpha(u)=\mathbb P (X_0 = u)$ for any state $u\in E.$

The linear drifting Markov model of order $1$ can be generalized to polynomial drifting Markov model of order $k$ and degree $d$. Let $\boldsymbol{\Pi}_{\frac{i}{d}} = (\boldsymbol{\Pi}_{\frac{i}{d}}(u_1, \dots, u_k, v))_{u_1, \dots, u_k,v \in E}$ be $d$ Markov transition matrices (of order $k$) over a state space $E.$

Definition 2 (polynomial drifting Markov chain of order $k$ and of length $n$). A sequence $X_0,$ $X_1,$ $\ldots,$ $X_n$ with state space $E=\mathbb \{1, 2, \ldots, s\}$ is said to be a polynomial drifting Markov chain of order k* and of length $n$ if the distribution of $X_t,$ $t = 1, \ldots, n,$ is defined by \[\label{eq_DriftPol1} \mathbb P(X_t=v \mid X_{t-1} = u_k, X_{t-2}=u_{k-1}, \ldots ) = \boldsymbol{\Pi}_{\frac{t}{n}}(u_1, \dots, u_k, v), \; u_1, \dots, u_k, v \in E,\] where \[\label{eq_DriftPol2} \boldsymbol{\Pi}_{\frac{t}{n}}(u_1, \dots, u_k, v) = \displaystyle\sum_{i=0}^d A_i(t)\boldsymbol{\Pi}_{\frac{i}{d}}(u_1, \dots, u_k,v), \; u_1, \dots, u_k, v \in E,\] with $A_i$ polynomials of degree $d$ such as, for any $i, j \in \{0, 1, \ldots, d\},$ $A_i\left(\frac{nj}{d}\right)=\mathbb{1}_{\{i=j\}}.$*

We would like to stress that the coherence between notations implies the choice of the notation $\boldsymbol{\Pi}_{\frac{i}{d}}$ and that, in fact, $A_i$ are Lagrange polynomials; see Vergne (2008) for more details on these two points.

2.2 Reliability of drifting Markov models

In order to undertake a reliability analysis of a system modeled by a DMM, let us assume that the state space of the system is partitioned into working and failure states, $E=U \cup D,$ with $U \cap D=\emptyset,$ where $U=\{1,\ldots,s_1\}$ represents the working states and $D=\{s_1+1,\ldots,s\}$ the failure states of the system. According to this partition of the state space we partition any matrix of vector we are working with and we denote the corresponding partitions accordingly (e.g., $\boldsymbol{\Pi}_0^{UU},$ $\boldsymbol{\Pi}_0^{DU},$ $\boldsymbol{\alpha}^{U}$ etc.).

For a linear drifting Markov chain of order $1$ $(X_t)_{0 \leq t \leq n},$ the reliability at time $l,$ $l \in \mathbb N,$ is given by \[\label{eq_RelDrift} R(l) = \boldsymbol{\alpha}^{U} \; \prod_{t=1}^l \left( \left(1 - \frac{t}{n}\right) \boldsymbol{\Pi}_0^{UU} + \frac{t}{n} \boldsymbol{\Pi}_1^{UU} \right) \; \mathbb{1}^U, \textrm{ where } \mathbb{1}^U = (\underbrace{1, \cdots, 1}_{s_1})^{\top}.\]

For a linear drifting Markov chain of order $1$ $(X_t)_{0 \leq t \leq n},$ the pointwise (or instantaneous) availability at time $l,$ $l \in \mathbb N,$ is given by \[\label{eq_AvDrift} A(l) = \boldsymbol{\alpha} \; \prod_{t=1}^l \left( \left(1 - \frac{t}{n}\right) \boldsymbol{\Pi}_0 + \frac{t}{n} \boldsymbol{\Pi}_1 \right) \; \mathbb{1}^{E,U}, \textrm{ where } \mathbb{1}^{E,U} = (\underbrace{1, \cdots, 1}_{s_1}, \underbrace{0, \cdots, 0}_{s-s_1} )^{\top}.\]

For a linear drifting Markov chain of order $1$ $(X_t)_{0 \leq t \leq n},$ the maintainability at time $l,$ $l \in \mathbb N,$ is given by \[\label{eq_MaintDrift_gen} M(l) = 1 - \boldsymbol{\alpha}^{D} \; \prod_{t=1}^l \left( \left(1 - \frac{t}{n}\right) \boldsymbol{\Pi}_0^{DD} + \frac{t}{n} \boldsymbol{\Pi}_1^{DD} \right) \; \mathbb{1}^D, \textrm{ where } \mathbb{1}^D = (\underbrace{1, \cdots, 1}_{s-s_1})^{\top}.\]

For a linear drifting Markov chain of order $1$ $(X_t)_{0 \leq t \leq n},$ the BMP-failure rate at time $l at time$l,$$l N,$is given by % \[\begin{equation} \label{eq_BMPDrift} \lambda(l) \left\{ \begin{array}{ll} 1- \frac{\mu_0^U \ \prod_{t=1}^{l}( \ (1-\frac{t}{n}) \pi_0^{UU} + (\frac{t}{n}) \pi_1^{UU}) \ \mathbb{1}^U}{\mu_0^U \ \prod_{t=1}^{l-1}( \ (1-\frac{t}{n}) \pi_0^{UU} + (\frac{t}{n}) \pi_1^{UU}) \ \mathbb{1}^U} \ \ , \ \text{si R(l-1) != 0 }\\ 0 \ , \ otherwise \\ \end{array} \right. \end{equation}\]

\end{proposition}

\begin{proposition} For a linear drifting Markov chain of order$1$$(X_t)\_0t n,$ the RG-failure rate at time$l at time $l,$ $l \in \mathbb N,$ is given by \[\label{eq_RGDrift} r(l) \left\{ \begin{array}{ll} -\ln \frac{\mu_0^U \ \prod_{t=1}^{l}( \ (1-\frac{t}{n}) \pi_0^{UU} + (\frac{t}{n}) \pi_1^{UU}) \ \mathbb{1}^U}{\mu_0^U \ \prod_{t=1}^{l-1}( \ (1-\frac{t}{n}) \pi_0^{UU} + (\frac{t}{n}) \pi_1^{UU}) \ \mathbb{1}^U} \ \ , \ if \ l \ge 1 \ ,\\ -\ln R(0) \ , \ if \ \ l = 0 \\ \end{array} \right.\]

See Barbu and Vergne (2019) for more details and Barbu and Limnios (2004); Barbu and Limnios (2008) for similar questions in a discrete-time semi-Markov framework.

2.3 Estimation of drifting Markov models

In this section we will consider different types of data for which the estimators of the characteristics of a drifting Markov chain and of the associated reliability indicators will be derived.

One can observe one sample path, that will be denoted by ${\mathcal H}(m,n):= (X_0,X_1, \ldots,X_{m}),$ where $m$ denotes the length of the sample path and $n$ the length of the drifting Markov chain. Two cases can be considered: (a1) $m=n$ (a complete sample path); (a2) $m < n$ (an incomplete sample path).

One can also observe $H$ i.i.d. sample paths, ${\mathcal H}_i(m_i,n_i), i=1, \ldots, H.$ Four cases can be considered here: (b1) $m_i=n_i=n$ for all $i=1, \ldots, H$ (complete sample paths of drifting Markov chains of the same length) ; (b2) $n_i=n$ for all $i=1, \ldots, H$ (incomplete sample paths of drifting Markov chains of the same length); (b3) $m_i=n_i$ for all $i=1, \ldots, H$ (complete sample paths of drifting Markov chains of different lengths) ; (b4) $m_i \leq n_i$ for all $i=1, \ldots, H$ (incomplete sample paths of drifting Markov chains of different lengths).

We have developed cf. Vergne (2008); Barbu and Vergne (2019) mean square estimators starting from data under these frameworks; we present here only the case of a linear drifting Markov chain of order $1$ under the sample framework (b1) (complete sample paths of drifting Markov chains of the same length).

Under the setting (b1), starting with $H$ complete sample paths of drifting Markov chains of the same length of a linear drifting Markov chain between two Markov transition matrices (of order 1) $\boldsymbol{\Pi}_0$ and $\boldsymbol{\Pi}_1,$ for any states $u, v \in E,$ the estimators of $\Pi_0(u,v)$ and $\Pi_1(u,v)$ are given by:

\[\begin{aligned} &&\widehat{\Pi}_{0; (n,H)} (u,v) = \frac{P_1(H, \boldsymbol{m}, \boldsymbol{n}) P_2(H, \boldsymbol{m}, \boldsymbol{n}) - P_3(H, \boldsymbol{m}, \boldsymbol{n})P_4(H, \boldsymbol{m}, \boldsymbol{n})} {P_5(H, \boldsymbol{m}, \boldsymbol{n}) P_1(H, \boldsymbol{m}, \boldsymbol{n}) - P_3(H, \boldsymbol{m}, \boldsymbol{n})^2}\\ && \widehat{\Pi}_{1; (n,H)} (u,v)=\frac{P_5(H, \boldsymbol{m}, \boldsymbol{n}) P_4(H, \boldsymbol{m}, \boldsymbol{n}) - P_3(H, \boldsymbol{m}, \boldsymbol{n}) P_2(H, \boldsymbol{m}, \boldsymbol{n})} {P_5(H, \boldsymbol{m}, \boldsymbol{n})P_1(H, \boldsymbol{m}, \boldsymbol{n}) - P_3(H, \boldsymbol{m}, \boldsymbol{n})^2},\end{aligned}\]

where we have introduced the following notation:

$\displaystyle P_1(H, \boldsymbol{m}, \boldsymbol{n}) = \sum_{t=1}^{m}\sum_{h=1}^H \mathbb{1}_{\{X_{t-1}^h=u \}} \left(\frac{t}{n}\right)^2,$ $\displaystyle P_2(H, \boldsymbol{m}, \boldsymbol{n}) = \sum_{t=1}^{m} \sum_{h=1}^H\mathbb{1}_{\{X_{t-1}^h=u, X_t^h=v \}} \left(1-\frac{t}{n}\right),$
$\displaystyle P_3(H, \boldsymbol{m}, \boldsymbol{n}) = \sum_{t=1}^{m}\sum_{h=1}^H\mathbb{1}_{\{X_{t-1}^h=u \}} \left(1-\frac{t}{n}\right)\left(\frac{t}{n}\right),$ $\displaystyle P_4(H, \boldsymbol{m}, \boldsymbol{n}) = \sum_{t=1}^{m}\sum_{h=1}^H\mathbb{1}_{\{X_{t-1}^h=u, X_t^h=v \}}\left(\frac{t}{n}\right),$
$\displaystyle P_5(H, \boldsymbol{m}, \boldsymbol{n}) = \displaystyle\sum_{t=1}^{m}\sum_{h=1}^H\mathbb{1}_{\{X_{t-1}^h=u \}}\left(1-\frac{t}{n}\right)^2,$ with $\boldsymbol{m}=(m_1, \ldots, m_H)$ and $\boldsymbol{n}=(n_1, \ldots, n_H).$

Note we can adapt the estimation procedures that we have previously obtained in order to get estimators of the drifting Markov models in the other cases; one can see Barbu and Vergne (2019) for more details. Using the expression of the reliability indicators of a drifting Markov chain previously obtained and the estimators of the characteristics of a drifting Markov chain, one immediately obtains the associated plug-in estimators of the reliability metrics.

3 The package

The package is principally devoted to the simulation and estimation of drifting Markov models, as well as to the estimation of associated reliability measures and to the computation of the probabilities of a word occurrence under a given model. All the different possibilities of the package are illustrated in Figure 2.

3.1 Estimation of drifting Markov models

The estimation of DMMs is carried out by the functions fitdmm, when starting from one (several) sample path (and obtain LSE). We will describe this function in the sequel.

The function fitdmm

The different arguments of this function are:

sequences: A list of character vector(s) representing one (several) sequence(s) from which the estimation is carried out
order: Order of the Markov chain
degree: Degree of the polynomials (e.g., linear drifting if degree=1, etc.)
states: Vector of states space of length s > 1
init.estim: Method used to estimate the initial law.

Here we have an example of an estimation of a drifting Markov model of order 1 and degree 1 using the function fitdmm, starting from a DNA sequence called lambda.

data(lambda, package = "drimmR")
states <- c("a","c","g","t")
dmm <- fitdmm(lambda, 1, 1, states, init.estim="freq", fit.method="sum")
## Le chargement a nécessité le package : foreach
## Le chargement a nécessité le package : iterators
## Le chargement a nécessité le package : parallel
dmm
## $states
## [1] "a" "c" "g" "t"
## 
## $order
## [1] 1
## 
## $degree
## [1] 1
## 
## $Polynomials
##     t^0 t^1
## A_0   1  -1
## A_1   0   1
## 
## $length
## [1] 48501
## 
## $matrices
## $matrices$Pi0
##           a         c         g         t
## a 0.2548330 0.2495270 0.2593907 0.2362493
## c 0.2353305 0.2465529 0.3360221 0.1820945
## g 0.2313610 0.3008693 0.2879023 0.1798674
## t 0.1356776 0.2174431 0.4129533 0.2339260
## 
## $matrices$Pi1
##           a         c         g         t
## a 0.3398007 0.1715523 0.1870335 0.3016135
## c 0.3339285 0.1911412 0.2075945 0.2673358
## g 0.2801369 0.2600487 0.2021395 0.2576749
## t 0.2218131 0.2283645 0.2305177 0.3193047
## 
## 
## $init.estim
## 
##         a         c         g         t 
## 0.2543400 0.2342171 0.2642778 0.2471651 
## 
## attr(,"class")
## [1] "dmm"

3.1.1 Estimation of corresponding model characteristics

Once a Markov drifting model has been estimated/constructed as described in the previous subsection, various characteristics of the model can be computed. These are: the log-likelihood of one or several sequences, the AIC and BIC information criteria, the stationary distribution on an entire sequence or only on a part of it, the distribution of the chain on the entire sequence or only on a part of it.

The functions loglik, aic and bic have the following arguments:

x: Object of class dmm
sequences: A list of character vector(s) representing one (several) sequence(s)

They return the numerical values of the log-likelihood, the AIC, the BIC of the sequence, respectively.

data(lambda, package = "drimmR")
sequence <- c("a","g","g","t","c","g","a","t","a","a","a")
dmm <-fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq")

loglik(dmm, sequence)
## [[1]]
## [1] -16.08786

aic(dmm, sequence)
## [[1]]
## [1] 80.17572

bic(dmm, sequence)
## [[1]]
## [1] 89.72521

The function getTransitionMatrix evaluates the transition matrix at a given position. The function has the following arguments :

x: Object of class dmm
pos: position along the sequence (integer)

data(lambda)
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'),init.estim = "freq")
t <- 10
getTransitionMatrix(dmm,pos=t)
##           a         c         g         t
## a 0.2548505 0.2495109 0.2593757 0.2362628
## c 0.2353509 0.2465415 0.3359956 0.1821121
## g 0.2313710 0.3008609 0.2878846 0.1798834
## t 0.1356954 0.2174453 0.4129157 0.2339436

The function getStationaryLaw evaluates the stationary law at a given position or for all positions along the list of sequence(s).

x: Object of class dmm
pos: position along the sequence (integer)
all.pos: FALSE (default, evaluation at pos index) ; TRUE (evaluation for all pos index)
internal: FALSE (default) ; TRUE (for internal use of fitdmm initial law)

data(lambda, package = "drimmR")
sequence <- sample(lambda, 30, replace=TRUE )
dmm <- fitdmm(sequence, 1, 1, c('a','c','g','t'), init.estim = "freq")
t <- 10
getStationaryLaw(dmm,pos=t)
##          a          c          g          t 
## 0.05329869 0.31418371 0.30487669 0.32764091
getStationaryLaw(dmm,all.pos=TRUE)
##                 a         c         g         t
## pos 1  0.03029875 0.3260096 0.2149455 0.4287462
## pos 2  0.03312811 0.3265113 0.2255174 0.4148431
## pos 3  0.03590952 0.3264993 0.2357585 0.4018327
## pos 4  0.03863421 0.3260104 0.2457567 0.3895987
## pos 5  0.04129360 0.3250730 0.2555927 0.3780407
## pos 6  0.04387909 0.3237083 0.2653414 0.3670712
## pos 7  0.04638187 0.3219310 0.2750746 0.3566125
## pos 8  0.04879275 0.3197500 0.2848618 0.3465955
## pos 9  0.05110190 0.3171683 0.2947723 0.3369575
## pos 10 0.05329869 0.3141837 0.3048767 0.3276409
## pos 11 0.05537143 0.3107884 0.3152481 0.3185921
## pos 12 0.05730709 0.3069687 0.3259641 0.3097601
## pos 13 0.05909095 0.3027051 0.3371083 0.3010957
## pos 14 0.06070628 0.2979708 0.3487727 0.2925502
## pos 15 0.06213378 0.2927319 0.3610597 0.2840746
## pos 16 0.06335102 0.2869451 0.3740857 0.2756182
## pos 17 0.06433169 0.2805567 0.3879841 0.2671276
## pos 18 0.06504463 0.2735004 0.4029103 0.2585447
## pos 19 0.06545256 0.2656943 0.4190477 0.2498055
## pos 20 0.06551052 0.2570370 0.4366152 0.2408372
## pos 21 0.06516360 0.2474028 0.4558779 0.2315556
## pos 22 0.06434406 0.2366345 0.4771606 0.2218609
## pos 23 0.06296727 0.2245336 0.5008668 0.2116324
## pos 24 0.06092606 0.2108471 0.5275058 0.2007210
## pos 25 0.05808260 0.1952480 0.5577304 0.1889390
## pos 26 0.05425659 0.1773069 0.5923925 0.1760440
## pos 27 0.04920750 0.1564495 0.6326267 0.1617162
## pos 28 0.04260709 0.1318905 0.6799796 0.1455228
## pos 29 0.03399540 0.1025281 0.7366160 0.1268606

3.1.2 Analysis of a model

Once a Markov drifting model has been estimated/constructed as described at the beginning of this section, the obtained model can be analyzed by means of several functions.

The computation of word probabilities according to an estimated/constructed DMM is carried out by means of the functions word_probability, word_probabilities and words_probabilities.

The function word_probability computes the probability of a word at a given position; it has the following arguments:

word: A word, i.e., a subsequence string of characters
pos: A position of the word along the sequence (numeric)
x: An object of class dmm
output_file: A file containing the probability
internal: FALSE (default) ; TRUE (for internal use of word applications)

It returns the numerical value of the probability.

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq")
word_probability("cgt",10,dmm,output_file=tempfile("drimmR"))
##        cgt 
## 0.01563115

The function word_probabilities computes the probabilities of a word at given positions. It has the following arguments:

word: A word, i.e., a subsequence string of characters
pos: Vector of integer positions of the word along the sequence; it is given by a start and end point
x: An object of class dmm
output_file: A file containing the probabilities
internal: FALSE (default) ; TRUE (for internal use of word applications)
plot: display a figure plot of word probabilities along the positions if TRUE

It returns a numerical vector of probabilities computed at the positions given in pos.

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq", fit.method="sum")
res <- word_probabilities("cgt", c(100, length(lambda)-3), dmm, output_file=tempfile("drimmR"), plot=TRUE)
head(res[[1]], n=10)
##    position probability
## 1       100  0.01562696
## 2       101  0.01562692
## 3       102  0.01562687
## 4       103  0.01562682
## 5       104  0.01562678
## 6       105  0.01562673
## 7       106  0.01562668
## 8       107  0.01562664
## 9       108  0.01562659
## 10      109  0.01562654
res[2]
## [[1]]

The function words_probabilities computes the probabilities of a word at given positions; it has the following arguments:

words: A vector of characters containing words
pos: Vector of integer positions of the word along the sequence; it is given by a start and end point
x: An object of class dmm
output_file: A file containing the matrix of probabilities
plot: display a figure plot of word probabilities along the positions if TRUE

It returns a data frame of probabilities computed for each element of words at the positions given in pos.

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq")
res <- words_probabilities(c("atcgattc", "taggct", "ggatcgg"), c(100, length(lambda)-8), dmm, output_file=tempfile("drimmR"), plot=TRUE)
head(res[[1]], n=10)
res[2]

The function lengthWord_probabilities computes the probabilities of occurrence of the observed word of given size in a sequence at several positions. This function has the following arguments:

n: Integer, the given length of the word
sequence: Vector of characters representing the sequence
pos: Vector of integer positions of the word along the sequence. It it is given by a start and end point
x: An object of class dmm
output_file A file containing the vector of probabilities
plot display several figure plots of probabilities of appearance of words along the positions if TRUE

It returns a data frame of probability by position

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq")
res <- lengthWord_probabilities(2, lambda, c(100, length(lambda)-2), dmm, output_file=tempfile("drimmR"), plot=TRUE)
res[2]
## [[1]]
## [[1]][[1]]

## 
## [[1]][[2]]
## Ignoring unknown labels:
## • fill : "Package :"

## 
## [[1]][[3]]
## Ignoring unknown labels:
## • fill : "Package :"

## 
## [[1]][[4]]
## Ignoring unknown labels:
## • fill : "Package :"

## 
## [[1]][[5]]
## Ignoring unknown labels:
## • fill : "Package :"

3.2 Simulation of drifting Markov models

The simulation of drifting Markov models is carried out by the function simulation that has the following arguments:

DRIMM: An object of class dmm
output_file: File containing the simulated sequence
model_size: Integer, the size of the model

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c('a','c','g','t'), init.estim = "freq")
# simulate a sequence of length 20 000 from dmm
simulate(dmm, tempfile("drimmR"), 20000)
## [1] "Write a simulated file from the model"
##     [1] "g" "t" "g" "c" "t" "g" "c" "c" "a" "c" "c" "g" "c" "c" "a" "g" "c" "c"
##    [19] "c" "a" "a" "g" "t" "g" "t" "c" "c" "c" "g" "a" "a" "t" "c" "c" "g" "a"
##    [37] "c" "g" "c" "a" "c" "g" "t" "t" "t" "a" "g" "t" "g" "a" "c" "a" "c" "g"
##    [55] "g" "g" "c" "c" "a" "c" "c" "g" "c" "g" "c" "g" "a" "t" "a" "c" "c" "c"
##    [73] "t" "c" "t" "t" "t" "a" "a" "c" "c" "c" "a" "t" "t" "t" "c" "a" "c" "c"
##    [91] "a" "g" "a" "c" "a" "a" "g" "c" "g" "c" "c" "g" "t" "t" "a" "t" "g" "t"
##   [109] "a" "g" "t" "c" "g" "g" "t" "a" "a" "a" "t" "c" "t" "a" "a" "g" "a" "g"
##   [127] "a" "g" "c" "a" "t" "t" "c" "a" "a" "c" "c" "a" "c" "t" "g" "g" "c" "a"
##   [145] "t" "c" "g" "g" "a" "c" "g" "c" "c" "a" "g" "c" "a" "c" "a" "g" "g" "c"
##   [163] "g" "g" "c" "g" "t" "c" "g" "a" "a" "c" "c" "g" "g" "g" "c" "g" "a" "g"
##   [181] "a" "t" "a" "c" "c" "t" "g" "g" "c" "c" "g" "c" "g" "a" "g" "a" "t" "c"
##   [199] "t" "g" "a" "c" "a" "g" "c" "a" "g" "c" "g" "t" "t" "t" "c" "c" "g" "t"
##   [217] "a" "c" "a" "g" "c" "g" "t" "c" "a" "a" "c" "g" "c" "c" "t" "c" "g" "c"
##   [235] "t" "t" "t" "c" "t" "t" "a" "a" "t" "c" "g" "a" "g" "c" "g" "g" "a" "t"
##   [253] "g" "a" "c" "g" "c" "g" "c" "g" "c" "g" "g" "g" "a" "a" "c" "a" "g" "a"
##   [271] "c" "g" "t" "c" "c" "c" "a" "c" "g" "t" "c" "c" "a" "c" "c" "c" "c" "a"
##   [289] "c" "c" "g" "c" "c" "a" "a" "c" "g" "a" "g" "c" "t" "c" "c" "c" "g" "t"
##   [307] "g" "c" "c" "t" "g" "c" "c" "g" "c" "g" "g" "c" "g" "g" "c" "a" "t" "g"
##   [325] "a" "c" "c" "g" "c" "a" "c" "c" "g" "c" "c" "g" "a" "t" "c" "c" "g" "c"
##   [343] "a" "c" "g" "a" "a" "c" "t" "g" "c" "t" "g" "g" "a" "c" "a" "g" "t" "c"
##   [361] "c" "t" "c" "a" "t" "a" "g" "c" "g" "g" "g" "t" "g" "g" "a" "c" "g" "g"
##   [379] "t" "g" "t" "g" "c" "g" "c" "c" "t" "t" "t" "t" "c" "a" "c" "a" "t" "c"
##   [397] "a" "g" "g" "g" "g" "t" "a" "t" "g" "a" "c" "t" "t" "g" "g" "g" "a" "g"
##   [415] "a" "g" "g" "a" "c" "c" "g" "a" "t" "t" "c" "g" "g" "g" "g" "g" "g" "g"
##   [433] "a" "c" "c" "c" "c" "c" "a" "c" "g" "a" "t" "g" "t" "g" "c" "g" "c" "c"
##   [451] "g" "a" "t" "g" "t" "c" "a" "g" "g" "g" "g" "t" "g" "c" "g" "c" "a" "t"
##   [469] "g" "a" "g" "g" "t" "g" "c" "c" "c" "c" "g" "c" "t" "a" "g" "a" "g" "c"
##   [487] "a" "c" "c" "t" "t" "t" "g" "a" "c" "t" "a" "t" "c" "c" "c" "a" "t" "a"
##   [505] "g" "c" "c" "g" "c" "a" "g" "g" "c" "g" "a" "t" "t" "c" "g" "a" "t" "c"
##   [523] "g" "t" "c" "g" "c" "t" "a" "c" "g" "t" "g" "a" "t" "g" "c" "c" "a" "c"
##   [541] "c" "a" "a" "t" "a" "c" "a" "a" "g" "a" "g" "t" "t" "g" "t" "t" "a" "g"
##   [559] "c" "a" "t" "c" "g" "a" "c" "g" "a" "a" "a" "g" "c" "t" "t" "t" "t" "g"
##   [577] "g" "c" "t" "c" "t" "g" "g" "a" "g" "c" "g" "c" "c" "g" "t" "g" "c" "c"
##   [595] "t" "c" "a" "c" "c" "a" "c" "g" "a" "t" "t" "a" "a" "g" "g" "t" "c" "g"
##   [613] "a" "a" "a" "c" "g" "g" "t" "c" "c" "c" "c" "g" "a" "a" "c" "a" "g" "a"
##   [631] "t" "g" "c" "t" "t" "a" "g" "c" "a" "a" "t" "g" "g" "g" "c" "a" "a" "c"
##   [649] "t" "a" "c" "c" "c" "a" "g" "g" "t" "g" "c" "c" "c" "a" "c" "g" "g" "c"
##   [667] "t" "g" "c" "g" "g" "a" "a" "g" "g" "g" "a" "c" "g" "a" "g" "a" "a" "t"
##   [685] "t" "c" "c" "a" "g" "c" "a" "a" "g" "a" "c" "a" "c" "g" "g" "t" "a" "a"
##   [703] "t" "c" "a" "g" "c" "c" "a" "g" "c" "g" "c" "a" "g" "c" "a" "c" "t" "t"
##   [721] "t" "c" "g" "t" "g" "g" "t" "c" "c" "g" "g" "a" "c" "g" "c" "a" "g" "g"
##   [739] "a" "c" "a" "t" "c" "c" "g" "c" "g" "t" "c" "g" "a" "c" "g" "a" "a" "c"
##   [757] "t" "g" "t" "t" "c" "g" "g" "g" "c" "g" "t" "t" "g" "g" "a" "g" "a" "g"
##   [775] "a" "c" "c" "g" "a" "c" "t" "g" "t" "g" "c" "g" "c" "c" "g" "t" "t" "g"
##   [793] "t" "c" "g" "a" "c" "g" "c" "a" "t" "a" "t" "t" "a" "c" "t" "g" "c" "a"
##   [811] "a" "c" "a" "g" "t" "g" "a" "t" "t" "c" "g" "g" "g" "a" "g" "a" "c" "a"
##   [829] "g" "c" "t" "c" "g" "c" "c" "c" "a" "a" "c" "g" "t" "c" "c" "c" "a" "c"
##   [847] "c" "a" "t" "c" "g" "a" "c" "c" "g" "c" "a" "c" "c" "g" "c" "c" "a" "g"
##   [865] "t" "g" "g" "g" "a" "a" "t" "g" "c" "g" "a" "c" "a" "a" "g" "c" "c" "a"
##   [883] "c" "t" "a" "g" "c" "t" "c" "a" "a" "a" "c" "g" "t" "c" "g" "a" "c" "t"
##   [901] "a" "c" "a" "a" "g" "c" "g" "a" "g" "g" "c" "c" "c" "a" "g" "c" "g" "g"
##   [919] "c" "t" "g" "g" "t" "g" "g" "g" "c" "c" "c" "c" "g" "t" "g" "c" "a" "g"
##   [937] "a" "g" "c" "a" "g" "a" "c" "a" "g" "g" "c" "a" "g" "t" "g" "a" "a" "g"
##   [955] "t" "a" "t" "a" "c" "g" "a" "g" "g" "c" "t" "t" "a" "t" "c" "g" "g" "c"
##   [973] "g" "a" "t" "g" "g" "g" "a" "g" "c" "g" "g" "a" "g" "g" "g" "g" "c" "g"
##   [991] "t" "c" "c" "c" "t" "g" "t" "g" "g" "c" "g" "c" "g" "g" "t" "t" "c" "a"
##  [1009] "c" "a" "g" "c" "a" "a" "g" "a" "g" "a" "g" "a" "t" "g" "c" "t" "a" "c"
##  [1027] "a" "g" "g" "a" "a" "c" "c" "t" "g" "g" "t" "a" "g" "c" "t" "c" "c" "g"
##  [1045] "t" "a" "g" "c" "g" "c" "c" "c" "g" "g" "a" "c" "t" "c" "t" "a" "t" "a"
##  [1063] "a" "c" "g" "t" "t" "g" "a" "a" "a" "a" "g" "c" "a" "g" "t" "t" "g" "c"
##  [1081] "g" "a" "c" "g" "g" "g" "g" "t" "t" "a" "t" "t" "g" "g" "c" "g" "g" "g"
##  [1099] "g" "c" "a" "g" "g" "g" "a" "g" "a" "a" "c" "c" "g" "a" "g" "t" "g" "t"
##  [1117] "a" "c" "c" "t" "c" "c" "g" "c" "a" "g" "t" "c" "g" "c" "c" "t" "g" "g"
##  [1135] "g" "g" "t" "a" "t" "g" "t" "g" "c" "a" "g" "a" "a" "c" "c" "c" "c" "c"
##  [1153] "c" "g" "g" "t" "t" "g" "a" "a" "a" "g" "a" "a" "t" "t" "c" "t" "t" "t"
##  [1171] "c" "t" "t" "t" "g" "c" "a" "g" "g" "g" "a" "g" "c" "c" "c" "t" "c" "t"
##  [1189] "a" "c" "a" "c" "a" "a" "g" "c" "g" "c" "a" "a" "g" "a" "t" "g" "g" "t"
##  [1207] "g" "c" "c" "a" "c" "g" "c" "a" "t" "t" "g" "g" "c" "t" "c" "c" "c" "a"
##  [1225] "a" "c" "c" "a" "a" "a" "g" "t" "a" "g" "t" "g" "a" "a" "g" "t" "a" "c"
##  [1243] "g" "c" "g" "c" "g" "t" "g" "t" "a" "c" "t" "g" "a" "a" "a" "c" "a" "c"
##  [1261] "a" "g" "c" "c" "g" "g" "g" "g" "a" "c" "t" "c" "a" "a" "c" "g" "g" "g"
##  [1279] "g" "t" "t" "c" "a" "c" "c" "c" "c" "g" "c" "a" "t" "c" "a" "c" "c" "g"
##  [1297] "g" "g" "c" "g" "g" "g" "g" "g" "a" "a" "g" "g" "t" "g" "g" "c" "a" "g"
##  [1315] "a" "c" "c" "c" "a" "c" "t" "a" "g" "a" "t" "c" "c" "a" "c" "t" "t" "t"
##  [1333] "g" "a" "t" "g" "g" "c" "c" "g" "c" "a" "g" "a" "g" "c" "g" "a" "t" "a"
##  [1351] "a" "g" "a" "g" "g" "c" "a" "c" "c" "a" "g" "a" "g" "a" "c" "g" "g" "t"
##  [1369] "t" "t" "c" "t" "g" "c" "t" "t" "g" "t" "c" "c" "g" "g" "c" "g" "c" "c"
##  [1387] "c" "a" "c" "a" "a" "g" "g" "t" "t" "t" "t" "g" "a" "t" "a" "a" "c" "g"
##  [1405] "a" "g" "g" "g" "g" "t" "a" "t" "t" "c" "a" "t" "t" "t" "c" "g" "a" "a"
##  [1423] "t" "g" "c" "t" "c" "c" "t" "g" "g" "c" "g" "c" "c" "a" "g" "g" "c" "c"
##  [1441] "g" "t" "g" "t" "a" "g" "c" "a" "c" "c" "c" "g" "a" "t" "g" "t" "c" "g"
##  [1459] "g" "t" "g" "t" "a" "t" "a" "g" "c" "a" "g" "a" "g" "a" "g" "c" "c" "g"
##  [1477] "t" "a" "g" "t" "c" "t" "g" "a" "a" "g" "g" "t" "c" "t" "g" "c" "c" "c"
##  [1495] "c" "a" "c" "a" "c" "g" "t" "a" "t" "g" "g" "c" "a" "g" "c" "a" "c" "a"
##  [1513] "t" "g" "c" "a" "a" "g" "g" "a" "a" "a" "g" "g" "a" "t" "c" "g" "a" "t"
##  [1531] "c" "c" "c" "c" "c" "g" "g" "a" "g" "t" "a" "c" "g" "c" "t" "a" "t" "g"
##  [1549] "a" "g" "t" "t" "c" "g" "t" "t" "c" "c" "g" "g" "t" "a" "t" "c" "a" "c"
##  [1567] "a" "c" "g" "c" "g" "c" "a" "t" "c" "c" "t" "c" "c" "a" "c" "c" "c" "a"
##  [1585] "t" "g" "g" "a" "c" "t" "a" "a" "t" "a" "g" "g" "a" "c" "c" "t" "c" "g"
##  [1603] "c" "t" "c" "a" "c" "a" "a" "a" "t" "a" "t" "t" "c" "g" "c" "c" "c" "g"
##  [1621] "a" "t" "a" "a" "g" "g" "c" "a" "a" "g" "g" "a" "c" "a" "a" "t" "g" "g"
##  [1639] "c" "c" "t" "t" "t" "c" "t" "a" "c" "a" "g" "c" "a" "a" "c" "t" "c" "g"
##  [1657] "a" "a" "t" "a" "c" "c" "c" "c" "a" "c" "c" "g" "c" "g" "a" "g" "t" "c"
##  [1675] "t" "a" "g" "g" "c" "c" "c" "g" "c" "c" "c" "a" "g" "g" "g" "c" "g" "c"
##  [1693] "g" "c" "a" "t" "g" "c" "c" "t" "t" "c" "c" "g" "c" "c" "t" "c" "c" "g"
##  [1711] "c" "t" "t" "g" "t" "t" "g" "c" "t" "a" "t" "g" "t" "t" "c" "c" "t" "c"
##  [1729] "g" "g" "g" "a" "a" "g" "c" "a" "g" "g" "g" "g" "g" "g" "c" "c" "c" "c"
##  [1747] "c" "t" "c" "a" "a" "g" "a" "g" "t" "g" "a" "t" "g" "c" "t" "g" "c" "t"
##  [1765] "c" "c" "g" "t" "c" "a" "g" "t" "g" "a" "t" "c" "c" "c" "c" "g" "g" "g"
##  [1783] "g" "c" "g" "t" "t" "t" "t" "g" "a" "c" "c" "c" "g" "t" "c" "c" "c" "g"
##  [1801] "t" "c" "a" "g" "g" "c" "g" "a" "g" "a" "g" "a" "c" "c" "g" "t" "g" "t"
##  [1819] "g" "c" "c" "t" "c" "c" "c" "g" "g" "a" "a" "g" "t" "a" "g" "a" "g" "a"
##  [1837] "t" "a" "a" "a" "a" "c" "c" "c" "t" "g" "c" "g" "c" "g" "g" "c" "a" "a"
##  [1855] "t" "t" "t" "c" "a" "t" "t" "g" "g" "c" "c" "a" "t" "g" "c" "g" "g" "t"
##  [1873] "g" "g" "a" "t" "g" "a" "c" "t" "a" "c" "g" "t" "a" "a" "g" "a" "a" "c"
##  [1891] "c" "c" "c" "t" "t" "c" "c" "t" "c" "c" "c" "c" "g" "a" "a" "c" "a" "t"
##  [1909] "a" "c" "c" "g" "t" "c" "a" "a" "g" "g" "g" "t" "t" "g" "c" "t" "g" "t"
##  [1927] "t" "g" "t" "c" "c" "c" "g" "c" "c" "a" "c" "c" "t" "t" "a" "a" "g" "t"
##  [1945] "g" "g" "g" "c" "a" "g" "a" "t" "c" "t" "a" "a" "t" "c" "c" "t" "a" "g"
##  [1963] "g" "g" "a" "g" "c" "g" "c" "a" "t" "c" "g" "a" "g" "g" "a" "g" "c" "t"
##  [1981] "a" "t" "c" "t" "a" "c" "c" "t" "g" "g" "a" "t" "a" "t" "a" "a" "c" "t"
##  [1999] "c" "a" "a" "c" "g" "c" "t" "g" "c" "t" "c" "a" "a" "t" "a" "c" "a" "c"
##  [2017] "a" "c" "g" "c" "g" "c" "c" "g" "c" "t" "t" "t" "g" "c" "t" "a" "c" "c"
##  [2035] "a" "t" "c" "c" "a" "c" "t" "c" "c" "t" "a" "a" "a" "g" "t" "a" "g" "c"
##  [2053] "c" "t" "g" "t" "g" "g" "g" "t" "c" "a" "c" "c" "a" "g" "g" "g" "t" "a"
##  [2071] "t" "a" "g" "g" "g" "c" "a" "t" "c" "c" "a" "a" "g" "g" "c" "a" "t" "a"
##  [2089] "g" "t" "t" "c" "g" "g" "t" "g" "g" "a" "a" "a" "g" "a" "g" "c" "t" "c"
##  [2107] "c" "t" "g" "g" "t" "a" "a" "c" "c" "c" "t" "c" "a" "a" "g" "g" "g" "c"
##  [2125] "a" "g" "c" "c" "c" "c" "a" "t" "t" "a" "g" "c" "g" "t" "t" "g" "c" "a"
##  [2143] "g" "a" "t" "g" "g" "c" "c" "a" "a" "a" "c" "c" "c" "c" "c" "t" "t" "a"
##  [2161] "c" "g" "a" "a" "g" "g" "c" "t" "t" "c" "g" "g" "a" "g" "a" "t" "c" "g"
##  [2179] "c" "g" "c" "c" "c" "c" "c" "g" "g" "t" "g" "g" "t" "t" "c" "a" "c" "g"
##  [2197] "a" "c" "c" "c" "g" "t" "g" "a" "g" "t" "g" "c" "a" "g" "c" "t" "t" "g"
##  [2215] "g" "g" "c" "a" "a" "g" "g" "g" "g" "t" "c" "c" "t" "t" "g" "g" "c" "c"
##  [2233] "a" "g" "t" "c" "t" "t" "a" "g" "g" "t" "c" "g" "g" "c" "c" "a" "t" "c"
##  [2251] "a" "c" "c" "c" "g" "c" "g" "t" "a" "g" "t" "c" "a" "c" "c" "g" "t" "g"
##  [2269] "a" "a" "a" "c" "c" "t" "a" "t" "c" "c" "g" "a" "g" "g" "g" "c" "g" "c"
##  [2287] "t" "t" "a" "g" "g" "c" "c" "c" "c" "a" "g" "c" "g" "g" "c" "t" "t" "c"
##  [2305] "c" "a" "c" "a" "g" "t" "g" "c" "g" "a" "c" "g" "t" "t" "a" "c" "t" "t"
##  [2323] "t" "a" "g" "g" "t" "c" "a" "g" "a" "a" "g" "t" "g" "c" "a" "g" "a" "g"
##  [2341] "g" "c" "t" "g" "g" "c" "g" "c" "a" "g" "t" "a" "c" "g" "c" "a" "c" "g"
##  [2359] "c" "c" "c" "g" "a" "g" "a" "a" "a" "c" "g" "g" "c" "g" "t" "t" "a" "g"
##  [2377] "t" "a" "t" "g" "t" "c" "a" "c" "t" "g" "c" "g" "c" "t" "a" "t" "t" "a"
##  [2395] "a" "a" "c" "c" "c" "a" "c" "t" "g" "a" "a" "c" "a" "c" "c" "t" "c" "c"
##  [2413] "c" "a" "t" "a" "g" "c" "g" "c" "t" "g" "c" "c" "c" "g" "c" "c" "a" "c"
##  [2431] "c" "t" "g" "c" "t" "g" "t" "c" "a" "c" "c" "g" "g" "g" "g" "g" "c" "c"
##  [2449] "a" "t" "t" "t" "a" "a" "c" "a" "g" "c" "g" "g" "g" "t" "t" "t" "t" "t"
##  [2467] "c" "t" "g" "g" "c" "t" "t" "a" "g" "t" "c" "a" "c" "t" "a" "c" "g" "a"
##  [2485] "g" "g" "a" "a" "c" "c" "a" "g" "g" "a" "a" "a" "g" "t" "a" "a" "t" "g"
##  [2503] "a" "c" "t" "a" "a" "t" "a" "c" "a" "t" "g" "a" "c" "t" "g" "g" "g" "a"
##  [2521] "c" "a" "a" "t" "g" "g" "t" "c" "g" "g" "c" "a" "c" "a" "g" "c" "t" "c"
##  [2539] "t" "a" "c" "t" "a" "a" "a" "a" "a" "g" "g" "a" "t" "a" "g" "a" "t" "t"
##  [2557] "t" "g" "a" "g" "c" "a" "g" "a" "g" "a" "g" "c" "g" "t" "g" "c" "g" "a"
##  [2575] "g" "c" "c" "a" "c" "a" "c" "g" "c" "g" "c" "g" "a" "a" "a" "c" "g" "c"
##  [2593] "a" "t" "c" "g" "c" "a" "a" "g" "a" "g" "a" "c" "t" "a" "a" "g" "t" "c"
##  [2611] "g" "a" "t" "a" "c" "c" "g" "t" "a" "t" "a" "c" "c" "c" "c" "c" "g" "g"
##  [2629] "c" "c" "c" "a" "c" "c" "a" "a" "g" "t" "a" "g" "t" "c" "g" "t" "g" "g"
##  [2647] "a" "a" "c" "g" "t" "a" "a" "t" "c" "g" "c" "g" "t" "a" "t" "c" "a" "t"
##  [2665] "a" "t" "a" "c" "g" "c" "a" "t" "g" "g" "a" "a" "g" "a" "a" "a" "c" "t"
##  [2683] "a" "c" "c" "g" "g" "a" "c" "a" "t" "g" "a" "g" "t" "a" "t" "t" "a" "c"
##  [2701] "c" "c" "t" "a" "c" "c" "g" "a" "t" "c" "c" "a" "c" "g" "g" "a" "t" "g"
##  [2719] "g" "t" "c" "g" "a" "g" "a" "c" "c" "t" "a" "a" "a" "g" "t" "g" "a" "a"
##  [2737] "a" "a" "c" "c" "a" "g" "c" "g" "g" "g" "c" "t" "g" "g" "c" "a" "g" "a"
##  [2755] "t" "t" "c" "c" "g" "a" "t" "a" "c" "g" "g" "c" "t" "c" "g" "a" "c" "t"
##  [2773] "c" "c" "a" "c" "a" "t" "t" "c" "a" "c" "a" "t" "c" "c" "a" "g" "a" "a"
##  [2791] "t" "a" "g" "g" "a" "g" "c" "a" "c" "g" "g" "a" "c" "g" "t" "a" "c" "t"
##  [2809] "t" "c" "t" "a" "g" "g" "c" "c" "c" "c" "c" "a" "g" "g" "g" "g" "c" "c"
##  [2827] "a" "g" "t" "c" "g" "g" "a" "g" "a" "g" "t" "c" "a" "c" "c" "c" "c" "g"
##  [2845] "c" "c" "g" "a" "c" "g" "g" "g" "a" "c" "a" "c" "a" "c" "c" "g" "g" "g"
##  [2863] "g" "t" "c" "g" "c" "c" "g" "t" "c" "t" "g" "t" "t" "a" "c" "g" "c" "c"
##  [2881] "c" "c" "g" "a" "c" "c" "g" "c" "a" "t" "t" "c" "c" "g" "a" "a" "g" "c"
##  [2899] "c" "a" "c" "g" "g" "g" "t" "c" "t" "g" "t" "g" "t" "c" "g" "t" "g" "t"
##  [2917] "t" "c" "t" "a" "a" "t" "g" "g" "g" "c" "c" "c" "a" "g" "g" "a" "a" "g"
##  [2935] "a" "c" "c" "c" "t" "a" "g" "c" "g" "t" "c" "c" "g" "a" "a" "c" "a" "g"
##  [2953] "g" "t" "t" "c" "g" "c" "t" "a" "g" "a" "t" "a" "c" "a" "g" "t" "a" "g"
##  [2971] "g" "c" "g" "g" "c" "a" "c" "t" "c" "c" "g" "g" "c" "g" "c" "g" "a" "c"
##  [2989] "t" "t" "t" "t" "t" "t" "g" "g" "g" "t" "a" "t" "c" "g" "c" "a" "c" "c"
##  [3007] "g" "a" "g" "t" "c" "g" "t" "c" "c" "a" "a" "a" "g" "g" "t" "g" "t" "c"
##  [3025] "t" "g" "a" "g" "g" "a" "a" "g" "c" "g" "g" "t" "g" "c" "a" "c" "c" "g"
##  [3043] "c" "a" "g" "t" "g" "g" "c" "c" "t" "t" "g" "c" "t" "a" "c" "t" "g" "c"
##  [3061] "c" "g" "t" "t" "c" "g" "c" "a" "g" "t" "a" "a" "a" "c" "a" "c" "t" "a"
##  [3079] "c" "g" "g" "g" "g" "a" "g" "c" "c" "t" "a" "t" "g" "t" "a" "c" "g" "g"
##  [3097] "g" "c" "t" "t" "g" "c" "t" "t" "c" "a" "t" "a" "g" "a" "g" "a" "g" "g"
##  [3115] "t" "g" "t" "c" "t" "c" "c" "c" "a" "t" "g" "c" "c" "a" "g" "a" "g" "t"
##  [3133] "a" "g" "a" "a" "c" "g" "c" "g" "t" "a" "c" "c" "a" "c" "g" "g" "a" "c"
##  [3151] "g" "c" "g" "c" "g" "a" "t" "c" "c" "c" "c" "a" "a" "g" "t" "c" "c" "c"
##  [3169] "a" "c" "g" "c" "c" "c" "g" "g" "c" "a" "a" "c" "g" "t" "c" "c" "g" "t"
##  [3187] "g" "g" "a" "g" "t" "t" "g" "g" "g" "t" "a" "t" "c" "c" "g" "t" "c" "a"
##  [3205] "g" "c" "g" "a" "a" "t" "t" "g" "a" "g" "t" "a" "a" "t" "g" "g" "c" "g"
##  [3223] "g" "a" "t" "t" "g" "c" "g" "t" "c" "g" "c" "a" "a" "c" "g" "a" "c" "c"
##  [3241] "g" "c" "c" "g" "g" "c" "c" "t" "g" "a" "t" "g" "a" "g" "g" "t" "g" "g"
##  [3259] "t" "g" "c" "t" "a" "a" "c" "g" "a" "t" "t" "c" "a" "g" "t" "g" "c" "g"
##  [3277] "t" "t" "t" "g" "c" "t" "a" "g" "c" "c" "t" "t" "t" "t" "t" "c" "t" "a"
##  [3295] "a" "a" "c" "c" "t" "c" "a" "c" "a" "g" "a" "c" "g" "g" "g" "g" "a" "g"
##  [3313] "c" "c" "c" "t" "g" "g" "a" "a" "c" "g" "c" "g" "t" "g" "g" "t" "g" "t"
##  [3331] "t" "c" "a" "g" "c" "a" "g" "c" "a" "g" "c" "c" "t" "c" "t" "c" "c" "g"
##  [3349] "c" "c" "t" "g" "a" "g" "t" "g" "c" "g" "c" "c" "c" "g" "a" "t" "g" "c"
##  [3367] "g" "c" "c" "a" "c" "c" "a" "a" "c" "t" "a" "c" "t" "a" "a" "c" "t" "c"
##  [3385] "g" "t" "c" "g" "g" "c" "a" "c" "c" "c" "c" "c" "a" "g" "t" "a" "a" "a"
##  [3403] "g" "g" "g" "g" "c" "c" "t" "c" "a" "c" "t" "g" "t" "a" "g" "t" "c" "g"
##  [3421] "g" "t" "t" "c" "a" "g" "t" "c" "g" "c" "a" "g" "a" "t" "t" "a" "t" "t"
##  [3439] "c" "c" "t" "g" "a" "g" "g" "a" "c" "c" "t" "t" "g" "c" "g" "g" "c" "t"
##  [3457] "t" "g" "t" "t" "a" "a" "g" "g" "a" "g" "c" "a" "c" "g" "c" "t" "t" "c"
##  [3475] "g" "g" "c" "a" "a" "c" "c" "g" "t" "c" "c" "a" "a" "t" "c" "g" "a" "g"
##  [3493] "a" "g" "c" "t" "a" "c" "t" "c" "a" "g" "g" "c" "g" "c" "a" "c" "t" "t"
##  [3511] "g" "a" "c" "c" "t" "g" "a" "c" "a" "a" "c" "a" "g" "c" "g" "c" "t" "a"
##  [3529] "a" "a" "t" "g" "t" "a" "c" "g" "a" "g" "c" "g" "c" "c" "c" "c" "c" "c"
##  [3547] "a" "a" "c" "c" "t" "a" "c" "c" "a" "a" "a" "t" "g" "t" "c" "t" "c" "c"
##  [3565] "a" "c" "c" "g" "a" "c" "c" "t" "g" "g" "c" "g" "c" "a" "a" "t" "a" "c"
##  [3583] "a" "g" "c" "g" "t" "c" "t" "g" "g" "g" "a" "c" "g" "c" "g" "g" "c" "a"
##  [3601] "g" "t" "a" "c" "c" "c" "c" "a" "g" "g" "c" "a" "c" "a" "a" "a" "g" "a"
##  [3619] "g" "a" "c" "g" "c" "a" "c" "a" "t" "g" "g" "c" "g" "a" "g" "g" "a" "c"
##  [3637] "c" "g" "g" "c" "a" "g" "c" "g" "a" "t" "c" "c" "t" "c" "c" "a" "a" "g"
##  [3655] "g" "c" "c" "c" "c" "c" "a" "a" "a" "c" "g" "c" "g" "t" "a" "t" "c" "g"
##  [3673] "c" "a" "g" "g" "a" "t" "c" "g" "g" "t" "a" "t" "g" "t" "a" "a" "g" "g"
##  [3691] "c" "c" "g" "c" "c" "c" "g" "t" "c" "a" "c" "t" "c" "c" "a" "g" "a" "t"
##  [3709] "t" "g" "a" "t" "g" "c" "t" "c" "g" "a" "g" "g" "t" "g" "a" "c" "a" "a"
##  [3727] "t" "g" "c" "c" "g" "a" "g" "t" "a" "g" "t" "c" "g" "c" "g" "c" "a" "g"
##  [3745] "c" "g" "a" "g" "c" "a" "t" "g" "c" "a" "a" "t" "a" "t" "c" "t" "a" "a"
##  [3763] "a" "a" "a" "a" "a" "a" "t" "a" "c" "g" "c" "c" "a" "c" "a" "g" "g" "t"
##  [3781] "g" "a" "c" "c" "a" "t" "t" "a" "c" "a" "g" "g" "g" "a" "t" "c" "g" "g"
##  [3799] "t" "c" "g" "a" "g" "c" "g" "a" "t" "a" "t" "g" "c" "g" "a" "a" "a" "c"
##  [3817] "c" "g" "a" "g" "t" "g" "c" "a" "g" "c" "g" "c" "c" "g" "c" "g" "c" "c"
##  [3835] "t" "a" "c" "g" "a" "t" "a" "c" "a" "c" "t" "g" "t" "c" "c" "t" "a" "a"
##  [3853] "a" "g" "a" "a" "t" "c" "g" "g" "c" "c" "g" "a" "c" "t" "g" "t" "c" "a"
##  [3871] "t" "t" "t" "a" "a" "t" "t" "g" "g" "a" "c" "g" "a" "c" "c" "t" "g" "a"
##  [3889] "a" "a" "c" "a" "c" "c" "c" "c" "t" "c" "c" "g" "t" "g" "t" "t" "g" "g"
##  [3907] "c" "g" "a" "a" "g" "t" "c" "c" "a" "t" "c" "g" "a" "g" "g" "g" "c" "t"
##  [3925] "g" "c" "g" "a" "g" "c" "a" "g" "t" "a" "c" "a" "c" "g" "a" "a" "a" "a"
##  [3943] "g" "g" "g" "c" "g" "c" "a" "g" "c" "g" "a" "g" "g" "a" "t" "g" "c" "c"
##  [3961] "g" "g" "t" "a" "g" "c" "c" "c" "c" "a" "c" "a" "c" "a" "g" "a" "t" "a"
##  [3979] "t" "g" "t" "c" "c" "t" "g" "c" "c" "c" "g" "c" "g" "g" "t" "a" "a" "c"
##  [3997] "c" "c" "a" "c" "c" "t" "t" "c" "g" "c" "c" "c" "a" "t" "t" "c" "c" "t"
##  [4015] "a" "c" "t" "c" "c" "c" "c" "t" "c" "a" "c" "a" "g" "g" "g" "a" "g" "t"
##  [4033] "a" "c" "c" "g" "a" "c" "c" "a" "c" "a" "g" "g" "t" "a" "g" "c" "a" "t"
##  [4051] "c" "t" "g" "t" "g" "g" "g" "t" "a" "t" "g" "a" "a" "t" "a" "a" "t" "c"
##  [4069] "g" "g" "a" "c" "t" "a" "c" "c" "c" "c" "t" "g" "c" "c" "g" "a" "g" "c"
##  [4087] "a" "c" "a" "a" "g" "g" "a" "a" "c" "c" "t" "c" "g" "t" "c" "t" "t" "c"
##  [4105] "t" "a" "a" "g" "c" "g" "g" "c" "g" "a" "g" "c" "g" "a" "g" "a" "c" "a"
##  [4123] "a" "c" "c" "g" "a" "t" "a" "g" "t" "c" "g" "a" "g" "t" "g" "c" "g" "a"
##  [4141] "g" "a" "a" "c" "t" "c" "a" "c" "g" "g" "a" "g" "t" "c" "c" "a" "a" "t"
##  [4159] "c" "a" "g" "c" "a" "a" "a" "g" "g" "c" "t" "a" "t" "a" "a" "a" "c" "g"
##  [4177] "a" "g" "t" "c" "t" "a" "g" "c" "g" "g" "a" "g" "g" "a" "c" "c" "c" "t"
##  [4195] "g" "g" "t" "t" "t" "t" "c" "g" "a" "c" "g" "g" "c" "a" "a" "c" "t" "c"
##  [4213] "a" "g" "a" "a" "g" "c" "c" "a" "c" "t" "g" "a" "c" "a" "g" "c" "c" "c"
##  [4231] "a" "g" "a" "c" "a" "a" "c" "c" "g" "t" "c" "g" "c" "a" "g" "a" "t" "t"
##  [4249] "a" "g" "t" "t" "g" "g" "a" "g" "c" "c" "c" "t" "a" "g" "a" "c" "g" "t"
##  [4267] "a" "g" "c" "a" "t" "a" "a" "a" "t" "t" "t" "t" "g" "a" "g" "c" "c" "t"
##  [4285] "t" "t" "c" "c" "g" "t" "c" "t" "g" "c" "c" "c" "c" "t" "t" "c" "c" "t"
##  [4303] "g" "t" "g" "c" "t" "a" "t" "g" "a" "t" "c" "a" "c" "c" "c" "a" "c" "g"
##  [4321] "c" "c" "t" "g" "t" "a" "t" "a" "a" "a" "c" "g" "a" "t" "a" "c" "t" "a"
##  [4339] "a" "g" "c" "c" "a" "g" "c" "t" "g" "a" "c" "a" "c" "a" "t" "g" "g" "c"
##  [4357] "c" "g" "a" "t" "a" "a" "g" "g" "c" "t" "c" "t" "a" "c" "c" "c" "g" "a"
##  [4375] "g" "a" "c" "g" "t" "g" "g" "g" "g" "c" "a" "a" "a" "a" "c" "g" "g" "a"
##  [4393] "t" "g" "c" "c" "a" "a" "a" "a" "t" "c" "c" "a" "t" "a" "c" "t" "a" "g"
##  [4411] "a" "c" "c" "g" "c" "c" "g" "a" "c" "c" "c" "g" "t" "t" "g" "t" "t" "c"
##  [4429] "g" "a" "g" "c" "c" "g" "a" "t" "g" "c" "g" "c" "c" "g" "c" "t" "g" "a"
##  [4447] "g" "a" "a" "a" "c" "a" "c" "c" "a" "c" "g" "c" "g" "c" "g" "t" "g" "a"
##  [4465] "c" "t" "g" "a" "a" "t" "a" "t" "t" "g" "g" "t" "t" "g" "a" "a" "g" "c"
##  [4483] "c" "t" "a" "a" "a" "a" "a" "c" "g" "a" "g" "c" "g" "g" "g" "g" "t" "c"
##  [4501] "t" "g" "c" "c" "c" "a" "g" "t" "g" "g" "a" "g" "c" "t" "c" "c" "g" "t"
##  [4519] "g" "t" "c" "t" "g" "g" "c" "t" "c" "c" "g" "g" "a" "t" "a" "a" "t" "t"
##  [4537] "c" "g" "c" "a" "g" "c" "a" "t" "c" "g" "t" "t" "a" "a" "t" "g" "a" "t"
##  [4555] "t" "c" "t" "c" "g" "t" "c" "a" "g" "a" "t" "c" "c" "g" "a" "t" "t" "c"
##  [4573] "g" "g" "a" "t" "a" "g" "g" "a" "a" "c" "c" "a" "t" "t" "g" "a" "c" "t"
##  [4591] "g" "g" "c" "c" "c" "t" "g" "g" "a" "g" "a" "c" "c" "g" "c" "t" "c" "a"
##  [4609] "a" "a" "g" "g" "g" "a" "g" "g" "c" "a" "a" "g" "c" "a" "g" "g" "g" "a"
##  [4627] "t" "c" "g" "c" "c" "a" "g" "g" "c" "c" "g" "g" "c" "g" "a" "c" "g" "a"
##  [4645] "g" "g" "g" "c" "a" "c" "a" "g" "t" "t" "a" "c" "g" "c" "g" "c" "a" "g"
##  [4663] "a" "a" "c" "c" "t" "t" "c" "c" "g" "g" "a" "a" "c" "g" "g" "g" "g" "c"
##  [4681] "t" "a" "c" "t" "c" "c" "c" "c" "c" "c" "g" "c" "g" "a" "c" "c" "g" "t"
##  [4699] "t" "t" "a" "c" "c" "a" "g" "a" "c" "g" "c" "g" "a" "g" "c" "c" "c" "c"
##  [4717] "g" "g" "a" "g" "g" "c" "t" "a" "c" "c" "g" "c" "t" "g" "t" "a" "a" "g"
##  [4735] "a" "a" "c" "g" "c" "g" "a" "g" "g" "t" "t" "g" "a" "a" "g" "c" "g" "g"
##  [4753] "t" "c" "c" "a" "a" "g" "g" "g" "g" "a" "g" "a" "c" "g" "t" "g" "t" "c"
##  [4771] "c" "g" "c" "g" "a" "g" "a" "c" "a" "c" "t" "c" "t" "g" "a" "t" "c" "g"
##  [4789] "t" "a" "g" "t" "g" "g" "c" "a" "c" "g" "c" "c" "g" "g" "g" "t" "g" "t"
##  [4807] "g" "c" "a" "a" "g" "a" "g" "g" "t" "t" "g" "t" "c" "c" "g" "c" "g" "c"
##  [4825] "t" "c" "c" "t" "a" "a" "c" "a" "t" "g" "a" "c" "t" "c" "c" "g" "g" "g"
##  [4843] "c" "a" "c" "t" "c" "c" "t" "g" "t" "a" "c" "c" "c" "g" "g" "a" "a" "g"
##  [4861] "c" "c" "a" "c" "a" "g" "c" "t" "t" "a" "c" "c" "a" "g" "g" "c" "a" "g"
##  [4879] "c" "c" "c" "g" "a" "a" "g" "g" "c" "a" "c" "c" "c" "c" "g" "c" "c" "a"
##  [4897] "a" "t" "a" "g" "g" "c" "a" "c" "c" "a" "t" "t" "a" "c" "c" "t" "g" "g"
##  [4915] "t" "c" "t" "c" "g" "g" "c" "g" "t" "a" "a" "a" "c" "a" "a" "t" "a" "g"
##  [4933] "g" "c" "c" "g" "g" "t" "a" "t" "c" "g" "a" "a" "a" "t" "g" "a" "t" "c"
##  [4951] "a" "a" "t" "t" "t" "a" "c" "a" "g" "c" "c" "c" "t" "g" "c" "a" "c" "g"
##  [4969] "c" "c" "c" "t" "a" "g" "g" "c" "c" "g" "a" "g" "t" "t" "g" "g" "a" "c"
##  [4987] "a" "a" "a" "g" "g" "c" "c" "g" "t" "a" "c" "t" "c" "a" "g" "g" "g" "a"
##  [5005] "g" "g" "g" "t" "c" "t" "c" "a" "t" "c" "a" "c" "c" "a" "c" "c" "t" "t"
##  [5023] "g" "t" "a" "c" "a" "a" "t" "a" "c" "a" "t" "g" "c" "a" "g" "a" "c" "c"
##  [5041] "g" "a" "g" "a" "t" "c" "c" "t" "g" "c" "g" "a" "t" "c" "g" "t" "a" "c"
##  [5059] "g" "g" "a" "c" "g" "a" "c" "c" "g" "c" "t" "c" "g" "g" "g" "c" "a" "g"
##  [5077] "a" "t" "t" "a" "c" "a" "c" "g" "c" "g" "a" "g" "g" "g" "g" "c" "a" "a"
##  [5095] "g" "c" "c" "g" "t" "g" "c" "c" "g" "t" "c" "c" "g" "a" "g" "t" "a" "c"
##  [5113] "t" "a" "c" "c" "g" "a" "g" "c" "a" "c" "c" "c" "a" "c" "g" "c" "c" "g"
##  [5131] "c" "g" "c" "c" "c" "g" "a" "a" "c" "t" "g" "g" "g" "g" "a" "c" "c" "c"
##  [5149] "t" "g" "g" "a" "t" "a" "g" "g" "t" "c" "a" "a" "g" "g" "t" "t" "t" "c"
##  [5167] "g" "c" "a" "c" "c" "a" "c" "a" "g" "c" "a" "c" "c" "c" "c" "c" "t" "g"
##  [5185] "c" "t" "g" "a" "t" "a" "a" "g" "a" "g" "a" "t" "g" "g" "t" "t" "a" "c"
##  [5203] "a" "a" "c" "t" "c" "a" "g" "c" "g" "c" "c" "g" "a" "t" "a" "a" "a" "c"
##  [5221] "a" "t" "a" "a" "t" "c" "t" "c" "g" "c" "c" "c" "a" "t" "t" "t" "a" "t"
##  [5239] "c" "c" "a" "a" "a" "g" "t" "c" "t" "a" "c" "c" "g" "g" "g" "t" "g" "g"
##  [5257] "c" "g" "c" "g" "a" "a" "a" "c" "a" "a" "g" "t" "g" "c" "c" "g" "g" "a"
##  [5275] "a" "g" "c" "g" "g" "a" "t" "c" "c" "c" "a" "t" "g" "c" "a" "c" "t" "g"
##  [5293] "c" "g" "g" "a" "t" "g" "a" "g" "a" "c" "a" "g" "g" "a" "t" "c" "t" "t"
##  [5311] "a" "g" "g" "c" "t" "c" "c" "a" "g" "t" "c" "t" "t" "t" "t" "a" "g" "a"
##  [5329] "g" "g" "g" "c" "a" "t" "g" "t" "t" "g" "g" "g" "t" "g" "c" "g" "c" "t"
##  [5347] "c" "g" "c" "t" "g" "g" "a" "c" "a" "c" "c" "c" "c" "a" "t" "g" "c" "c"
##  [5365] "c" "a" "c" "g" "t" "c" "a" "c" "c" "c" "a" "a" "c" "a" "g" "a" "c" "g"
##  [5383] "a" "g" "t" "g" "a" "g" "c" "t" "c" "c" "t" "c" "g" "g" "g" "c" "c" "t"
##  [5401] "a" "g" "g" "c" "c" "c" "c" "g" "a" "g" "g" "t" "a" "g" "t" "c" "c" "a"
##  [5419] "a" "a" "t" "t" "a" "c" "g" "g" "g" "g" "g" "c" "a" "a" "a" "g" "g" "c"
##  [5437] "a" "g" "a" "g" "g" "t" "g" "a" "g" "a" "g" "a" "g" "c" "t" "a" "a" "a"
##  [5455] "a" "g" "c" "c" "g" "g" "g" "t" "t" "g" "g" "g" "c" "g" "g" "a" "a" "a"
##  [5473] "a" "g" "g" "a" "c" "a" "a" "g" "a" "g" "c" "c" "g" "c" "c" "c" "g" "t"
##  [5491] "c" "c" "a" "g" "a" "a" "c" "t" "a" "a" "c" "g" "g" "c" "a" "t" "t" "a"
##  [5509] "c" "g" "t" "g" "a" "g" "c" "t" "a" "c" "g" "g" "t" "a" "g" "t" "g" "c"
##  [5527] "c" "g" "g" "t" "a" "t" "a" "g" "c" "t" "c" "a" "t" "g" "g" "t" "t" "t"
##  [5545] "c" "g" "a" "c" "a" "a" "g" "g" "a" "g" "c" "c" "g" "g" "g" "c" "g" "c"
##  [5563] "c" "a" "a" "c" "c" "a" "t" "g" "a" "g" "g" "t" "g" "a" "g" "g" "g" "c"
##  [5581] "g" "t" "t" "a" "g" "c" "a" "c" "a" "c" "c" "c" "a" "c" "g" "c" "t" "a"
##  [5599] "c" "a" "c" "c" "a" "g" "a" "a" "c" "g" "t" "t" "t" "g" "a" "c" "c" "g"
##  [5617] "a" "a" "c" "t" "g" "g" "t" "c" "t" "c" "a" "c" "a" "g" "a" "t" "c" "g"
##  [5635] "c" "a" "a" "g" "a" "c" "t" "c" "g" "g" "c" "c" "g" "g" "g" "t" "a" "g"
##  [5653] "g" "g" "g" "c" "t" "c" "t" "c" "g" "c" "a" "t" "a" "t" "g" "g" "t" "g"
##  [5671] "a" "t" "c" "g" "g" "t" "c" "a" "g" "a" "t" "g" "a" "t" "t" "a" "a" "t"
##  [5689] "c" "t" "g" "g" "a" "g" "g" "c" "a" "a" "a" "c" "g" "t" "c" "g" "g" "c"
##  [5707] "c" "t" "t" "c" "c" "c" "g" "c" "t" "t" "g" "a" "g" "c" "a" "c" "g" "g"
##  [5725] "c" "a" "g" "c" "c" "c" "t" "a" "g" "g" "a" "g" "a" "c" "a" "g" "g" "c"
##  [5743] "c" "a" "c" "a" "a" "g" "c" "c" "g" "c" "g" "g" "t" "a" "g" "g" "t" "g"
##  [5761] "t" "t" "c" "a" "c" "g" "a" "c" "c" "a" "t" "c" "a" "c" "c" "a" "c" "g"
##  [5779] "c" "g" "g" "g" "c" "c" "c" "a" "c" "a" "a" "t" "t" "c" "c" "t" "g" "t"
##  [5797] "g" "c" "c" "g" "t" "c" "g" "g" "g" "c" "c" "g" "c" "a" "g" "t" "c" "a"
##  [5815] "t" "c" "c" "a" "c" "t" "t" "a" "c" "g" "c" "a" "t" "g" "g" "a" "g" "c"
##  [5833] "c" "a" "a" "t" "a" "g" "g" "g" "c" "c" "t" "a" "g" "g" "a" "g" "t" "g"
##  [5851] "g" "a" "a" "a" "c" "t" "a" "g" "c" "a" "t" "g" "c" "t" "g" "c" "a" "t"
##  [5869] "c" "g" "g" "c" "c" "g" "c" "a" "t" "t" "g" "c" "a" "t" "c" "t" "c" "c"
##  [5887] "g" "c" "c" "t" "a" "a" "t" "g" "a" "c" "c" "g" "c" "t" "c" "c" "a" "t"
##  [5905] "g" "c" "g" "c" "a" "a" "t" "g" "c" "t" "a" "g" "g" "t" "t" "g" "t" "c"
##  [5923] "c" "c" "a" "g" "a" "c" "a" "g" "a" "t" "a" "c" "g" "t" "t" "a" "c" "t"
##  [5941] "c" "a" "g" "a" "g" "a" "c" "a" "t" "g" "g" "c" "g" "a" "g" "a" "c" "g"
##  [5959] "a" "g" "a" "g" "c" "g" "a" "a" "t" "g" "a" "c" "t" "c" "a" "c" "a" "c"
##  [5977] "t" "t" "g" "g" "t" "g" "a" "a" "c" "c" "g" "g" "t" "t" "c" "c" "c" "c"
##  [5995] "c" "c" "a" "g" "c" "t" "g" "c" "c" "t" "a" "a" "t" "c" "a" "a" "g" "a"
##  [6013] "c" "t" "c" "t" "a" "g" "c" "a" "a" "g" "g" "c" "t" "t" "a" "g" "c" "a"
##  [6031] "a" "c" "c" "a" "g" "a" "a" "g" "a" "g" "c" "g" "c" "c" "g" "c" "a" "t"
##  [6049] "t" "g" "t" "c" "a" "t" "c" "c" "g" "c" "t" "t" "c" "a" "t" "c" "a" "c"
##  [6067] "c" "g" "t" "c" "c" "g" "t" "a" "g" "a" "c" "g" "c" "g" "c" "a" "c" "a"
##  [6085] "c" "g" "a" "a" "t" "t" "g" "g" "c" "t" "a" "t" "g" "g" "a" "g" "a" "c"
##  [6103] "c" "c" "t" "c" "a" "c" "t" "c" "t" "g" "a" "g" "g" "a" "t" "a" "c" "a"
##  [6121] "a" "c" "c" "c" "a" "g" "a" "c" "g" "t" "c" "a" "a" "c" "g" "g" "a" "c"
##  [6139] "a" "a" "c" "c" "c" "g" "c" "g" "a" "c" "a" "c" "a" "c" "g" "c" "t" "c"
##  [6157] "t" "c" "c" "g" "c" "c" "a" "t" "c" "c" "g" "a" "t" "c" "g" "g" "t" "a"
##  [6175] "c" "g" "c" "c" "c" "a" "c" "c" "c" "g" "g" "g" "c" "c" "c" "t" "c" "c"
##  [6193] "a" "g" "g" "a" "c" "c" "a" "c" "c" "c" "a" "t" "c" "c" "g" "c" "c" "c"
##  [6211] "g" "t" "a" "a" "c" "c" "g" "g" "g" "a" "t" "a" "t" "g" "t" "c" "g" "a"
##  [6229] "g" "g" "c" "a" "c" "g" "c" "g" "t" "g" "t" "c" "c" "a" "g" "c" "t" "a"
##  [6247] "a" "g" "c" "c" "t" "c" "t" "g" "c" "g" "c" "c" "t" "a" "t" "c" "t" "g"
##  [6265] "c" "c" "t" "g" "a" "g" "c" "t" "t" "g" "t" "a" "a" "g" "g" "a" "g" "t"
##  [6283] "a" "c" "c" "t" "g" "a" "t" "g" "c" "c" "g" "g" "a" "g" "g" "c" "c" "t"
##  [6301] "a" "c" "t" "g" "c" "c" "a" "a" "t" "c" "g" "t" "t" "a" "t" "c" "c" "a"
##  [6319] "g" "t" "c" "c" "c" "g" "c" "c" "t" "c" "g" "a" "g" "g" "a" "a" "g" "t"
##  [6337] "c" "g" "c" "a" "a" "c" "a" "c" "t" "c" "c" "c" "t" "c" "a" "a" "g" "a"
##  [6355] "g" "a" "t" "a" "t" "a" "t" "a" "g" "t" "a" "g" "c" "c" "g" "c" "c" "c"
##  [6373] "g" "t" "g" "c" "g" "g" "t" "t" "t" "g" "g" "c" "g" "a" "c" "c" "g" "g"
##  [6391] "c" "c" "c" "a" "g" "c" "t" "a" "a" "g" "t" "a" "t" "a" "c" "a" "g" "t"
##  [6409] "g" "a" "t" "a" "g" "c" "t" "a" "a" "g" "t" "g" "a" "g" "g" "c" "c" "t"
##  [6427] "a" "t" "g" "c" "c" "c" "c" "c" "a" "a" "c" "g" "t" "c" "g" "g" "a" "g"
##  [6445] "g" "t" "c" "t" "a" "t" "a" "c" "a" "t" "c" "a" "c" "g" "g" "c" "g" "t"
##  [6463] "a" "c" "g" "c" "g" "c" "g" "g" "t" "c" "g" "g" "t" "t" "c" "c" "a" "a"
##  [6481] "a" "g" "t" "g" "g" "c" "g" "t" "a" "c" "c" "a" "a" "g" "g" "g" "c" "t"
##  [6499] "g" "c" "c" "g" "g" "c" "a" "a" "g" "g" "t" "c" "a" "g" "c" "c" "a" "g"
##  [6517] "g" "g" "t" "t" "c" "c" "g" "a" "a" "a" "t" "a" "g" "a" "g" "c" "g" "g"
##  [6535] "c" "c" "t" "g" "g" "g" "c" "g" "a" "t" "g" "g" "g" "a" "g" "t" "t" "g"
##  [6553] "c" "a" "t" "g" "a" "t" "g" "c" "g" "c" "a" "a" "t" "c" "a" "c" "g" "t"
##  [6571] "c" "t" "g" "g" "t" "t" "g" "g" "g" "a" "t" "a" "a" "c" "c" "a" "c" "g"
##  [6589] "a" "t" "g" "g" "g" "c" "t" "t" "t" "a" "c" "c" "g" "c" "c" "t" "a" "a"
##  [6607] "t" "g" "g" "c" "a" "g" "c" "a" "c" "c" "a" "c" "g" "c" "c" "c" "g" "a"
##  [6625] "c" "t" "c" "t" "t" "a" "a" "c" "c" "g" "g" "g" "a" "g" "t" "a" "a" "t"
##  [6643] "a" "t" "c" "t" "g" "g" "c" "c" "a" "c" "t" "a" "g" "a" "a" "t" "t" "c"
##  [6661] "t" "t" "c" "g" "t" "g" "t" "t" "a" "a" "t" "c" "a" "c" "t" "t" "c" "c"
##  [6679] "a" "c" "c" "a" "g" "g" "c" "a" "c" "g" "g" "c" "a" "c" "c" "g" "g" "c"
##  [6697] "a" "a" "g" "t" "a" "a" "c" "a" "t" "t" "a" "c" "g" "g" "t" "t" "g" "a"
##  [6715] "c" "t" "g" "a" "g" "g" "a" "c" "g" "g" "g" "g" "g" "t" "g" "g" "t" "c"
##  [6733] "g" "g" "t" "t" "c" "g" "g" "t" "t" "a" "t" "c" "a" "c" "g" "t" "c" "a"
##  [6751] "c" "a" "a" "a" "a" "c" "a" "t" "a" "t" "t" "c" "t" "g" "g" "t" "c" "g"
##  [6769] "t" "t" "g" "c" "g" "a" "c" "a" "t" "a" "c" "t" "t" "a" "c" "t" "g" "a"
##  [6787] "t" "g" "g" "c" "c" "a" "t" "c" "a" "t" "a" "c" "c" "c" "g" "a" "a" "c"
##  [6805] "g" "c" "g" "t" "c" "c" "a" "g" "t" "c" "a" "g" "a" "g" "g" "c" "g" "g"
##  [6823] "c" "a" "c" "t" "g" "c" "a" "g" "a" "c" "g" "a" "c" "a" "a" "g" "g" "t"
##  [6841] "a" "c" "t" "t" "c" "a" "c" "c" "c" "c" "g" "g" "g" "c" "t" "t" "c" "g"
##  [6859] "c" "t" "c" "t" "c" "c" "a" "g" "c" "c" "t" "c" "g" "g" "t" "t" "c" "a"
##  [6877] "g" "c" "a" "c" "g" "g" "c" "a" "c" "a" "c" "c" "g" "c" "g" "a" "g" "g"
##  [6895] "a" "g" "t" "c" "c" "a" "a" "g" "g" "c" "g" "c" "a" "g" "g" "g" "c" "t"
##  [6913] "c" "t" "c" "g" "t" "g" "t" "c" "c" "c" "g" "t" "c" "c" "a" "g" "a" "c"
##  [6931] "g" "a" "g" "t" "t" "g" "a" "c" "t" "c" "a" "a" "g" "a" "a" "a" "t" "g"
##  [6949] "g" "g" "a" "t" "c" "a" "g" "g" "c" "a" "a" "t" "t" "a" "c" "a" "t" "c"
##  [6967] "a" "t" "a" "t" "a" "c" "t" "g" "g" "c" "a" "c" "g" "c" "c" "g" "g" "a"
##  [6985] "g" "t" "t" "c" "c" "t" "c" "c" "t" "c" "c" "a" "g" "a" "g" "a" "g" "t"
##  [7003] "c" "g" "c" "c" "t" "t" "t" "c" "g" "t" "g" "a" "c" "t" "a" "a" "c" "a"
##  [7021] "g" "c" "g" "g" "t" "t" "a" "g" "g" "t" "c" "a" "c" "a" "t" "g" "t" "c"
##  [7039] "a" "c" "t" "c" "a" "t" "c" "a" "c" "t" "a" "t" "a" "g" "g" "c" "c" "c"
##  [7057] "t" "c" "c" "c" "c" "g" "g" "g" "t" "t" "g" "a" "g" "a" "t" "a" "a" "g"
##  [7075] "g" "t" "t" "c" "g" "t" "g" "g" "a" "a" "a" "t" "g" "g" "c" "g" "g" "a"
##  [7093] "c" "c" "t" "c" "a" "a" "a" "g" "a" "t" "g" "c" "g" "c" "a" "a" "t" "a"
##  [7111] "c" "t" "g" "g" "c" "c" "c" "t" "a" "t" "g" "g" "c" "g" "a" "a" "a" "c"
##  [7129] "c" "c" "c" "c" "a" "a" "a" "a" "g" "t" "t" "g" "c" "c" "a" "c" "a" "t"
##  [7147] "c" "g" "g" "t" "t" "c" "a" "g" "g" "c" "c" "a" "g" "t" "t" "a" "c" "t"
##  [7165] "g" "a" "g" "c" "t" "g" "g" "g" "c" "t" "c" "g" "g" "g" "c" "a" "g" "a"
##  [7183] "a" "c" "g" "c" "a" "g" "g" "a" "c" "c" "g" "a" "c" "c" "g" "c" "c" "a"
##  [7201] "c" "c" "a" "g" "t" "c" "c" "g" "c" "c" "c" "c" "a" "t" "a" "c" "t" "t"
##  [7219] "g" "g" "c" "c" "c" "a" "g" "t" "g" "a" "a" "a" "g" "g" "g" "g" "a" "g"
##  [7237] "a" "t" "g" "c" "g" "g" "a" "a" "a" "t" "a" "t" "c" "a" "a" "c" "g" "g"
##  [7255] "c" "c" "a" "a" "t" "g" "c" "t" "a" "c" "c" "a" "g" "g" "c" "c" "g" "a"
##  [7273] "c" "t" "c" "a" "a" "g" "c" "a" "g" "g" "a" "t" "a" "g" "a" "g" "g" "t"
##  [7291] "c" "c" "c" "c" "c" "c" "c" "c" "c" "c" "a" "c" "t" "a" "t" "g" "g" "g"
##  [7309] "c" "c" "a" "g" "a" "a" "g" "c" "c" "g" "g" "c" "t" "g" "c" "t" "t" "g"
##  [7327] "g" "c" "t" "c" "g" "g" "a" "g" "a" "g" "a" "a" "g" "t" "c" "a" "a" "c"
##  [7345] "t" "g" "c" "g" "t" "c" "t" "g" "a" "c" "t" "a" "t" "c" "g" "t" "g" "a"
##  [7363] "a" "g" "a" "a" "g" "c" "c" "a" "c" "a" "c" "t" "a" "c" "g" "c" "a" "c"
##  [7381] "g" "a" "g" "g" "a" "t" "t" "c" "t" "g" "t" "t" "t" "g" "c" "a" "c" "a"
##  [7399] "c" "g" "c" "g" "g" "a" "t" "c" "a" "g" "g" "g" "g" "t" "c" "g" "c" "c"
##  [7417] "g" "a" "t" "g" "a" "g" "c" "t" "c" "c" "a" "a" "c" "c" "c" "t" "g" "g"
##  [7435] "a" "t" "t" "a" "a" "g" "t" "c" "c" "c" "g" "g" "t" "a" "c" "a" "t" "a"
##  [7453] "a" "g" "a" "a" "c" "c" "a" "t" "t" "t" "c" "a" "t" "c" "g" "a" "a" "t"
##  [7471] "c" "g" "c" "c" "a" "c" "g" "g" "c" "a" "c" "g" "a" "t" "t" "g" "t" "g"
##  [7489] "c" "a" "g" "g" "a" "a" "c" "c" "c" "t" "g" "t" "g" "g" "c" "a" "a" "a"
##  [7507] "c" "t" "g" "g" "c" "a" "c" "a" "a" "g" "c" "c" "g" "g" "t" "c" "c" "c"
##  [7525] "a" "g" "g" "t" "a" "a" "t" "c" "g" "c" "t" "g" "c" "c" "c" "t" "t" "c"
##  [7543] "c" "c" "g" "g" "g" "g" "g" "c" "t" "c" "a" "t" "g" "g" "a" "t" "a" "a"
##  [7561] "t" "c" "g" "a" "a" "g" "t" "c" "a" "c" "c" "c" "c" "a" "a" "t" "g" "a"
##  [7579] "t" "c" "g" "c" "a" "a" "a" "t" "c" "a" "g" "g" "t" "a" "c" "a" "g" "c"
##  [7597] "t" "g" "g" "a" "t" "t" "a" "c" "g" "a" "t" "t" "c" "t" "t" "a" "c" "a"
##  [7615] "c" "c" "a" "g" "c" "t" "a" "c" "t" "g" "a" "c" "g" "c" "g" "a" "g" "c"
##  [7633] "g" "t" "t" "c" "t" "a" "a" "c" "a" "a" "a" "g" "a" "g" "t" "c" "c" "t"
##  [7651] "a" "c" "g" "t" "c" "t" "t" "c" "c" "c" "c" "t" "a" "c" "c" "t" "a" "g"
##  [7669] "c" "a" "a" "c" "c" "a" "c" "a" "g" "g" "g" "c" "t" "a" "c" "g" "g" "c"
##  [7687] "t" "a" "t" "c" "c" "c" "a" "t" "c" "c" "a" "g" "a" "t" "g" "c" "t" "c"
##  [7705] "t" "a" "c" "t" "g" "a" "g" "t" "a" "g" "t" "c" "c" "a" "a" "g" "c" "c"
##  [7723] "a" "g" "a" "c" "t" "g" "a" "g" "g" "t" "c" "a" "c" "t" "a" "a" "a" "c"
##  [7741] "g" "a" "g" "c" "a" "a" "c" "c" "c" "a" "c" "t" "c" "c" "g" "g" "c" "t"
##  [7759] "g" "g" "t" "c" "g" "a" "g" "t" "g" "a" "c" "c" "a" "t" "c" "g" "c" "a"
##  [7777] "c" "a" "g" "t" "g" "a" "a" "g" "g" "c" "t" "c" "c" "a" "a" "c" "a" "a"
##  [7795] "c" "c" "c" "g" "a" "g" "c" "g" "c" "a" "g" "c" "c" "a" "c" "a" "c" "c"
##  [7813] "t" "c" "c" "t" "t" "g" "c" "a" "c" "g" "t" "a" "c" "g" "a" "a" "t" "a"
##  [7831] "g" "a" "a" "c" "c" "t" "a" "g" "c" "c" "a" "g" "g" "t" "a" "a" "g" "c"
##  [7849] "a" "c" "a" "a" "t" "g" "t" "c" "g" "a" "t" "c" "g" "g" "g" "a" "c" "c"
##  [7867] "a" "c" "c" "a" "t" "g" "a" "g" "g" "a" "c" "t" "c" "t" "g" "g" "a" "g"
##  [7885] "c" "c" "a" "c" "c" "a" "a" "a" "a" "c" "c" "g" "a" "a" "t" "g" "c" "c"
##  [7903] "a" "a" "c" "c" "g" "a" "t" "g" "c" "a" "a" "a" "c" "a" "t" "c" "t" "a"
##  [7921] "g" "t" "g" "g" "a" "t" "c" "a" "c" "c" "g" "c" "c" "c" "g" "c" "g" "a"
##  [7939] "t" "c" "c" "g" "c" "c" "c" "g" "c" "a" "g" "t" "g" "c" "t" "g" "g" "g"
##  [7957] "t" "t" "a" "g" "a" "t" "a" "a" "t" "g" "t" "g" "g" "a" "g" "c" "c" "g"
##  [7975] "a" "c" "g" "a" "t" "t" "c" "g" "g" "c" "t" "g" "g" "g" "g" "t" "a" "c"
##  [7993] "a" "c" "a" "c" "g" "c" "c" "c" "g" "c" "c" "c" "c" "g" "t" "c" "a" "a"
##  [8011] "c" "t" "a" "g" "a" "t" "t" "g" "a" "g" "c" "c" "g" "a" "c" "t" "c" "c"
##  [8029] "a" "a" "t" "g" "c" "c" "t" "t" "c" "g" "a" "c" "a" "c" "g" "a" "a" "t"
##  [8047] "g" "a" "t" "c" "g" "t" "c" "t" "a" "g" "a" "t" "t" "a" "t" "a" "t" "g"
##  [8065] "c" "a" "t" "a" "c" "g" "c" "a" "a" "t" "c" "t" "g" "t" "g" "t" "a" "c"
##  [8083] "g" "g" "c" "t" "a" "c" "t" "c" "g" "g" "a" "c" "g" "a" "a" "a" "a" "c"
##  [8101] "a" "a" "a" "a" "g" "t" "a" "g" "c" "c" "g" "t" "a" "t" "t" "g" "g" "g"
##  [8119] "a" "c" "a" "c" "c" "c" "g" "c" "g" "c" "a" "a" "c" "g" "a" "c" "g" "a"
##  [8137] "a" "g" "c" "t" "g" "a" "t" "g" "a" "g" "g" "t" "c" "a" "g" "t" "a" "a"
##  [8155] "a" "g" "g" "a" "c" "a" "c" "c" "t" "a" "c" "a" "c" "t" "g" "c" "t" "g"
##  [8173] "g" "t" "g" "c" "c" "a" "t" "t" "c" "g" "g" "t" "t" "t" "c" "t" "t" "a"
##  [8191] "c" "a" "g" "g" "a" "g" "c" "a" "c" "g" "a" "a" "c" "a" "c" "a" "c" "a"
##  [8209] "t" "c" "a" "t" "t" "a" "c" "a" "c" "t" "a" "c" "g" "t" "t" "c" "a" "t"
##  [8227] "t" "c" "a" "c" "g" "a" "c" "c" "g" "a" "g" "a" "a" "a" "g" "t" "a" "g"
##  [8245] "c" "t" "a" "a" "t" "t" "a" "c" "g" "c" "c" "t" "g" "g" "g" "a" "a" "c"
##  [8263] "t" "t" "c" "t" "c" "a" "a" "t" "c" "c" "t" "a" "a" "a" "g" "c" "a" "g"
##  [8281] "t" "a" "t" "t" "g" "g" "a" "c" "g" "c" "a" "g" "a" "t" "a" "c" "c" "g"
##  [8299] "t" "c" "g" "c" "g" "c" "g" "c" "g" "c" "g" "t" "g" "c" "g" "g" "a" "t"
##  [8317] "c" "t" "c" "g" "t" "c" "c" "c" "t" "g" "g" "t" "a" "c" "g" "c" "c" "a"
##  [8335] "g" "a" "c" "c" "g" "g" "t" "a" "t" "g" "a" "c" "t" "a" "c" "g" "c" "t"
##  [8353] "c" "g" "a" "t" "c" "g" "c" "a" "c" "t" "c" "c" "a" "g" "c" "c" "t" "c"
##  [8371] "c" "a" "g" "c" "t" "g" "a" "g" "g" "t" "g" "c" "c" "c" "c" "t" "t" "t"
##  [8389] "a" "g" "t" "t" "a" "c" "a" "t" "c" "a" "a" "c" "a" "c" "c" "t" "c" "g"
##  [8407] "c" "g" "g" "a" "g" "g" "t" "g" "c" "a" "c" "g" "g" "g" "a" "a" "g" "c"
##  [8425] "t" "t" "c" "t" "c" "a" "g" "c" "g" "c" "g" "t" "t" "a" "c" "c" "t" "g"
##  [8443] "g" "c" "a" "a" "c" "g" "a" "g" "g" "c" "c" "g" "t" "a" "g" "a" "g" "g"
##  [8461] "t" "t" "c" "c" "c" "t" "c" "g" "g" "a" "a" "c" "c" "t" "t" "t" "g" "g"
##  [8479] "a" "c" "t" "g" "a" "a" "g" "g" "c" "g" "c" "c" "t" "t" "g" "c" "c" "c"
##  [8497] "t" "c" "g" "a" "c" "t" "t" "c" "a" "c" "c" "c" "g" "t" "c" "c" "a" "t"
##  [8515] "c" "t" "a" "g" "a" "g" "c" "a" "t" "t" "c" "g" "g" "c" "c" "c" "t" "c"
##  [8533] "a" "a" "c" "g" "g" "c" "g" "a" "t" "g" "c" "a" "a" "a" "g" "a" "a" "c"
##  [8551] "a" "t" "t" "a" "a" "c" "g" "g" "t" "a" "c" "c" "c" "t" "g" "t" "c" "c"
##  [8569] "c" "c" "g" "c" "a" "g" "g" "c" "g" "a" "a" "a" "a" "t" "c" "c" "g" "g"
##  [8587] "g" "t" "c" "c" "g" "c" "g" "g" "c" "t" "g" "g" "g" "a" "c" "a" "c" "a"
##  [8605] "c" "a" "c" "g" "g" "g" "g" "t" "a" "g" "t" "g" "a" "t" "t" "t" "c" "c"
##  [8623] "a" "g" "g" "g" "t" "c" "a" "a" "a" "t" "g" "t" "t" "g" "g" "t" "t" "c"
##  [8641] "a" "c" "g" "a" "g" "a" "g" "a" "g" "c" "c" "g" "c" "c" "g" "t" "c" "c"
##  [8659] "t" "c" "a" "a" "g" "a" "c" "g" "c" "c" "a" "a" "a" "a" "a" "t" "a" "c"
##  [8677] "t" "c" "g" "t" "g" "t" "a" "g" "t" "c" "a" "a" "g" "t" "g" "t" "c" "g"
##  [8695] "t" "g" "c" "a" "c" "g" "c" "c" "a" "g" "t" "c" "g" "a" "t" "t" "c" "c"
##  [8713] "g" "t" "g" "c" "g" "c" "t" "t" "t" "c" "t" "c" "t" "c" "g" "a" "g" "a"
##  [8731] "a" "c" "c" "t" "t" "a" "t" "t" "c" "c" "a" "g" "c" "c" "c" "c" "g" "g"
##  [8749] "g" "t" "t" "t" "c" "c" "t" "g" "g" "c" "g" "g" "c" "t" "t" "c" "c" "g"
##  [8767] "c" "a" "g" "c" "t" "a" "a" "a" "g" "a" "g" "t" "g" "c" "a" "c" "a" "c"
##  [8785] "g" "a" "g" "g" "c" "c" "g" "c" "g" "c" "g" "a" "t" "a" "c" "t" "t" "c"
##  [8803] "a" "g" "t" "t" "a" "g" "a" "g" "a" "a" "t" "a" "c" "t" "c" "g" "t" "c"
##  [8821] "c" "c" "t" "t" "a" "t" "a" "g" "t" "t" "g" "g" "c" "g" "g" "t" "a" "c"
##  [8839] "g" "c" "a" "a" "g" "c" "t" "g" "c" "g" "t" "g" "c" "a" "c" "c" "a" "g"
##  [8857] "a" "g" "t" "a" "a" "t" "c" "c" "a" "g" "t" "g" "a" "t" "t" "t" "g" "c"
##  [8875] "g" "c" "c" "g" "a" "a" "c" "c" "t" "a" "g" "c" "a" "a" "a" "t" "t" "c"
##  [8893] "c" "a" "g" "g" "g" "a" "a" "c" "g" "c" "c" "a" "t" "a" "a" "c" "c" "c"
##  [8911] "a" "g" "c" "a" "c" "g" "g" "t" "t" "t" "t" "g" "g" "a" "c" "c" "g" "c"
##  [8929] "g" "c" "c" "t" "t" "g" "a" "t" "c" "c" "g" "a" "g" "a" "g" "t" "g" "g"
##  [8947] "c" "a" "c" "c" "g" "g" "g" "t" "t" "g" "t" "g" "g" "c" "c" "g" "c" "a"
##  [8965] "a" "c" "t" "g" "t" "g" "t" "g" "g" "a" "t" "g" "c" "c" "a" "c" "c" "a"
##  [8983] "g" "g" "t" "g" "t" "c" "c" "g" "a" "a" "c" "a" "c" "t" "c" "c" "t" "t"
##  [9001] "c" "c" "a" "a" "t" "a" "c" "g" "t" "a" "g" "c" "c" "a" "a" "g" "c" "c"
##  [9019] "g" "a" "g" "t" "c" "c" "c" "c" "c" "g" "c" "a" "g" "a" "g" "a" "a" "g"
##  [9037] "c" "t" "c" "c" "c" "g" "t" "g" "c" "t" "a" "c" "a" "c" "c" "t" "c" "g"
##  [9055] "g" "g" "a" "c" "c" "t" "g" "g" "c" "g" "t" "a" "c" "a" "g" "a" "c" "c"
##  [9073] "g" "c" "t" "t" "c" "t" "c" "t" "g" "g" "t" "t" "g" "g" "g" "g" "t" "c"
##  [9091] "c" "c" "a" "t" "c" "c" "a" "a" "a" "a" "c" "g" "a" "c" "t" "c" "g" "t"
##  [9109] "g" "g" "c" "c" "c" "a" "g" "g" "c" "g" "a" "c" "t" "g" "a" "a" "a" "a"
##  [9127] "t" "t" "c" "g" "c" "a" "a" "g" "c" "g" "a" "t" "a" "c" "t" "a" "a" "a"
##  [9145] "t" "a" "a" "a" "g" "g" "t" "t" "t" "c" "t" "a" "a" "g" "t" "c" "t" "a"
##  [9163] "g" "a" "a" "t" "a" "c" "c" "g" "g" "g" "c" "g" "t" "a" "c" "g" "c" "a"
##  [9181] "c" "c" "t" "g" "a" "a" "g" "t" "a" "a" "a" "t" "c" "t" "t" "t" "c" "g"
##  [9199] "g" "c" "a" "t" "a" "g" "a" "t" "c" "c" "t" "a" "g" "g" "c" "g" "g" "t"
##  [9217] "g" "c" "g" "t" "c" "c" "a" "t" "a" "g" "c" "c" "c" "g" "c" "a" "c" "t"
##  [9235] "a" "a" "g" "g" "c" "t" "a" "t" "g" "a" "c" "t" "a" "c" "t" "t" "g" "a"
##  [9253] "a" "a" "a" "c" "a" "g" "g" "a" "t" "c" "c" "g" "a" "c" "t" "g" "c" "g"
##  [9271] "t" "c" "c" "t" "a" "a" "a" "a" "g" "t" "a" "g" "c" "a" "a" "g" "c" "a"
##  [9289] "c" "c" "c" "g" "a" "t" "g" "a" "c" "a" "g" "a" "g" "c" "t" "g" "t" "g"
##  [9307] "g" "t" "g" "a" "g" "t" "a" "c" "t" "g" "c" "a" "t" "t" "a" "c" "t" "a"
##  [9325] "g" "t" "t" "t" "c" "t" "t" "g" "a" "a" "g" "a" "g" "a" "g" "c" "t" "c"
##  [9343] "a" "c" "c" "t" "g" "a" "c" "a" "g" "g" "c" "t" "c" "g" "c" "a" "a" "t"
##  [9361] "c" "a" "c" "t" "a" "g" "g" "c" "g" "c" "a" "c" "g" "g" "c" "t" "g" "g"
##  [9379] "g" "c" "g" "g" "a" "t" "g" "g" "t" "g" "a" "c" "c" "a" "g" "c" "c" "c"
##  [9397] "t" "g" "a" "g" "c" "a" "a" "a" "a" "t" "g" "a" "t" "t" "c" "g" "c" "t"
##  [9415] "a" "c" "t" "g" "g" "a" "a" "c" "t" "t" "a" "g" "a" "g" "t" "c" "t" "c"
##  [9433] "c" "t" "a" "a" "g" "c" "c" "a" "c" "g" "c" "c" "g" "c" "c" "g" "a" "g"
##  [9451] "g" "c" "g" "t" "t" "c" "a" "t" "g" "c" "a" "c" "g" "g" "c" "a" "c" "t"
##  [9469] "t" "c" "c" "t" "a" "a" "c" "a" "t" "a" "a" "a" "c" "g" "g" "c" "c" "c"
##  [9487] "c" "a" "c" "g" "c" "a" "g" "c" "c" "g" "a" "t" "a" "c" "c" "a" "c" "c"
##  [9505] "a" "g" "t" "g" "c" "a" "a" "a" "g" "a" "c" "t" "g" "a" "g" "t" "c" "c"
##  [9523] "a" "t" "g" "g" "a" "g" "c" "a" "a" "t" "c" "a" "c" "c" "g" "g" "g" "g"
##  [9541] "c" "t" "t" "g" "a" "t" "c" "a" "t" "a" "c" "a" "c" "a" "g" "t" "g" "g"
##  [9559] "c" "a" "t" "c" "t" "g" "c" "t" "c" "a" "t" "c" "t" "g" "a" "a" "c" "a"
##  [9577] "t" "c" "t" "c" "g" "c" "c" "g" "a" "t" "c" "t" "t" "t" "t" "a" "t" "g"
##  [9595] "c" "c" "t" "a" "c" "a" "a" "c" "g" "g" "c" "c" "a" "a" "t" "a" "a" "a"
##  [9613] "t" "c" "g" "a" "c" "t" "t" "a" "c" "c" "c" "g" "g" "t" "t" "a" "a" "c"
##  [9631] "a" "c" "c" "c" "c" "c" "a" "g" "g" "g" "t" "a" "t" "a" "a" "c" "a" "g"
##  [9649] "t" "c" "c" "t" "g" "g" "t" "g" "c" "a" "a" "g" "a" "t" "a" "t" "g" "a"
##  [9667] "c" "a" "a" "c" "a" "c" "c" "t" "t" "g" "a" "c" "c" "c" "t" "a" "c" "g"
##  [9685] "t" "a" "c" "c" "g" "t" "c" "t" "a" "a" "t" "c" "t" "a" "a" "c" "a" "g"
##  [9703] "t" "a" "c" "a" "c" "c" "t" "a" "g" "t" "c" "a" "a" "g" "a" "c" "g" "g"
##  [9721] "g" "t" "a" "g" "c" "c" "a" "g" "g" "g" "t" "g" "c" "c" "a" "g" "c" "a"
##  [9739] "t" "t" "g" "c" "g" "g" "c" "a" "g" "a" "a" "g" "c" "a" "t" "a" "t" "c"
##  [9757] "t" "c" "t" "g" "g" "a" "t" "c" "g" "a" "t" "a" "t" "a" "t" "g" "c" "g"
##  [9775] "a" "c" "t" "g" "c" "c" "a" "t" "c" "g" "a" "a" "c" "c" "g" "c" "g" "c"
##  [9793] "c" "a" "a" "g" "g" "g" "c" "a" "g" "g" "g" "c" "g" "g" "g" "a" "g" "a"
##  [9811] "a" "t" "g" "c" "c" "c" "c" "c" "a" "a" "t" "a" "t" "g" "g" "g" "c" "g"
##  [9829] "a" "a" "a" "g" "g" "t" "c" "c" "c" "g" "a" "g" "a" "a" "g" "t" "g" "c"
##  [9847] "g" "c" "t" "a" "t" "a" "t" "a" "a" "t" "a" "g" "t" "g" "c" "t" "g" "a"
##  [9865] "g" "t" "t" "g" "c" "c" "g" "c" "g" "c" "t" "a" "a" "a" "c" "c" "a" "a"
##  [9883] "t" "t" "c" "c" "c" "g" "a" "g" "c" "c" "g" "g" "c" "g" "c" "t" "c" "a"
##  [9901] "g" "c" "a" "a" "c" "t" "g" "g" "a" "c" "g" "a" "g" "c" "c" "t" "a" "c"
##  [9919] "c" "g" "c" "a" "g" "g" "g" "c" "g" "c" "g" "g" "g" "t" "c" "c" "t" "t"
##  [9937] "c" "g" "c" "c" "c" "c" "g" "c" "g" "a" "g" "c" "c" "g" "a" "a" "c" "a"
##  [9955] "c" "g" "c" "c" "g" "t" "g" "a" "a" "g" "g" "c" "t" "c" "t" "c" "t" "g"
##  [9973] "g" "g" "t" "a" "c" "g" "c" "t" "a" "c" "g" "a" "a" "c" "c" "t" "a" "t"
##  [9991] "c" "a" "a" "g" "g" "c" "g" "c" "g" "g" "t" "g" "g" "c" "g" "a" "c" "a"
## [10009] "g" "g" "c" "c" "c" "g" "a" "a" "g" "c" "c" "g" "a" "a" "t" "g" "c" "g"
## [10027] "a" "c" "c" "c" "t" "a" "t" "c" "a" "g" "a" "c" "g" "t" "a" "t" "a" "c"
## [10045] "c" "t" "g" "c" "c" "a" "c" "g" "g" "g" "a" "a" "a" "g" "t" "a" "g" "a"
## [10063] "g" "g" "a" "c" "a" "t" "c" "c" "c" "g" "t" "g" "a" "c" "g" "g" "a" "c"
## [10081] "c" "g" "g" "g" "t" "g" "c" "g" "a" "c" "g" "t" "a" "g" "a" "c" "g" "a"
## [10099] "a" "a" "g" "t" "g" "t" "c" "g" "a" "g" "c" "c" "a" "a" "c" "c" "t" "t"
## [10117] "a" "a" "a" "c" "a" "c" "a" "t" "t" "g" "a" "c" "c" "a" "g" "a" "g" "a"
## [10135] "g" "g" "a" "t" "c" "c" "g" "a" "g" "c" "g" "t" "c" "c" "a" "c" "t" "g"
## [10153] "a" "a" "g" "a" "g" "g" "g" "c" "a" "g" "a" "a" "a" "a" "c" "g" "c" "g"
## [10171] "c" "c" "c" "t" "a" "g" "g" "t" "a" "a" "c" "a" "c" "g" "c" "g" "t" "g"
## [10189] "t" "a" "a" "c" "c" "a" "a" "a" "c" "g" "c" "g" "g" "t" "c" "a" "g" "a"
## [10207] "t" "a" "a" "c" "c" "a" "a" "c" "c" "g" "t" "c" "t" "t" "a" "a" "t" "c"
## [10225] "a" "c" "t" "c" "t" "t" "c" "a" "c" "c" "a" "c" "a" "c" "a" "c" "g" "g"
## [10243] "t" "t" "g" "t" "g" "t" "a" "g" "t" "g" "t" "g" "c" "c" "c" "a" "c" "t"
## [10261] "c" "c" "a" "g" "c" "g" "t" "c" "c" "a" "c" "c" "c" "a" "t" "g" "t" "a"
## [10279] "a" "g" "g" "a" "g" "c" "c" "t" "g" "c" "t" "c" "t" "a" "g" "t" "c" "t"
## [10297] "a" "g" "g" "g" "c" "c" "a" "c" "c" "t" "g" "a" "c" "c" "t" "c" "a" "g"
## [10315] "t" "g" "c" "a" "a" "a" "c" "g" "a" "t" "g" "a" "g" "c" "t" "t" "t" "g"
## [10333] "c" "g" "t" "g" "c" "a" "g" "g" "a" "c" "a" "a" "g" "c" "c" "t" "c" "c"
## [10351] "c" "g" "g" "a" "t" "c" "c" "a" "g" "c" "c" "c" "t" "c" "a" "t" "c" "c"
## [10369] "a" "c" "c" "a" "a" "c" "a" "c" "t" "c" "a" "a" "g" "c" "c" "g" "c" "c"
## [10387] "g" "c" "g" "a" "c" "c" "c" "a" "c" "c" "c" "a" "g" "t" "c" "a" "c" "g"
## [10405] "t" "a" "g" "a" "g" "t" "t" "t" "a" "a" "a" "g" "c" "a" "g" "a" "t" "a"
## [10423] "g" "a" "c" "t" "g" "t" "g" "a" "g" "a" "t" "a" "c" "g" "a" "c" "c" "g"
## [10441] "g" "g" "a" "t" "g" "t" "a" "c" "a" "t" "g" "t" "g" "g" "t" "a" "g" "a"
## [10459] "t" "c" "c" "g" "a" "g" "a" "a" "t" "c" "c" "a" "t" "t" "c" "t" "t" "t"
## [10477] "g" "t" "c" "a" "c" "c" "a" "t" "t" "c" "c" "c" "t" "c" "c" "c" "g" "g"
## [10495] "c" "t" "g" "t" "a" "a" "g" "a" "t" "c" "t" "g" "t" "c" "g" "c" "t" "t"
## [10513] "g" "c" "c" "g" "c" "c" "c" "g" "g" "a" "c" "t" "t" "a" "t" "c" "c" "g"
## [10531] "c" "g" "a" "a" "a" "a" "t" "a" "t" "t" "a" "c" "a" "a" "c" "a" "g" "c"
## [10549] "a" "c" "t" "t" "a" "g" "a" "t" "g" "t" "a" "t" "a" "a" "c" "g" "g" "a"
## [10567] "c" "c" "g" "a" "c" "a" "g" "g" "g" "a" "c" "a" "t" "t" "a" "a" "a" "c"
## [10585] "t" "c" "t" "t" "g" "c" "c" "c" "c" "t" "c" "a" "g" "t" "c" "g" "c" "a"
## [10603] "g" "a" "a" "g" "t" "c" "c" "c" "g" "g" "g" "t" "g" "c" "t" "a" "c" "t"
## [10621] "g" "g" "g" "g" "g" "c" "t" "c" "c" "g" "t" "a" "t" "c" "t" "c" "a" "g"
## [10639] "g" "g" "c" "t" "c" "t" "a" "a" "g" "g" "t" "c" "t" "c" "a" "t" "g" "t"
## [10657] "a" "c" "t" "g" "a" "g" "c" "a" "c" "a" "g" "c" "g" "c" "c" "c" "g" "t"
## [10675] "t" "c" "a" "t" "g" "g" "c" "t" "c" "a" "g" "c" "t" "g" "t" "g" "t" "t"
## [10693] "c" "a" "t" "t" "c" "g" "t" "t" "c" "c" "g" "c" "a" "c" "t" "g" "c" "g"
## [10711] "g" "a" "t" "a" "a" "c" "t" "t" "c" "t" "t" "t" "c" "t" "g" "t" "a" "c"
## [10729] "a" "a" "c" "c" "g" "g" "c" "g" "a" "t" "a" "a" "t" "g" "a" "t" "a" "a"
## [10747] "c" "t" "g" "g" "a" "a" "g" "c" "a" "g" "a" "a" "a" "c" "t" "g" "t" "t"
## [10765] "c" "a" "a" "c" "c" "g" "g" "c" "c" "t" "c" "a" "c" "a" "g" "c" "a" "a"
## [10783] "g" "c" "c" "c" "t" "a" "c" "g" "c" "g" "c" "a" "t" "t" "a" "g" "a" "c"
## [10801] "a" "g" "g" "c" "c" "a" "g" "g" "g" "t" "a" "a" "c" "g" "t" "t" "c" "a"
## [10819] "g" "c" "g" "t" "g" "c" "g" "c" "t" "g" "g" "g" "t" "g" "c" "c" "g" "g"
## [10837] "t" "g" "a" "c" "g" "g" "a" "g" "c" "g" "a" "t" "c" "c" "g" "t" "g" "g"
## [10855] "c" "t" "t" "c" "c" "c" "g" "c" "g" "g" "a" "c" "g" "a" "g" "c" "a" "t"
## [10873] "c" "t" "c" "c" "a" "a" "a" "g" "c" "c" "a" "a" "a" "c" "c" "t" "t" "g"
## [10891] "g" "c" "a" "g" "t" "c" "g" "a" "t" "g" "c" "t" "g" "a" "t" "a" "c" "t"
## [10909] "g" "a" "t" "g" "c" "a" "c" "c" "g" "a" "g" "c" "c" "c" "g" "a" "c" "g"
## [10927] "c" "c" "a" "g" "g" "g" "a" "t" "t" "c" "g" "g" "a" "a" "c" "c" "c" "c"
## [10945] "c" "c" "a" "g" "t" "g" "t" "a" "a" "a" "a" "t" "a" "g" "c" "g" "t" "a"
## [10963] "a" "c" "t" "g" "g" "a" "c" "c" "a" "t" "c" "g" "c" "c" "c" "g" "a" "c"
## [10981] "t" "t" "c" "t" "c" "t" "g" "t" "g" "c" "a" "a" "g" "g" "t" "c" "a" "a"
## [10999] "t" "t" "c" "c" "g" "c" "g" "t" "c" "g" "t" "t" "c" "a" "c" "c" "c" "a"
## [11017] "g" "t" "c" "g" "c" "t" "c" "g" "t" "t" "g" "g" "t" "g" "a" "g" "t" "a"
## [11035] "t" "a" "c" "c" "c" "c" "c" "c" "c" "g" "a" "a" "c" "c" "t" "a" "c" "a"
## [11053] "a" "c" "c" "a" "c" "t" "t" "g" "a" "c" "c" "c" "t" "a" "a" "c" "a" "t"
## [11071] "a" "c" "t" "a" "c" "a" "t" "g" "t" "g" "t" "a" "a" "c" "c" "t" "a" "a"
## [11089] "a" "c" "t" "c" "g" "t" "a" "c" "g" "c" "c" "g" "c" "g" "g" "a" "g" "t"
## [11107] "c" "a" "g" "a" "g" "g" "c" "g" "g" "c" "a" "g" "c" "a" "c" "c" "g" "t"
## [11125] "c" "c" "g" "g" "a" "g" "g" "c" "c" "g" "c" "g" "t" "c" "t" "c" "c" "t"
## [11143] "t" "c" "t" "c" "a" "c" "a" "a" "c" "t" "c" "t" "a" "a" "c" "c" "t" "a"
## [11161] "g" "c" "c" "c" "t" "g" "a" "t" "c" "a" "c" "g" "c" "a" "c" "g" "a" "g"
## [11179] "g" "c" "a" "c" "a" "a" "g" "t" "a" "g" "t" "c" "c" "t" "a" "c" "c" "a"
## [11197] "c" "g" "c" "c" "g" "t" "c" "t" "t" "a" "c" "g" "t" "t" "c" "g" "t" "a"
## [11215] "c" "a" "g" "t" "c" "g" "c" "t" "c" "c" "t" "a" "t" "t" "c" "t" "a" "c"
## [11233] "g" "c" "a" "c" "c" "t" "c" "g" "g" "t" "a" "c" "t" "a" "c" "g" "g" "a"
## [11251] "a" "g" "a" "a" "t" "c" "t" "g" "t" "c" "t" "t" "g" "c" "g" "a" "g" "a"
## [11269] "g" "c" "c" "t" "t" "c" "c" "t" "t" "a" "t" "g" "a" "t" "a" "a" "c" "t"
## [11287] "t" "a" "a" "c" "c" "t" "g" "g" "t" "g" "g" "g" "c" "a" "c" "a" "g" "a"
## [11305] "c" "g" "a" "g" "c" "t" "g" "g" "t" "a" "t" "a" "g" "t" "c" "a" "a" "t"
## [11323] "a" "t" "c" "c" "c" "c" "t" "a" "g" "c" "a" "t" "g" "a" "c" "a" "c" "a"
## [11341] "g" "t" "g" "a" "a" "t" "c" "c" "a" "a" "c" "g" "c" "c" "t" "g" "c" "a"
## [11359] "a" "c" "c" "g" "c" "c" "g" "c" "c" "t" "c" "t" "t" "g" "c" "g" "c" "g"
## [11377] "c" "a" "c" "g" "c" "t" "g" "t" "g" "g" "a" "g" "t" "a" "a" "t" "g" "g"
## [11395] "c" "t" "a" "g" "g" "c" "c" "c" "g" "g" "a" "g" "a" "t" "g" "t" "a" "a"
## [11413] "g" "g" "c" "t" "c" "g" "a" "g" "c" "c" "t" "g" "g" "c" "c" "t" "a" "g"
## [11431] "c" "a" "g" "c" "c" "g" "c" "a" "t" "t" "g" "c" "g" "t" "a" "g" "t" "g"
## [11449] "g" "g" "c" "g" "g" "g" "a" "c" "g" "c" "t" "c" "g" "c" "a" "g" "t" "c"
## [11467] "g" "t" "a" "t" "g" "a" "g" "a" "t" "a" "c" "a" "a" "c" "g" "t" "a" "c"
## [11485] "t" "c" "t" "a" "c" "a" "g" "a" "c" "t" "g" "c" "g" "c" "g" "c" "c" "c"
## [11503] "g" "a" "g" "t" "c" "t" "g" "g" "a" "c" "a" "a" "t" "c" "g" "t" "c" "c"
## [11521] "a" "a" "t" "t" "a" "t" "t" "g" "a" "t" "a" "c" "t" "g" "a" "t" "t" "a"
## [11539] "t" "t" "g" "a" "a" "g" "g" "c" "a" "c" "c" "a" "a" "a" "g" "g" "a" "g"
## [11557] "t" "t" "c" "g" "g" "g" "a" "t" "g" "t" "a" "c" "g" "c" "t" "g" "t" "c"
## [11575] "c" "c" "c" "c" "a" "t" "g" "a" "a" "t" "g" "t" "c" "g" "a" "a" "c" "c"
## [11593] "g" "a" "c" "g" "a" "a" "t" "t" "a" "t" "t" "t" "c" "g" "t" "g" "c" "g"
## [11611] "c" "a" "c" "a" "c" "g" "a" "g" "t" "c" "a" "t" "c" "a" "a" "g" "g" "a"
## [11629] "a" "t" "t" "g" "t" "c" "g" "a" "g" "c" "t" "g" "g" "c" "g" "g" "a" "g"
## [11647] "c" "g" "t" "c" "a" "c" "g" "a" "c" "c" "g" "c" "g" "t" "c" "a" "a" "g"
## [11665] "g" "c" "c" "a" "a" "t" "g" "a" "g" "c" "c" "c" "t" "a" "g" "t" "t" "t"
## [11683] "g" "a" "a" "g" "c" "t" "g" "t" "a" "a" "t" "t" "t" "a" "c" "t" "c" "g"
## [11701] "g" "c" "c" "a" "t" "t" "t" "c" "c" "c" "a" "a" "g" "c" "g" "a" "a" "a"
## [11719] "a" "g" "g" "t" "a" "g" "c" "c" "t" "g" "t" "t" "g" "t" "a" "a" "g" "a"
## [11737] "g" "a" "g" "g" "a" "g" "c" "a" "c" "c" "a" "g" "g" "c" "t" "a" "a" "g"
## [11755] "a" "g" "c" "a" "g" "a" "c" "a" "g" "c" "a" "g" "c" "a" "g" "c" "a" "a"
## [11773] "g" "g" "t" "g" "c" "t" "c" "g" "a" "c" "a" "g" "a" "t" "g" "a" "a" "g"
## [11791] "c" "a" "c" "t" "g" "g" "g" "g" "g" "c" "g" "c" "g" "c" "c" "g" "c" "c"
## [11809] "c" "c" "g" "c" "c" "g" "t" "t" "a" "c" "c" "c" "a" "g" "c" "a" "t" "g"
## [11827] "c" "a" "g" "a" "a" "c" "a" "g" "g" "g" "t" "t" "g" "g" "c" "c" "a" "t"
## [11845] "a" "t" "a" "g" "t" "g" "c" "a" "a" "a" "c" "a" "t" "g" "g" "t" "g" "c"
## [11863] "a" "a" "c" "c" "c" "t" "c" "c" "g" "c" "c" "c" "g" "a" "t" "c" "c" "g"
## [11881] "g" "t" "g" "a" "g" "c" "t" "t" "c" "c" "g" "t" "t" "a" "t" "a" "t" "c"
## [11899] "c" "t" "a" "a" "a" "c" "t" "g" "c" "a" "a" "g" "g" "a" "a" "c" "g" "c"
## [11917] "a" "c" "c" "t" "a" "g" "g" "t" "a" "t" "c" "a" "a" "c" "a" "c" "c" "c"
## [11935] "t" "t" "c" "a" "g" "t" "t" "c" "c" "a" "g" "c" "g" "a" "c" "c" "c" "g"
## [11953] "t" "t" "c" "c" "a" "c" "a" "g" "c" "g" "a" "a" "c" "a" "c" "a" "a" "c"
## [11971] "g" "t" "a" "a" "g" "g" "g" "c" "g" "c" "g" "a" "c" "t" "c" "g" "a" "c"
## [11989] "c" "c" "a" "a" "t" "a" "g" "t" "g" "g" "a" "a" "g" "a" "g" "c" "a" "t"
## [12007] "g" "g" "t" "t" "t" "c" "g" "a" "g" "g" "g" "t" "c" "a" "t" "c" "t" "g"
## [12025] "g" "t" "a" "c" "c" "g" "t" "t" "c" "g" "t" "c" "c" "g" "c" "c" "g" "g"
## [12043] "c" "a" "a" "t" "c" "t" "g" "g" "c" "c" "a" "a" "t" "g" "g" "a" "g" "g"
## [12061] "c" "c" "c" "a" "g" "t" "a" "a" "c" "a" "a" "g" "g" "c" "a" "g" "c" "g"
## [12079] "g" "c" "c" "g" "c" "t" "a" "a" "g" "g" "c" "g" "t" "c" "a" "c" "g" "g"
## [12097] "g" "t" "a" "c" "g" "a" "a" "c" "t" "g" "g" "a" "a" "t" "a" "a" "c" "t"
## [12115] "c" "t" "c" "g" "a" "t" "c" "a" "g" "g" "a" "c" "g" "a" "c" "a" "c" "t"
## [12133] "c" "t" "t" "g" "g" "g" "c" "g" "a" "g" "c" "g" "c" "t" "a" "t" "t" "g"
## [12151] "g" "c" "g" "c" "a" "c" "t" "t" "c" "c" "c" "c" "a" "c" "a" "a" "t" "a"
## [12169] "c" "g" "c" "c" "g" "g" "c" "g" "g" "g" "a" "c" "a" "a" "g" "a" "g" "c"
## [12187] "g" "a" "c" "a" "g" "c" "a" "c" "g" "t" "g" "c" "t" "t" "g" "t" "t" "g"
## [12205] "g" "c" "g" "g" "t" "a" "g" "a" "c" "a" "a" "c" "c" "c" "c" "g" "a" "g"
## [12223] "g" "c" "g" "g" "g" "c" "c" "c" "g" "a" "a" "c" "c" "c" "c" "a" "c" "a"
## [12241] "g" "c" "a" "g" "g" "t" "c" "a" "g" "a" "g" "a" "g" "g" "t" "g" "a" "c"
## [12259] "c" "c" "a" "c" "c" "a" "t" "g" "t" "a" "a" "g" "t" "c" "c" "g" "g" "c"
## [12277] "c" "g" "t" "g" "a" "c" "g" "g" "a" "a" "g" "g" "c" "c" "a" "t" "c" "c"
## [12295] "g" "g" "c" "c" "c" "t" "c" "t" "c" "t" "g" "c" "a" "t" "g" "g" "t" "a"
## [12313] "a" "c" "t" "t" "a" "c" "a" "t" "g" "t" "c" "a" "c" "g" "a" "g" "g" "c"
## [12331] "a" "t" "a" "g" "c" "c" "g" "t" "t" "g" "t" "c" "a" "a" "g" "t" "c" "g"
## [12349] "g" "a" "g" "c" "g" "c" "c" "a" "a" "g" "c" "c" "t" "g" "a" "t" "c" "c"
## [12367] "a" "g" "t" "t" "t" "g" "g" "c" "t" "a" "t" "c" "g" "a" "t" "g" "c" "a"
## [12385] "t" "c" "t" "c" "t" "c" "g" "g" "c" "g" "c" "a" "c" "g" "a" "g" "a" "g"
## [12403] "t" "g" "c" "g" "c" "g" "t" "c" "t" "a" "a" "a" "c" "g" "t" "c" "a" "t"
## [12421] "g" "t" "c" "c" "g" "t" "g" "a" "g" "c" "c" "t" "a" "a" "a" "g" "a" "a"
## [12439] "g" "c" "g" "t" "c" "g" "a" "g" "c" "c" "g" "g" "t" "t" "c" "a" "g" "g"
## [12457] "g" "c" "t" "c" "a" "t" "c" "g" "g" "c" "a" "a" "c" "c" "c" "a" "g" "c"
## [12475] "c" "t" "c" "a" "a" "c" "a" "a" "c" "c" "c" "t" "a" "c" "t" "t" "c" "c"
## [12493] "a" "c" "g" "g" "t" "c" "c" "g" "c" "t" "c" "g" "a" "t" "t" "a" "c" "g"
## [12511] "c" "g" "t" "a" "c" "t" "c" "a" "c" "a" "a" "g" "c" "a" "a" "g" "a" "g"
## [12529] "g" "g" "t" "g" "t" "t" "a" "g" "g" "t" "a" "c" "a" "c" "a" "c" "c" "g"
## [12547] "c" "c" "t" "a" "g" "g" "t" "t" "a" "g" "a" "g" "g" "g" "c" "c" "a" "a"
## [12565] "t" "g" "a" "a" "g" "t" "c" "c" "a" "a" "g" "t" "c" "a" "c" "t" "a" "c"
## [12583] "g" "a" "g" "g" "c" "g" "c" "a" "t" "t" "a" "c" "c" "g" "c" "a" "t" "t"
## [12601] "t" "a" "c" "a" "a" "a" "g" "c" "a" "g" "a" "t" "g" "a" "g" "g" "c" "g"
## [12619] "t" "g" "g" "c" "a" "g" "c" "c" "c" "a" "a" "c" "c" "g" "g" "c" "g" "c"
## [12637] "c" "t" "g" "a" "a" "a" "g" "t" "g" "a" "g" "t" "a" "a" "a" "t" "g" "a"
## [12655] "g" "t" "t" "a" "a" "a" "t" "g" "g" "g" "t" "g" "g" "a" "g" "c" "g" "a"
## [12673] "a" "t" "c" "t" "c" "g" "g" "g" "g" "g" "a" "t" "t" "c" "a" "a" "t" "a"
## [12691] "c" "a" "t" "a" "c" "g" "a" "t" "g" "c" "a" "g" "a" "a" "a" "t" "a" "g"
## [12709] "a" "a" "a" "g" "c" "c" "t" "g" "a" "c" "t" "g" "a" "g" "a" "a" "c" "g"
## [12727] "g" "t" "a" "t" "g" "g" "g" "g" "c" "c" "g" "a" "c" "g" "c" "g" "a" "t"
## [12745] "g" "c" "g" "t" "g" "t" "a" "g" "a" "c" "g" "c" "t" "g" "c" "t" "a" "c"
## [12763] "c" "c" "a" "g" "c" "c" "t" "c" "g" "g" "t" "t" "c" "g" "g" "c" "a" "t"
## [12781] "t" "c" "t" "c" "t" "a" "t" "t" "c" "g" "c" "g" "c" "g" "t" "t" "g" "a"
## [12799] "c" "c" "a" "g" "a" "g" "c" "g" "g" "g" "g" "t" "g" "t" "a" "t" "t" "c"
## [12817] "c" "c" "t" "t" "c" "t" "g" "a" "g" "a" "g" "g" "g" "c" "c" "g" "c" "a"
## [12835] "c" "t" "a" "g" "g" "a" "c" "t" "t" "c" "c" "t" "a" "c" "a" "t" "g" "t"
## [12853] "c" "c" "a" "a" "t" "g" "a" "g" "g" "t" "a" "a" "a" "a" "t" "g" "c" "c"
## [12871] "a" "a" "a" "c" "c" "t" "t" "g" "a" "c" "g" "a" "c" "a" "c" "c" "c" "g"
## [12889] "a" "a" "a" "t" "c" "g" "a" "g" "t" "a" "a" "a" "t" "c" "c" "g" "a" "c"
## [12907] "g" "c" "c" "g" "c" "g" "t" "g" "a" "t" "g" "c" "a" "g" "a" "c" "c" "g"
## [12925] "c" "g" "a" "c" "t" "c" "c" "a" "g" "g" "t" "g" "a" "a" "t" "g" "c" "c"
## [12943] "t" "c" "g" "a" "t" "c" "c" "t" "c" "c" "g" "g" "c" "a" "a" "c" "g" "c"
## [12961] "a" "c" "a" "c" "c" "c" "a" "a" "a" "g" "a" "c" "g" "c" "c" "t" "t" "t"
## [12979] "a" "g" "c" "g" "a" "a" "c" "t" "g" "a" "a" "a" "a" "a" "a" "a" "a" "t"
## [12997] "t" "g" "c" "g" "c" "c" "a" "g" "g" "a" "a" "a" "a" "c" "c" "c" "g" "g"
## [13015] "c" "t" "g" "g" "c" "a" "c" "c" "t" "t" "a" "a" "a" "c" "g" "g" "t" "c"
## [13033] "t" "c" "g" "g" "a" "c" "a" "g" "c" "a" "c" "t" "t" "c" "t" "c" "g" "t"
## [13051] "g" "g" "g" "t" "t" "g" "c" "a" "c" "a" "g" "c" "c" "c" "c" "g" "a" "t"
## [13069] "a" "a" "c" "c" "a" "g" "g" "g" "g" "c" "a" "t" "g" "a" "a" "a" "g" "a"
## [13087] "t" "g" "t" "t" "t" "t" "a" "g" "g" "c" "a" "c" "t" "g" "g" "g" "g" "g"
## [13105] "g" "a" "g" "c" "g" "a" "t" "g" "t" "t" "c" "c" "c" "c" "g" "t" "a" "c"
## [13123] "g" "c" "a" "t" "c" "a" "c" "c" "c" "c" "c" "a" "c" "c" "a" "a" "g" "t"
## [13141] "g" "t" "t" "g" "c" "g" "t" "g" "g" "a" "t" "c" "a" "t" "g" "g" "c" "c"
## [13159] "a" "a" "g" "a" "c" "g" "a" "c" "c" "c" "g" "g" "t" "t" "t" "c" "g" "c"
## [13177] "g" "g" "c" "g" "t" "g" "a" "c" "t" "g" "a" "a" "g" "g" "g" "a" "g" "c"
## [13195] "a" "c" "c" "c" "a" "c" "t" "t" "c" "c" "c" "t" "c" "c" "g" "a" "g" "g"
## [13213] "g" "t" "a" "a" "t" "t" "t" "g" "g" "c" "a" "g" "g" "a" "g" "c" "g" "a"
## [13231] "a" "g" "a" "a" "c" "c" "a" "g" "c" "g" "c" "c" "g" "c" "g" "g" "g" "t"
## [13249] "a" "a" "c" "c" "a" "t" "a" "g" "a" "g" "a" "g" "t" "g" "a" "t" "g" "a"
## [13267] "a" "c" "t" "t" "g" "g" "g" "a" "g" "t" "a" "g" "c" "g" "a" "c" "t" "c"
## [13285] "c" "t" "c" "c" "c" "a" "t" "a" "g" "t" "g" "c" "g" "t" "a" "c" "c" "a"
## [13303] "c" "a" "a" "c" "t" "t" "g" "a" "c" "a" "g" "a" "a" "c" "g" "c" "g" "g"
## [13321] "c" "c" "a" "a" "g" "t" "c" "a" "c" "a" "c" "t" "a" "a" "a" "c" "t" "a"
## [13339] "t" "t" "t" "t" "c" "g" "c" "g" "t" "c" "a" "c" "g" "c" "g" "g" "a" "a"
## [13357] "g" "a" "a" "a" "t" "g" "a" "c" "c" "t" "a" "t" "c" "g" "c" "a" "t" "g"
## [13375] "c" "a" "t" "c" "t" "g" "g" "c" "t" "g" "g" "c" "c" "a" "t" "g" "t" "g"
## [13393] "g" "a" "g" "g" "c" "g" "g" "t" "c" "g" "t" "t" "a" "g" "a" "t" "g" "c"
## [13411] "a" "a" "c" "c" "t" "c" "g" "c" "t" "c" "g" "g" "t" "a" "t" "c" "t" "c"
## [13429] "t" "t" "t" "t" "c" "c" "g" "g" "t" "t" "c" "c" "g" "a" "t" "c" "c" "t"
## [13447] "c" "c" "g" "a" "t" "c" "c" "c" "c" "g" "a" "c" "c" "a" "g" "c" "t" "c"
## [13465] "c" "t" "g" "c" "c" "t" "c" "g" "c" "t" "a" "a" "a" "g" "a" "g" "c" "t"
## [13483] "g" "g" "g" "t" "t" "a" "t" "c" "a" "t" "g" "c" "g" "t" "c" "g" "c" "g"
## [13501] "g" "g" "c" "g" "a" "g" "t" "g" "c" "c" "g" "c" "a" "c" "c" "c" "g" "t"
## [13519] "g" "c" "t" "t" "g" "g" "g" "a" "a" "g" "t" "g" "g" "t" "a" "c" "c" "a"
## [13537] "g" "c" "a" "c" "a" "a" "c" "c" "g" "g" "c" "a" "g" "a" "g" "a" "c" "a"
## [13555] "g" "t" "g" "c" "a" "g" "a" "g" "t" "a" "c" "g" "c" "c" "g" "c" "c" "t"
## [13573] "c" "a" "g" "a" "t" "c" "c" "t" "g" "g" "c" "t" "c" "a" "g" "t" "g" "c"
## [13591] "c" "c" "t" "c" "a" "t" "g" "g" "g" "g" "c" "a" "t" "c" "c" "c" "t" "c"
## [13609] "g" "c" "t" "a" "g" "g" "g" "t" "c" "c" "c" "g" "t" "a" "c" "c" "c" "a"
## [13627] "c" "g" "c" "g" "a" "c" "c" "g" "c" "c" "a" "t" "c" "c" "a" "c" "c" "a"
## [13645] "c" "t" "g" "a" "a" "a" "c" "t" "c" "a" "c" "a" "c" "g" "a" "c" "c" "g"
## [13663] "c" "a" "t" "g" "t" "g" "c" "t" "c" "t" "g" "a" "g" "c" "c" "g" "c" "g"
## [13681] "t" "a" "c" "g" "g" "c" "t" "c" "a" "c" "c" "g" "a" "c" "c" "a" "c" "c"
## [13699] "a" "t" "t" "g" "c" "c" "a" "a" "c" "c" "c" "a" "t" "a" "t" "a" "t" "a"
## [13717] "a" "t" "c" "t" "t" "a" "t" "g" "g" "c" "t" "g" "g" "t" "g" "c" "c" "g"
## [13735] "a" "g" "t" "c" "c" "c" "a" "g" "c" "c" "c" "c" "g" "g" "a" "a" "t" "g"
## [13753] "g" "a" "c" "t" "g" "g" "g" "a" "t" "g" "a" "a" "t" "g" "t" "g" "t" "t"
## [13771] "a" "t" "a" "t" "c" "g" "g" "g" "t" "c" "c" "c" "a" "c" "c" "g" "c" "t"
## [13789] "c" "c" "g" "c" "g" "a" "c" "a" "a" "g" "g" "c" "g" "g" "g" "c" "a" "g"
## [13807] "a" "t" "a" "g" "g" "c" "a" "c" "a" "a" "c" "g" "t" "c" "c" "c" "c" "g"
## [13825] "g" "c" "t" "a" "a" "c" "a" "t" "a" "a" "c" "c" "a" "t" "c" "c" "c" "t"
## [13843] "a" "c" "c" "c" "a" "a" "a" "t" "c" "a" "a" "a" "c" "c" "g" "t" "a" "g"
## [13861] "t" "t" "g" "c" "g" "g" "t" "g" "c" "c" "c" "a" "c" "a" "a" "c" "a" "c"
## [13879] "a" "g" "a" "g" "a" "g" "g" "a" "c" "a" "c" "t" "g" "c" "c" "a" "g" "g"
## [13897] "c" "g" "c" "g" "c" "t" "t" "g" "g" "g" "c" "t" "g" "t" "a" "c" "c" "g"
## [13915] "c" "a" "a" "c" "a" "g" "a" "t" "t" "c" "c" "t" "g" "a" "t" "t" "c" "c"
## [13933] "t" "g" "c" "g" "c" "g" "g" "a" "c" "a" "a" "a" "t" "t" "a" "t" "g" "g"
## [13951] "a" "t" "g" "a" "g" "a" "g" "g" "t" "g" "c" "t" "c" "c" "c" "t" "a" "a"
## [13969] "g" "a" "t" "t" "c" "a" "g" "g" "c" "g" "g" "a" "g" "t" "g" "g" "g" "c"
## [13987] "c" "t" "c" "g" "a" "g" "c" "a" "g" "c" "c" "a" "a" "c" "t" "t" "g" "g"
## [14005] "c" "c" "t" "t" "t" "a" "a" "t" "g" "a" "a" "c" "c" "a" "c" "a" "a" "g"
## [14023] "g" "t" "g" "c" "t" "g" "c" "a" "c" "c" "a" "c" "c" "t" "a" "g" "c" "c"
## [14041] "a" "a" "g" "c" "c" "c" "g" "g" "c" "c" "g" "t" "c" "c" "g" "c" "a" "c"
## [14059] "t" "t" "a" "g" "a" "g" "c" "a" "c" "t" "g" "c" "t" "c" "g" "t" "a" "a"
## [14077] "g" "a" "t" "t" "g" "t" "g" "t" "a" "g" "g" "c" "a" "a" "c" "a" "c" "c"
## [14095] "g" "a" "t" "g" "g" "t" "a" "c" "a" "t" "g" "t" "g" "c" "g" "c" "c" "g"
## [14113] "a" "a" "c" "g" "g" "a" "t" "c" "g" "t" "c" "t" "g" "t" "g" "g" "c" "a"
## [14131] "g" "g" "c" "g" "g" "g" "a" "g" "t" "c" "a" "c" "t" "c" "a" "c" "a" "c"
## [14149] "t" "t" "a" "a" "t" "a" "t" "g" "a" "a" "t" "t" "g" "t" "g" "t" "a" "c"
## [14167] "c" "c" "a" "g" "t" "g" "t" "a" "t" "c" "g" "g" "g" "a" "g" "c" "t" "a"
## [14185] "c" "a" "g" "c" "a" "g" "c" "c" "c" "g" "g" "t" "g" "g" "c" "a" "a" "a"
## [14203] "a" "g" "a" "a" "a" "t" "a" "g" "g" "c" "g" "g" "c" "c" "c" "t" "c" "a"
## [14221] "g" "c" "g" "c" "g" "a" "a" "a" "g" "g" "c" "g" "c" "a" "a" "c" "t" "c"
## [14239] "c" "a" "a" "g" "g" "c" "c" "t" "c" "g" "g" "a" "t" "g" "c" "a" "t" "g"
## [14257] "t" "t" "g" "a" "c" "t" "c" "g" "a" "g" "g" "c" "c" "t" "t" "a" "g" "c"
## [14275] "c" "a" "c" "c" "a" "a" "c" "g" "t" "a" "g" "c" "t" "g" "c" "g" "c" "a"
## [14293] "g" "t" "c" "c" "c" "g" "t" "t" "c" "t" "c" "g" "a" "a" "g" "a" "t" "a"
## [14311] "g" "c" "a" "c" "g" "g" "g" "t" "c" "c" "c" "a" "a" "a" "a" "a" "g" "c"
## [14329] "a" "g" "g" "c" "g" "c" "c" "a" "a" "a" "g" "a" "a" "a" "a" "c" "g" "c"
## [14347] "c" "a" "g" "a" "c" "t" "g" "g" "c" "c" "g" "c" "t" "c" "c" "g" "c" "a"
## [14365] "c" "g" "t" "t" "a" "c" "t" "t" "c" "a" "c" "a" "t" "c" "c" "g" "t" "t"
## [14383] "g" "g" "c" "a" "t" "t" "g" "t" "t" "c" "g" "g" "a" "a" "g" "c" "g" "c"
## [14401] "g" "g" "c" "t" "g" "g" "t" "g" "t" "c" "t" "c" "g" "g" "g" "c" "a" "a"
## [14419] "a" "t" "t" "c" "t" "g" "a" "g" "c" "t" "c" "t" "a" "c" "g" "c" "c" "t"
## [14437] "g" "c" "t" "c" "c" "a" "a" "t" "a" "a" "g" "t" "g" "g" "c" "t" "c" "a"
## [14455] "g" "t" "a" "t" "c" "a" "g" "a" "c" "g" "t" "c" "a" "a" "c" "t" "a" "c"
## [14473] "c" "c" "g" "a" "c" "c" "g" "a" "c" "g" "t" "g" "g" "a" "g" "t" "a" "t"
## [14491] "c" "g" "g" "c" "a" "g" "a" "c" "a" "g" "c" "g" "a" "a" "a" "c" "t" "c"
## [14509] "c" "a" "g" "a" "a" "g" "t" "c" "a" "g" "g" "t" "c" "c" "g" "g" "g" "t"
## [14527] "c" "c" "a" "g" "g" "g" "c" "a" "a" "c" "t" "c" "g" "a" "g" "a" "a" "c"
## [14545] "t" "a" "t" "c" "c" "g" "t" "g" "g" "a" "g" "c" "t" "g" "c" "g" "c" "a"
## [14563] "g" "t" "t" "g" "t" "c" "a" "t" "a" "t" "t" "c" "a" "c" "a" "g" "c" "t"
## [14581] "g" "t" "a" "a" "t" "g" "c" "c" "a" "t" "c" "t" "c" "a" "t" "g" "t" "g"
## [14599] "c" "c" "c" "g" "a" "c" "t" "c" "a" "a" "g" "c" "t" "c" "a" "a" "t" "g"
## [14617] "g" "c" "a" "t" "c" "a" "a" "c" "a" "c" "t" "a" "c" "t" "g" "t" "c" "t"
## [14635] "c" "c" "c" "c" "t" "t" "g" "a" "g" "c" "g" "a" "g" "c" "t" "a" "a" "a"
## [14653] "g" "c" "g" "a" "g" "a" "a" "a" "a" "g" "g" "g" "t" "t" "t" "a" "t" "g"
## [14671] "g" "g" "c" "c" "a" "g" "c" "c" "t" "g" "a" "a" "c" "c" "t" "c" "c" "g"
## [14689] "a" "a" "t" "c" "t" "a" "a" "g" "c" "a" "c" "g" "a" "c" "g" "t" "a" "a"
## [14707] "c" "c" "a" "t" "a" "g" "c" "t" "c" "a" "a" "c" "g" "c" "c" "g" "t" "c"
## [14725] "c" "t" "a" "a" "c" "g" "g" "g" "g" "g" "a" "c" "t" "t" "c" "a" "c" "g"
## [14743] "a" "c" "c" "g" "c" "c" "a" "a" "t" "a" "c" "a" "c" "g" "a" "a" "t" "g"
## [14761] "c" "t" "c" "c" "g" "g" "g" "g" "a" "c" "g" "a" "a" "t" "a" "c" "t" "t"
## [14779] "g" "c" "c" "c" "t" "c" "c" "a" "c" "c" "c" "g" "c" "a" "a" "a" "a" "t"
## [14797] "g" "a" "a" "g" "t" "g" "c" "c" "g" "c" "a" "c" "g" "a" "c" "a" "g" "a"
## [14815] "t" "g" "g" "g" "c" "a" "t" "g" "t" "t" "c" "c" "t" "a" "g" "a" "c" "c"
## [14833] "g" "c" "t" "c" "g" "g" "t" "a" "a" "a" "t" "a" "g" "a" "a" "g" "a" "g"
## [14851] "a" "g" "g" "c" "c" "c" "g" "a" "a" "t" "g" "t" "t" "g" "g" "c" "a" "a"
## [14869] "a" "c" "c" "g" "c" "c" "g" "g" "g" "a" "c" "t" "c" "c" "t" "c" "t" "t"
## [14887] "a" "a" "c" "c" "c" "g" "a" "t" "a" "g" "g" "g" "t" "g" "c" "g" "a" "a"
## [14905] "a" "c" "t" "c" "g" "t" "t" "a" "c" "g" "c" "g" "a" "g" "a" "c" "g" "g"
## [14923] "a" "c" "g" "t" "g" "c" "a" "c" "c" "g" "g" "a" "a" "a" "t" "t" "c" "a"
## [14941] "g" "c" "c" "t" "t" "t" "g" "a" "g" "a" "t" "t" "g" "t" "a" "g" "g" "g"
## [14959] "t" "t" "t" "g" "a" "a" "c" "t" "a" "c" "g" "t" "g" "a" "a" "t" "t" "g"
## [14977] "a" "c" "g" "c" "t" "c" "a" "t" "a" "c" "g" "t" "a" "a" "c" "a" "a" "a"
## [14995] "c" "a" "c" "a" "a" "g" "a" "g" "g" "t" "c" "a" "c" "c" "t" "a" "c" "c"
## [15013] "g" "g" "t" "a" "t" "g" "c" "c" "a" "c" "c" "c" "g" "a" "a" "c" "g" "a"
## [15031] "c" "t" "g" "c" "a" "g" "t" "t" "t" "a" "c" "g" "c" "g" "t" "c" "c" "c"
## [15049] "c" "t" "t" "g" "a" "c" "t" "c" "g" "g" "a" "c" "t" "g" "a" "g" "g" "c"
## [15067] "a" "g" "t" "a" "g" "c" "g" "t" "g" "c" "a" "g" "c" "g" "a" "a" "g" "g"
## [15085] "g" "a" "c" "c" "t" "g" "a" "g" "g" "g" "t" "c" "a" "c" "a" "a" "a" "a"
## [15103] "t" "a" "g" "t" "t" "c" "a" "a" "t" "g" "t" "g" "c" "a" "t" "t" "c" "g"
## [15121] "c" "t" "a" "c" "g" "t" "g" "g" "t" "a" "t" "t" "a" "c" "a" "c" "g" "a"
## [15139] "a" "a" "t" "t" "g" "c" "c" "a" "t" "c" "c" "c" "g" "t" "g" "a" "g" "a"
## [15157] "g" "a" "a" "c" "a" "c" "c" "a" "t" "a" "g" "a" "a" "g" "a" "g" "g" "a"
## [15175] "c" "g" "c" "c" "a" "g" "t" "g" "a" "g" "a" "a" "c" "c" "g" "g" "g" "t"
## [15193] "g" "g" "g" "g" "c" "t" "a" "g" "c" "a" "a" "g" "c" "a" "g" "g" "g" "a"
## [15211] "g" "t" "c" "g" "t" "t" "t" "c" "a" "t" "t" "c" "c" "g" "t" "g" "a" "c"
## [15229] "g" "a" "a" "t" "c" "g" "c" "c" "t" "g" "c" "g" "c" "c" "g" "c" "t" "a"
## [15247] "t" "g" "g" "a" "c" "a" "c" "g" "a" "g" "t" "t" "g" "c" "t" "g" "g" "t"
## [15265] "c" "c" "t" "g" "g" "t" "t" "t" "t" "c" "c" "a" "t" "c" "a" "a" "t" "t"
## [15283] "t" "g" "t" "g" "c" "g" "a" "a" "g" "c" "g" "a" "c" "t" "c" "c" "a" "a"
## [15301] "t" "g" "c" "g" "t" "c" "g" "g" "c" "t" "a" "g" "a" "c" "g" "a" "c" "c"
## [15319] "t" "c" "t" "c" "g" "t" "g" "c" "g" "a" "c" "a" "c" "t" "g" "a" "c" "a"
## [15337] "a" "t" "t" "g" "g" "t" "a" "g" "t" "c" "c" "a" "g" "a" "c" "a" "g" "t"
## [15355] "g" "c" "g" "c" "c" "t" "c" "g" "a" "t" "c" "t" "t" "t" "t" "g" "c" "a"
## [15373] "t" "a" "c" "t" "g" "t" "c" "c" "a" "c" "g" "g" "t" "c" "g" "c" "a" "c"
## [15391] "c" "t" "c" "c" "t" "c" "t" "t" "a" "g" "t" "g" "t" "t" "t" "g" "t" "g"
## [15409] "c" "g" "t" "g" "t" "t" "g" "c" "a" "c" "g" "g" "g" "a" "c" "c" "c" "g"
## [15427] "c" "c" "a" "c" "g" "t" "t" "c" "g" "g" "t" "a" "a" "a" "g" "t" "g" "c"
## [15445] "t" "a" "c" "a" "t" "t" "a" "a" "c" "g" "t" "t" "t" "g" "c" "g" "t" "a"
## [15463] "c" "c" "a" "c" "a" "c" "t" "a" "a" "t" "t" "c" "t" "a" "t" "a" "a" "a"
## [15481] "c" "a" "t" "t" "g" "c" "t" "a" "c" "c" "c" "t" "t" "a" "c" "a" "c" "a"
## [15499] "g" "g" "a" "c" "c" "a" "a" "a" "c" "g" "c" "c" "a" "a" "g" "g" "t" "c"
## [15517] "g" "t" "g" "a" "g" "c" "t" "t" "t" "c" "a" "c" "a" "g" "c" "c" "c" "c"
## [15535] "c" "a" "c" "g" "t" "g" "c" "g" "a" "t" "a" "c" "t" "t" "a" "t" "c" "t"
## [15553] "g" "c" "c" "g" "g" "a" "g" "t" "g" "c" "t" "t" "t" "g" "c" "g" "c" "c"
## [15571] "g" "a" "a" "c" "g" "g" "a" "g" "a" "c" "a" "t" "t" "a" "a" "a" "g" "c"
## [15589] "a" "g" "t" "t" "c" "t" "c" "g" "g" "t" "g" "c" "t" "c" "g" "a" "g" "g"
## [15607] "a" "t" "g" "t" "c" "c" "a" "c" "t" "t" "a" "t" "g" "g" "a" "a" "a" "c"
## [15625] "c" "g" "g" "g" "a" "a" "c" "t" "c" "t" "c" "t" "c" "g" "g" "a" "c" "g"
## [15643] "c" "t" "t" "g" "c" "c" "a" "g" "g" "g" "a" "t" "g" "c" "c" "g" "a" "a"
## [15661] "t" "g" "g" "a" "c" "c" "c" "g" "t" "c" "c" "c" "g" "g" "c" "t" "a" "g"
## [15679] "a" "g" "a" "c" "t" "c" "c" "c" "g" "a" "t" "t" "t" "t" "a" "c" "c" "t"
## [15697] "g" "t" "c" "g" "c" "a" "c" "t" "g" "t" "g" "a" "g" "g" "c" "a" "a" "c"
## [15715] "g" "g" "c" "g" "c" "c" "g" "c" "a" "t" "a" "a" "t" "g" "a" "c" "g" "g"
## [15733] "t" "t" "g" "a" "g" "a" "a" "c" "g" "g" "c" "a" "g" "t" "t" "g" "a" "c"
## [15751] "g" "t" "c" "a" "a" "a" "t" "c" "t" "c" "g" "c" "c" "g" "t" "a" "g" "a"
## [15769] "g" "c" "t" "g" "g" "a" "t" "c" "g" "g" "t" "a" "g" "c" "c" "a" "c" "c"
## [15787] "c" "t" "t" "c" "g" "a" "c" "c" "a" "c" "c" "t" "a" "c" "g" "c" "c" "g"
## [15805] "c" "t" "t" "c" "a" "t" "t" "c" "c" "g" "a" "a" "g" "t" "c" "a" "c" "g"
## [15823] "t" "a" "t" "a" "c" "c" "g" "c" "c" "g" "g" "c" "t" "g" "t" "c" "t" "t"
## [15841] "t" "c" "g" "g" "t" "g" "a" "t" "g" "c" "t" "g" "g" "t" "t" "g" "a" "c"
## [15859] "a" "t" "a" "c" "g" "c" "g" "g" "t" "c" "c" "c" "c" "a" "t" "g" "a" "t"
## [15877] "a" "c" "a" "g" "t" "a" "g" "a" "g" "g" "g" "t" "c" "t" "a" "g" "c" "t"
## [15895] "a" "c" "c" "g" "t" "c" "g" "t" "c" "c" "t" "g" "a" "t" "c" "c" "a" "t"
## [15913] "t" "g" "t" "c" "a" "a" "t" "g" "a" "c" "a" "t" "g" "c" "g" "g" "c" "g"
## [15931] "a" "g" "c" "t" "a" "c" "a" "g" "t" "t" "a" "t" "c" "c" "t" "a" "c" "c"
## [15949] "t" "t" "t" "t" "g" "a" "t" "a" "c" "a" "c" "c" "a" "g" "a" "c" "c" "g"
## [15967] "a" "t" "a" "g" "t" "a" "a" "c" "t" "c" "g" "a" "a" "g" "a" "c" "a" "t"
## [15985] "c" "a" "c" "g" "c" "a" "g" "g" "t" "c" "g" "c" "t" "c" "c" "g" "c" "t"
## [16003] "a" "g" "g" "g" "g" "t" "c" "g" "t" "g" "g" "c" "a" "c" "t" "c" "t" "t"
## [16021] "g" "c" "a" "g" "c" "c" "t" "a" "t" "c" "c" "c" "t" "g" "g" "g" "t" "g"
## [16039] "a" "a" "c" "g" "g" "t" "a" "t" "c" "c" "c" "g" "g" "a" "t" "a" "a" "a"
## [16057] "t" "t" "g" "t" "t" "a" "g" "a" "t" "c" "a" "g" "c" "g" "c" "g" "g" "t"
## [16075] "t" "t" "t" "t" "g" "a" "a" "g" "t" "c" "a" "c" "g" "a" "a" "t" "t" "g"
## [16093] "c" "g" "a" "g" "g" "g" "c" "a" "c" "a" "a" "t" "c" "g" "g" "c" "t" "c"
## [16111] "g" "g" "c" "g" "t" "c" "a" "c" "g" "g" "a" "g" "c" "g" "c" "c" "a" "t"
## [16129] "t" "c" "a" "c" "g" "t" "g" "g" "a" "c" "t" "c" "a" "t" "g" "a" "c" "t"
## [16147] "a" "a" "c" "t" "g" "a" "t" "c" "a" "c" "c" "t" "g" "g" "g" "t" "c" "g"
## [16165] "c" "a" "t" "c" "g" "t" "c" "a" "a" "a" "c" "g" "g" "a" "t" "t" "c" "g"
## [16183] "g" "c" "c" "c" "a" "t" "a" "t" "c" "a" "g" "g" "g" "g" "t" "g" "t" "g"
## [16201] "c" "a" "c" "a" "t" "c" "g" "g" "a" "t" "g" "t" "c" "a" "a" "c" "a" "t"
## [16219] "c" "g" "a" "t" "c" "t" "c" "g" "c" "g" "a" "t" "a" "c" "a" "g" "c" "c"
## [16237] "t" "g" "t" "a" "c" "t" "c" "t" "g" "c" "t" "t" "g" "g" "t" "c" "a" "a"
## [16255] "a" "c" "g" "g" "a" "c" "t" "c" "g" "t" "t" "g" "c" "c" "t" "c" "c" "t"
## [16273] "c" "g" "c" "t" "t" "c" "t" "g" "a" "c" "g" "c" "a" "t" "g" "c" "g" "a"
## [16291] "t" "c" "a" "t" "a" "g" "c" "a" "t" "c" "t" "a" "g" "c" "g" "c" "c" "g"
## [16309] "a" "g" "t" "t" "a" "a" "t" "a" "c" "t" "a" "a" "c" "a" "g" "c" "a" "g"
## [16327] "t" "a" "c" "c" "t" "t" "a" "c" "a" "g" "t" "a" "t" "c" "g" "c" "g" "g"
## [16345] "c" "t" "t" "c" "g" "g" "c" "c" "c" "g" "t" "a" "t" "a" "c" "t" "a" "c"
## [16363] "c" "a" "g" "g" "g" "t" "c" "a" "c" "c" "t" "t" "c" "c" "t" "c" "g" "a"
## [16381] "a" "g" "c" "t" "t" "c" "t" "a" "a" "c" "c" "c" "c" "t" "t" "g" "g" "g"
## [16399] "g" "g" "g" "a" "a" "c" "g" "g" "t" "g" "a" "t" "g" "g" "t" "g" "a" "a"
## [16417] "c" "c" "c" "t" "a" "g" "g" "c" "g" "t" "a" "g" "g" "c" "c" "t" "g" "c"
## [16435] "t" "c" "c" "c" "c" "t" "a" "c" "t" "c" "c" "t" "t" "t" "c" "g" "c" "g"
## [16453] "c" "c" "a" "a" "t" "c" "a" "c" "g" "a" "a" "c" "a" "t" "c" "t" "t" "a"
## [16471] "c" "c" "t" "a" "g" "a" "a" "t" "c" "g" "c" "a" "c" "c" "c" "g" "c" "c"
## [16489] "c" "c" "t" "a" "a" "a" "a" "c" "c" "g" "g" "g" "g" "g" "c" "g" "g" "c"
## [16507] "a" "g" "t" "c" "t" "g" "g" "c" "g" "c" "c" "c" "c" "t" "t" "c" "c" "a"
## [16525] "a" "g" "c" "g" "t" "t" "g" "c" "g" "g" "a" "a" "t" "c" "c" "t" "a" "c"
## [16543] "g" "c" "c" "t" "a" "t" "c" "t" "c" "t" "t" "a" "c" "g" "g" "g" "g" "c"
## [16561] "t" "a" "c" "a" "g" "c" "a" "a" "a" "c" "g" "t" "c" "a" "c" "c" "a" "g"
## [16579] "c" "a" "c" "a" "c" "g" "a" "t" "g" "a" "c" "a" "t" "c" "g" "c" "g" "c"
## [16597] "g" "g" "a" "a" "t" "c" "g" "t" "c" "g" "a" "a" "c" "g" "c" "a" "c" "a"
## [16615] "c" "t" "c" "a" "a" "t" "c" "a" "t" "t" "g" "g" "a" "c" "a" "g" "t" "a"
## [16633] "t" "g" "c" "t" "c" "g" "t" "a" "g" "a" "c" "a" "a" "c" "t" "t" "t" "c"
## [16651] "g" "a" "a" "g" "g" "g" "g" "g" "t" "g" "c" "c" "t" "g" "a" "c" "c" "a"
## [16669] "g" "g" "c" "g" "a" "c" "c" "a" "t" "c" "c" "a" "a" "g" "a" "c" "a" "c"
## [16687] "t" "g" "g" "t" "g" "g" "g" "g" "a" "g" "c" "a" "t" "t" "t" "c" "c" "t"
## [16705] "g" "a" "t" "g" "g" "g" "g" "g" "t" "g" "t" "c" "t" "c" "a" "t" "t" "c"
## [16723] "g" "t" "t" "c" "a" "g" "a" "c" "c" "a" "a" "a" "t" "g" "a" "g" "g" "g"
## [16741] "a" "a" "c" "a" "c" "g" "g" "a" "t" "c" "c" "t" "c" "a" "c" "g" "c" "t"
## [16759] "g" "c" "a" "t" "a" "a" "c" "g" "c" "g" "t" "g" "t" "a" "g" "g" "t" "g"
## [16777] "a" "a" "c" "g" "g" "t" "c" "t" "t" "a" "a" "g" "g" "g" "a" "g" "t" "t"
## [16795] "t" "g" "g" "t" "a" "c" "g" "t" "t" "c" "g" "a" "a" "a" "a" "a" "a" "c"
## [16813] "a" "t" "g" "g" "t" "a" "c" "t" "c" "t" "c" "a" "g" "g" "a" "g" "a" "g"
## [16831] "t" "c" "g" "a" "t" "g" "g" "g" "c" "t" "a" "a" "a" "c" "g" "t" "a" "c"
## [16849] "t" "t" "a" "a" "a" "c" "c" "g" "c" "g" "c" "a" "a" "c" "c" "g" "a" "a"
## [16867] "t" "g" "t" "a" "g" "g" "a" "t" "c" "a" "c" "a" "c" "a" "a" "c" "g" "g"
## [16885] "g" "a" "g" "a" "g" "a" "a" "a" "g" "c" "c" "a" "g" "a" "c" "c" "t" "t"
## [16903] "g" "t" "a" "t" "a" "t" "g" "g" "g" "t" "c" "c" "c" "t" "g" "g" "g" "c"
## [16921] "t" "a" "a" "c" "c" "t" "g" "g" "a" "g" "g" "a" "g" "c" "g" "t" "t" "a"
## [16939] "a" "t" "g" "c" "c" "t" "t" "g" "c" "c" "c" "a" "t" "a" "c" "a" "g" "g"
## [16957] "g" "c" "c" "t" "g" "g" "g" "c" "t" "c" "g" "g" "c" "c" "g" "a" "a" "c"
## [16975] "a" "c" "g" "a" "g" "g" "g" "a" "c" "g" "c" "c" "a" "g" "c" "g" "t" "c"
## [16993] "c" "a" "c" "g" "g" "t" "a" "a" "a" "g" "g" "a" "g" "t" "g" "g" "t" "a"
## [17011] "c" "g" "a" "c" "g" "t" "c" "g" "t" "g" "c" "g" "c" "t" "a" "c" "a" "a"
## [17029] "t" "t" "c" "c" "c" "g" "a" "t" "t" "t" "g" "c" "g" "t" "a" "t" "c" "g"
## [17047] "t" "t" "g" "t" "g" "c" "c" "a" "c" "a" "a" "a" "t" "t" "t" "t" "a" "a"
## [17065] "c" "a" "g" "g" "t" "c" "g" "c" "a" "t" "c" "a" "g" "c" "c" "t" "a" "c"
## [17083] "a" "g" "t" "c" "c" "c" "t" "g" "g" "a" "a" "a" "a" "g" "g" "c" "c" "g"
## [17101] "g" "a" "c" "c" "c" "t" "c" "g" "g" "g" "g" "a" "c" "a" "t" "a" "c" "c"
## [17119] "c" "a" "c" "a" "c" "a" "a" "g" "c" "g" "a" "a" "g" "t" "a" "a" "c" "a"
## [17137] "a" "c" "g" "t" "g" "g" "c" "c" "t" "a" "c" "a" "g" "g" "c" "c" "a" "g"
## [17155] "c" "c" "g" "a" "c" "g" "a" "g" "g" "c" "c" "c" "c" "a" "t" "t" "g" "t"
## [17173] "a" "c" "a" "c" "c" "c" "t" "t" "t" "t" "g" "a" "g" "a" "g" "c" "c" "t"
## [17191] "a" "t" "g" "c" "c" "c" "c" "g" "c" "a" "g" "c" "g" "a" "g" "g" "c" "a"
## [17209] "c" "c" "c" "c" "g" "a" "c" "g" "g" "t" "t" "c" "c" "t" "c" "g" "g" "g"
## [17227] "c" "c" "t" "a" "c" "c" "c" "c" "c" "t" "c" "c" "c" "t" "c" "a" "a" "t"
## [17245] "t" "a" "c" "c" "t" "g" "t" "a" "g" "t" "c" "g" "g" "t" "t" "c" "g" "c"
## [17263] "a" "a" "g" "a" "t" "a" "c" "t" "g" "t" "c" "c" "g" "t" "a" "g" "a" "c"
## [17281] "g" "a" "t" "g" "c" "g" "c" "c" "a" "c" "g" "a" "c" "a" "c" "g" "a" "g"
## [17299] "a" "g" "g" "g" "g" "t" "a" "g" "t" "g" "g" "a" "t" "c" "g" "a" "t" "a"
## [17317] "c" "a" "c" "t" "g" "c" "c" "g" "t" "g" "g" "g" "a" "c" "g" "c" "t" "a"
## [17335] "g" "c" "a" "g" "a" "t" "c" "t" "c" "t" "g" "t" "g" "g" "a" "a" "a" "t"
## [17353] "g" "c" "g" "a" "a" "c" "t" "c" "a" "c" "t" "t" "a" "c" "a" "c" "g" "t"
## [17371] "a" "c" "a" "a" "c" "a" "g" "g" "c" "c" "c" "a" "a" "c" "g" "a" "t" "t"
## [17389] "a" "g" "g" "t" "c" "a" "a" "c" "g" "a" "t" "a" "c" "a" "c" "g" "t" "c"
## [17407] "a" "a" "c" "a" "c" "g" "c" "c" "a" "c" "g" "c" "a" "c" "t" "g" "c" "c"
## [17425] "c" "t" "a" "t" "t" "t" "g" "c" "t" "a" "t" "g" "c" "c" "c" "t" "c" "c"
## [17443] "a" "c" "g" "g" "g" "a" "t" "a" "g" "t" "a" "t" "c" "a" "a" "c" "c" "c"
## [17461] "a" "g" "a" "a" "g" "a" "c" "t" "g" "g" "g" "c" "c" "t" "a" "c" "g" "g"
## [17479] "c" "g" "a" "g" "c" "t" "t" "t" "g" "g" "a" "t" "t" "t" "g" "t" "c" "c"
## [17497] "a" "a" "t" "t" "c" "c" "a" "t" "t" "t" "a" "c" "a" "t" "c" "t" "a" "a"
## [17515] "c" "g" "t" "a" "a" "g" "g" "c" "t" "a" "g" "a" "a" "c" "g" "a" "c" "c"
## [17533] "g" "c" "g" "a" "g" "c" "t" "t" "t" "g" "c" "g" "c" "g" "c" "c" "t" "a"
## [17551] "a" "t" "t" "a" "g" "c" "g" "a" "t" "c" "g" "g" "a" "c" "a" "t" "a" "a"
## [17569] "t" "c" "g" "g" "t" "g" "c" "c" "a" "g" "c" "t" "a" "a" "g" "g" "g" "g"
## [17587] "g" "a" "t" "a" "c" "t" "t" "a" "a" "c" "c" "a" "c" "a" "c" "a" "a" "g"
## [17605] "c" "g" "c" "c" "c" "g" "g" "a" "t" "c" "c" "c" "g" "t" "t" "c" "c" "c"
## [17623] "g" "g" "g" "g" "a" "g" "c" "g" "g" "a" "g" "c" "g" "a" "a" "t" "t" "t"
## [17641] "g" "c" "t" "c" "a" "a" "g" "c" "a" "c" "c" "t" "g" "c" "t" "a" "c" "c"
## [17659] "a" "c" "a" "a" "t" "t" "a" "c" "t" "c" "c" "c" "c" "t" "t" "t" "g" "t"
## [17677] "t" "t" "t" "g" "a" "t" "g" "a" "t" "c" "t" "t" "g" "g" "c" "g" "g" "a"
## [17695] "c" "t" "a" "c" "t" "t" "g" "a" "g" "a" "c" "t" "g" "t" "g" "a" "g" "g"
## [17713] "g" "t" "c" "t" "t" "t" "t" "a" "g" "t" "c" "a" "g" "g" "t" "t" "c" "c"
## [17731] "t" "c" "c" "g" "g" "g" "g" "g" "a" "c" "a" "g" "a" "a" "c" "t" "a" "t"
## [17749] "c" "c" "a" "g" "c" "a" "t" "c" "g" "c" "c" "a" "t" "a" "t" "a" "a" "a"
## [17767] "t" "a" "a" "t" "c" "c" "c" "g" "g" "c" "g" "a" "g" "a" "g" "c" "a" "a"
## [17785] "t" "g" "g" "a" "g" "c" "c" "a" "c" "t" "a" "a" "c" "c" "g" "a" "a" "t"
## [17803] "c" "c" "a" "a" "g" "g" "a" "c" "g" "t" "a" "g" "g" "c" "a" "c" "g" "c"
## [17821] "t" "t" "g" "g" "t" "c" "a" "c" "g" "t" "t" "g" "a" "t" "c" "g" "c" "g"
## [17839] "t" "a" "g" "c" "t" "g" "a" "t" "t" "a" "c" "c" "c" "c" "g" "g" "t" "a"
## [17857] "a" "a" "c" "a" "c" "a" "g" "c" "c" "g" "a" "t" "a" "a" "a" "a" "g" "g"
## [17875] "g" "a" "g" "c" "c" "g" "g" "t" "a" "g" "a" "a" "a" "t" "a" "a" "c" "g"
## [17893] "g" "t" "g" "t" "c" "c" "c" "t" "c" "a" "g" "c" "t" "c" "c" "t" "t" "a"
## [17911] "c" "c" "g" "t" "c" "g" "g" "g" "a" "t" "t" "c" "a" "c" "c" "a" "t" "g"
## [17929] "t" "g" "c" "c" "c" "a" "t" "g" "a" "g" "t" "g" "t" "c" "a" "c" "g" "c"
## [17947] "c" "g" "c" "c" "c" "t" "c" "c" "t" "g" "g" "a" "g" "g" "c" "a" "t" "a"
## [17965] "c" "a" "t" "a" "c" "c" "c" "t" "g" "g" "g" "a" "c" "g" "a" "g" "g" "a"
## [17983] "a" "t" "a" "t" "g" "c" "c" "g" "c" "t" "a" "t" "g" "a" "g" "g" "t" "g"
## [18001] "a" "c" "c" "a" "g" "g" "g" "a" "c" "c" "t" "a" "a" "a" "c" "t" "a" "c"
## [18019] "g" "t" "g" "c" "a" "g" "a" "a" "g" "g" "c" "t" "a" "a" "a" "a" "t" "a"
## [18037] "c" "a" "c" "a" "a" "c" "g" "a" "c" "a" "g" "t" "c" "g" "g" "c" "a" "t"
## [18055] "t" "a" "a" "c" "t" "a" "a" "c" "g" "c" "t" "t" "c" "t" "c" "c" "g" "a"
## [18073] "c" "t" "c" "c" "g" "g" "g" "t" "a" "g" "c" "c" "c" "t" "c" "c" "a" "g"
## [18091] "a" "a" "t" "t" "c" "g" "g" "a" "a" "a" "c" "t" "c" "a" "g" "t" "a" "g"
## [18109] "c" "g" "c" "c" "g" "c" "c" "t" "g" "t" "c" "c" "c" "c" "a" "a" "c" "c"
## [18127] "c" "g" "c" "g" "a" "g" "g" "a" "g" "t" "c" "c" "t" "t" "g" "c" "c" "g"
## [18145] "t" "g" "g" "g" "g" "t" "c" "t" "a" "a" "g" "c" "c" "g" "a" "t" "a" "a"
## [18163] "g" "t" "a" "a" "t" "t" "c" "c" "t" "a" "g" "c" "a" "g" "a" "a" "c" "t"
## [18181] "g" "c" "c" "c" "g" "c" "c" "g" "a" "g" "a" "a" "g" "t" "a" "a" "c" "g"
## [18199] "c" "a" "c" "g" "a" "g" "t" "t" "a" "t" "c" "g" "a" "a" "a" "t" "c" "a"
## [18217] "g" "c" "g" "c" "t" "c" "g" "t" "g" "c" "g" "t" "a" "g" "a" "c" "t" "t"
## [18235] "c" "c" "a" "t" "c" "t" "a" "c" "g" "c" "a" "c" "c" "g" "a" "c" "c" "t"
## [18253] "a" "c" "c" "a" "a" "a" "g" "c" "a" "a" "g" "t" "t" "g" "c" "a" "g" "c"
## [18271] "g" "a" "t" "c" "g" "t" "t" "c" "t" "a" "a" "c" "g" "g" "t" "a" "a" "g"
## [18289] "a" "a" "c" "c" "a" "g" "a" "c" "t" "a" "c" "c" "a" "c" "c" "g" "t" "a"
## [18307] "g" "a" "t" "c" "a" "t" "c" "t" "t" "c" "g" "c" "g" "g" "c" "t" "g" "t"
## [18325] "g" "c" "a" "c" "c" "t" "c" "c" "t" "a" "a" "t" "t" "g" "c" "a" "g" "a"
## [18343] "c" "c" "g" "a" "a" "t" "g" "g" "g" "g" "c" "a" "a" "g" "g" "t" "g" "c"
## [18361] "g" "c" "g" "g" "t" "g" "c" "g" "g" "t" "a" "g" "g" "c" "t" "c" "c" "g"
## [18379] "g" "c" "t" "c" "g" "c" "g" "g" "t" "t" "a" "t" "t" "t" "t" "t" "a" "c"
## [18397] "g" "c" "c" "t" "g" "a" "a" "c" "g" "c" "g" "c" "t" "a" "g" "g" "c" "c"
## [18415] "c" "c" "t" "a" "g" "a" "c" "g" "g" "g" "t" "c" "c" "g" "g" "t" "a" "c"
## [18433] "c" "a" "a" "c" "g" "c" "c" "a" "c" "a" "t" "c" "a" "t" "c" "a" "c" "g"
## [18451] "g" "a" "a" "c" "c" "c" "c" "c" "c" "a" "c" "t" "g" "g" "a" "t" "g" "g"
## [18469] "c" "c" "t" "t" "t" "c" "t" "g" "g" "t" "g" "a" "c" "g" "a" "c" "t" "t"
## [18487] "c" "a" "g" "g" "c" "c" "t" "c" "c" "g" "c" "g" "t" "a" "a" "c" "g" "a"
## [18505] "a" "g" "g" "a" "g" "c" "g" "t" "g" "c" "a" "t" "g" "g" "c" "g" "a" "t"
## [18523] "a" "c" "g" "t" "g" "a" "a" "c" "g" "g" "c" "g" "g" "t" "a" "t" "a" "c"
## [18541] "g" "t" "g" "c" "a" "t" "c" "g" "g" "c" "a" "c" "g" "c" "g" "t" "t" "t"
## [18559] "a" "a" "a" "c" "g" "c" "g" "g" "c" "a" "t" "t" "g" "a" "g" "a" "g" "g"
## [18577] "g" "a" "c" "g" "c" "t" "g" "t" "c" "c" "a" "t" "a" "a" "c" "c" "t" "a"
## [18595] "t" "t" "t" "a" "a" "g" "t" "c" "c" "t" "a" "c" "g" "c" "t" "g" "t" "c"
## [18613] "c" "a" "t" "g" "g" "a" "g" "g" "g" "a" "a" "t" "c" "c" "a" "a" "g" "c"
## [18631] "a" "g" "t" "c" "t" "g" "t" "a" "c" "c" "g" "t" "g" "g" "c" "t" "c" "g"
## [18649] "a" "t" "a" "t" "g" "c" "t" "t" "c" "t" "c" "t" "a" "t" "t" "a" "c" "t"
## [18667] "a" "a" "t" "t" "c" "c" "g" "g" "a" "t" "t" "t" "a" "a" "c" "a" "t" "g"
## [18685] "a" "c" "a" "c" "g" "c" "g" "a" "a" "g" "a" "a" "c" "a" "c" "a" "g" "g"
## [18703] "c" "g" "c" "g" "c" "g" "c" "c" "a" "a" "a" "a" "t" "t" "c" "g" "t" "t"
## [18721] "t" "c" "g" "t" "g" "a" "g" "a" "c" "g" "t" "g" "g" "a" "a" "g" "t" "c"
## [18739] "g" "t" "t" "a" "t" "c" "c" "a" "g" "t" "g" "g" "c" "g" "a" "a" "c" "a"
## [18757] "g" "g" "a" "g" "c" "a" "c" "t" "t" "t" "t" "g" "a" "a" "c" "g" "c" "g"
## [18775] "t" "a" "c" "a" "g" "a" "g" "t" "c" "t" "c" "t" "t" "g" "c" "t" "g" "t"
## [18793] "a" "c" "c" "a" "c" "g" "c" "c" "t" "a" "g" "g" "g" "t" "a" "a" "c" "a"
## [18811] "t" "g" "g" "c" "a" "g" "g" "a" "t" "c" "g" "c" "a" "c" "g" "t" "c" "a"
## [18829] "c" "c" "g" "c" "t" "a" "a" "a" "g" "c" "c" "g" "c" "t" "a" "g" "t" "t"
## [18847] "a" "g" "g" "c" "t" "t" "g" "c" "c" "a" "t" "c" "g" "c" "a" "g" "t" "c"
## [18865] "c" "g" "g" "g" "a" "t" "g" "g" "c" "g" "c" "c" "g" "t" "t" "g" "a" "g"
## [18883] "c" "t" "t" "t" "c" "a" "g" "c" "c" "c" "t" "t" "t" "c" "g" "g" "a" "g"
## [18901] "g" "g" "g" "c" "t" "t" "g" "g" "a" "c" "a" "a" "c" "c" "c" "g" "t" "c"
## [18919] "c" "a" "c" "t" "a" "a" "a" "a" "t" "a" "c" "a" "g" "g" "a" "t" "a" "t"
## [18937] "t" "a" "g" "t" "c" "g" "g" "a" "t" "c" "a" "g" "t" "c" "a" "t" "a" "g"
## [18955] "t" "c" "g" "g" "c" "a" "t" "g" "c" "g" "a" "t" "t" "c" "a" "g" "t" "c"
## [18973] "t" "c" "c" "t" "a" "a" "g" "g" "t" "c" "t" "g" "c" "a" "g" "a" "c" "g"
## [18991] "t" "g" "t" "a" "c" "c" "g" "a" "c" "a" "g" "t" "a" "c" "t" "a" "g" "t"
## [19009] "t" "g" "a" "g" "c" "g" "c" "a" "a" "c" "c" "c" "g" "a" "g" "a" "g" "a"
## [19027] "t" "a" "c" "c" "a" "g" "g" "a" "c" "c" "g" "a" "t" "g" "t" "c" "a" "a"
## [19045] "c" "c" "t" "g" "t" "a" "t" "c" "t" "t" "t" "c" "c" "g" "g" "t" "g" "a"
## [19063] "t" "g" "a" "a" "a" "t" "c" "t" "c" "g" "g" "g" "g" "g" "a" "t" "a" "c"
## [19081] "t" "c" "g" "a" "a" "a" "g" "c" "g" "g" "c" "c" "c" "t" "g" "a" "t" "t"
## [19099] "g" "a" "c" "g" "g" "g" "c" "g" "c" "g" "t" "c" "a" "c" "a" "t" "a" "c"
## [19117] "c" "a" "g" "g" "c" "t" "a" "g" "a" "a" "g" "a" "a" "g" "t" "c" "a" "t"
## [19135] "g" "t" "c" "g" "t" "c" "t" "a" "a" "a" "a" "c" "c" "g" "g" "a" "a" "c"
## [19153] "c" "c" "t" "c" "a" "c" "c" "g" "g" "g" "c" "g" "c" "g" "g" "c" "c" "t"
## [19171] "g" "a" "c" "c" "g" "a" "g" "g" "g" "c" "t" "g" "g" "g" "c" "a" "c" "c"
## [19189] "t" "g" "a" "a" "g" "a" "g" "c" "t" "a" "t" "a" "c" "a" "g" "a" "g" "g"
## [19207] "t" "a" "g" "c" "a" "t" "c" "a" "a" "a" "t" "c" "g" "c" "a" "g" "a" "c"
## [19225] "t" "t" "t" "c" "g" "a" "g" "g" "c" "g" "t" "t" "g" "g" "t" "t" "a" "c"
## [19243] "c" "a" "g" "t" "a" "c" "g" "g" "c" "c" "g" "a" "c" "a" "a" "c" "g" "g"
## [19261] "c" "g" "c" "c" "t" "a" "t" "t" "a" "c" "a" "t" "t" "a" "a" "a" "c" "c"
## [19279] "g" "t" "c" "g" "g" "g" "c" "g" "g" "a" "g" "a" "a" "t" "c" "c" "g" "t"
## [19297] "t" "a" "c" "a" "t" "g" "g" "a" "a" "g" "c" "g" "g" "t" "c" "c" "a" "c"
## [19315] "c" "c" "c" "a" "c" "c" "g" "g" "g" "g" "c" "t" "g" "t" "a" "a" "g" "a"
## [19333] "c" "t" "a" "g" "g" "a" "c" "c" "c" "c" "t" "c" "t" "c" "c" "c" "c" "t"
## [19351] "a" "c" "t" "a" "a" "g" "t" "t" "c" "t" "a" "t" "t" "a" "c" "a" "g" "a"
## [19369] "t" "a" "t" "g" "a" "c" "c" "t" "a" "g" "a" "a" "c" "a" "c" "g" "g" "t"
## [19387] "a" "c" "g" "a" "c" "a" "g" "t" "c" "t" "g" "t" "c" "t" "g" "c" "a" "c"
## [19405] "c" "c" "c" "g" "a" "g" "t" "a" "t" "a" "t" "t" "t" "a" "g" "a" "a" "a"
## [19423] "g" "c" "g" "t" "c" "t" "c" "t" "c" "c" "c" "a" "t" "a" "a" "a" "t" "g"
## [19441] "a" "a" "a" "a" "c" "t" "c" "a" "g" "g" "g" "c" "t" "g" "g" "a" "g" "a"
## [19459] "a" "c" "t" "t" "a" "c" "g" "c" "c" "g" "a" "c" "t" "g" "g" "a" "a" "g"
## [19477] "g" "t" "a" "g" "a" "c" "a" "a" "t" "g" "a" "c" "g" "c" "a" "g" "c" "c"
## [19495] "g" "t" "a" "t" "a" "t" "g" "g" "t" "a" "c" "a" "g" "c" "c" "g" "t" "g"
## [19513] "c" "t" "c" "t" "a" "a" "a" "c" "g" "c" "t" "a" "t" "c" "a" "g" "g" "a"
## [19531] "c" "g" "c" "c" "g" "g" "c" "a" "c" "a" "t" "c" "a" "c" "c" "a" "g" "a"
## [19549] "c" "t" "c" "a" "c" "t" "a" "c" "c" "t" "t" "g" "g" "g" "t" "c" "t" "c"
## [19567] "a" "g" "a" "t" "g" "t" "c" "a" "t" "c" "a" "c" "a" "a" "c" "g" "a" "a"
## [19585] "c" "g" "a" "a" "t" "c" "c" "t" "a" "g" "a" "c" "g" "a" "g" "c" "a" "c"
## [19603] "c" "a" "c" "c" "t" "t" "c" "t" "c" "c" "a" "g" "c" "t" "c" "g" "a" "g"
## [19621] "a" "g" "t" "t" "t" "t" "a" "g" "a" "c" "a" "t" "a" "a" "g" "a" "g" "t"
## [19639] "t" "a" "c" "t" "g" "g" "a" "a" "g" "a" "a" "g" "t" "a" "g" "t" "t" "g"
## [19657] "t" "g" "t" "a" "c" "g" "g" "g" "a" "c" "g" "g" "g" "a" "t" "c" "g" "a"
## [19675] "t" "c" "g" "t" "a" "a" "g" "a" "a" "c" "a" "t" "t" "t" "a" "c" "t" "t"
## [19693] "c" "a" "c" "a" "t" "a" "c" "a" "c" "g" "t" "a" "a" "t" "t" "g" "c" "c"
## [19711] "t" "c" "t" "c" "c" "t" "c" "g" "a" "g" "c" "t" "g" "t" "c" "a" "g" "a"
## [19729] "c" "a" "g" "c" "a" "a" "g" "t" "a" "t" "g" "a" "a" "c" "t" "c" "t" "a"
## [19747] "c" "g" "a" "t" "c" "c" "a" "g" "a" "c" "c" "c" "a" "g" "a" "g" "c" "a"
## [19765] "a" "c" "t" "a" "a" "a" "a" "a" "a" "c" "a" "g" "t" "a" "c" "a" "g" "g"
## [19783] "t" "t" "a" "a" "c" "g" "a" "g" "c" "a" "c" "c" "g" "a" "g" "c" "c" "c"
## [19801] "a" "t" "c" "a" "g" "t" "a" "t" "c" "c" "g" "g" "c" "t" "c" "c" "g" "c"
## [19819] "c" "g" "g" "t" "a" "t" "t" "g" "c" "c" "c" "a" "a" "c" "c" "c" "t" "g"
## [19837] "a" "g" "t" "c" "c" "g" "c" "a" "t" "g" "a" "g" "a" "g" "a" "g" "a" "a"
## [19855] "g" "g" "g" "c" "t" "t" "a" "g" "g" "c" "c" "a" "t" "g" "c" "a" "c" "t"
## [19873] "g" "g" "a" "c" "c" "c" "c" "g" "g" "t" "t" "a" "c" "a" "a" "g" "t" "g"
## [19891] "c" "a" "c" "a" "a" "t" "a" "a" "a" "c" "c" "g" "a" "t" "t" "a" "t" "a"
## [19909] "t" "g" "a" "t" "c" "a" "t" "g" "c" "c" "c" "t" "t" "a" "a" "g" "g" "a"
## [19927] "g" "a" "c" "t" "a" "g" "c" "t" "g" "a" "c" "t" "a" "a" "g" "c" "a" "g"
## [19945] "c" "a" "g" "a" "g" "c" "c" "t" "c" "c" "g" "g" "a" "c" "g" "g" "g" "c"
## [19963] "a" "g" "g" "t" "g" "a" "c" "t" "g" "a" "g" "c" "t" "a" "a" "t" "t" "t"
## [19981] "c" "g" "t" "g" "t" "a" "a" "a" "a" "a" "a" "a" "g" "c" "c" "c" "c" "t"
## [19999] "c" "a"

3.3 Reliability of drifting Markov models

The reliability measures are computed by means of the following functions: A (availability), R (reliability,survival function), M (maintainability), errorRate (the classical failure rate, called BMP-failure rate and a more recent failure rate adapted to discrete data, called RG-failure rate). For more details on these reliability measures one can see Barbu and Limnios (2004); Barbu and Limnios (2008). All these functions have the following arguments:

x: An object of class dmm
k1: An integer, start position for the computation of the corresponding reliability measure
k2: An integer, end position for the computation of the corresponding reliability measure
s1: Character vector of the subspace working states (up-states) among the state space vector s.t. s1 < s
output_file: File containing the estimated/computed values of the corresponding reliability measure at each position of the selected frame
plot: display figure plot of the corresponding reliability measure along the positions if TRUE

The function errorRate enables to select the type of evaluated failure rate with an additional argument :

error.rate: Default="BMP", then BMP-failure-rate is the method used to estimate the error rate. If error.rate= "RG", then RG-failure rate is the method used to estimate the error rate.

These functions return vectors with the values of the corresponding reliability measure.

3.3.1 Estimation of reliability

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c("a","c","g","t"), init.estim ="freq", fit.method="sum")
reliability(dmm, k1=1, k2=10, upstates=c("a","c"),output_file = tempfile("drimmR"), plot=TRUE)
## [[1]]
##       positions  reliability
##  [1,]         1 0.4885571729
##  [2,]         2 0.2411445408
##  [3,]         3 0.1188995038
##  [4,]         4 0.0586239943
##  [5,]         5 0.0289048794
##  [6,]         6 0.0142517252
##  [7,]         7 0.0070269076
##  [8,]         8 0.0034646676
##  [9,]         9 0.0017082814
## [10,]        10 0.0008422825
## 
## [[2]]

3.3.2 Estimation of availability

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c("a", "c", "g", "t"), init.estim = "freq", fit.method="sum")
availability(dmm, k1=1, k2=10, upstates=c("a","c"), output_file = tempfile("drimmR"), plot=TRUE)
## [[1]]
##       positions availability
##  [1,]         1    0.4690808
##  [2,]         2    0.4761674
##  [3,]         3    0.4766474
##  [4,]         4    0.4766457
##  [5,]         5    0.4766422
##  [6,]         6    0.4766424
##  [7,]         7    0.4766428
##  [8,]         8    0.4766432
##  [9,]         9    0.4766436
## [10,]        10    0.4766440
## 
## [[2]]

3.3.3 Estimation of maintainability

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c("a","c","g","t"), init.estim="freq", fit.method="sum")
maintainability(dmm, k1=1, k2=10, upstates=c("a","c"), output_file = tempfile("drimmR"), plot=TRUE)
## [[1]]
##       positions maintainability
##  [1,]         1       0.4885572
##  [2,]         2       0.7164935
##  [3,]         3       0.8485148
##  [4,]         4       0.9189865
##  [5,]         5       0.9566753
##  [6,]         6       0.9768307
##  [7,]         7       0.9876094
##  [8,]         8       0.9933737
##  [9,]         9       0.9964564
## [10,]        10       0.9981049
## 
## [[2]]

3.3.4 Estimation of failure rates

data(lambda, package = "drimmR")
dmm <- fitdmm(lambda, 1, 1, c("a","c","g","t"), init.estim="freq", fit.method="sum")
failureRate(dmm, k1=1, k2=10, upstates=c("a","c"), failure.rate="BMP", output_file=tempfile("drimmR"), plot=TRUE)
## [[1]]
##       positions       BMP
##  [1,]         1 0.5064149
##  [2,]         2 0.5069368
##  [3,]         3 0.5069450
##  [4,]         4 0.5069446
##  [5,]         5 0.5069440
##  [6,]         6 0.5069434
##  [7,]         7 0.5069428
##  [8,]         8 0.5069422
##  [9,]         9 0.5069416
## [10,]        10 0.5069410
## 
## [[2]]

4 Concluding remarks

To conclude, in this paper we have presented drimmR, an R package for simulation, estimation and reliability and survival analysis of drifting Markov models. These are versatile stochastic models of Markov type capable of taking into account a time non-homogeneity of a known, controlled shape. For this reason these models can represent interesting modeling alternatives to classical models (like Markov models, semi-Markov models, etc.) and can be useful for researchers, practitioners and engineers in various fields.

The package addresses several items important from practical point of view, when carrying out the estimation: we consider one or several samples, the sample paths can be complete or not, can come from models of the same length or of different lengths.

The fields of application of the DMMs can be numerous; we only want to point out three of them.

Survival analysis and reliability theory: note that important indicators like reliability or survival function, availability, maintainability and failure rates are computed in our package;
Bioinformatics: note that important quantities for OMICS data in general, DNA analysis in particular, are computed by the proposed package. Thus we have the computation of the probabilities and expectations of appearance of a given word along the sequence.

Flowchart of the drimmR package.

Barbu, Vlad Stefan, Caroline Bérard, Dominique Cellier, Mathilde Sautreuil, and Nicolas Vergne. 2018. “SMM: An r Package for Estimation and Simulation of Discrete-Time Semi-Markov Models.” R Journal 10 (2).

Barbu, Vlad Stefan, Geoffray Brelurut, Annthomy Gilles, Arnaud Lefebvre, Victor Mataigne, and Nicolas Vergne. 2018. “DRIMM’r and WebDRIMM: An r Package and a Web Interface for Drifting Markov Model Estimation and Reliability.” JDS (Paris, France). https://normandie-univ.hal.science/hal-02337261.

Barbu, Vlad Stefan, and Nikolaos Limnios. 2004. “Discrete Time Semi-Markov Processes for Reliability and Survival Analysis - a Nonparametric Estimation Approach.” In Parametric and Semiparametric Models with Applications to Reliability, Survival Analysis and Quality of Life, edited by M. Nikulin, N. Balakrishnan, M. Mesbah, and N. Limnios. Statistics for Industry and Technology. Birkhäuser, Boston. https://doi.org/10.1007/978-0-8176-8206-4_30.

Barbu, Vlad Stefan, and Nikolaos Limnios. 2008. “Reliability of Semi-Markov Systems in Discrete Time: Modeling and Estimation.” In Handbook on Performability Engineering, edited by K. B. Misra. Springer. https://doi.org/10.1007/978-1-84800-131-2_24.

Barbu, Vlad Stefan, and Nicolas Vergne. 2019. “Reliability and Survival Analysis for Drifting Markov Models: Modeling and Estimation.” Methodology and Computing in Applied Probability 21 (4): 1407–29.

Bulla, Jan, Ingo Bulla, and Oleg Nenadić. 2010. “Hsmm — An R Package for Analyzing Hidden Semi-Markov Models.” Computational Statistics & Data Analysis, Second Special Issue on Statistical Algorithms and Software, vol. 54 (3): 611–19. https://doi.org/10.1016/j.csda.2008.08.025.

Król, Agnieszka, and Philippe Saint-Pierre. 2015. “SemiMarkov : An r Package for Parametric Estimation in Multi-State Semi-Markov Models.” Journal of Statistical Software 66 (6). https://doi.org/10.18637/jss.v066.i06.

O’Connell, Jared, and Søren Højsgaard. 2011. “Hidden Semi Markov Models for Multiple Observation Sequences: The Mhsmm Package for r.” Journal of Statistical Software 39 (4). https://doi.org/10.18637/jss.v039.i04.

Vergne, Nicolas. 2008. “Drifting Markov Models with Polynomial Drift and Applications to DNA Sequences.” Statistical Applications in Genetics and Molecular Biology 7 (1). https://doi.org/10.2202/1544-6115.1326.

drimmR : R Package for Estimation, Simulation and Reliability of drifting Markov Models