 and finding properties of $\aleph$.]{} From assumption (A1), if a nucleotide is proposed as a candidate for joining the end of one of the strands, then its complement is automatically a candidate for joining the beginning of the other strand. Therefore, the probability of a nucleotide being selected as a candidate for joining one strand is the same as the probability of its complement being selected for the other strand, that is,\
[**($H1$)**]{} $\mu_A=\mu_T$ and $\mu_C=\mu_G$. In other words, $\mu$ takes the form $\mu(m)=(m,0.5-m,0.5-m,m)$, where $0\leq m \leq
0.5$.\
Now, since the matrix $\aleph$ is assumed to be the same for all genomes, any sequence produced by the model is a realization of a Markov chain belonging to a family of Markov chains parameterized by $m$ and $t$, where $0< m< 0.5$ and $0\leq t\leq 1$.

Further, from (A2) we have that the probability of a candidate nucleotide of type $j$ being accepted to follow a nucleotide of type $i$ at the end of one strand is equal to the probability of a nucleotide of type $\alpha(i)$ preceding a nucleotide of type $\alpha(j)$ at the beginning of the other strand, that is, the matrix $\aleph$ has the form\
[**($H2$)**]{} $a_{ij}=a_{\alpha(j)\alpha(i)}$.\
We use the [*lsqcurvefit*]{} function of the optimization toolbox in [Matlab]{} to construct estimators for the matrices $\aleph(n)$, vectors $\mu(n)$ and values $t(n)$. The [*lsqcurvefit*]{} function solves (in the least-squares sense) the nonlinear problem of determining the parameters for which the right side of equation most closely approximates $P(n)$. Firstly, we carried out free estimation without the [*a priori*]{} assumption of the properties ($H1$) and ($H2$).

The solution to the optimization problem is sensitive to the value of $t(n)$ chosen to initialize the optimization algorithm. Since previous tests suggest that neither strand is significantly favored over the other during the construction of the sequence, we used $t(n)=0.5
and finding properties of $ \aleph$. ] { } From assumption (A1), if a nucleotide is proposed as a campaigner for join the end of one of the strands, then its complement is mechanically a candidate for join the beginning of the other strand. consequently, the probability of a nucleotide being selected as a candidate for joining one chain is the same as the probability of its complement being selected for the early strand, that is,\ 
 [ * * ($ H1 $) * * ] { } $ \mu_A=\mu_T$ and $ \mu_C=\mu_G$. In other words, $ \mu$ takes the form $ \mu(m)=(m,0.5 - m,0.5 - m, m)$, where $ 0\leq m \leq 
 0.5$.\ 
 nowadays, since the matrix   $ \aleph$ is assumed to be the same for all genomes, any sequence produced by the exemplar is a realization of a Markov chain belonging to a family of Markov chain parameterized by   $ m$ and   $ t$, where $ 0 < m < 0.5 $ and $ 0\leq t\leq 1$. 

 Further, from (A2) we have that the probability of a candidate nucleotide of type   $ j$ being accept to follow a nucleotide of type   $ i$ at the end of one strand is equal to the probability of a nucleotide of character $ \alpha(i)$ preceding a nucleotide of type $ \alpha(j)$ at the beginning of the other strand, that is, the matrix $ \aleph$ has the form\ 
 [ * * ($ H2 $) * * ] { } $ a_{ij}=a_{\alpha(j)\alpha(i)}$.\ 
 We use the [ * lsqcurvefit * ] { } function of the optimization toolbox in [ Matlab ] { } to construct calculator for the matrix $ \aleph(n)$, vectors   $ \mu(n)$ and values $ t(n)$. The [ * lsqcurvefit * ] { } function solves (in the least - squares sense) the nonlinear problem of determining the parameters for which the right slope of equation most close approximate   $ P(n)$. first, we carried out free estimation without the [ * a priori * ] { } assumption of the properties ($ H1 $) and ($ H2 $). 

 The solution to the optimization problem is sensitive to the value of   $ t(n)$ chosen to format the optimization algorithm. Since previous tests suggest that neither strand is significantly favored over the other during the construction of the sequence, we used $ t(n)=0.5
 anf finding properties of $\xleph$.]{} From assumption (E1), if a hucleotiae is proposed as a candidatx foe joibing the end of one of the strajds, then its xomplement is automaticampy a randidate for jpining the beginning of dhd lther strand. Therefore, the probabilyty of s jucleotide beigg stlested as a candidate for joining one stdand is the same as yhe probability of its comolemfnt being selected for the otyer frrand, that ir,\
[**($H1$)**]{} $\mu_A=\mu_T$ and $\mu_C=\mu_G$. Ih other words, $\mu$ takes the form $\mu(m)=(m,0.5-k,0.5-m,m)$, where $0\oew m \neq
0.5$.\
Now, sincx the iatrix $\aleph$ lx assuked to ne the same fov all gebomes, any sequence prmduced by the modej is a redlnzation of a Markov cyaun benonghng gi a fajiky of Magkot chains padameterized by $m$ and $t$, where $0< m< 0.5$ age $0\leq t\leq 1$.

Fudther, srjm (A2) we have that the probability of a bandjdate nucleotide of typw $j$ being accepted to vollow a gucleotide of type $i$ at the end of one strand is exual vo tht ivobacullty of a nucleotide of type $\alpha(i)$ preceding w nickeotide of typc $\alpha(j)$ at the bebijnogg of the othgr straus, fhat is, the matrix $\aleph$ ras tye form\
[**($H2$)**]{} $w_{ij}=a_{\slpha(j)\alpha(i)}$.\
We use the [*lsqcyrvefit*]{} funcnion of the optimizatiln toolbox nn [Matkab]{} tp construct estimators yor ths matrices $\wleph(n)$, vedgors $\mu(n)$ and valuds $n(n)$. Tve [*lsqcurvefit*]{} function sojves (in tie lecst-squards sgnse) thq nonlineag pronnem of determining the karameders for wjich the right side of equation most closely sp[roqimates $P(n)$. Firsbly, we carried jut free estimction wichout ghe [*a priogi*]{} assumpvion of the [roperties ($H1$) djd ($H2$).

The soluvion to tre optimuzation oroblem is senxitive to the value of $t(n)$ chosen to inltialkae the optimizacnob algorithm. Sinve oredilux easts suggest thag ndother strand is wkgnigicantly favored oves ths other during the cjnstructuon of tre sequence, wr used $t(n)=0.5
and finding properties of $\aleph$.]{} From assumption a is proposed a candidate for of strands, then its is automatically a for joining the beginning of the strand. Therefore, the probability of a nucleotide being selected as a candidate for one strand is the same as the probability of its complement being selected the strand, is,\ $\mu_A=\mu_T$ and $\mu_C=\mu_G$. In other words, $\mu$ takes the form $\mu(m)=(m,0.5-m,0.5-m,m)$, where $0\leq m \leq 0.5$.\ since the matrix $\aleph$ is assumed to be same for all genomes, sequence produced by the model a of a chain to family of Markov parameterized by $m$ and $t$, where $0< m< 0.5$ and $0\leq t\leq 1$. Further, from (A2) we that the a candidate of $j$ accepted to follow of type $i$ at the end is equal to the probability of a nucleotide type $\alpha(i)$ a nucleotide of type $\alpha(j)$ at beginning of the other strand, that is, the $\aleph$ has the form\ [**($H2$)**]{} $a_{ij}=a_{\alpha(j)\alpha(i)}$.\ We use the [*lsqcurvefit*]{} function of the optimization toolbox to construct estimators for matrices $\aleph(n)$, vectors and $t(n)$. [*lsqcurvefit*]{} solves (in least-squares sense) the nonlinear problem of determining the parameters for which right side of equation most closely approximates $P(n)$. Firstly, we free without the [*a assumption of the properties and The solution to the is to $t(n)$ to the optimization algorithm. Since tests suggest that neither strand significantly favored over the the sequence, we used $t(n)=0.5
 and finding properties of $\alePh$.]{} From assuMptioN (A1), iF a nUcLeotIde iS proposed as a caNDidaTe for joining the end of onE of thE sTRandS, ThEn its ComplemENt IS AutOmAtIcaLlY A cAndidAte For joinIng the begiNniNg Of the other stRAnD. Therefore, The Probability oF a nUcleotIdE beINg selEctEd as a CandidATe for jOining one StRAnd is tHE same as THE pRobaBility of its compleMEnT Being selected fOr the oThER sTRAnd, ThaT is,\
[**($H1$)**]{} $\mu_A=\mu_T$ AnD $\mu_C=\mU_g$. In otheR WoRDS, $\Mu$ tAKes the form $\mu(m)=(M,0.5-m,0.5-m,m)$, where $0\leQ M \leQ
0.5$.\
Now, siNcE thE Matrix $\Aleph$ Is ASsuMed to be the sAme fOr all genoMes, any SEquence PRoduced By the mOdeL is A reaLIzAtIon Of A marKOv ChaIN beLonging tO a FaMily oF MarKOV CHainS paRameTerizEd by $m$ and $t$, wherE $0< m< 0.5$ aNd $0\leQ T\leQ 1$.

FurtHer, frOm (A2) wE hAve thAt the pRobabIlIty of a candidate NuclEotide of tYpe $J$ bEinG aCceptED to folLow A nuCleotidE of type $I$ At tHe END Of One strand is equal to ThE PRoBability Of a nucLEoTiDE of type $\aLpHa(i)$ PrecEDIng a nUcleOTiDe of type $\Alpha(j)$ AT tHe BeginniNg Of the oThEr sTraNd, thaT Is, thE matriX $\aleph$ haS the fORm\
[**($H2$)**]{} $a_{ij}=a_{\alpha(j)\aLPha(i)}$.\
We use the [*lSQcURVeFIt*]{} fuNctIon of the optImizATion ToolBOx In [MATlab]{} tO consTrUCt EStimators for the matrIcEs $\alepH(n)$, vecTors $\mu(n)$ and valUes $t(n)$. The [*lsQCURvefit*]{} fuNctiON sOLves (in the least-SquarEs sense) the NOnlinear ProblEm of deteRmining thE PArameterS foR whIch The RIGhT side of equatiON Most ClOsely apProXimates $p(n)$. FIrsTly, We cArRied out frEe estimaTiOn WiThOut The [*a pRIori*]{} assuMpTioN oF thE propERties ($H1$) And ($H2$).

THe soLuTiON to The optiMIzATIon pRoBlEm is SenSiTive tO the VAluE of $t(n)$ chOsen to iniTiaLIze tHe OpTimizatIon algorithm. SInCe previous TeSts SuggesT THat neithEr strand is significantly FAvored oVer The otHer dUring the cOnsTructiOn oF The seqUence, wE used $T(n)=0.5
 and finding properties of  $\aleph$. ]{} F rom  as su mpti on ( A1), if a nucl e otid e is proposed as a can didat ef or j o in ing t he endo fo n e o fth e s tr a nd s, th enits com plement is  au to matically ac an didate for  jo ining the be gin ning o fthe other  st rand.  There f ore, t he probab il i ty ofa  nucleo t i de  bei ng selected as ac an d idate for join ing on es tr a n d i s t he same as  t he pr o babilit y  o f i tsc omplement bei ng selected for  the o th ers trand,  that  i s ,\[**($H1$)** ]{}$\mu_A=\m u_T$ a n d $\mu_ C =\mu_G$ . In o the r w ords ,  $ \m u$ta k est he  fo r m $ \mu(m)=( m, 0. 5-m,0 .5-m , m ) $ , wh ere  $0\ leq m  \leq
0.5$.\Now , si n cethe m atrix  $\a le ph$ i s assu med t obe the same for  all  genomes,  an yseq ue nce p r oduced  by  th e model  is a r e ali za t i o nof a Markov chainbe l o ng ing to a  famil y  o fM arkov ch ai nspara m e teriz ed b y  $ m$ and $ t$, wh e re  $ 0< m< 0 .5 $ and$0 \le q t \leq1 $.
Furthe r, from(A2)w e have that th e  probabilityo fa ca n dida tenucleotideof t y pe $ j$ b e in g a c cepte d tofo l lo w  a nucleotide of ty pe  $i$ a t the  end of one s trand is e q u a l to the  pro b ab i lity of a nucl eotid e of type$ \alpha(i )$ pr ecedinga nucleot i d e of typ e $ \al pha (j) $ at  the beginnin g of t he  otherstr and, th atis,  th e m at rix $\ale ph$ hasth efo rm \
[ **($H 2 $)**]{}$a _{i j} =a_ {\alp h a(j)\a lpha( i)}$ .\ 
W e  us e the [ * ls q c urve fi t* ]{}fun ct ion o f th e  op timizat ion toolb oxi n [M at la b]{} to  construct es ti mators for  t hematric e s  $\aleph (n)$, vectors $\mu(n)$a nd valu es$t(n) $. T he [*lsqc urv efit*] {}f unctio n solv es (i nthe l east- s q ua res  s ense) then o nli nearpr oble m of de termining the para m ete rs for whichthe  rig h t  s ide of equ at i onm o st closely appr oximates $ P( n )$ . Firstly, weca rried o ut free  esti m ation w ithout th e [*a pri or i*]{ } ass umption of  the pro perties ( $ H1$)a nd  ($H2 $). 

Theso lut ion t o theo pti mizat ion pr ob lem is  sens it ive to t he value of $t(n)$ chos en toiniti ali ze the op tim i zat ion algor ithm . Since pr evi ous  test s s u ggest  tha t  n eit h er st rand is signif i ca ntl y fa vored overt h e  ot her d uri n g thecons truction of the s e quence, we use d $t ( n )=0 .5
 and_finding properties_of $\aleph$.]{} From assumption_(A1), if_a_nucleotide is_proposed_as a candidate_for joining the_end of one of_the strands, then_its_complement is automatically a candidate for joining the beginning of the other strand. Therefore,_the_probability of_a_nucleotide_being selected as a candidate_for joining one strand is_the same_as the probability of its complement being selected_for_the other strand,_that is,\
[**($H1$)**]{} $\mu_A=\mu_T$ and $\mu_C=\mu_G$. In other words, $\mu$_takes the form $\mu(m)=(m,0.5-m,0.5-m,m)$, where $0\leq_m \leq
0.5$.\
Now, since_the_matrix $\aleph$_is assumed to be_the same for all genomes, any_sequence produced by the model is_a realization of a Markov chain belonging_to a family of Markov chains_parameterized by $m$ and $t$, where $0<_m< 0.5$_and $0\leq t\leq 1$.

Further, from_(A2) we have_that the_probability of a_candidate nucleotide of type $j$ being accepted_to follow a_nucleotide of type $i$ at the end_of_one strand is_equal_to_the probability_of a nucleotide_of_type $\alpha(i)$_preceding_a nucleotide of type $\alpha(j)$ at_the_beginning of the other strand, that is,_the matrix $\aleph$ has_the_form\
[**($H2$)**]{} $a_{ij}=a_{\alpha(j)\alpha(i)}$.\
We use the_[*lsqcurvefit*]{} function of the optimization_toolbox in [Matlab]{} to construct estimators_for the_matrices $\aleph(n)$,_vectors $\mu(n)$ and values $t(n)$. The [*lsqcurvefit*]{} function solves (in the least-squares_sense) the nonlinear problem of determining_the parameters for which_the right_side_of equation most_closely_approximates $P(n)$. Firstly,_we carried out free estimation without the_[*a priori*]{}_assumption of the properties ($H1$) and_($H2$).

The solution to the_optimization_problem is sensitive to the value_of $t(n)$ chosen to initialize the optimization_algorithm. Since previous tests suggest_that_neither_strand is significantly favored over_the other during the construction of_the sequence, we_used $t(n)=0.5