<title>Introduction</title> <p>A fundamental task in sequence analysis is calculating the probability of a multiple alignment given a phylogenetic tree relating the sequences and an evolutionary model describing how sequences change over time. This quantity is already at the heart of phylogenetic tree inference by either maximum likelihood ##REF##7288891##[1]## or Bayesian approaches ##UREF##0##[2]##–##UREF##1##[5]##. It is also desirable to integrate evolutionary modeling into the probabilistic models widely used for other sequence analysis problems such as HMMs and SCFGs for homology search and genefinding ##UREF##2##[6]##. However, whereas HMMs, SCFGs, and other sequence analysis models account for insertion and deletion events in their probabilistic framework, many widely used phylogenetic models only account for residue substitution events.</p> <p>The general approach of modeling residue substitution as a continuous-time Markov process was introduced by the Jukes-Cantor model of nucleotide substitution in DNA ##UREF##3##[7]## and the Dayhoff pam model of amino acid substitution in proteins ##UREF##4##[8]##. It has since been extensively developed both for nucleotides ##REF##7463489##[9]##–##REF##8583911##[14]## and amino acids ##REF##8913301##[15]##–##REF##15064085##[18]##, and extended to models of more than a single residue, such as codon to codon substitutions ##REF##8583885##[19]##,##REF##10790415##[20]## and RNA basepair to basepair substitutions ##REF##10383470##[21]##–##REF##15261152##[23]##.</p> <p>Given a substitution model and a tree, one can efficiently calculate the probability of an <italic>ungapped</italic> multiple alignment using Felsenstein's peeling algorithm ##REF##7288891##[1]##. The Felsenstein algorithm scales linearly with the length of the alignment and the number of sequences. Because of its economy, it is the basis of maximum likelihood methods in many practical phylogenetic inference tools, including <sc>phylip</sc>
##UREF##7##[24]##, <sc>paup</sc>* ##UREF##8##[25]##, and others ##UREF##9##[26]##–##REF##15608047##[33]##. The Felsenstein algorithm is readily integrated with other probabilistic models for ungapped alignment analysis, including HMMs ##REF##8583911##[14]##, ##REF##7713447##[34]##–##REF##15285899##[38]## and SCFGs ##REF##10383470##[21]##. However, when this approach is applied to gapped multiple sequence alignments, gap characters are typically treated as missing data (an unknown residue), effectively equivalent to ignoring them.</p> <p>A variety of more formal approaches for treating insertions and deletions exist ##UREF##13##[39]##–##REF##17110370##[49]##. The canonical model in this active area of research is the Thorne-Kishino-Felsenstein model (<sc>tkf</sc>91) ##REF##1920447##[50]##. <sc>tkf</sc>91 treats insertion and deletion events as a continuous-time process governed by explicit insertion and deletion rate parameters, allowing multiple insertions and deletions to accumulate at the same “site” over long times. <sc>tkf</sc>91 improved, for example, upon methods that parsimoniously assume no more than one change per site per branch, including the pioneering Bishop and Thompson pairwise alignment likelihood model that preceded <sc>tkf</sc>91 ##REF##3641921##[51]##. Many extensions of <sc>tkf</sc>91 have appeared ##REF##1556741##52##–##REF##14694074##56##, including practical applications for pairwise alignment ##REF##11590097##[57]##,##REF##14529629##[58]## and multiple alignment ##UREF##15##[59]##–##REF##16085574##[62]##. Several approaches based on <sc>tkf</sc>91 or related models have addressed the problem of simultaneously aligning and inferring the phylogeny of a group of related sequences ##REF##10964574##[63]##–##REF##15804354##[68]##. In general all algorithms in this class have difficult time complexities, worst-case exponential in the number of sequences, and at least in the case of parsimony models akin to <sc>tkf</sc>91 , the problem of inferring the optimal insertion and deletion history given a tree has been formally shown to be NP-complete ##REF##16960972##[47]##. This complexity is inherent to the problem. Any implementation of evolutionary models that allow insertions and deletions, <sc>tkf</sc>91-based or otherwise, seeks to make approximations that make calculations tractable.</p> <p>Here we are specifically concerned with the problem of calculating the probability of a <italic>given</italic> multiple alignment and phylogeny. A <sc>tkf</sc>91-based approach for this problem ##UREF##19##[69]## used a clever approach of using ungapped columns in the alignment to constrain and subdivide the solution space (a so-called “homology structure”), but (because <sc>tkf</sc>91 is not invariant under column rearrangement) a sum over all subalignments compatible with a given homology structure is still required. As a result, though time complexity exponential in the number of sequences can be avoided in the average case where not many gaps occur, the approach remains expensive in absolute terms; a Bayesian Markov chain Monte Carlo phylogenetic inference for an alignment of 10 globins was reported to require 3 CPU hours on a 1.25 GHz G4 Apple Macintosh ##UREF##19##[69]##. Although <sc>tkf</sc>91 has many desirable and realistic properties as a model of evolution, it would be advantageous to have even more computationally efficient approaches, particularly for problems where a gain in efficiency might outweigh sacrificing some of the realism of the model.</p> <p>An alternative is to use a continuous-time Markov process for insertions and deletions that remains fully compatible with the Felsenstein algorithm ##REF##11264399##[41]##,##REF##15780137##[46]##. This requires an unrealistic assumption of column independence, but nonetheless, it may be better than ignoring gaps altogether. Here we further explore such models. We propose a non-reversible generative model based on a Markov process that we readily incorporated into an existing phylogenetic inference application, resulting in a gain in its accuracy.</p>
<title>Materials and Methods</title> <p>The C source code for the modified <sc>phylip</sc> 3.66 package ##REF##8583911##[14]## that contains the program <sc>dnaml</sc>
<italic>ε</italic> , the C source code for evolving sequences with the generative model (<italic>ε</italic>
<sc>rate</sc> ), the modified <sc>rose</sc> package (version 1.3) ##REF##9545448##[76]##, as well as all the Perl scripts and datasets used to generate the results presented in this paper are provided as a tarball in ##SUPPL##0##Dataset S1##. The program <sc>dnaml</sc>
<italic>ε</italic> uses the <sc>easel</sc> sequence analysis library (SRE, unpublished) which is also provided.</p>
<title>Results</title> <title>Premises of the Model</title> <p>The originating idea is to use an extended (<italic>K</italic>+1)×(<italic>K</italic>+1) rate matrix for <italic>K</italic> residues (4 nucleotides or 20 amino acids) plus the gap character to describe the rate of change of a residue to a residue (substitution), a residue to a gap (deletion), and a gap to a residue (insertion) ##UREF##2##[6]##,##REF##11264399##[41]##,##REF##15780137##[46]##. From this Markov process, we construct a generative model of sequence evolution that includes insertions and deletions. Several consequences flow from this, which here we discuss informally by way of introduction to the rest of the paper.</p> <p>This model describes the evolution of single residues in one column of a multiple alignment, given a phylogenetic tree. The total probability of the alignment is then assumed to be an independent product of each column probability. For each column, a variant of the Felsenstein peeling algorithm recursively infers the probability of ancestral characters at each tree node, where an ancestral character is either a residue or a gap. Assuming that insertion and deletion events happen one residue at a time necessarily implies a linear “gap cost”. This is a much less satisfactory model of insertion and deletion processes than models that can assume an affine or arbitrary gap cost.</p> <p>It has generally been thought that models based on a gap-extended rate matrix must be conceptually flawed, because it appears necessary to assume that all ancestral and descendant sequences fit in a fixed number of columns. This fundamentally conflicts with allowing any number of insertions and deletions to occur, and it produces a so-called “memory effect” artifact ##UREF##20##[70]## in which descendants “remember” how many gap characters were present in ancestral sequences, allowing insertions up to that length and precluding longer ones. A related conceptual flaw would be treating gaps like residues, assigning a probability to a gap/gap alignment (as one would do for any residue/residue alignment), rather than recognizing that a gap/gap alignment may represent no evolutionary event at all from the standpoint of just the descendant and ancestral sequence; rather, gap/gap alignments are imposed by events that occurred in <italic>other</italic> sequences. For a model to be at all satisfactory, one must be able to describe a generative evolutionary model <italic>unconditional</italic> on any fixed sequence length, in terms of substitution, insertion and deletion events that evolve one (unaligned) sequence to another, and show how that generative process relates uniquely to the column-by-column inference algorithm that one will apply to a given multiple alignment.</p> <p>Here we will develop such a generative model, by borrowing terms from a Markov process for a rate matrix extended for the gap character. The rate matrix includes the gap character, but the subsequent generative evolutionary model does not treat gaps as an extra residue. Rather, it describes evolutionary insertion and deletion (birth-death) events ##REF##16589686##[71]##, where the evolved sequences form alignments of arbitrary length. The existence of a generative evolutionary model for unaligned sequences, and the mapping of its events to the column-by-column inference procedure, is the crucial point of differentiation between our work and previous work on efficient column-based phylogenetic inference with gap-extended rate matrices ##UREF##2##[6]##,##REF##11264399##[41]##.</p> <p>Similarly, one must have a consistent way of dealing with columns that are unobserved in the alignment of extant sequences – that is, places where ancestral residues have been inserted and deleted, where alignment columns would exist if all the ancestral sequences were known in addition to just the extant sequences. Therefore the likelihood we calculate for an extant multiple alignment will be its marginal likelihood, marginalized over all possible ancestral sequences including unobserved alignment columns that left no trace in the observed alignment.</p> <p>Another conceptual problem of column-based models arises if one adopts the usual practice of making the substitution process reversible (in the sense that the probability of an ancestor/descendant sequence alignment is independent of the direction of time along the branch that connects them; this is mathematically convenient for applying the “pulley principle” and using the Felsenstein peeling algorithm on unrooted trees). If one assumes reversibility for a substitution process that includes gaps, one necessarily imposes a frequency of gap characters that is constant with respect to divergence time ##REF##11264399##[41]##; but obviously in the limit of zero divergence time, there are no gap characters in an alignment of homologous sequences. Moreover, reversibility clearly cannot hold if insertion and deletion rates are free parameters. For example, for a nonzero insertion rate and zero deletion rate, the shorter of two homologous sequences automatically must be the ancestor. A reversible insertion/deletion model would also imply that all homologous sequences have the same expected length at all divergence times. We overcome these problems by adopting a non-reversible model that states an explicit prior length distribution for unaligned ancestral sequences (as opposed to assuming that all <italic>K</italic>+1 characters including gaps occur in the ancestor at the stationary frequencies of the Markov process). Because our model is non-reversible, we must always work with rooted phylogenies.</p> <p>Because we map generative evolutionary model events onto an alignment column, the column is assumed to be correctly aligned phylogenetically – all aligned residues are assumed to be homologous and related only by substitution events. This means that no more than one insertion event may occur in any given column. Enforcing this assumption requires modification of the usual Felsenstein peeling algorithm to include some extra bookkeeping.</p> <p>Given alignments are unlikely to be phylogenetically correct, because humans and alignment programs tend to produce aesthetically pleasing alignments that compress columns containing few residues. Importantly, for any arbitrary tree topology and any arrangement of observed residues and gap characters in an alignment column, there exists at least one possible assignment of characters to ancestral nodes that makes all extant aligned residues homologous. Therefore the problem with using phylogenetically incorrect alignments is not that the algorithm will fail altogether (as would happen if some combinations of alignments and trees were impossible), but rather that we can expect its inference ability to be degraded by forced inference of incorrect histories in phylogenetically incorrectly aligned columns. How much the overall inference is degraded by this and by the other assumptions described above is a matter for empirical testing, which we describe in the second half of the paper, after we describe the model itself.</p> <title>Solving the Markov Process for a Gap-Extended Rate Matrix Model</title> <p>First, we start by solving the Markov process associated with a rate matrix extended to include a gap character. For an alphabet of <italic>K</italic> residues, probabilistic substitution models are defined by a <italic>K</italic>×<italic>K</italic> rate matrix <italic>R</italic> such that the matrix of conditional probabilities <italic>Q<sub>t</sub></italic>(<italic>i</italic>,<italic>j</italic>) ≡ <italic>P</italic>(<italic>j</italic>|<italic>i</italic>,<italic>t</italic>) is given by
</p> <p>We extend this to include the gap character by augmenting the rate matrix to a (<italic>K</italic>+1)×(<italic>K</italic>+1) matrix <italic>R<sup>ε</sup></italic> that depends on arbitrary rates of deletions <italic>μ</italic>≥0 and insertions <italic>λ</italic>≥0:where <italic>p</italic> = (<italic>p</italic>
<sub>1</sub>,…,<italic>p<sub>K</sub></italic>) is the distribution of inserted residues , and <italic>δ<sub>ij</sub></italic> stands for the Kronecker delta in the <italic>K</italic>×<italic>K</italic> subspace (valued one if <italic>i</italic> = <italic>j</italic> and zero otherwise). Since the extended rate matrix <italic>R<sup>ε</sup></italic> has the property that each row adds up to zero, we can construct a model of evolution for the extended rate matrix defined asAt zero divergence, the probabilities of any insertion or deletion of a residue or any substitution of a residue to a different residue are all zero.</p> <p>Generally, the extended conditionals  can be cast into the form,with the conditions  and , for each row <italic>i</italic>, and where <italic>M<sub>t</sub></italic> is the <italic>K</italic>×<italic>K</italic> conditional substitution matrix (to be defined later).</p> <p>In particular, for a reversible <italic>K</italic>×<italic>K</italic> rate matrix <italic>R</italic>, if we assume that the distribution of inserted residues in equation (2) is the stationary distribution associated to the reversible rate <italic>R</italic> (<italic>π<sub>i</sub>R</italic>(<italic>i</italic>,<italic>j</italic>) = <italic>π<sub>j</sub>R</italic>(<italic>j</italic>,<italic>i</italic>)) then one can derive the analytic expression for  in terms of the solution for the <italic>K</italic>×<italic>K</italic> conditional matrix <italic>Q<sub>t</sub></italic> = <italic>e<sup>tR</sup></italic> and the rates of insertion and deletion. This particular solution for the extended conditional probabilities is given bywhere the gap-specific functions <italic>γ<sub>t</sub></italic> and <italic>ξ<sub>t</sub></italic> are given by
when at least one of the two rates is positive. For the particular case of no insertions and deletions (<italic>λ</italic> = <italic>μ</italic> = 0) both functions are defined as identically zero.</p> <p>Finally, we have to describe how to obtain the <italic>K</italic>×<italic>K</italic> conditional substitution matrix subspace <italic>M<sub>t</sub></italic> in  from a rate matrix <italic>R</italic> and the insertion and deletion rates. Generally, any biologically relevant <italic>K</italic>×<italic>K</italic> substitution rate matrix <italic>R</italic> must have zero as a non-degenerate eigenvalue, and at least one other negative eigenvalue. To express <italic>M<sub>t</sub></italic> in general form, let (0,−<italic>e</italic>
<sub>1</sub>,…,−<italic>e<sub>A</sub></italic>) represent the eigenvalues of <italic>R</italic>, for 1≤<italic>A</italic>≤<italic>K</italic>−1, with <italic>e<sub>a</sub></italic>&gt;0 for 1≤<italic>a</italic>≤<italic>A</italic>. Then the standard conditional substitution probabilities could be expressed as:where <italic>O<sub>a</sub></italic> is a <italic>K</italic>×<italic>K</italic> matrix of real numbers specific for each unique nonzero eigenvalue. ##SUPPL##1##Text S1## shows how to derive the <italic>O<sub>a</sub></italic> matrices from the similarity transformation that relates R with its diagonal form.</p> <p>Now the extended conditional matrix in the substitution-only subspace <italic>M<sub>t</sub></italic> is expressed as:If <italic>λ</italic> = <italic>μ</italic> = 0, <italic>M<sub>t</sub></italic>(<italic>i</italic>,<italic>j</italic>) is defined as the standard substitution conditional matrix <italic>Q<sub>t</sub></italic>(<italic>i</italic>,<italic>j</italic>) in Equation 8.</p> <p>For example, the F84 model ##REF##8583911##[14]## is defined by the substitution rate <italic>R</italic>(<italic>i</italic>≠<italic>j</italic>) = <italic>βπ<sub>j</sub></italic>+<italic>α</italic>Δ<italic><sub>ij</sub></italic>, where , and where the function <italic>ε<sub>ij</sub></italic> is 1 if <italic>i</italic>,<italic>j</italic> are both either purines or pyrimidines, and 0 otherwise. The F84 rate depends on two non-negative parameters <italic>α</italic> and <italic>β</italic>, and on a stationary residue distribution <italic>π</italic>. The F84 rate matrix has non-zero eigenvalues <italic>e</italic>
<sub>1</sub> = <italic>β</italic>, and <italic>e</italic>
<sub>2</sub> = <italic>α</italic>+<italic>β</italic>. The matrix <italic>O</italic>
<sub>1</sub> is given by <italic>O</italic>
<sub>1</sub>(<italic>i</italic>,<italic>j</italic>) = Δ<italic><sub>ij</sub></italic>−<italic>π<sub>j</sub></italic>, and the matrix <italic>O</italic>
<sub>2</sub> is given by <italic>O</italic>
<sub>2</sub>(<italic>i</italic>,<italic>j</italic>) = <italic>δ<sub>ij</sub></italic>−Δ<italic><sub>ij</sub></italic>. For the gap-extended F84 model, we obtainwhich for the particular case <italic>λ</italic> = <italic>μ</italic> = 0 is defined as the original F84 model for substitutions,
</p> <p>The conditional model of McGuire <italic>et al.</italic>
##REF##11264399##[41]## describes a particular solution of the extended F84 model presented here, in which the rates of insertions and deletions are constrained to satisfy the conditions <italic>λ</italic> = <italic>β</italic> and <italic>μ</italic> ∝ <italic>β</italic>.</p> <p>A more detailed description of the characteristic differential equations for an arbitrary gap-augmented rate matrix and the particular solution described above is provided in ##SUPPL##1##Text S1##. Under different assumptions, such as a reversible substitution rate matrix but using a distribution of inserted residues other than the stationary distribution for <italic>R</italic>, or a non-reversible substitution rate matrix,  could still be obtained numerically ##UREF##21##[72]##.</p> <p>A note about reversibility. For a reversible substitution rate matrix, the particular Markov process including gaps solved here is also reversible, as can be seen by using the marginal frequencies <italic>π<sub>i</sub>λ</italic>/(<italic>λ</italic>+<italic>μ</italic>) for a residue, and <italic>μ</italic>/(<italic>λ</italic>+<italic>μ</italic>) for the gap character. However, one has to distinguish between the reversibility of a Markov process and the reversibility of an evolutionary process constructed using that Markov chain. A Markov chain is said to be reversible if and only if there exists a marginal distribution that satisfies the reversibility condition ##UREF##22##[73]##. However, regardless of whether the Markov process is reversible in the strict sense, when constructing an evolutionary process from it one may specify an ancestral marginal distribution other than what reversibility requires. In evolutionary models, “reversibility” is generally taken to mean that the joint probability of an ancestral and a descendant sequence is invariant regardless of which sequence is used as the ancestor and which is used as the descendant. In probabilistic inference, this latter definition of reversibility is the most relevant (for instance to invoke the pulley principle ##REF##7288891##[1]##). From here on in this paper we use the term reversibility in this latter (broader) sense. Thus, the evolutionary model that we construct in the next section is not reversible in this (usual) sense though the Markov chain that it is based on is reversible in the strict sense.</p> <p>Up to this point, this is essentially McGuire's model ##REF##11264399##[41]## (with a minor generalization). That model has the conceptual problems we described in the preamble, thought to be inherent to approaches based on gap-extended rate matrices. In the next section, we show how these problems may be circumvented.</p> <title>The Generative Model</title> <p>We construct the generative model as an independent product of single-event (substitution, deletion or insertion) contributions. What we will do to construct the model is to borrow terms from the conditional Markov process introduced in the previous section. We describe the probability of an insertion as proportional to <italic>ξ<sub>t</sub>π<sub>i</sub></italic> in Equation 7, the probability of a deletion as proportional to <italic>γ<sub>t</sub></italic> in Equation 6, and that of a substitution as proportional to <italic>M<sub>t</sub></italic>(<italic>i</italic>,<italic>j</italic>) in Equation 9, but we ignore the gap to gap transition of the Markov process. This allows us to derive a generative model that describes insertions and deletions not as mere “gap character” replacements on a fixed length alignment but as true evolutionary events (births and deaths of residues).</p> <p>A given sequence <bold>x</bold> = {<italic>x</italic>
<sub>1</sub>…<italic>x<sub>l</sub></italic>} that evolves to another sequence <bold>y</bold> = {<italic>y</italic>
<sub>1</sub>…<italic>y<sub>l</sub></italic>
<sub>′</sub>} corresponds unambiguously to a pairwise alignment in which a substitution is represented by a conserved or mismatched column, a deletion by an ancestral residue aligned to a gap in the descendant sequence, and an insertion is represented by an ancestral gap aligned to a residue. Let <bold>x</bold>ˆ = <italic>xˆ</italic>
<sub>1</sub>‥<italic>xˆ</italic>
<sub>L</sub> and <bold>y</bold>ˆ = <italic>yˆ</italic>
<sub>1</sub>‥<italic>yˆ</italic>
<sub>L</sub> mean the aligned sequences <bold>x</bold> and <bold>y</bold>. Specifically, the probability that <bold>y</bold> was generated from <bold>x</bold> after time <italic>t</italic> with pairwise alignment <bold>x</bold>ˆ<bold>y</bold>ˆ that includes <italic>s</italic>≤<italic>l</italic> substitutions, (<italic>l</italic>−<italic>s</italic>) deletions, and (<italic>l</italic>′−<italic>s</italic>) insertions {I<sub>1</sub>…I<italic><sub>l</sub></italic>
<sub>′−<italic>s</italic></sub>} (where the subset of residues  from <bold>x</bold> are substituted by  in sequence <bold>y</bold>) is constructed asThe functions <italic>γ<sub>t</sub></italic>, <italic>ξ<sub>t</sub></italic> and  are given by the Markov model solutions (Equations 6, 7, an 9). The residue distribution <italic>π</italic>, is set to the stationary distribution of the substitutions rate matrix. In that way, for a reversible substitution rate matrix, this generative model is quasi-reversible (<italic>i.e.</italic> reversible in the substitutions subspace).</p> <p>In equation (12), everything but the term (1−<italic>ξ<sub>t</sub></italic>)<italic><sup>l</sup></italic>
<sup>+1</sup> is borrowed from the Markov process. This extra term is responsible for having a normalized distribution, such that the sum of the contributions of all sequences <bold>y</bold> of all possible lengths and all possible alignments is one. The extra term accounts for the fact given a sequence <bold>x</bold> of length <italic>l</italic>, insertions can occur at (<italic>l</italic>+1) places. That is, at each of the (<italic>l</italic>+1) places that an insertion in <bold>x</bold> could occur, an insertion of length <italic>z</italic> occurs with probability  (a normalized geometric distribution).</p> <p>Using the generative probability distribution (Equation 12), one can calculate the expected length (in residues without gaps) of descendant sequences originated from an ancestral sequence of length <italic>l</italic> after summing to all possible patterns of substitutions, insertions and deletions, that is given by (full derivation in ##SUPPL##2##Text S2##):
</p> <p>Crucially, this model does not assume that all observed sequences are generated from a preset number of aligned columns. Thus there is no “memory effect” ##UREF##20##[70]##. Our model is a generative probabilistic model with a close correspondence to the birth-death description of the <sc>tkf</sc>91 model ##REF##1920447##[50]##, (see ##SUPPL##2##Text S2##, for more detail).</p> <p>We have implemented this generative model in a computer program for evolving sequences named <italic>ε</italic>
<sc>rate</sc>.</p> <title>The Joint Probability of a Pairwise Ancestor/Descendant Alignment</title> <p>Given the <italic>conditional</italic> probabilities of the generative model, one wants to calculate the <italic>joint</italic> probability of a multiple alignment of sequences generated with the model, given a phylogenetic tree. Before doing that, let us consider the joint probability of a pairwise ancestor/descendant alignment, as a building block for a full-fledged algorithm for phylogenetic inference on multiple alignments.</p> <p>Using the same notation as in the previous section, we can calculate the joint probability of the pairwise ancestor/descendant alignment (<bold>x</bold>ˆ,<bold>y</bold>ˆ) after divergence time <italic>t</italic> using the expressionwhere <italic>P<sup>ε</sup></italic>(<bold>x</bold>) is the probability of the ancestral sequence. This distribution is a prior, not determined by the generative model. In order to have this prior factorize effectively into columns, we will assume the length of ancestral sequences follows a geometric distribution (1−<italic>p</italic>)<italic>p<sup>l</sup></italic> for a sequence of length <italic>l</italic> with arbitrary Bernoulli frequency parameter 0&lt;<italic>p</italic>&lt;1.</p> <p>The joint probability of a pairwise ancestor/descendant alignment can then be factorized as a product of terms over alignment columns as:where the length of the alignment <sc>l</sc> = <italic>s</italic>+<italic>d</italic>+<italic>i</italic>. Normalization of this joint probability requires including the extra term <italic>P<sup>ε</sup></italic>(★|<italic>t</italic>) = (1−<italic>p</italic>)(1−<italic>ξ<sub>t</sub></italic>) that cannot be associated to any observed column. Intuitively this term can be viewed as an extra (terminating) column (★). An analogous “extra column” term will appear in the multiple alignment case with the Felsenstein pruning algorithm.</p> <p>In ##SUPPL##3##Text S3##, we use the joint probabilities in Equation 15 to calculate other related length distributions of the model such as the length distribution for descendant sequences and the length distribution of alignment length, by summing over the other variables (marginalization). Because the model is non-reversible, the distribution of descendant sequences is different from that of ancestral sequences, and depends on the divergence time (non-stationary).</p> <p>From the joint alignment probability distribution, we can also calculate the expected frequencies of insertions  and deletions  in pairwise ancestor/descendant alignments (see ##SUPPL##3##Text S3## for derivation). These quantities illustrate some important properties of the model. They will also be useful in making comparisons to the properties of other reversible models (see <xref ref-type="sec" rid="s3">Discussion</xref>). They are,

</p> <p>Because the model is non-reversible,  and  are different (note that we avoid the traditional term “indel” in this paper, because we are not treating them reversibly). When the insertion rate is zero (and the deletion rate is positive), the expected frequency of insertions is zero for all divergence times but the expected frequency of deletions is not, as one would expect. Similarly, when the deletion rate is zero (and the insertion rate is positive), the expected frequency of deletions is zero for all divergence times while the expected frequency of insertions is not. For more detail, see ##SUPPL##3##Text S3##.</p> <title>Felsenstein's Peeling Algorithm Extended to Gaps</title> <p>Felsenstein's peeling algorithm is an efficient algorithm for calculating the probability of a multiple alignment given a tree and a Markov substitution model. Extending the Felsenstein algorithm to include insertion and deletion events with the model we propose here requires four modifications. Those are: extra bookkeeping in the Felsenstein recursions to enforce that no more than one insertion occurs per column, so that all aligned residues are homologous; including a term from the prior ancestral sequence length distribution in the calculation of each individual column likelihood; including in the overall alignment likelihood the extra normalization terms collected in the “extra column” (★); and finally, marginalizing the contributions of possible ancestral residues that have left no trace in extant sequences. Otherwise, the substitution model assumed by the Felsenstein peeling algorithm is simply replaced by the generative model described in the previous section by Equation 12 which includes substitutions, insertions and deletions.</p> <p>Using the notation of ##UREF##2##[6]##, for a given position <italic>u</italic> in the alignment, let <italic>P<sub>u</sub></italic>(<italic>L<sub>k</sub></italic>,<italic>i</italic>) be the probability up to node <italic>k</italic> given that the character (residue or gap) at node <italic>k</italic> is <italic>i</italic>. For a residue <italic>i</italic>, <italic>P<sub>u</sub></italic>(<italic>L<sub>k</sub></italic>,<italic>i</italic>) cannot contain any insertion only substitutions and deletions. For a gap, <italic>P<sub>u</sub></italic>(<italic>L<sub>k</sub></italic>,–) has to include one and only one insertion. Thus, <italic>P<sub>u</sub></italic>(<italic>L<sub>k</sub></italic>,–) is defined to be zero if all the nodes under <italic>k</italic> are gaps for that position. These probabilities are calculated recursively starting from the leaves of the binary tree as,</p> <p>If node <italic>k</italic> is a leaf, for a residue <italic>i</italic>,for a gap,
</p> <p>If node <italic>k</italic> is not a leaf, for a residue <italic>i</italic>,for a gap,where  and  are the two daughters of node <italic>k</italic>, , and  are the distances from node <italic>k</italic> to its left and right child respectively, and where the probabilities for the daughter nodes have already been calculated by the recursion. <italic>u<sub>k</sub></italic> stands for the subset of leaves under node <italic>k</italic> for column <italic>u</italic>, and <italic>u<sub>k</sub></italic> = – indicates that all leaves under node <italic>k</italic> are gaps for column <italic>u</italic>. The single-event conditional probabilities are dictated by the generative model in Equation 12 as,

for 1≤<italic>i</italic>,<italic>j</italic>≤<italic>K</italic>, where the functions <italic>γ<sub>t</sub></italic>, <italic>ξ<sub>t</sub></italic> and  are given by the Markov model solutions (Equations 6, 7, and 9). ##FIG##0##Figure 1## shows a graphical interpretation of the Felsenstein recursions described in Equations 20 and 21.</p> <p>The second modification is due to the existence of a length distribution for ancestral sequences for this model. In order to calculate the total probability of a given aligned column <italic>P</italic>(<italic>u</italic>|<italic>T</italic>,<italic>R<sup>ε</sup></italic>,<italic>p</italic>), we need to use the probability of the sequence at the root. The total probability of site <italic>u</italic> is given bywhere <italic>p</italic> is the parameter of the geometric distribution of ancestral sequences, and <italic>π</italic> are the residue frequencies at the root.</p> <p>Third, the factorization in columns of the unconditional in length alignment distribution leaves some normalization terms that we gather together into what we think of as an “extra column” (★) contribution. Thus, when calculating the total probability of a multiple alignment as the product of <italic>l</italic> individual columns, there is an additional term in the equation:The term <italic>P</italic>(★|<italic>T</italic>,<italic>R<sup>ε</sup></italic>,<italic>p</italic>) = (1−<italic>p</italic>)<italic>P</italic>
<sub>★</sub>(<italic>L<sub>root</sub></italic>) is calculated by a similar peeling algorithm,where .</p> <p>Lastly, one could delete an ancestral residue leaving no trace in extant sequences—a column could exist but not be present in the observed alignment—therefore we are interested in calculating <italic>P</italic>(<italic>L</italic>
<sub>0</sub>) the probability of an alignment with <italic>L</italic>
<sub>0</sub> observed columns after marginalizing the unobserved columns. The probability of an alignment with <italic>L</italic>
<sub>0</sub> observed columns and <italic>k</italic> unobserved columns is given bywhere the probability of an all-gaps column <italic>u</italic>
<sub>gap</sub> is a particular case of Equation. 25. Then, for a given alignment of <italic>L</italic>
<sub>0</sub> observed columns, we marginalize to all possible unobserved columns as,This is the final expression for the probability of a given alignment of <italic>L</italic>
<sub>0</sub> columns, after marginalizing all possible unobserved ancestral residues.</p> <p>This algorithm will now reproduce the results of the original Felsenstein algorithm for ungapped alignments when the parameters <italic>λ</italic> = <italic>μ</italic> = 0 (except for a geometric term which is the same for all trees, thus will not affect the maximum likelihood estimation of a tree). In addition, <italic>λ</italic> = <italic>μ</italic> = 0 are the optimal parameter choices for any ungapped multiple alignment for any given tree. Notice that for the previous statement to be true, it is crucial to have the extra term (1−<italic>ξ<sub>t</sub></italic>) in the generative process of a substitution. Once we set <italic>μ</italic> = 0, the residue-residue substitution process in the presence of gaps (described by Equation 9) reaches the same asymptotic value (described by Equation 8) for any value of <italic>λ</italic>, including both the limit <italic>λ</italic> = 0 and the limit <italic>λ</italic> = ∞. It is the extra term that renders the total probability of any finite-length alignment to zero in the case <italic>λ</italic> = ∞, and makes <italic>λ</italic> = 0 optimal on ungapped alignments.</p> <p>The extended peeling algorithm has worse-case time complexity  for an alphabet of size <italic>K</italic> and a multiple alignment of <sc>l</sc> columns and <italic>n</italic> sequences.</p> <title>An Implementation in phylip: dnaml<italic>ε</italic>
</title> <p>To test an application of our model, we modified the program <sc>dnaml</sc> (from the <sc>phylip</sc> package version 3.66, 4 August 2006 ##UREF##7##[24]##) to use our generative model. Given a nucleotide multiple alignment, the <sc>dnaml</sc> program infers a phylogenetic tree for the sequences by maximum likelihood under an F84 rate matrix model ##REF##8583911##[14]##. Our modified program <sc>dnaml</sc>
<italic>ε</italic> uses the extended F84 rate matrix given by Equations 5, 6, 7, and 10, and implements the extended Felsenstein algorithm described in the previous section.</p> <p>
<sc>dnaml</sc>
<italic>ε</italic> uses the same core algorithms for maximum likelihood tree inference as the original <sc>dnaml</sc>. For every <sc>dnaml</sc> function we implemented a <sc>dnaml</sc>
<italic>ε</italic> counterpart. In addition, <sc>dnaml</sc>
<italic>ε</italic> optimizes the gap parameters (<italic>λ</italic>, <italic>μ</italic>). After each new branch is added to the tree using <sc>dnaml</sc>'s Newton-Raphson method, <sc>dnaml</sc>
<italic>ε</italic> midpoint roots the tree, and then jointly optimizes both branch lengths and gap parameters using conjugate gradient descent. Midpoint rooting is an easy but simplistic rooting method; we made no attempt in <sc>dnaml</sc>
<italic>ε</italic> to optimize the placement of the root, though this would be possible without a significant efficiency cost, using methods described in ##REF##17060197##[74]##, for example.</p> <p>For simplicity, in this paper's results, we approximate the average length of ancestral sequences () by the average length of sequences in the given alignment. This is in line with common practice in maximum likelihood methods for phylogenetic inference for substitutions, which approximate the prior residue distribution at the root by the observed residue frequencies in the data ##UREF##7##[24]##. In Bayesian methods, the root residue distribution and other parameters of the rate matrix are determined as part of the inference process usually in combination with Markov-chain Monte Carlo (MCMC) methods as for instance in ##UREF##1##[5]##; parameter <italic>p</italic> could be estimated using similar techniques.</p> <p>Given an ungapped alignment, <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> are expected to produce identical results (same unrooted tree with same likelihood, same branch lengths, and confidence limits). Given a gapped alignment, comparison of the two implementations allows us to ask how much performance improves when gaps are treated by our extended model as opposed to treating gaps as missing data.</p> <title>Benchmarking</title> <p>We compared <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> in three types of benchmarking experiments. First, using simulated ungapped alignment data, we confirmed that the two programs give essentially identical results when no insertions and deletions are present. Second, we assessed the ability of the two programs to accurately infer phylogenetic tree topologies for simulated alignment data with insertions and deletions from arbitrarily generated trees. Third, we assessed the programs' ability to correctly infer phylogenetic trees on real ribosomal RNA data.</p> <p>The problem with benchmarking phylogenetic inference methods is knowing the ground truth. Evolutionary trees are only experimentally known in a few unusual cases for short time scales and rapidly evolving organisms; all other trees have been inferred. The strength of simulated data experiments is that the tree is known, but their weakness is that the simulation must assume an evolutionary model that may be biologically unrealistic and/or too similar to the evolutionary model assumed by the inference method to be tested. On the other hand, the strength of testing on real biological alignments is that the data are realistic, and the weakness is that the true tree is unknown. To evaluate inferences on real data, we have developed a “concordance test” that does not rely on the true tree being known. Testing on both simulated and real data should help compensate for weaknesses of either approach.</p> <title>Tests on simulated ungapped alignments</title> <p>We first made sure that <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> produce essentially identical inferences when no insertions or deletions are present. This is a control experiment, making sure that <sc>dnaml</sc>
<italic>ε</italic> correctly infers that optimal deletion and insertion rates are 0 on ungapped alignments, and that the likelihood calculation correctly reduces to the original ungapped version.</p> <p>We generated simulated 8-taxon rooted trees using the algorithm of Kuhn and Felsenstein ##REF##8015439##[75]##, which samples a variety of branch lengths and topologies. Each tree was rescaled to a chosen average branch length.</p> <p>We used the program <sc>rose</sc> to generate simulated alignments from these sampled phylogenetic trees ##REF##9545448##[76]##. We use <sc>rose</sc> because it allows us to implement reasonably realistic models of insertion and deletion, which we will describe and use in the next section; in these initial ungapped control experiments, we turn off the insertion and deletion parameters. We modified <sc>rose</sc> to use the F84 rate matrix model for residue substitution, the same model used in <sc>dnaml</sc> and <sc>dnaml</sc>
<italic>ε</italic>. We used a uniform stationary distribution of 25% for each residue, and a transition/transversion ratio of 2.0 (the <sc>dnaml</sc> default).</p> <p>We sampled trees for 9 different average branch lengths, ranging from 0.005 to 2.0 substitutions/site and branch in roughly 2–3× multiplicative steps (corresponding to average pairwise identities ranging from 98% down to a fully saturated 25% in the 8 sequences). (The units of evolutionary time are in principle arbitrary. Substitution rate matrices are traditionally normalized to units of substitutions/site. <sc>dnaml</sc>
<italic>ε</italic> reports time in units of <italic>changes</italic> per site, where changes include substitutions, insertions, and deletions.) For each choice of average branch length, we sample 100 different trees. For each tree, we generate different alignment lengths ranging from 50 to 1000 in steps of 5. A total of 900 trees and 171900 alignments were generated (100 samples each for 9 choices of branch length and 191 choices of alignment length). For each alignment, we infer a maximum likelihood tree topology using <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc>.</p> <p>To evaluate the correctness of inferred tree topologies, we computed three standard measures: the fraction of correct topologies (true positives, TP), the symmetric difference distance (SDD) ##UREF##23##[77]##, and the branch score distance (BSD) ##REF##8015439##[75]##. TP simply counts an inferred topology as right or wrong (high TP is better). SDD counts the number of non-identical internal branches in the two (unrooted) trees, where “identical” means splitting the two trees into the same disjoint sets of taxa; for a comparison of trees of 8 taxa, SDD ranges from 0 for identical trees to 10 for maximally dissimilar topologies. BSD evaluates not just the topology but also the correctness of inferred branch lengths, by summing the cost assigned to each internal branch of the square of the difference to the identical branch in the other tree (if there is no identical branch in the other tree, the cost is the whole branch square). Because BSD depends on total tree branch length, in order to compare results across different branch lengths, we calculate a normalized BSD (nBSD) on trees rescaled to an average branch length of one. Both SDD and BSD were calculated using the <sc>phylip</sc> program <sc>treedist</sc>. Better inferences are indicated by larger TP and by smaller SDD and nBSD.</p> <p>The results are shown in ##FIG##1##Figure 2##. As expected, the accuracy of tree topologies inferred by the two methods is essentially identical by the TP and SDD measures. As in any phylogenetic inference method, accuracy improves with alignment length (more data is better), and shows an optimum average branch length of about 0.05, corresponding to average pairwise sequence identities of about 82% (more similar sequences have fewer substitutions and less signal, and less similar sequences are more saturated). In the best cases (0.05 branch length, alignments longer than 800 nt) both methods infer the correct tree about 78% of the time (##FIG##2##Figure 3A##). nBSD values are also essentially identical for most choices of average branch length. We do see significant differences in branch length estimation (either by the nBSD test, or by plotting average inferred branch length) at large, saturating choices of average branch lengths of 1.0 or 2.0, corresponding to almost uncorrelated random sequences. In this extreme (and not biologically relevant) regime, branch lengths make little difference in likelihood so long as they are large (and accordingly, likelihoods assigned by the two methods are not significantly different, despite the different inferred branch lengths), and the inferred branch lengths become sensitive to details of the optimization method. Thus these extreme cases appear to be identifying a biologically irrelevant weakness in <sc>dnaml</sc>'s optimizer, which infers overly large branch lengths (5.7±0.5 instead of 2.0 for the largest alignments), whereas <sc>dnaml</sc>
<italic>ε</italic> results, produced by a different implementation of optimizer, are closer to the correct value (1.55±0.05). Overall, with the exception of this minor difference in numerical optimization, these results indicate that <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> do produce essentially identical results on ungapped alignments, as expected.</p> <title>Tests on simulated alignments with insertions and deletions</title> <p>Using the same protocol as above, we then asked whether <sc>dnaml</sc>
<italic>ε</italic> would produce better inferences on gapped alignments than <sc>dnaml</sc>. <sc>rose</sc> simulates insertions and deletions by inserting and deleting a block of residues at a random position according to a chosen length distribution, and controlled by insertion/deletion probability parameters <italic>p<sub>ins</sub></italic> and <italic>p<sub>del</sub></italic> per unit branch length. Note that because <sc>rose</sc> generates multi-residue insertions and deletions, it simulates a different and more biologically realistic process than the column-independent insertion/deletion process our model assumes. We therefore expect our model to suffer in the simulation from its assumption that gaps are uncorrelated, much as it would in real data; the question is whether this is outweighed by having at least some model of insertion and deletion, compared to <sc>dnaml</sc>'s treatment of gaps as missing data.</p> <title>Tests using a Poisson insertion and deletion length distribution</title> <p>We had <sc>rose</sc> sample insertion and deletion lengths <italic>l</italic> from a Poisson distribution  with <italic>λ</italic> = 0.5. This particular distribution allows with probability 0.3935 insertions or deletions of more than one nucleotide at the time. This is still a weak violation of the column independence assumptions, and in the next section we will test the limits of our model regarding this issue. We used a range of 14 different gap probability parameters <italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic>, ranging from <italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic> = 0 (no gaps) to <italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic> = 0.1.</p> <p>Results of an illustrative comparison of two choices of gap probability and average branch length are shown in ##FIG##2##Figure 3##. These two datasets have roughly equal frequencies of observed substitution and therefore show roughly equal performance with <sc>dnaml</sc> (left side of ##FIG##2##Figure 3##): in the first (red points) <italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic> = 0 and average substitution branch length of 0.06 produced ungapped alignment datasets in which all changes were substitutions (pairwise comparisons averaged 21% substitutions and 79% overall identity), and in the second (blue points), setting <italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic> = 0.001 and an average substitution branch length of 0.07 produced alignment datasets with substitution events that averaged 20% in pairwise comparisons, and insertion/deletion events that averaged 13% in pairwise comparisons. These average percentages of substitutions and gaps in the blue dataset are reasonably realistic, comparable to those observed in our real rRNA alignments; see ##TAB##3##Table 4##. The right side of ##FIG##2##Figure 3## then shows the effect of <sc>dnaml</sc>
<italic>ε</italic> being able to use the information in gaps in addition to substitutions: for <sc>dnaml</sc>
<italic>ε</italic> inferences, the blue dataset (the one with gaps) is more informative, and accuracy increases.</p> <p>
##TAB##0##Table 1## shows the relative performance of <sc>dnaml</sc>
<italic>ε</italic> versus <sc>dnaml</sc> for a large range of alignments with gaps. In order to summarize with one number curves similar to those presented in ##FIG##2##Figure 3##, we calculate the area under the curve (<sc>auc</sc>) in analogy to ROC curves. In particular, we introduce the relative <sc>auc</sc> difference, . For a given measure <italic>f</italic>, the quantity  describes (as a percentage) the difference between the <sc>auc</sc> for <sc>dnaml</sc>
<italic>ε</italic> versus that of <sc>dnaml</sc> relative to the larger of the two <sc>auc</sc>'s. <sc>dnaml</sc>
<italic>ε</italic> shows consistently better performance than <sc>dnaml</sc>, and the improvement increases with the percentage of gaps in the alignment.</p> <p>
##FIG##2##Figure 3## and ##TAB##0##Table 1## also show the mean branch lengths inferred by the two methods. For the gapped alignments (blue), <sc>dnaml</sc>
<italic>ε</italic> infers longer branch lengths than <sc>dnaml</sc> (in units of changes/site) because it counts insertions and deletions as changes, in addition to substitutions.</p> <p>
##TAB##0##Table 1## also shows results (in gray) when instead of using the phylogenetically correct alignments, we allow the sequences to be realigned using a standard (and not phylogenetically aware) alignment algorithm. In particular for the experiments in this paper, we used <sc>clustalw</sc> version 1.83. The more noticeable effect of realigning is that the columns of the alignment get compressed. The number of substitutions increases and the number of gaps decreases respect to the phylogenetically correct alignments. As expected, performance of <sc>dnaml</sc>
<italic>ε</italic> declines in terms of predicting the phylogenetic topology correctly, (although the two methods become almost equivalent in terms of topological measures as the number of gaps in the realignments increases). However, as seen in ##TAB##0##Table 1##, <sc>dnaml</sc>
<italic>ε</italic> still shows improvement in all cases in terms of the nBSD measure which takes into account the combined effect of both the accuracy of the topology and the accuracy of the branch lengths. These results show that the loss due to the lack of a phylogenetically correct alignment seems to be compensated by the contribution to branch lengths introduced by treating insertions and deletions as real evolutionary events.</p> <p>The time performance of <sc>dnaml</sc>
<italic>ε</italic> is similar to that of <sc>dnaml</sc>. Both scale linearly with alignment length (and with number of taxa, though this is not shown by these experiments, which use a fixed number of 8 taxa). For both methods, the scaling factor increases as the number of substitutions and the number of gaps increases (see ##TAB##0##Table 1##). <sc>dnaml</sc>
<italic>ε</italic> shows (for trees of eight taxa and alignments of 1,000 residues) about a three-fold increase in average run time compared to <sc>dnaml</sc>, whether the alignment has gaps or not.</p> <title>Tests sampling from the generative model</title> <p>As a matter of completeness, we have also sampled from the generative model described in this work using the program <italic>ε</italic>
<sc>rate</sc>. The model inserts residues following a geometric distribution controlled by the time-dependent parameter <italic>ξ<sub>t</sub></italic> in Equation 7. We applied the protocol described in the previous section. Results for gapped alignments with a large range of degrees divergence (comparable to those presented in the previous section) are given in ##TAB##1##Table 2##. As expected, this is a relatively easier test than that posed by alignments with Poisson distributed gaps, but in general, similar trends are observed in the results (compare ##TAB##0##Table 1## with ##TAB##1##Table 2##, see also ##TAB##2##Table 3##).</p> <title>Test of the column independence assumption</title> <p>In order to test the limitations of considering individual gaps as independent events, we also had <sc>rose</sc> sample insertions and deletions from a distribution that simulates a DNA alignment of a protein-coding region. To that purpose, we used a nucleotide version of the empirically-derived coding gap distribution introduced in <sc>simprot</sc>
##REF##16188037##[78]## that allows gaps to occur only as multiple of three. We also added an option in <sc>rose</sc> in order to allow insertions and deletions to occur randomly but only at multiple of three positions with respect to the start of the alignment.</p> <p>The <sc>simprot</sc> gap length distribution depends on one parameter <italic>c</italic>/<italic>t</italic>. We selected three different values, in increasing order of divergence: <italic>c</italic>/<italic>t</italic> = 100 for which <italic>p</italic>(3 nts) = 0.7238 and <italic>p</italic>(24 nts) = 0.0001; a second distribution <italic>c</italic>/<italic>t</italic> = 6 which is a good approximation to a empirically-determined distribution trained on protein sequences with less than 100 PAM sequence divergence ##REF##15276848##[79]##; and finally <italic>c</italic>/<italic>t</italic> = 1 which corresponds exactly to a distribution empirically determined from proteins sharing no more than 25% sequence identity ##REF##11536366##[80]##. We limited the maximum number of inserted or deleted nucleotides to 100. For the <italic>c</italic>/<italic>t</italic> = 6 distribution, one has <italic>p</italic>(3 nts) = 0.4510 and <italic>p</italic>(99 nts) = 0.0009. For the <italic>c</italic>/<italic>t</italic> = 1 distribution, one has <italic>p</italic>(3 nts) = 0.2938 and <italic>p</italic>(99 nts) = 0.0034. For each of the three protein-coding gap distributions, we applied the same protocol as for the Poisson distribution used previously.</p> <p>As we observe in ##TAB##2##Table 3##, the <italic>c</italic>/<italic>t</italic> = 100 coding gap length distribution <sc>dnaml</sc>
<italic>ε</italic> still shows an improvement respect to <sc>dnaml</sc> at least for the two topological measures, this is despite the fact that this distribution breaks the column independence assumption more strongly than the Poisson distribution used in the previous section (<italic>p</italic>(<italic>l</italic>≥6 nts) is 0.2762 for this distribution, compared to 0.0002 for the Poisson distribution). However, for the two more divergent coding gap length distributions, the size of the insertion/deletion blocks is so large that by taking gaps into account, we cannot even reproduce the topology of the original substitutions-only tree.</p> <title>Test for long-branch attraction</title> <p>To make sure that <sc>dnaml</sc>
<italic>ε</italic> does not suffer for some unexpected reason from significantly more systematic long-branch attraction than <sc>dnaml</sc> itself, we used 4-taxon trees depending on two parameters ##UREF##24##[81]##. One parameter is the branch length of the internal branch and two opposite leaves (the three-branch parameter, <italic>t</italic>
<sub>3</sub>), the other parameter is the branch length of the other two opposite leaves (the two-branch parameter, <italic>t</italic>
<sub>2</sub>). Using the same protocol as before, we assessed the programs performance on alignments generated according to 4-taxon trees in which the two-branch parameter is 10 times longer than the three-branch parameter, for all gap length distributions used above.</p> <p>Results are presented in ##TAB##3##Table 4##. For the phylogenetically correct alignments, we do not observe any long branch attraction effect. In fact, in most cases (except for the two most extreme coding distributions and the most divergent cases) <sc>dnaml</sc>
<italic>ε</italic> tends to perform better than <sc>dnaml</sc>, and the improvement is larger for the case prone to long branch attraction (<italic>t</italic>
<sub>3</sub> = 10<italic>t</italic>
<sub>2</sub> or <italic>r</italic> = 10) than for the balanced case in which all five branches are equal (<italic>t</italic>
<sub>3</sub> = <italic>t</italic>
<sub>2</sub> or <italic>r</italic> = 1). For the realigned tests, both <sc>dnaml</sc> and <sc>dnaml</sc>
<italic>ε</italic> show a certain long-branch attraction tendency when the amount of substitutions nearly triples that of gaps. For the more extreme coding distribution, the effect seems to be more severe for <sc>dnaml</sc>
<italic>ε</italic>. This effect can be attributed to the column independence assumption. These results show that while there is some systematic long-branch attraction for <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> that effect occurs in some extreme situations which we do not see reproduced in real alignment such as those of rRNA.</p> <title>Concordance test on ribosomal RNA alignments</title> <p>We developed a test to evaluate a phylogenetic inference method on real data, called the “concordance test”. We split an alignment randomly into two disjoint sets of columns, infer a tree on each, and ask if the two trees are identical. The frequency that the same tree is obtained for different subsets of columns from the same alignment should be correlated with phylogenetic inference accuracy, because in general, we only expect to obtain the same tree if it is the correct tree. The concordance test should be well suited for evaluating whether a modified method extracts more information from a given alignment, which is the question at hand here. The test may be less suited for evaluating the absolute accuracy of an entirely new inference program, because some systematic errors, such as long branch attraction, could result in agreement on the wrong tree; simulations could be used to detect such systematic problems, though.</p> <p>We found it is important to choose columns randomly, rather than simply splitting an alignment in half, because real alignments often contain preferentially 5′- or 3′-truncated sequence fragments.</p> <p>In order to evaluate the performance of <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> on real alignments, we applied the concordance test to a large number of ribosomal RNA alignments. We obtained curated rRNA alignments from the Comparative RNA Web Site (CRW; <ext-link ext-link-type="uri" xlink:href="http://www.rna.ccbb.utexas.edu">http://www.rna.ccbb.utexas.edu</ext-link>) ##UREF##25##[82]## for different domains of life – five small subunit (SSU) alignments and four large subunit (LSU) alignments – and filtered out sequences more than 95% identical to another. These datasets are summarized in ##TAB##4##Table 5##. We randomly sampled a large number of eight-taxon subalignments from these datasets and applied the concordance test with the <sc>dnaml</sc>
<italic>ε</italic> and <sc>dnaml</sc> methods.</p> <p>Results are summarized in ##FIG##3##Figure 4##. Overall, <sc>dnaml</sc>
<italic>ε</italic> shows tree concordance of 27.9% for SSU and 46.6% for LSU, while <sc>dnaml</sc> shows tree concordance in 16.9% for SSU and 35.7% for LSU. The error estimate for all these results is about 0.5–0.6%, which indicates that the improvement obtained by <sc>dnaml</sc>
<italic>ε</italic> is significant. LSU alignments are longer than SSU (4205±1179 versus 1959±579), probably explaining the better performance. For alignments with few gaps, the two methods produce similar results. The improvement of <sc>dnaml</sc>
<italic>ε</italic> over <sc>dnaml</sc> increases with the frequency of gaps in the alignments.</p> <p>With respect to computational time, ##FIG##3##Figure 4D## shows that both methods scale similarly, with <sc>dnaml</sc>
<italic>ε</italic> taking about two- to three-fold longer.</p>
<title>Discussion</title> <p>We have presented a non-reversible probabilistic model of sequence evolution accounting for substitutions, insertions, and deletions that is based on a continuous-time Markov process. This model does not assume a pre-stated number of columns. Rather, it describes a generative evolutionary model of substitution, insertion, and deletion events, starting from an explicit prior distribution over ancestral sequences of any length. This avoids the conceptual flaws that easily arise in column-independent models based on a Markov process that includes the gap character. The model remains compatible with efficient post-order transversal algorithms to calculate the likelihood of a phylogenetic tree. The model can also be used to calculate the probability of ancestral sequences with an arbitrary number of residues that leave no trace in the observed alignment.</p> <p>To do this, the model assumes column independence, an assumption that is problematic because insertions and deletions typically involve multiple residues at the same time. Although we can indeed produce synthetic examples for which our model breaks down due to the column independence assumption, our results for real rRNA alignments show that the gap-extended model is able to produce better trees than the standard model, indicating that the cost of the assumption is outweighed by the gain in modeling gaps instead of ignoring them as missing data. McGuire <italic>et al.</italic> have already made the same observation ##REF##11264399##[41]##. There are ways in which we might relax the column independence assumption; for example, one could extend the ideas of context-dependent substitution rate matrices and context-dependent residue distributions ##REF##15292512##[83]##.</p> <p>The model introduces as free parameters the rates of insertions and deletions, and the geometric probability parameter for the distribution of ancestral sequences. In the results presented here, for simplicity, we have used the F84 substitution rate matrix, identical deletion rates for all residues, and assumed that inserted residues have the same probability distribution as the stationary distribution of the substitution process. Our model does not require these simplifications. Our results may be generalized to arbitrary substitution rate matrices, deletion rates, and inserted residue probability distributions.</p> <p>It would be straightforward to integrate this method into any probabilistic model of sequence alignment, including phylo-HMMs ##REF##15285899##[38]##,##REF##16024819##[84]## and phylo-SCFGs ##REF##10383470##[21]##, in order to account for insertion and deletion events as well as substitution events, without changing the algorithmic complexity of the existing algorithms, and where the insertion and deletion parameters can be inferred from the data ##REF##15780137##[46]##.</p> <title>Comparison to Reversible Models</title> <p>The <sc>tkf</sc>91 model and McGuire's model are designed to be reversible. Our model is non-reversible with respect to insertions and deletions. Reversibility is a mathematically convenient assumption in phylogenetic inference (allowing inference on unrooted trees, for example), and for residue substitution events, usually seem reasonable. However, in the case of insertion/deletion processes, reversibility seem less easy to justify, if one expects the insertion rate <italic>λ</italic> and the deletion rate <italic>μ</italic> to behave as independent parameters.</p> <p>For a model to be reversible with respect to insertions/deletions, obviously constraints must be imposed on the insertion rate <italic>λ</italic> and deletion rate <italic>μ</italic>; moreover, the constraints imposed by reversibility can be counterintuitive. Consider what happens in a case of zero deletion rate and a positive insertion rate. For McGuire's model, reversibility requires that the frequency of gaps is a constant (<italic>π</italic>
<sub>–</sub>), and that it is related to the rates of insertions and deletions by the condition <italic>π</italic>
<sub>–</sub> = <italic>μ</italic>/(<italic>λ</italic>+<italic>μ</italic>). For the <sc>tkf</sc>91 model, reversibility requires a length distribution for evolved sequences identical to that of ancestral sequences which is geometric with parameter <italic>λ</italic>/<italic>μ</italic>
##UREF##26##[85]##. So, in McGuire's model, imposing <italic>μ</italic> = 0 automatically implies that the gap frequency is zero, effectively converting the model into a substitution-only model, regardless of the insertion rate. In <sc>tkf</sc>91, the length distribution breaks down entirely in the <italic>μ</italic>&lt;<italic>λ</italic> regime, and gives arbitrarily large joint “probabilities”. Our model remains valid for any arbitrary (positive or zero) values for insertion and deletion rate.</p> <p>Another consequence of imposing reversibility is that the expected frequencies of insertions and deletions in a pairwise alignment must be identical. For McGuire's model, the expected frequencies of insertions and deletions in a pairwise alignment are given by,
</p> <p>where <italic>π</italic>
<sub>–</sub> is the constant frequency of a gap, and <italic>λ</italic>+<italic>μ</italic> = <italic>β</italic>(1−<italic>π</italic>
<sub>–</sub>), where <italic>β</italic> is one of the two parameters of the F84 substitution model ##REF##8583911##[14]##. For the <sc>tkf</sc>91 model, the expected frequencies of insertions and deletions in a pairwise alignment are given by (see ##SUPPL##4##Text S4##):
</p> <p>Therefore in <sc>tkf</sc>91 for small divergence times, the expected frequency of <italic>insertions</italic> is solely dependent on the <italic>deletion</italic> rate <italic>μ</italic>; indeed, Equation 31 does not depend on the insertion rate at all. If one desires insertion and deletion rates to be independent parameters, it is more logical to expect that the observed frequencies of deletions and insertions should be in general different, and should depend on both deletion and insertion rates. This is the case in our model, shown in Equations 16 and 17. At one extreme, for <italic>λ</italic> = 0 (zero insertion rate),  and . At another extreme, for <italic>λ</italic> = ∞, ,  for any positive time <italic>t</italic>.</p> <p>However, this is not to say that our model is free of its own problems. Most significantly, our model's assumption of column independence is problematic. Rather, in return for this simplification, our model's advantage is that it allows computationally efficient likelihood inference while using a birth-death generative model allowing arbitrary rates of insertion and deletion.</p>
<p>Conceived and designed the experiments: ER SRE. Performed the experiments: ER. Analyzed the data: ER. Wrote the paper: ER SRE.</p> <p>A fundamental task in sequence analysis is to calculate the probability of a multiple alignment given a phylogenetic tree relating the sequences and an evolutionary model describing how sequences change over time. However, the most widely used phylogenetic models only account for residue substitution events. We describe a probabilistic model of a multiple sequence alignment that accounts for insertion and deletion events in addition to substitutions, given a phylogenetic tree, using a rate matrix augmented by the gap character. Starting from a continuous Markov process, we construct a non-reversible generative (birth–death) evolutionary model for insertions and deletions. The model assumes that insertion and deletion events occur one residue at a time. We apply this model to phylogenetic tree inference by extending the program <sc>dnaml</sc> in <sc>phylip</sc>. Using standard benchmarking methods on simulated data and a new “concordance test” benchmark on real ribosomal RNA alignments, we show that the extended program <sc>dnaml</sc>
<italic>ε</italic> improves accuracy relative to the usual approach of ignoring gaps, while retaining the computational efficiency of the Felsenstein peeling algorithm.</p> <title>Author Summary</title> <p>We describe a computationally efficient method to use insertion and deletion events, in addition to substitutions, in phylogenetic inference. To date, many evolutionary models in probabilistic phylogenetic inference methods have only accounted for substitution events, not for insertions and deletions. As a result, not only do tree inference methods use less sequence information than they could, but also it has remained difficult to integrate phylogenetic modeling into sequence alignment methods (such as profiles and profile-hidden Markov models) that inherently require a model of insertion and deletion events. Therefore an important goal in the field has been to develop tractable evolutionary models of insertion/deletion events over time of sufficient accuracy to increase the resolution of phylogenetic inference methods and to increase the power of profile-based sequence homology searches. Our model offers a partial answer to this problem. We show that our model generally improves inference power in both simulated and real data and that it is easily implemented in the framework of standard inference packages with little effect on computational efficiency (we extended <sc>dnaml</sc>, in Felsenstein's popular <sc>phylip</sc> package).</p>
<title>Supporting Information</title>
<p>This work was conceived and partially implemented at the Benasque Center for Science, Spain, during the workshop, “Computational Methods for Nonconding RNAs,” July, 2006. Thanks to Eric Nawrocki for a Stockholm-formatted version of the rRNA alignments from the Comparative RNA Web Site.</p>
<fig id="pcbi-1000172-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.g001</object-id><label>Figure 1</label><caption><title>Felsenstein's peeling algorithm extended to gaps.</title><p>Graphical description of the extended Felsenstein algorithm to calculate the probability of an alignment given a phylogenetic tree using the generative model for gaps presented in this work, given by Equations 20 and 21. For a given column in a multiple alignment, recursion A corresponds to the probability of a tree up to node <italic>k</italic> with a residue at the node (Equation 20). Recursion B corresponds to the probability of a tree up to node <italic>k</italic> with a gap at the node (Equation. 21). The algorithm needs to consider only evolutionarily correct events, thus recursion A includes substitutions and deletions but no insertions, and recursion B has to include one and only one insertion. In the case of having a residue <italic>i</italic> at node <italic>k</italic> (recursion A), the A1 term corresponds to the original substitution-only Felsenstein algorithm, where <italic>q</italic> and <italic>s</italic> are residues, and ∑<italic>q</italic> and ∑<italic>s</italic> stand for the sum to all possible residue substitutions. The terms A2 and A3 represent a deletion occurring for the right and left child respectively (and a substitution for the other child node). Because no insertion can occur in this recursion, there are no evolutionary events happening for descendants for the child that suffers the deletion. Thus, for a A2 or A3 term to have a contribution all nodes down to the leaves for the child that suffers the deletion have to be gaps. (We represent with gray the situation of no evolutionary event happening, and with an empty gray triangle with gaps at the bottom the situation of gaps at all nodes including the leaves for a given subtree.) The A4 term corresponds to a deletion for both children nodes, and no evolutionary event from there on. A A4 term contributes only if with all the leaves under node <italic>k</italic> are gaps. In the case of having a gap at node <italic>k</italic> (recursion B), the insertion might occur for the left or right child, which is represented by graphs B1 and B2 respectively. The insertion might instead be delayed to a node under the left or right child, which is represented by graphs B3 and B4 respectively. In all four cases, the child node for which no insertion occurs does not have any evolutionary events, and has to include all gaps at the internal nodes and leaves. The substitution, deletion and insertion probabilities are given by the generative model as per Equations 22–24, respectively.</p></caption></fig> <fig id="pcbi-1000172-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.g002</object-id><label>Figure 2</label><caption><title>Comparison of <sc>dnaml</sc> versus <sc>dnaml</sc>
<italic>ε</italic> for ungapped alignments.</title><p>Tree reconstruction for ungapped alignments generated according to a F84 substitution model (2.0 transition to transversion ratio and equiprobable residues), for nine different time divergences, ranging from 0.005 to 2.0 substitutions per site and branch. For a given divergence value, 100 random trees with eight taxa were used. For each tree, single alignments were generated with lengths ranging from 50 to 1000 residues in 5 residue increments. For each alignment, a tree was inferred using the programs <sc>dnaml</sc> and <sc>dnaml</sc>
<italic>ε</italic>. Results are displayed as a function of the length of the alignments. (A) Fraction of trees which topology was correctly inferred as a function of the alignment length. The best performance occurs for alignments that contain about 18% pairwise substitutions on average (0.05 substitutions per site and branch). In this case, detectability seems to asymptote to approximately 78% for alignments of at least 800 residues. (B,C) Corresponding results when using the SDD and nBSD measures respectively. (D) Average mean branch length for each length bin. Overall, the two methods show similar performance for ungapped alignment. We mark with an arrow some extreme cases in which the two methods perform differently when inferring the tree branch lengths. (E) Comparison of computational time performance.</p></caption></fig> <fig id="pcbi-1000172-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.g003</object-id><label>Figure 3</label><caption><title>Comparison of <sc>dnaml</sc> versus <sc>dnaml</sc>
<italic>ε</italic> for synthetic alignments with gaps.</title><p>Tree reconstruction test for simulated alignments with gaps generated with the program <sc>rose</sc>
##REF##9545448##[76]## according to 100 random eight-taxon trees, using the F84 model for substitutions and a Poisson gap length distribution (, for <italic>λ</italic> = 0.5). Alignments lengths range from 50 to 1,000 residues. Here we compare the performance of <sc>dnaml</sc> (left) versus <sc>dnaml</sc>
<italic>ε</italic> (right) for two different evolutionary situations: in red alignments with no gaps and 0.06 substitutions per site and branch, which produce alignments with 21%±6% pairwise average substitutions; in blue alignments with gaps (<italic>p<sub>ins</sub></italic> = <italic>p<sub>del</sub></italic> = 0.001 and 0.07 substitutions per site and branch) with a similar percentage of pairwise substitutions (20%±5%) as the ungapped alignments. The alignments with gaps have a percentage of pairwise gaps of 13%±4%. Results are presented as a function of the geometric mean of sequences lengths. The tree reconstruction test assesses the similarity between the inferred tree and the original tree. Three measures of tree similarity are displayed in (A), (B), and (C), respectively: a binary count of whether the trees are topologically identical or not (TP), the Symmetric Difference Distance (SDD) and the normalized Branch Scoring Distance (nBSD). (D) Mean branch length of the inferred trees, and (E) the running time required for the different inferences.</p></caption></fig> <fig id="pcbi-1000172-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.g004</object-id><label>Figure 4</label><caption><title>Comparison of <sc>dnaml</sc> versus <sc>dnaml</sc>
<italic>ε</italic> using the “tree concordance test” on ribosomal RNA alignments.</title><p>Tree concordance test for SSU (left) and LSU (right) rRNA alignments displayed as a function of the total fraction of gaps present in the alignment. We used five SSU and four LSU alignments described in ##TAB##4##Table 5##. For each alignment, we randomly selected a large number of eight taxa alignments (4,000 for the Archaea and Chloroplasts alignments, and 10,000 for the Eukarya, Bacteria and Mitochondria alignments). Each eight taxa alignment was first shuffled and then split in two halves. The tree concordance test assesses the similarity between the two trees inferred for the two sections of the alignment. Three measures of tree similarity are displayed: a binary count of whether the trees are topologically identical or not (TP), the Symmetric Difference Distance (SDD) and the normalized Branch Scoring Distance (nBSD). Results for all SSU (LSU) tests have been summarized together. (A) A histogram of total alignments, as well as the number of TPs for <sc>dnaml</sc> (magenta) and <sc>dnaml</sc>
<italic>ε</italic> (cyan) as a function of the total fraction of gaps in the alignment. (B,C,D) Results for the fraction of TPs, the SDD, and the nBSD respectively. Overall tree concordance for SSU rRNA is 27.9% (10,589/38,000) for <sc>dnaml</sc>
<italic>ε</italic>, versus for 16.9% (6,418/38,000) <sc>dnaml</sc>. Overall tree concordance for LSU rRNA is 46.6% (13,048/28,000) for <sc>dnaml</sc>
<italic>ε</italic> , versus 35.7% (10,002/28,000) for <sc>dnaml</sc>. (E) shows a comparison of time performance. <sc>dnaml</sc>
<italic>ε</italic> shows on average a two to three fold time increase respect to <sc>dnaml</sc> for eight taxa alignments.</p></caption></fig>
<table-wrap id="pcbi-1000172-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.t001</object-id><label>Table 1</label><caption><title>Poisson <italic>λ</italic> = 0.5 gap distribution for 8-taxon alignments.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Ave Subs per Site and Branch</td><td align="left" rowspan="1" colspan="1">
<sc>rose</sc> Gap Parameter</td><td align="left" rowspan="1" colspan="1">Pairwise % ID</td><td align="left" rowspan="1" colspan="1">Pairwise % SUBS</td><td align="left" rowspan="1" colspan="1">Pairwise % GAPS</td><td align="left" rowspan="1" colspan="1">
 (TP)</td><td align="left" rowspan="1" colspan="1">
 (SDD)</td><td align="left" rowspan="1" colspan="1">
(nBSD)</td><td colspan="2" align="left" rowspan="1">Mean MBL</td><td colspan="2" align="left" rowspan="1">Time (L = 1000) (s)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>23.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>18.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>17.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.007</td><td align="left" rowspan="1" colspan="1">2.9±1.7</td><td align="left" rowspan="1" colspan="1">4.7±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>3.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>5.0</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>16.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">0.008</td><td align="left" rowspan="1" colspan="1">2.5±1.1</td><td align="left" rowspan="1" colspan="1">4.6±0.6</td></tr><tr><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">94±2</td><td align="left" rowspan="1" colspan="1">4±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>14.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>16.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>7.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">2.4±0.8</td><td align="left" rowspan="1" colspan="1">5.3±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±2</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">−4.8</td><td align="left" rowspan="1" colspan="1">2.3</td><td align="left" rowspan="1" colspan="1">−<bold>6.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.013</td><td align="left" rowspan="1" colspan="1">2.4±0.9</td><td align="left" rowspan="1" colspan="1">5.2±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">92±3</td><td align="left" rowspan="1" colspan="1">4±1</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">
<bold>22.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>26.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>15.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">2.4±0.6</td><td align="left" rowspan="1" colspan="1">5.9±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">91±3</td><td align="left" rowspan="1" colspan="1">5±2</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">−4.9</td><td align="left" rowspan="1" colspan="1">4.3</td><td align="left" rowspan="1" colspan="1">−<bold>14.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">0.018</td><td align="left" rowspan="1" colspan="1">2.4±0.4</td><td align="left" rowspan="1" colspan="1">5.8±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">88±4</td><td align="left" rowspan="1" colspan="1">8±2</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">
<bold>12.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>19.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>6.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.024</td><td align="left" rowspan="1" colspan="1">2.4±0.2</td><td align="left" rowspan="1" colspan="1">6.0±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">87±4</td><td align="left" rowspan="1" colspan="1">9±3</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">−8.7</td><td align="left" rowspan="1" colspan="1">12.3</td><td align="left" rowspan="1" colspan="1">−<bold>5.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.023</td><td align="left" rowspan="1" colspan="1">0.027</td><td align="left" rowspan="1" colspan="1">2.4±0.2</td><td align="left" rowspan="1" colspan="1">5.9±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">85±5</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">
<bold>20.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>32.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>12.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">2.8±0.6</td><td align="left" rowspan="1" colspan="1">7.0±1.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">83±6</td><td align="left" rowspan="1" colspan="1">10±4</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">−5.2</td><td align="left" rowspan="1" colspan="1">7.1</td><td align="left" rowspan="1" colspan="1">−<bold>11.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.039</td><td align="left" rowspan="1" colspan="1">2.9±0.6</td><td align="left" rowspan="1" colspan="1">7.1±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">86±4</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">
<bold>7.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>14.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>2.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">2.5±0.7</td><td align="left" rowspan="1" colspan="1">6.0±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">85±4</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">−12.8</td><td align="left" rowspan="1" colspan="1">15.9</td><td align="left" rowspan="1" colspan="1">−<bold>1.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">2.7±0.6</td><td align="left" rowspan="1" colspan="1">6.1±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">83±5</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">
<bold>14.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>22.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>6.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">2.7±0.5</td><td align="left" rowspan="1" colspan="1">7.0±1.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">81±5</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">5±2</td><td align="left" rowspan="1" colspan="1">−10.5</td><td align="left" rowspan="1" colspan="1">13.3</td><td align="left" rowspan="1" colspan="1">−<bold>4.8</bold>
</td><td align="left" rowspan="1" colspan="1">0.036</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">2.8±0.5</td><td align="left" rowspan="1" colspan="1">6.9±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">78±6</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">11±4</td><td align="left" rowspan="1" colspan="1">
<bold>18.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>35.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>12.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">2.9±0.6</td><td align="left" rowspan="1" colspan="1">8.2±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">75±8</td><td align="left" rowspan="1" colspan="1">16±5</td><td align="left" rowspan="1" colspan="1">9±3</td><td align="left" rowspan="1" colspan="1">−5.5</td><td align="left" rowspan="1" colspan="1">5.8</td><td align="left" rowspan="1" colspan="1">−<bold>9.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.062</td><td align="left" rowspan="1" colspan="1">3.1±0.4</td><td align="left" rowspan="1" colspan="1">8.6±1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">82±5</td><td align="left" rowspan="1" colspan="1">14±4</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">
<bold>9.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>17.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>2.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">2.7±0.4</td><td align="left" rowspan="1" colspan="1">6.5±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">81±5</td><td align="left" rowspan="1" colspan="1">15±4</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">−16.4</td><td align="left" rowspan="1" colspan="1">19.3</td><td align="left" rowspan="1" colspan="1">−<bold>0.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.043</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">2.7±0.3</td><td align="left" rowspan="1" colspan="1">6.7±0.6</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">78±6</td><td align="left" rowspan="1" colspan="1">14±4</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">
<bold>15.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>26.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>6.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.045</td><td align="left" rowspan="1" colspan="1">3.0±0.5</td><td align="left" rowspan="1" colspan="1">8.0±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">76±7</td><td align="left" rowspan="1" colspan="1">17±5</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">−10.1</td><td align="left" rowspan="1" colspan="1">12.3</td><td align="left" rowspan="1" colspan="1">−<bold>3.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.051</td><td align="left" rowspan="1" colspan="1">0.058</td><td align="left" rowspan="1" colspan="1">3.2±0.8</td><td align="left" rowspan="1" colspan="1">7.9±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0015</td><td align="left" rowspan="1" colspan="1">75±7</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">
<bold>16.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>33.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>8.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">3.1±0.4</td><td align="left" rowspan="1" colspan="1">8.9±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">72±8</td><td align="left" rowspan="1" colspan="1">19±6</td><td align="left" rowspan="1" colspan="1">9±3</td><td align="left" rowspan="1" colspan="1">−7.1</td><td align="left" rowspan="1" colspan="1">7.8</td><td align="left" rowspan="1" colspan="1">−<bold>6.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.061</td><td align="left" rowspan="1" colspan="1">0.071</td><td align="left" rowspan="1" colspan="1">3.4±0.6</td><td align="left" rowspan="1" colspan="1">9.1±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">73±7</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">15±4</td><td align="left" rowspan="1" colspan="1">
<bold>20.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>37.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>10.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.052</td><td align="left" rowspan="1" colspan="1">3.4±0.5</td><td align="left" rowspan="1" colspan="1">9.4±1.3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">68±9</td><td align="left" rowspan="1" colspan="1">21±7</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">−4.7</td><td align="left" rowspan="1" colspan="1">4.6</td><td align="left" rowspan="1" colspan="1">−<bold>8.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.072</td><td align="left" rowspan="1" colspan="1">0.085</td><td align="left" rowspan="1" colspan="1">3.7±0.4</td><td align="left" rowspan="1" colspan="1">10.8±1.3</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">71±7</td><td align="left" rowspan="1" colspan="1">22±5</td><td align="left" rowspan="1" colspan="1">7±3</td><td align="left" rowspan="1" colspan="1">
<bold>12.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>21.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>3.0</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.071</td><td align="left" rowspan="1" colspan="1">3.3±0.4</td><td align="left" rowspan="1" colspan="1">9.6±1.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">70±8</td><td align="left" rowspan="1" colspan="1">24±6</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">−11.6</td><td align="left" rowspan="1" colspan="1">15.0</td><td align="left" rowspan="1" colspan="1">−<bold>0.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.080</td><td align="left" rowspan="1" colspan="1">0.085</td><td align="left" rowspan="1" colspan="1">3.4±0.5</td><td align="left" rowspan="1" colspan="1">9.3±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">66±8</td><td align="left" rowspan="1" colspan="1">20±5</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">
<bold>16.0</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>31.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>6.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.074</td><td align="left" rowspan="1" colspan="1">3.7±0.4</td><td align="left" rowspan="1" colspan="1">10.8±1.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">62±10</td><td align="left" rowspan="1" colspan="1">28±8</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">−9.2</td><td align="left" rowspan="1" colspan="1">9.3</td><td align="left" rowspan="1" colspan="1">−<bold>3.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.101</td><td align="left" rowspan="1" colspan="1">0.109</td><td align="left" rowspan="1" colspan="1">3.9±0.3</td><td align="left" rowspan="1" colspan="1">11.6±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0015</td><td align="left" rowspan="1" colspan="1">62±9</td><td align="left" rowspan="1" colspan="1">19±5</td><td align="left" rowspan="1" colspan="1">19±5</td><td align="left" rowspan="1" colspan="1">
<bold>20.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>39.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>8.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.078</td><td align="left" rowspan="1" colspan="1">4.0±0.6</td><td align="left" rowspan="1" colspan="1">12.5±1.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">57±11</td><td align="left" rowspan="1" colspan="1">31±8</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">−7.1</td><td align="left" rowspan="1" colspan="1">5.0</td><td align="left" rowspan="1" colspan="1">−<bold>9.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.126</td><td align="left" rowspan="1" colspan="1">0.135</td><td align="left" rowspan="1" colspan="1">4.7±0.6</td><td align="left" rowspan="1" colspan="1">14.6±1.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">58±10</td><td align="left" rowspan="1" colspan="1">18±4</td><td align="left" rowspan="1" colspan="1">24±7</td><td align="left" rowspan="1" colspan="1">
<bold>23.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>41.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>10.8</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.083</td><td align="left" rowspan="1" colspan="1">4.3±0.7</td><td align="left" rowspan="1" colspan="1">13.5±1.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">52±12</td><td align="left" rowspan="1" colspan="1">35±9</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">−4.4</td><td align="left" rowspan="1" colspan="1">3.1</td><td align="left" rowspan="1" colspan="1">−<bold>11.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.154</td><td align="left" rowspan="1" colspan="1">0.162</td><td align="left" rowspan="1" colspan="1">5.1±0.7</td><td align="left" rowspan="1" colspan="1">16.8±2.4</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">63±9</td><td align="left" rowspan="1" colspan="1">28±6</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">
<bold>17.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>29.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>3.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">4.0±0.7</td><td align="left" rowspan="1" colspan="1">12.1±1.3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">60±10</td><td align="left" rowspan="1" colspan="1">32±8</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">−11.7</td><td align="left" rowspan="1" colspan="1">12.7</td><td align="left" rowspan="1" colspan="1">−<bold>0.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.121</td><td align="left" rowspan="1" colspan="1">0.125</td><td align="left" rowspan="1" colspan="1">4.1±0.5</td><td align="left" rowspan="1" colspan="1">12.0±1.5</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">57±9</td><td align="left" rowspan="1" colspan="1">25±6</td><td align="left" rowspan="1" colspan="1">18±6</td><td align="left" rowspan="1" colspan="1">
<bold>20.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>40.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>6.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">4.5±0.7</td><td align="left" rowspan="1" colspan="1">14.5±1.3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">52±12</td><td align="left" rowspan="1" colspan="1">36±9</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">−7.9</td><td align="left" rowspan="1" colspan="1">6.6</td><td align="left" rowspan="1" colspan="1">−<bold>6.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.160</td><td align="left" rowspan="1" colspan="1">0.165</td><td align="left" rowspan="1" colspan="1">4.7±0.5</td><td align="left" rowspan="1" colspan="1">15.4±2.3</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0015</td><td align="left" rowspan="1" colspan="1">52±10</td><td align="left" rowspan="1" colspan="1">23±5</td><td align="left" rowspan="1" colspan="1">25±6</td><td align="left" rowspan="1" colspan="1">
<bold>25.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>44.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>8.8</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.105</td><td align="left" rowspan="1" colspan="1">4.9±0.6</td><td align="left" rowspan="1" colspan="1">16.6±2.2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">47±12</td><td align="left" rowspan="1" colspan="1">40±10</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">−4.0</td><td align="left" rowspan="1" colspan="1">2.0</td><td align="left" rowspan="1" colspan="1">−<bold>12.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.201</td><td align="left" rowspan="1" colspan="1">0.202</td><td align="left" rowspan="1" colspan="1">5.7±1.0</td><td align="left" rowspan="1" colspan="1">19.0±3.3</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">48±11</td><td align="left" rowspan="1" colspan="1">21±4</td><td align="left" rowspan="1" colspan="1">31±7</td><td align="left" rowspan="1" colspan="1">
<bold>27.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>50.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>11.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.110</td><td align="left" rowspan="1" colspan="1">5.4±0.7</td><td align="left" rowspan="1" colspan="1">18.3±2.2</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">43±12</td><td align="left" rowspan="1" colspan="1">43±10</td><td align="left" rowspan="1" colspan="1">14±4</td><td align="left" rowspan="1" colspan="1">
<bold>0.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>0.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>16.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.242</td><td align="left" rowspan="1" colspan="1">0.234</td><td align="left" rowspan="1" colspan="1">6.0±1.0</td><td align="left" rowspan="1" colspan="1">20.7±3.1</td></tr><tr><td colspan="12" align="left" rowspan="1">

</td></tr><tr><td colspan="12" align="left" rowspan="1">For a given method M, the area under the curve (<sc>auc</sc>): ,</td></tr><tr><td colspan="12" align="left" rowspan="1">where:</td></tr><tr><td colspan="12" align="left" rowspan="1">
, for alignments  with geometric mean of sequence length <sc>l</sc>, and Δ<sc>l</sc> = 5 nts, for these experiments.</td></tr><tr><td colspan="12" align="left" rowspan="1">
<bold>Nomenclature:</bold>
</td></tr><tr><td colspan="12" align="left" rowspan="1">TP = fraction of true positive trees.</td></tr><tr><td colspan="12" align="left" rowspan="1">SDD = Symmetric Difference Distance.</td></tr><tr><td colspan="12" align="left" rowspan="1">NBSD = normalized Branch Score Distance.</td></tr><tr><td colspan="12" align="left" rowspan="1">MBL = mean branch length.</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000172-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.t002</object-id><label>Table 2</label><caption><title>Geometric gap distribution for 8-taxon alignments.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Ave Subs per Site and Branch</td><td align="left" rowspan="1" colspan="1">
<italic>λ</italic> = <italic>μ</italic> Parameter</td><td align="left" rowspan="1" colspan="1">Pairwise % ID</td><td align="left" rowspan="1" colspan="1">Pairwise % SUBS</td><td align="left" rowspan="1" colspan="1">Pairwise % GAPS</td><td align="left" rowspan="1" colspan="1">
 (TP)</td><td align="left" rowspan="1" colspan="1">
 (SDD)</td><td align="left" rowspan="1" colspan="1">
 (nBSD)</td><td colspan="2" align="left" rowspan="1">Mean MBL</td><td colspan="2" align="left" rowspan="1">Time (L = 1000) (s)</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>28.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>27.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>24.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.007</td><td align="left" rowspan="1" colspan="1">2.8±1.5</td><td align="left" rowspan="1" colspan="1">4.4±0.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">95±2</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>7.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>7.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>22.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.006</td><td align="left" rowspan="1" colspan="1">0.009</td><td align="left" rowspan="1" colspan="1">2.2±0.5</td><td align="left" rowspan="1" colspan="1">4.5±0.6</td></tr><tr><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">94±2</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>19.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>21.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">2.2±0.7</td><td align="left" rowspan="1" colspan="1">5.0±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">93±2</td><td align="left" rowspan="1" colspan="1">5±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">−3.2</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">−<bold>11.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">2.2±0.2</td><td align="left" rowspan="1" colspan="1">5.0±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">92±3</td><td align="left" rowspan="1" colspan="1">4±1</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">
<bold>28.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>36.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>20.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.014</td><td align="left" rowspan="1" colspan="1">2.8±1.3</td><td align="left" rowspan="1" colspan="1">5.7±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">91±3</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">
<bold>0.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.2</td><td align="left" rowspan="1" colspan="1">−<bold>17.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.015</td><td align="left" rowspan="1" colspan="1">0.021</td><td align="left" rowspan="1" colspan="1">2.4±0.3</td><td align="left" rowspan="1" colspan="1">6.0±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">88±4</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">
<bold>16.0</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>27.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>11.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.024</td><td align="left" rowspan="1" colspan="1">2.5±0.4</td><td align="left" rowspan="1" colspan="1">5.9±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">87±4</td><td align="left" rowspan="1" colspan="1">9±3</td><td align="left" rowspan="1" colspan="1">4±1</td><td align="left" rowspan="1" colspan="1">−7.0</td><td align="left" rowspan="1" colspan="1">7.6</td><td align="left" rowspan="1" colspan="1">−<bold>8.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.025</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">2.5±0.3</td><td align="left" rowspan="1" colspan="1">6.1±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">85±4</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">22.3</td><td align="left" rowspan="1" colspan="1">−40.6</td><td align="left" rowspan="1" colspan="1">−18.0</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.028</td><td align="left" rowspan="1" colspan="1">2.7±0.4</td><td align="left" rowspan="1" colspan="1">6.7±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">86±4</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">0.0</td><td align="left" rowspan="1" colspan="1">−<bold>11.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.036</td><td align="left" rowspan="1" colspan="1">0.045</td><td align="left" rowspan="1" colspan="1">2.4±0.3</td><td align="left" rowspan="1" colspan="1">5.7±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.06</td><td align="left" rowspan="1" colspan="1">86±4</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>12.0</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>18.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>5.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">2.5±0.4</td><td align="left" rowspan="1" colspan="1">5.9±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">86±4</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">−9.3</td><td align="left" rowspan="1" colspan="1">11.3</td><td align="left" rowspan="1" colspan="1">−<bold>3.5</bold>
</td><td align="left" rowspan="1" colspan="1">0.032</td><td align="left" rowspan="1" colspan="1">0.053</td><td align="left" rowspan="1" colspan="1">2.4±0.3</td><td align="left" rowspan="1" colspan="1">5.7±0.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">83±5</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">
<bold>17.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>33.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>10.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.035</td><td align="left" rowspan="1" colspan="1">2.6±0.5</td><td align="left" rowspan="1" colspan="1">6.7±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">81±6</td><td align="left" rowspan="1" colspan="1">14±5</td><td align="left" rowspan="1" colspan="1">5±2</td><td align="left" rowspan="1" colspan="1">−5.0</td><td align="left" rowspan="1" colspan="1">6.7</td><td align="left" rowspan="1" colspan="1">−<bold>6.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.041</td><td align="left" rowspan="1" colspan="1">0.048</td><td align="left" rowspan="1" colspan="1">2.9±0.4</td><td align="left" rowspan="1" colspan="1">6.9±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">78±6</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">
<bold>22.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>45.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>16.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">3.0±0.4</td><td align="left" rowspan="1" colspan="1">8.2±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">73±8</td><td align="left" rowspan="1" colspan="1">19±6</td><td align="left" rowspan="1" colspan="1">8±2</td><td align="left" rowspan="1" colspan="1">
<bold>1.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>1.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>9.5</bold>
</td><td align="left" rowspan="1" colspan="1">0.060</td><td align="left" rowspan="1" colspan="1">0.071</td><td align="left" rowspan="1" colspan="1">3.3±0.3</td><td align="left" rowspan="1" colspan="1">9.0±1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.06</td><td align="left" rowspan="1" colspan="1">82±5</td><td align="left" rowspan="1" colspan="1">14±4</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">
<bold>11.0</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>22.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>5.0</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.042</td><td align="left" rowspan="1" colspan="1">2.7±0.3</td><td align="left" rowspan="1" colspan="1">6.5±0.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">81±5</td><td align="left" rowspan="1" colspan="1">16±5</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">−10.4</td><td align="left" rowspan="1" colspan="1">14.6</td><td align="left" rowspan="1" colspan="1">−<bold>2.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.044</td><td align="left" rowspan="1" colspan="1">0.047</td><td align="left" rowspan="1" colspan="1">2.7±0.4</td><td align="left" rowspan="1" colspan="1">6.5±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">78±6</td><td align="left" rowspan="1" colspan="1">14±4</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">
<bold>18.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>35.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>9.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.045</td><td align="left" rowspan="1" colspan="1">2.9±0.4</td><td align="left" rowspan="1" colspan="1">7.7±0.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">75±7</td><td align="left" rowspan="1" colspan="1">19±6</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">−4.1</td><td align="left" rowspan="1" colspan="1">4.4</td><td align="left" rowspan="1" colspan="1">−<bold>4.8</bold>
</td><td align="left" rowspan="1" colspan="1">0.058</td><td align="left" rowspan="1" colspan="1">0.066</td><td align="left" rowspan="1" colspan="1">3.1±0.4</td><td align="left" rowspan="1" colspan="1">8.2±1.0</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">75±7</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">
<bold>20.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>42.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.049</td><td align="left" rowspan="1" colspan="1">3.2±0.5</td><td align="left" rowspan="1" colspan="1">8.8±1.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">70±8</td><td align="left" rowspan="1" colspan="1">22±7</td><td align="left" rowspan="1" colspan="1">8±2</td><td align="left" rowspan="1" colspan="1">−0.6</td><td align="left" rowspan="1" colspan="1">0.6</td><td align="left" rowspan="1" colspan="1">−<bold>7.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.074</td><td align="left" rowspan="1" colspan="1">0.085</td><td align="left" rowspan="1" colspan="1">3.6±0.4</td><td align="left" rowspan="1" colspan="1">9.6±1.2</td></tr><tr><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">73±7</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">15±4</td><td align="left" rowspan="1" colspan="1">
<bold>23.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>46.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.9</bold>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.052</td><td align="left" rowspan="1" colspan="1">3.3±0.5</td><td align="left" rowspan="1" colspan="1">9.3±1.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">66±9</td><td align="left" rowspan="1" colspan="1">25±7</td><td align="left" rowspan="1" colspan="1">9±3</td><td align="left" rowspan="1" colspan="1">
<bold>2.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>1.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>8.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.085</td><td align="left" rowspan="1" colspan="1">0.097</td><td align="left" rowspan="1" colspan="1">4.0±0.5</td><td align="left" rowspan="1" colspan="1">11.0±1.6</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.06</td><td align="left" rowspan="1" colspan="1">72±7</td><td align="left" rowspan="1" colspan="1">22±5</td><td align="left" rowspan="1" colspan="1">7±3</td><td align="left" rowspan="1" colspan="1">
<bold>14.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>30.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>4.6</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.071</td><td align="left" rowspan="1" colspan="1">3.3±0.6</td><td align="left" rowspan="1" colspan="1">8.9±0.9</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">70±8</td><td align="left" rowspan="1" colspan="1">25±7</td><td align="left" rowspan="1" colspan="1">5±2</td><td align="left" rowspan="1" colspan="1">−7.5</td><td align="left" rowspan="1" colspan="1">10.4</td><td align="left" rowspan="1" colspan="1">−<bold>1.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.082</td><td align="left" rowspan="1" colspan="1">0.086</td><td align="left" rowspan="1" colspan="1">3.4±0.4</td><td align="left" rowspan="1" colspan="1">9.1±1.1</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">67±7</td><td align="left" rowspan="1" colspan="1">20±5</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">
<bold>23.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>42.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>9.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.074</td><td align="left" rowspan="1" colspan="1">3.8±0.7</td><td align="left" rowspan="1" colspan="1">11.3±1.3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">61±10</td><td align="left" rowspan="1" colspan="1">31±8</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">−3.2</td><td align="left" rowspan="1" colspan="1">2.9</td><td align="left" rowspan="1" colspan="1">−<bold>5.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.113</td><td align="left" rowspan="1" colspan="1">0.121</td><td align="left" rowspan="1" colspan="1">4.1±0.6</td><td align="left" rowspan="1" colspan="1">12.2±2.0</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">61±9</td><td align="left" rowspan="1" colspan="1">19±4</td><td align="left" rowspan="1" colspan="1">20±5</td><td align="left" rowspan="1" colspan="1">
<bold>25.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>47.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>11.5</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.079</td><td align="left" rowspan="1" colspan="1">4.2±0.7</td><td align="left" rowspan="1" colspan="1">12.8±1.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">55±11</td><td align="left" rowspan="1" colspan="1">35±9</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">−4.4</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">−<bold>7.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.147</td><td align="left" rowspan="1" colspan="1">0.154</td><td align="left" rowspan="1" colspan="1">4.8±0.8</td><td align="left" rowspan="1" colspan="1">15.0±2.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">58±9</td><td align="left" rowspan="1" colspan="1">18±4</td><td align="left" rowspan="1" colspan="1">24±6</td><td align="left" rowspan="1" colspan="1">
<bold>24.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>47.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.2</bold>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.083</td><td align="left" rowspan="1" colspan="1">4.4±0.5</td><td align="left" rowspan="1" colspan="1">13.9±1.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">51±11</td><td align="left" rowspan="1" colspan="1">38±9</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">−3.0</td><td align="left" rowspan="1" colspan="1">1.8</td><td align="left" rowspan="1" colspan="1">−<bold>10.1</bold>
</td><td align="left" rowspan="1" colspan="1">0.170</td><td align="left" rowspan="1" colspan="1">0.176</td><td align="left" rowspan="1" colspan="1">5.0±0.6</td><td align="left" rowspan="1" colspan="1">16.5±2.8</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.06</td><td align="left" rowspan="1" colspan="1">64±8</td><td align="left" rowspan="1" colspan="1">28±6</td><td align="left" rowspan="1" colspan="1">8±2</td><td align="left" rowspan="1" colspan="1">
<bold>16.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>34.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>5.5</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.097</td><td align="left" rowspan="1" colspan="1">3.8±0.6</td><td align="left" rowspan="1" colspan="1">11.2±1.5</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">61±10</td><td align="left" rowspan="1" colspan="1">33±8</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">−11.6</td><td align="left" rowspan="1" colspan="1">10.1</td><td align="left" rowspan="1" colspan="1">−<bold>0.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.123</td><td align="left" rowspan="1" colspan="1">0.127</td><td align="left" rowspan="1" colspan="1">3.9±0.4</td><td align="left" rowspan="1" colspan="1">11.7±1.7</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">57±10</td><td align="left" rowspan="1" colspan="1">25±6</td><td align="left" rowspan="1" colspan="1">18±5</td><td align="left" rowspan="1" colspan="1">
<bold>22.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>44.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>9.8</bold>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">4.5±0.6</td><td align="left" rowspan="1" colspan="1">14.8±4.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">51±11</td><td align="left" rowspan="1" colspan="1">39±9</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">−6.1</td><td align="left" rowspan="1" colspan="1">3.4</td><td align="left" rowspan="1" colspan="1">−<bold>7.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.175</td><td align="left" rowspan="1" colspan="1">0.179</td><td align="left" rowspan="1" colspan="1">5.0±0.5</td><td align="left" rowspan="1" colspan="1">16.1±3.2</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">51±10</td><td align="left" rowspan="1" colspan="1">22±5</td><td align="left" rowspan="1" colspan="1">26±6</td><td align="left" rowspan="1" colspan="1">
<bold>29.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>50.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>11.5</bold>
</td><td align="left" rowspan="1" colspan="1">0.101</td><td align="left" rowspan="1" colspan="1">0.105</td><td align="left" rowspan="1" colspan="1">5.3±0.8</td><td align="left" rowspan="1" colspan="1">17.0±2.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">46±12</td><td align="left" rowspan="1" colspan="1">43±10</td><td align="left" rowspan="1" colspan="1">11±3</td><td align="left" rowspan="1" colspan="1">
<bold>4.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>2.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>12.4</bold>
</td><td align="left" rowspan="1" colspan="1">0.223</td><td align="left" rowspan="1" colspan="1">0.220</td><td align="left" rowspan="1" colspan="1">5.6±0.7</td><td align="left" rowspan="1" colspan="1">19.1±4.1</td></tr><tr><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">47±11</td><td align="left" rowspan="1" colspan="1">21±4</td><td align="left" rowspan="1" colspan="1">32±7</td><td align="left" rowspan="1" colspan="1">
<bold>29.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>53.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.102</td><td align="left" rowspan="1" colspan="1">0.110</td><td align="left" rowspan="1" colspan="1">5.4±0.8</td><td align="left" rowspan="1" colspan="1">18.2±2.4</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">43±12</td><td align="left" rowspan="1" colspan="1">46±10</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">
<bold>4.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>3.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>15.3</bold>
</td><td align="left" rowspan="1" colspan="1">0.254</td><td align="left" rowspan="1" colspan="1">0.243</td><td align="left" rowspan="1" colspan="1">5.8±0.8</td><td align="left" rowspan="1" colspan="1">21.0±3.8</td></tr><tr><td colspan="12" align="left" rowspan="1">

</td></tr><tr><td colspan="12" align="left" rowspan="1">For a given method M, the area under the curve (<sc>auc</sc>): 
</td></tr><tr><td colspan="12" align="left" rowspan="1">where:</td></tr><tr><td colspan="12" align="left" rowspan="1">
, for alignments  with geometric mean of sequence length <sc>l</sc>, and Δ<sc>l</sc> = 5 nts, for these experiments.</td></tr><tr><td colspan="12" align="left" rowspan="1">
<bold>Nomenclature:</bold>
</td></tr><tr><td colspan="12" align="left" rowspan="1">TP = fraction of true positive trees.</td></tr><tr><td colspan="12" align="left" rowspan="1">SDD = Symmetric Difference Distance.</td></tr><tr><td colspan="12" align="left" rowspan="1">NBSD = normalized Branch Score Distance.</td></tr><tr><td colspan="12" align="left" rowspan="1">MBL = mean branch length.</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000172-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.t003</object-id><label>Table 3</label><caption><title>8-Taxon alignments.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col width="98*" align="left" span="1"/><col width="98*" align="center" span="1"/><col width="98*" align="center" span="1"/><col width="98*" align="center" span="1"/><col width="98*" align="center" span="1"/><col width="98*" align="center" span="1"/><col width="98*" align="center" span="1"/><col width="111*" align="center" span="1"/><col width="93*" align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gap Length Distribution</td><td align="left" rowspan="1" colspan="1">Ave Subs per Site and Branch</td><td align="left" rowspan="1" colspan="1">Gap Parameter</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>id</sc>
</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>subs</sc>
</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>gaps</sc>
</td><td align="left" rowspan="1" colspan="1">
 (TP)</td><td align="left" rowspan="1" colspan="1">

</td><td align="left" rowspan="1" colspan="1">

</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>28.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>27.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>24.9</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">
<bold>24.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>19.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>16.8</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">96±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±2</td><td align="left" rowspan="1" colspan="1">
<bold>4.8</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>2.5</bold>
</td><td align="left" rowspan="1" colspan="1">11.8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">96±4</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±4</td><td align="left" rowspan="1" colspan="1">−4.5</td><td align="left" rowspan="1" colspan="1">5.0</td><td align="left" rowspan="1" colspan="1">36.8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">96±4</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±4</td><td align="left" rowspan="1" colspan="1">−8.1</td><td align="left" rowspan="1" colspan="1">7.0</td><td align="left" rowspan="1" colspan="1">41.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">85±4</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">
<bold>22.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>40.6</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>18.0</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">85±5</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">
<bold>20.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>32.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>12.3</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">85±5</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">7±4</td><td align="left" rowspan="1" colspan="1">
<bold>2.9</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>0.5</bold>
</td><td align="left" rowspan="1" colspan="1">11.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">84±8</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">9±7</td><td align="left" rowspan="1" colspan="1">−11.6</td><td align="left" rowspan="1" colspan="1">19.9</td><td align="left" rowspan="1" colspan="1">43.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.020</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">85±8</td><td align="left" rowspan="1" colspan="1">7±2</td><td align="left" rowspan="1" colspan="1">7±8</td><td align="left" rowspan="1" colspan="1">−15.1</td><td align="left" rowspan="1" colspan="1">20.9</td><td align="left" rowspan="1" colspan="1">51.6</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">73±7</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">15±4</td><td align="left" rowspan="1" colspan="1">
<bold>23.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>46.5</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.9</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">73±7</td><td align="left" rowspan="1" colspan="1">12±3</td><td align="left" rowspan="1" colspan="1">15±4</td><td align="left" rowspan="1" colspan="1">
<bold>20.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>37.2</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>10.6</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">73±8</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">14±5</td><td align="left" rowspan="1" colspan="1">
<bold>0.7</bold>
</td><td align="left" rowspan="1" colspan="1">0.2</td><td align="left" rowspan="1" colspan="1">8.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">72±10</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">16±10</td><td align="left" rowspan="1" colspan="1">−21.9</td><td align="left" rowspan="1" colspan="1">31.3</td><td align="left" rowspan="1" colspan="1">40.4</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.040</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">73±11</td><td align="left" rowspan="1" colspan="1">13±4</td><td align="left" rowspan="1" colspan="1">14±11</td><td align="left" rowspan="1" colspan="1">−26.9</td><td align="left" rowspan="1" colspan="1">36.4</td><td align="left" rowspan="1" colspan="1">50.1</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">58±9</td><td align="left" rowspan="1" colspan="1">18±4</td><td align="left" rowspan="1" colspan="1">24±6</td><td align="left" rowspan="1" colspan="1">
<bold>24.1</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>47.3</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>13.2</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">58±10</td><td align="left" rowspan="1" colspan="1">18±4</td><td align="left" rowspan="1" colspan="1">24±7</td><td align="left" rowspan="1" colspan="1">
<bold>23.7</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>41.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>10.8</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">60±10</td><td align="left" rowspan="1" colspan="1">18±4</td><td align="left" rowspan="1" colspan="1">22±7</td><td align="left" rowspan="1" colspan="1">
<bold>2.4</bold>
</td><td align="left" rowspan="1" colspan="1">−<bold>2.4</bold>
</td><td align="left" rowspan="1" colspan="1">5.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">57±12</td><td align="left" rowspan="1" colspan="1">17±5</td><td align="left" rowspan="1" colspan="1">26±12</td><td align="left" rowspan="1" colspan="1">−28.6</td><td align="left" rowspan="1" colspan="1">34.5</td><td align="left" rowspan="1" colspan="1">38.0</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.070</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">59±13</td><td align="left" rowspan="1" colspan="1">18±5</td><td align="left" rowspan="1" colspan="1">23±13</td><td align="left" rowspan="1" colspan="1">−39.4</td><td align="left" rowspan="1" colspan="1">50.1</td><td align="left" rowspan="1" colspan="1">48.9</td></tr><tr><td colspan="9" align="left" rowspan="1">
<bold>Nomenclature:</bold>
</td></tr><tr><td colspan="9" align="left" rowspan="1">
<sc>auc</sc> = area under the curve</td></tr><tr><td colspan="9" align="left" rowspan="1">TP = fraction of true positive trees.</td></tr><tr><td colspan="9" align="left" rowspan="1">SDD = Symmetric Difference Distance.</td></tr><tr><td colspan="9" align="left" rowspan="1">NBSD = normalized Branch Score Distance.</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000172-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.t004</object-id><label>Table 4</label><caption><title>4-Taxon alignments.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">Gap Length Distribution</td><td align="left" rowspan="1" colspan="1">Ave Subs per Site and Branch</td><td align="left" rowspan="1" colspan="1">Gap Parameter</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>id</sc>
</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>subs</sc>
</td><td align="left" rowspan="1" colspan="1">Pairwise % <sc>gaps</sc>
</td><td align="left" rowspan="1" colspan="1">〈TP〉<sub>%</sub>
<italic>r</italic> = 1</td><td align="left" rowspan="1" colspan="1">〈TP〉<sub>%</sub>
<italic>r</italic> = 1</td><td align="left" rowspan="1" colspan="1">〈TP〉<sub>%</sub>
<italic>r</italic> = 10</td><td align="left" rowspan="1" colspan="1">〈TP〉<sub>%</sub>
<italic>r</italic> = 10</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
</td><td align="left" rowspan="1" colspan="1">
<sc>dnaml</sc>
<italic>ε</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">94.3</td><td align="left" rowspan="1" colspan="1">96.2</td><td align="left" rowspan="1" colspan="1">67.1</td><td align="left" rowspan="1" colspan="1">80.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">93.6</td><td align="left" rowspan="1" colspan="1">95.6</td><td align="left" rowspan="1" colspan="1">64.0</td><td align="left" rowspan="1" colspan="1">69.4</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">93.5</td><td align="left" rowspan="1" colspan="1">95.1</td><td align="left" rowspan="1" colspan="1">65.2</td><td align="left" rowspan="1" colspan="1">76.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">93.5</td><td align="left" rowspan="1" colspan="1">94.5</td><td align="left" rowspan="1" colspan="1">63.8</td><td align="left" rowspan="1" colspan="1">69.7</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">91.9</td><td align="left" rowspan="1" colspan="1">92.6</td><td align="left" rowspan="1" colspan="1">63.3</td><td align="left" rowspan="1" colspan="1">72.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">92.2</td><td align="left" rowspan="1" colspan="1">92.4</td><td align="left" rowspan="1" colspan="1">63.1</td><td align="left" rowspan="1" colspan="1">69.8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">97±2</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">2±3</td><td align="left" rowspan="1" colspan="1">91.4</td><td align="left" rowspan="1" colspan="1">91.6</td><td align="left" rowspan="1" colspan="1">62.9</td><td align="left" rowspan="1" colspan="1">96.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">97±3</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">2±3</td><td align="left" rowspan="1" colspan="1">91.1</td><td align="left" rowspan="1" colspan="1">91.0</td><td align="left" rowspan="1" colspan="1">61.5</td><td align="left" rowspan="1" colspan="1">68.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.005</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">97±4</td><td align="left" rowspan="1" colspan="1">1±1</td><td align="left" rowspan="1" colspan="1">1±3</td><td align="left" rowspan="1" colspan="1">90.7</td><td align="left" rowspan="1" colspan="1">91.4</td><td align="left" rowspan="1" colspan="1">62.7</td><td align="left" rowspan="1" colspan="1">69.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">97±4</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">1±4</td><td align="left" rowspan="1" colspan="1">91.3</td><td align="left" rowspan="1" colspan="1">90.4</td><td align="left" rowspan="1" colspan="1">62.3</td><td align="left" rowspan="1" colspan="1">69.0</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±2</td><td align="left" rowspan="1" colspan="1">97.4</td><td align="left" rowspan="1" colspan="1">98.3</td><td align="left" rowspan="1" colspan="1">78.2</td><td align="left" rowspan="1" colspan="1">89.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">4±2</td><td align="left" rowspan="1" colspan="1">2±1</td><td align="left" rowspan="1" colspan="1">97.1</td><td align="left" rowspan="1" colspan="1">97.1</td><td align="left" rowspan="1" colspan="1">71.3</td><td align="left" rowspan="1" colspan="1">70.0</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">97.6</td><td align="left" rowspan="1" colspan="1">98.6</td><td align="left" rowspan="1" colspan="1">77.4</td><td align="left" rowspan="1" colspan="1">86.4</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">97.2</td><td align="left" rowspan="1" colspan="1">98.1</td><td align="left" rowspan="1" colspan="1">74.0</td><td align="left" rowspan="1" colspan="1">72.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±2</td><td align="left" rowspan="1" colspan="1">96.9</td><td align="left" rowspan="1" colspan="1">97.3</td><td align="left" rowspan="1" colspan="1">76.0</td><td align="left" rowspan="1" colspan="1">81.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±3</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±2</td><td align="left" rowspan="1" colspan="1">96.8</td><td align="left" rowspan="1" colspan="1">96.8</td><td align="left" rowspan="1" colspan="1">75.4</td><td align="left" rowspan="1" colspan="1">74.4</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">94±5</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±4</td><td align="left" rowspan="1" colspan="1">96.5</td><td align="left" rowspan="1" colspan="1">97.0</td><td align="left" rowspan="1" colspan="1">75.6</td><td align="left" rowspan="1" colspan="1">78.4</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±5</td><td align="left" rowspan="1" colspan="1">3±2</td><td align="left" rowspan="1" colspan="1">3±4</td><td align="left" rowspan="1" colspan="1">96.6</td><td align="left" rowspan="1" colspan="1">96.2</td><td align="left" rowspan="1" colspan="1">75.3</td><td align="left" rowspan="1" colspan="1">75.8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.010</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">94±5</td><td align="left" rowspan="1" colspan="1">3±1</td><td align="left" rowspan="1" colspan="1">3±5</td><td align="left" rowspan="1" colspan="1">96.8</td><td align="left" rowspan="1" colspan="1">96.8</td><td align="left" rowspan="1" colspan="1">74.5</td><td align="left" rowspan="1" colspan="1">77.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">94±6</td><td align="left" rowspan="1" colspan="1">3±2</td><td align="left" rowspan="1" colspan="1">3±5</td><td align="left" rowspan="1" colspan="1">96.4</td><td align="left" rowspan="1" colspan="1">96.4</td><td align="left" rowspan="1" colspan="1">75.1</td><td align="left" rowspan="1" colspan="1">77.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">83±7</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">99.9</td><td align="left" rowspan="1" colspan="1">89.9</td><td align="left" rowspan="1" colspan="1">96.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">80±7</td><td align="left" rowspan="1" colspan="1">14±6</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">99.3</td><td align="left" rowspan="1" colspan="1">99.4</td><td align="left" rowspan="1" colspan="1">65.8</td><td align="left" rowspan="1" colspan="1">
<bold>43.6</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">84±7</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">8±4</td><td align="left" rowspan="1" colspan="1">99.4</td><td align="left" rowspan="1" colspan="1">99.8</td><td align="left" rowspan="1" colspan="1">89.8</td><td align="left" rowspan="1" colspan="1">95.4</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">83±7</td><td align="left" rowspan="1" colspan="1">11±5</td><td align="left" rowspan="1" colspan="1">7±3</td><td align="left" rowspan="1" colspan="1">99.2</td><td align="left" rowspan="1" colspan="1">99.4</td><td align="left" rowspan="1" colspan="1">76.3</td><td align="left" rowspan="1" colspan="1">55.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">84±7</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">8±4</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">90.0</td><td align="left" rowspan="1" colspan="1">90.4</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">83±7</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">7±4</td><td align="left" rowspan="1" colspan="1">99.3</td><td align="left" rowspan="1" colspan="1">98.9</td><td align="left" rowspan="1" colspan="1">85.1</td><td align="left" rowspan="1" colspan="1">66.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">83±9</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">10±8</td><td align="left" rowspan="1" colspan="1">99.2</td><td align="left" rowspan="1" colspan="1">98.6</td><td align="left" rowspan="1" colspan="1">88.4</td><td align="left" rowspan="1" colspan="1">85.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">83±9</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">9±7</td><td align="left" rowspan="1" colspan="1">99.0</td><td align="left" rowspan="1" colspan="1">96.8</td><td align="left" rowspan="1" colspan="1">85.4</td><td align="left" rowspan="1" colspan="1">67.9</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.030</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">83±10</td><td align="left" rowspan="1" colspan="1">8±3</td><td align="left" rowspan="1" colspan="1">8±9</td><td align="left" rowspan="1" colspan="1">98.2</td><td align="left" rowspan="1" colspan="1">97.2</td><td align="left" rowspan="1" colspan="1">89.2</td><td align="left" rowspan="1" colspan="1">85.3</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">83±10</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">8±8</td><td align="left" rowspan="1" colspan="1">98.0</td><td align="left" rowspan="1" colspan="1">95.4</td><td align="left" rowspan="1" colspan="1">86.2</td><td align="left" rowspan="1" colspan="1">74.5</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">75±10</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">14±6</td><td align="left" rowspan="1" colspan="1">99.6</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">92.86</td><td align="left" rowspan="1" colspan="1">98.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">72±10</td><td align="left" rowspan="1" colspan="1">23±8</td><td align="left" rowspan="1" colspan="1">8±2</td><td align="left" rowspan="1" colspan="1">99.3</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">
<bold>46.8</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>25.7</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">75±10</td><td align="left" rowspan="1" colspan="1">11±4</td><td align="left" rowspan="1" colspan="1">13±6</td><td align="left" rowspan="1" colspan="1">99.8</td><td align="left" rowspan="1" colspan="1">99.9</td><td align="left" rowspan="1" colspan="1">92.5</td><td align="left" rowspan="1" colspan="1">97.1</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">73±10</td><td align="left" rowspan="1" colspan="1">19±5</td><td align="left" rowspan="1" colspan="1">10±2</td><td align="left" rowspan="1" colspan="1">99.2</td><td align="left" rowspan="1" colspan="1">99.4</td><td align="left" rowspan="1" colspan="1">62.6</td><td align="left" rowspan="1" colspan="1">
<bold>36.1</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">76±10</td><td align="left" rowspan="1" colspan="1">11±4</td><td align="left" rowspan="1" colspan="1">12±6</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">99.6</td><td align="left" rowspan="1" colspan="1">92.2</td><td align="left" rowspan="1" colspan="1">92.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">74±9</td><td align="left" rowspan="1" colspan="1">15±6</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">99.4</td><td align="left" rowspan="1" colspan="1">98.8</td><td align="left" rowspan="1" colspan="1">82.7</td><td align="left" rowspan="1" colspan="1">52.3</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">74±12</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">15±10</td><td align="left" rowspan="1" colspan="1">99.3</td><td align="left" rowspan="1" colspan="1">98.7</td><td align="left" rowspan="1" colspan="1">90.8</td><td align="left" rowspan="1" colspan="1">84.7</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">72±12</td><td align="left" rowspan="1" colspan="1">14±6</td><td align="left" rowspan="1" colspan="1">13±8</td><td align="left" rowspan="1" colspan="1">98.4</td><td align="left" rowspan="1" colspan="1">95.6</td><td align="left" rowspan="1" colspan="1">82.9</td><td align="left" rowspan="1" colspan="1">55.8</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.050</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">74±12</td><td align="left" rowspan="1" colspan="1">12±4</td><td align="left" rowspan="1" colspan="1">13±11</td><td align="left" rowspan="1" colspan="1">99.0</td><td align="left" rowspan="1" colspan="1">97.6</td><td align="left" rowspan="1" colspan="1">91.4</td><td align="left" rowspan="1" colspan="1">83.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">74±13</td><td align="left" rowspan="1" colspan="1">14±6</td><td align="left" rowspan="1" colspan="1">12±10</td><td align="left" rowspan="1" colspan="1">98.2</td><td align="left" rowspan="1" colspan="1">94.7</td><td align="left" rowspan="1" colspan="1">86.0</td><td align="left" rowspan="1" colspan="1">63.6</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>ε</italic>
<sc>rate</sc>
</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.30</td><td align="left" rowspan="1" colspan="1">58±14</td><td align="left" rowspan="1" colspan="1">19±5</td><td align="left" rowspan="1" colspan="1">24±9</td><td align="left" rowspan="1" colspan="1">99.8</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">92.9</td><td align="left" rowspan="1" colspan="1">98.8</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">53±11</td><td align="left" rowspan="1" colspan="1">38±9</td><td align="left" rowspan="1" colspan="1">10±3</td><td align="left" rowspan="1" colspan="1">98.1</td><td align="left" rowspan="1" colspan="1">98.3</td><td align="left" rowspan="1" colspan="1">
<bold>34.9</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>32.4</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>poisson</sc>
<italic>λ</italic> = 0.5</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">55±14</td><td align="left" rowspan="1" colspan="1">19±5</td><td align="left" rowspan="1" colspan="1">24±9</td><td align="left" rowspan="1" colspan="1">99.6</td><td align="left" rowspan="1" colspan="1">100.0</td><td align="left" rowspan="1" colspan="1">93.7</td><td align="left" rowspan="1" colspan="1">98.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">54±11</td><td align="left" rowspan="1" colspan="1">35±10</td><td align="left" rowspan="1" colspan="1">11±4</td><td align="left" rowspan="1" colspan="1">98.3</td><td align="left" rowspan="1" colspan="1">98.5</td><td align="left" rowspan="1" colspan="1">
<bold>35.1</bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>33.8</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 100</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0020</td><td align="left" rowspan="1" colspan="1">59±14</td><td align="left" rowspan="1" colspan="1">17±5</td><td align="left" rowspan="1" colspan="1">22±9</td><td align="left" rowspan="1" colspan="1">99.8</td><td align="left" rowspan="1" colspan="1">99.7</td><td align="left" rowspan="1" colspan="1">94.0</td><td align="left" rowspan="1" colspan="1">91.5</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">53±11</td><td align="left" rowspan="1" colspan="1">29±9</td><td align="left" rowspan="1" colspan="1">13±5</td><td align="left" rowspan="1" colspan="1">99.0</td><td align="left" rowspan="1" colspan="1">98.0</td><td align="left" rowspan="1" colspan="1">51.4</td><td align="left" rowspan="1" colspan="1">
<bold>43.8</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 6</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0010</td><td align="left" rowspan="1" colspan="1">54±16</td><td align="left" rowspan="1" colspan="1">18±6</td><td align="left" rowspan="1" colspan="1">28±14</td><td align="left" rowspan="1" colspan="1">98.4</td><td align="left" rowspan="1" colspan="1">98.4</td><td align="left" rowspan="1" colspan="1">92.1</td><td align="left" rowspan="1" colspan="1">78.6</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">54±15</td><td align="left" rowspan="1" colspan="1">26±9</td><td align="left" rowspan="1" colspan="1">20±10</td><td align="left" rowspan="1" colspan="1">97.5</td><td align="left" rowspan="1" colspan="1">92.4</td><td align="left" rowspan="1" colspan="1">51.3</td><td align="left" rowspan="1" colspan="1">
<bold>38.9</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<sc>coding</sc>
<italic>c</italic>/<italic>t</italic> = 1</td><td align="left" rowspan="1" colspan="1">0.100</td><td align="left" rowspan="1" colspan="1">0.0005</td><td align="left" rowspan="1" colspan="1">56±16</td><td align="left" rowspan="1" colspan="1">19±6</td><td align="left" rowspan="1" colspan="1">24±15</td><td align="left" rowspan="1" colspan="1">99.0</td><td align="left" rowspan="1" colspan="1">97.4</td><td align="left" rowspan="1" colspan="1">92.5</td><td align="left" rowspan="1" colspan="1">77.0</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">56±15</td><td align="left" rowspan="1" colspan="1">25±9</td><td align="left" rowspan="1" colspan="1">19±11</td><td align="left" rowspan="1" colspan="1">97.9</td><td align="left" rowspan="1" colspan="1">91.4</td><td align="left" rowspan="1" colspan="1">63.8</td><td align="left" rowspan="1" colspan="1">
<bold>44.1</bold>
</td></tr><tr><td colspan="10" align="left" rowspan="1">
<bold>nomenclature:</bold>
</td></tr><tr><td colspan="10" align="left" rowspan="1">〈TP〉<sub>%</sub> = mean percentage of true positive trees.</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pcbi-1000172-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000172.t005</object-id><label>Table 5</label><caption><title>rRNA alignment statistics.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">SSU</td><td align="left" rowspan="1" colspan="1">Archaea</td><td align="left" rowspan="1" colspan="1">Chloroplasts</td><td align="left" rowspan="1" colspan="1">Bacteria</td><td align="left" rowspan="1" colspan="1">Eukarya</td><td align="left" rowspan="1" colspan="1">Mitochondria</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">No. seqs</td><td align="left" rowspan="1" colspan="1">74</td><td align="left" rowspan="1" colspan="1">75</td><td align="left" rowspan="1" colspan="1">1601</td><td align="left" rowspan="1" colspan="1">900</td><td align="left" rowspan="1" colspan="1">617</td></tr><tr><td align="left" rowspan="1" colspan="1">Alignment length</td><td align="left" rowspan="1" colspan="1">1756</td><td align="left" rowspan="1" colspan="1">2456</td><td align="left" rowspan="1" colspan="1">3094</td><td align="left" rowspan="1" colspan="1">7160</td><td align="left" rowspan="1" colspan="1">4716</td></tr><tr><td align="left" rowspan="1" colspan="1">Geometric mean seqs</td><td align="left" rowspan="1" colspan="1">1446</td><td align="left" rowspan="1" colspan="1">1446</td><td align="left" rowspan="1" colspan="1">1497</td><td align="left" rowspan="1" colspan="1">1788</td><td align="left" rowspan="1" colspan="1">1011</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % ID</td><td align="left" rowspan="1" colspan="1">71±7</td><td align="left" rowspan="1" colspan="1">74±11</td><td align="left" rowspan="1" colspan="1">68±6</td><td align="left" rowspan="1" colspan="1">56±14</td><td align="left" rowspan="1" colspan="1">49±20</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % SUBS</td><td align="left" rowspan="1" colspan="1">23±5</td><td align="left" rowspan="1" colspan="1">19±7</td><td align="left" rowspan="1" colspan="1">22±4</td><td align="left" rowspan="1" colspan="1">20±5</td><td align="left" rowspan="1" colspan="1">25±6</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % GAP</td><td align="left" rowspan="1" colspan="1">6±4</td><td align="left" rowspan="1" colspan="1">7±5</td><td align="left" rowspan="1" colspan="1">11±5</td><td align="left" rowspan="1" colspan="1">24±12</td><td align="left" rowspan="1" colspan="1">26±20</td></tr><tr><td align="left" rowspan="1" colspan="1">Total % gaps</td><td align="left" rowspan="1" colspan="1">17.6</td><td align="left" rowspan="1" colspan="1">39.0</td><td align="left" rowspan="1" colspan="1">53.2</td><td align="left" rowspan="1" colspan="1">74.8</td><td align="left" rowspan="1" colspan="1">78.0</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">LSU</td><td align="left" rowspan="1" colspan="1">Archaea</td><td align="left" rowspan="1" colspan="1">Chloroplasts</td><td align="left" rowspan="1" colspan="1">Bacteria</td><td align="left" rowspan="1" colspan="1">Eukarya</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">No. seqs</td><td align="left" rowspan="1" colspan="1">26</td><td align="left" rowspan="1" colspan="1">32</td><td align="left" rowspan="1" colspan="1">120</td><td align="left" rowspan="1" colspan="1">89</td></tr><tr><td align="left" rowspan="1" colspan="1">Alignment length</td><td align="left" rowspan="1" colspan="1">3346</td><td align="left" rowspan="1" colspan="1">4311</td><td align="left" rowspan="1" colspan="1">4555</td><td align="left" rowspan="1" colspan="1">9055</td></tr><tr><td align="left" rowspan="1" colspan="1">Geometric mean seqs</td><td align="left" rowspan="1" colspan="1">2994</td><td align="left" rowspan="1" colspan="1">2936</td><td align="left" rowspan="1" colspan="1">2918</td><td align="left" rowspan="1" colspan="1">3623</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % ID</td><td align="left" rowspan="1" colspan="1">65±9</td><td align="left" rowspan="1" colspan="1">70±10</td><td align="left" rowspan="1" colspan="1">67±6</td><td align="left" rowspan="1" colspan="1">50±12</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % SUBS</td><td align="left" rowspan="1" colspan="1">29±7</td><td align="left" rowspan="1" colspan="1">20±5</td><td align="left" rowspan="1" colspan="1">25±4</td><td align="left" rowspan="1" colspan="1">22±5</td></tr><tr><td align="left" rowspan="1" colspan="1">Pairwise % GAP</td><td align="left" rowspan="1" colspan="1">6±2</td><td align="left" rowspan="1" colspan="1">11±5</td><td align="left" rowspan="1" colspan="1">9±4</td><td align="left" rowspan="1" colspan="1">28±11</td></tr><tr><td align="left" rowspan="1" colspan="1">Total % gaps</td><td align="left" rowspan="1" colspan="1">10.5</td><td align="left" rowspan="1" colspan="1">31.8</td><td align="left" rowspan="1" colspan="1">35.9</td><td align="left" rowspan="1" colspan="1">59.6</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula><label>(1)</label></disp-formula> <disp-formula><label>(2)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(3)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(4)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(5)</label></disp-formula> <disp-formula><label>(6)</label></disp-formula> <disp-formula><label>(7)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(8)</label></disp-formula> <disp-formula><label>(9)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(10)</label></disp-formula> <disp-formula><label>(11)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(12)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(13)</label></disp-formula> <disp-formula><label>(14)</label></disp-formula> <disp-formula><label>(15)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(16)</label></disp-formula> <disp-formula><label>(17)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(18)</label></disp-formula> <disp-formula><label>(19)</label></disp-formula> <disp-formula><label>(20)</label></disp-formula> <disp-formula><label>(21)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(22)</label></disp-formula> <disp-formula><label>(23)</label></disp-formula> <disp-formula><label>(24)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(25)</label></disp-formula> <disp-formula><label>(26)</label></disp-formula> <disp-formula><label>(27)</label></disp-formula> <inline-formula></inline-formula> <disp-formula><label>(28)</label></disp-formula> <disp-formula><label>(29)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(30)</label></disp-formula> <disp-formula><label>(31)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<supplementary-material content-type="local-data" id="pcbi.1000172.s001"><label>Dataset S1</label><caption><p>Supplemental Material</p><p>(24.89 MB GZ)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000172.s002"><label>Text S1</label><caption><p>Appendix 1</p><p>(0.16 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000172.s003"><label>Text S2</label><caption><p>Appendix 2</p><p>(0.11 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000172.s004"><label>Text S3</label><caption><p>Appendix 3</p><p>(0.14 MB PDF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000172.s005"><label>Text S4</label><caption><p>Appendix 4</p><p>(0.13 MB PDF)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><p>For a large variety of synthetic alignments of 8 taxa with gaps generated using a Poisson (<italic>λ</italic> = 0.5) length distribution. We compare the performance of <sc>dnaml</sc>
<italic>ε</italic> respect to <sc>dnaml</sc> using three different measures. For each measure, we report the relative area under the curve (<sc>auc</sc>) difference , which for a given measure <italic>f</italic> describes (as a percentage) the difference between the <sc>auc</sc> for <sc>dnaml</sc>
<italic>ε</italic> versus that of <sc>dnaml</sc> relative to the larger of the two <sc>auc</sc>'s. In white, we show the results obtained with the phylogenetically correct alignments, in gray after realigning with <sc>clustalw</sc> 1.83. In bold when the result is favorable to <sc>dnaml</sc>
<italic>ε</italic>.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><p>Similarly to ##TAB##0##Table 1##, we show results for a large variety of synthetic alignments of 8 taxa with gaps generated using the generative model described in this paper. We compare the performance of <sc>dnaml</sc>
<italic>ε</italic>
<sc>respect</sc>
<sc>to</sc>
<sc>dnaml</sc> using three different measures. As in ##TAB##0##Table 1## and for each measure, we report the relative area under the curve (<sc>auc</sc>) difference  defined above. White rows indicate results for the phylogenetically correct alignments. Gray rows indicate results after realigning the evolved sequences with <sc>clustalw</sc> 1.83. In bold when the result is favorable to <sc>dnaml</sc>
<italic>ε</italic>.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><p>We assess the performance of <sc>dnaml</sc>
<italic>ε</italic> for five different gap length distributions using alignments of 8 taxa with similar average pairwise substitutions and gaps. The gap distributions used to create the synthetic alignments are ordered by increasing degree of violation of the column independence assumption. For a given measure, we report the relative <sc>auc</sc> difference  introduced in ##TAB##0##Table 1##. In bold when the result is favorable to <sc>dnaml</sc>
<italic>ε</italic>.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt104"><p>We test the possibility of spurious long branch attraction associated to <sc>dnaml</sc>
<italic>ε</italic> using 4-taxon trees depending on two parameters: <italic>t</italic>
<sub>3</sub> the length of the internal branch and two opposite leaves, and <italic>t</italic>
<sub>2</sub>, the length of the other two opposite leaves ##UREF##24##[81]##. We compare the performance of <sc>dnaml</sc>
<italic>ε</italic> respect to that of <sc>dnaml</sc> using alignments generated by trees with similar average substitutions branch length (abl) but such that in one case the four leaves are identical (<italic>t</italic>
<sub>3</sub> = <italic>t</italic>
<sub>2</sub> or <italic>r</italic> = 1) with another extreme case prone to long branch attraction in which two opposite leaves are 10 times longer than the other two leaves and the internal branch (<italic>t</italic>
<sub>2</sub> = 10-<italic>t</italic>
<sub>3</sub> or <italic>r</italic> = 10). For a given “abl” and a given parameter <italic>r</italic>, one has  and <italic>t</italic>
<sub>2</sub> = <italic>rt</italic>
<sub>3</sub>. For each method (<sc>dnaml</sc> and <sc>dnaml</sc>
<italic>ε</italic>) and tree configuration (<italic>r</italic> = 1 and <italic>r</italic> = 10), we report the mean percentage of true positives (TP). No long branch attraction is observed for the phylogenetically correct alignment (in white). For the <sc>clustalw</sc> 1.83 alignments (in gray), both methods have a tendency to infer the wrong tree as the divergence increases. In bold when more than 50% of the trees are incorrectly predicted.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt105"><p>Statistics of the rRNA alignments obtained from the Comparative RNA Web Site ##UREF##25##[82]## after sequences with more than 95% identity to each other have been removed from the alignments.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was funded by the Howard Hughes Medical Institute.</p></fn></fn-group>
<graphic xlink:href="pcbi.1000172.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e002.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e003.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e004.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e005.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e006.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e007.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e009.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e010.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e011.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e012.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e013.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e014.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e015.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e016.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e017.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e018.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e019.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e020.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e021.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e022.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e023.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e024.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e025.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e026.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e027.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e028.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e029.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e081.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e082.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e030.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e031.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e032.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e033.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e034.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e035.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e036.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e037.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e038.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e039.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e040.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e041.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e042.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e043.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.g001"/> <graphic xlink:href="pcbi.1000172.e044.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e045.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e046.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e047.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e048.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e049.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e050.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e051.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.g002"/> <inline-graphic xlink:href="pcbi.1000172.e052.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.g003"/> <inline-graphic xlink:href="pcbi.1000172.e053.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e054.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e055.jpg" mimetype="image"/> <graphic id="pcbi-1000172-t001-1" xlink:href="pcbi.1000172.t001"/> <inline-graphic xlink:href="pcbi.1000172.e056.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e057.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e058.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e059.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e060.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e061.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e062.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e063.jpg" mimetype="image"/> <graphic id="pcbi-1000172-t002-2" xlink:href="pcbi.1000172.t002"/> <inline-graphic xlink:href="pcbi.1000172.e064.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e065.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e066.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000172.e067"/> <inline-graphic xlink:href="pcbi.1000172.e068.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e069.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e070.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e071.jpg" mimetype="image"/> <graphic id="pcbi-1000172-t003-3" xlink:href="pcbi.1000172.t003"/> <inline-graphic xlink:href="pcbi.1000172.e072.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e073.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e074.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e075.jpg" mimetype="image"/> <graphic id="pcbi-1000172-t004-4" xlink:href="pcbi.1000172.t004"/> <inline-graphic xlink:href="pcbi.1000172.e076.jpg" mimetype="image"/> <graphic id="pcbi-1000172-t005-5" xlink:href="pcbi.1000172.t005"/> <graphic xlink:href="pcbi.1000172.g004"/> <graphic xlink:href="pcbi.1000172.e083.jpg" mimetype="image" position="float"/> <graphic xlink:href="pcbi.1000172.e084.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000172.e077.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e078.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e079.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000172.e080.jpg" mimetype="image"/>
<media xlink:href="pcbi.1000172.s001.gz"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000172.s002.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000172.s003.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000172.s004.pdf"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000172.s005.pdf"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["2"], "element-citation": ["\n"], "surname": ["Mau"], "given-names": ["B"], "year": ["1996"], "article-title": ["Bayesian phylogenetic inference via Markov chain Monte Carlo methods."], "comment": ["Ph.D. Dissertation, University of Wisconsin, Madison"]}, {"label": ["5"], "element-citation": ["\n"], "surname": ["Larget", "Simon"], "given-names": ["B", "D"], "year": ["1999"], "article-title": ["Markov chain Monte Carlo algorithms for the Bayesian analysis of phylogenetic trees."], "source": ["Mol Biol Evol"], "volume": ["16"], "fpage": ["750"], "lpage": ["759"]}, {"label": ["6"], "element-citation": ["\n"], "surname": ["Durbin", "Eddy", "Krogh", "Mitchison"], "given-names": ["R", "SR", "A", "GJ"], "year": ["1998"], "source": ["Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids"], "publisher-loc": ["Cambridge, UK"], "publisher-name": ["Cambridge University Press"]}, {"label": ["7"], "element-citation": ["\n"], "surname": ["Jukes", "Cantor"], "given-names": ["TH", "C"], "year": ["1965"], "article-title": ["Evolution of protein molecules."], "source": ["Mammalian Protein Metabolism"], "publisher-loc": ["New York"], "publisher-name": ["Academic Press"], "fpage": ["21"], "lpage": ["132"]}, {"label": ["8"], "element-citation": ["\n"], "surname": ["Dayhoff", "Schwartz", "Orcutt"], "given-names": ["M", "R", "B"], "year": ["1978"], "article-title": ["A model of evolutionary change in protein."], "source": ["Atlas of Protein Sequence Structure"], "volume": ["5"], "fpage": ["345"], "lpage": ["352"]}, {"label": ["10"], "element-citation": ["\n"], "surname": ["Hasegawa", "Kishino", "Yano"], "given-names": ["M", "H", "T"], "year": ["1985"], "article-title": ["Dating of the human-ape splitting by a molecular clock of mitochondrial DNA."], "source": ["J Mol Evol"], "volume": ["21"], "fpage": ["160"], "lpage": ["174"]}, {"label": ["11"], "element-citation": ["\n"], "surname": ["Tavar\u00e9"], "given-names": ["S"], "year": ["1986"], "article-title": ["Some probabilistic and statistical problems in the analysis of DNA sequences."], "source": ["Lect Math Life Sci"], "volume": ["17"], "fpage": ["57"], "lpage": ["86"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Felsenstein"], "given-names": ["J"], "year": ["2006"], "article-title": ["PHYLIP (Phylogeny Inference Package), version 3.66."], "comment": ["Distributed by the author. Department of Genome Sciences, University of Washington, Seattle"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Swofford"], "given-names": ["DL"], "year": ["2003"], "source": ["PAUP*. Phylogenetic analysis using parsimony (*and other methods,. version 4"], "publisher-loc": ["Sunderland, Massachusetts"], "publisher-name": ["Sinauer Associates"]}, {"label": ["26"], "element-citation": ["\n"], "surname": ["Adachi", "Hasegawa"], "given-names": ["J", "M"], "year": ["1995"], "source": ["MOLPHY programs for molecular phylogenetics, version 2.3"], "publisher-loc": ["Tokyo"], "publisher-name": ["Institute of Statistical Mathematics"]}, {"label": ["29"], "element-citation": ["\n"], "surname": ["Simon", "Larget"], "given-names": ["D", "B"], "year": ["2000"], "article-title": ["Bayesian analysis in molecular biology and evolution (BAMBE), version 2.03 beta."], "comment": ["Department of Mathematics and Computer Science, Duquesne University"]}, {"label": ["30"], "element-citation": ["\n"], "surname": ["Ronquist", "Huelsenbeck"], "given-names": ["F", "JP"], "year": ["2003"], "article-title": ["MRBAYES: Bayesian inference of phylogenetic trees."], "source": ["Bioinformatics"], "volume": ["17"], "fpage": ["754"], "lpage": ["7555"]}, {"label": ["36"], "element-citation": ["\n"], "surname": ["Coin", "Durbin"], "given-names": ["L", "R"], "year": ["2004"], "article-title": ["Improved techniques for the identification of pseudogenes."], "source": ["Bioinformatics Suppl"], "volume": ["1"], "fpage": ["I94"], "lpage": ["I100"]}, {"label": ["39"], "element-citation": ["\n"], "surname": ["Mitchison", "Durbin"], "given-names": ["GJ", "RM"], "year": ["1995"], "article-title": ["Tree-based maximal likelihood substitutions matrices and hidden Markov models."], "source": ["J Mol Evol"], "volume": ["41"], "fpage": ["1139"], "lpage": ["11351"]}, {"label": ["54"], "element-citation": ["\n"], "surname": ["Mikl\u00f3s", "Toroczkai", "Gascuel", "Moret"], "given-names": ["I", "Z", "O", "BME"], "year": ["2001"], "article-title": ["An improved model for statistical aligment."], "source": ["WABI 2001"], "publisher-loc": ["Berlin Heidelberg"], "publisher-name": ["Springer-Verlag"], "fpage": ["1"], "lpage": ["10"]}, {"label": ["59"], "element-citation": ["\n"], "surname": ["Holmes"], "given-names": ["I"], "year": ["2003"], "article-title": ["Using guide trees to construct multiple-sequence evolutionary HMMs."], "source": ["Bioinformatics Suppl"], "volume": ["1"], "fpage": ["147"], "lpage": ["157"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Steel", "Hein"], "given-names": ["M", "J"], "year": ["2001"], "article-title": ["Applying the Thorne-Kishino-Felsenstein model to sequence evolution on a star-shaped tree."], "source": ["Appl Math Lett"], "volume": ["14"], "fpage": ["679"], "lpage": ["684"]}, {"label": ["65"], "element-citation": ["\n"], "surname": ["Hein"], "given-names": ["J"], "year": ["2001"], "article-title": ["An algorithm ofr statistical alignment of sequences related by a binary tree."], "source": ["Pac Symp Biocomput"], "volume": ["6"], "fpage": ["179"], "lpage": ["190"]}, {"label": ["67"], "element-citation": ["\n"], "surname": ["Lunter", "Mikl\u00f3s", "Song", "Hein"], "given-names": ["G", "I", "YS", "J"], "year": ["2003"], "article-title": ["An efficient algorithm for statistical multiple alignment on arbitrary phylogenetic trees."], "source": ["J Mol Biol"], "volume": ["10"], "fpage": ["869"], "lpage": ["889"]}, {"label": ["69"], "element-citation": ["\n"], "surname": ["Lunter", "Mikl\u00f3s", "Drummond", "Jensen", "Hein"], "given-names": ["G", "I", "A", "J", "J"], "year": ["2003"], "article-title": ["Bayesian phylogenetic inference under a statistical insertion-deletion model. Proceedings of WABI'03."], "source": ["Lect Notes Bioinformatics"], "volume": ["2812"], "fpage": ["228"], "lpage": ["244"]}, {"label": ["70"], "element-citation": ["\n"], "surname": ["Felsenstein"], "given-names": ["J"], "year": ["2004"], "source": ["Inferring Phylogenies"], "publisher-loc": ["Sunderland, Massachusetts"], "publisher-name": ["Sinauer"], "fpage": ["508"]}, {"label": ["72"], "element-citation": ["\n"], "surname": ["Moler", "Loan"], "given-names": ["C", "CV"], "year": ["2003"], "article-title": ["Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later."], "source": ["SIAM Rev"], "volume": ["45"], "fpage": ["3"], "lpage": ["49"]}, {"label": ["73"], "element-citation": ["\n"], "surname": ["Yang"], "given-names": ["Z"], "year": ["2006"], "source": ["Computational molecular evolution"], "publisher-loc": ["London"], "publisher-name": ["Oxford University Press"], "fpage": ["33"]}, {"label": ["77"], "element-citation": ["\n"], "surname": ["Robinson", "Foulds"], "given-names": ["DF", "LR"], "year": ["1981"], "article-title": ["Comparison of phylogenetic trees."], "source": ["Math Biosci"], "volume": ["53"], "fpage": ["131"], "lpage": ["147"]}, {"label": ["81"], "element-citation": ["\n"], "surname": ["Huelsenbeck"], "given-names": ["JP"], "year": ["1995"], "article-title": ["The performance of phylogenetic methods in simulation."], "source": ["Syst Biol"], "volume": ["44"], "fpage": ["17"], "lpage": ["48"]}, {"label": ["82"], "element-citation": ["\n"], "surname": ["Cannone", "Subramanian", "Schnare", "Collett", "D'Souza"], "given-names": ["JJ", "S", "MN", "JR", "LM"], "year": ["2002"], "article-title": ["The Comparative RNA Web (CRW) Site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs. BMC Bioinformatics 3: 2."], "source": ["Correction: BMC Bioinformatics"], "volume": ["3"], "fpage": ["15"]}, {"label": ["85"], "element-citation": ["\n"], "surname": ["Nielsen"], "given-names": ["R"], "year": ["2005"], "source": ["Statistical Methods in Molecular Evolution"], "publisher-loc": ["New York"], "publisher-name": ["Springer-Verlag"], "fpage": ["378"]}]
PLoS Comput Biol. 2008 Sep 19; 4(9):e1000172