Skip to content

Instantly share code, notes, and snippets.

@spond
Last active August 13, 2018 10:55
Show Gist options
  • Save spond/c827d8c2f68cf1daa762 to your computer and use it in GitHub Desktop.
Save spond/c827d8c2f68cf1daa762 to your computer and use it in GitHub Desktop.
Root-to-tip for intra-host HIV-1 data

Root-to-tip regression for intrahost HIV-1 phylogenies.

In this example, we read in trees with branch lengths (FastTree2), based on haplotypes and counts reported in Lorenzo-Redondo et al, construct root-to-tip distances for time-stamped tips, and regress those on time to obtain crude evolutionary rate estimates.

The initial rooting of the tree is obtained by considering all time-point 0 sequences and selecting one that maximizes R2 (regardless of the slope). The user can perform arbitrary rerooting and see what the effect on the regression slope is.

Further, 1000 simulated perturbation of haplotype counts are considered for a given rooting: the estimated frequency of each haplotype (encoded in the tip name) is multiplied by a random number from [0.05-10] distributed uniformly. The proportion of perturbations with positive inferred regression slopes is shown on the chart.

<!DOCTYPE html>
<html lang='en'>
<head>
<meta charset="utf-8">
<meta http-equiv="Content-Security-Policy" content="upgrade-insecure-requests">
<link href="//veg.github.io/phylotree.js/phylotree.css" rel="stylesheet">
<link rel="stylesheet" href="//netdna.bootstrapcdn.com/bootstrap/3.3.5/css/bootstrap.min.css">
<link rel="stylesheet" href="//netdna.bootstrapcdn.com/bootstrap/3.3.5/css/bootstrap-theme.min.css">
<script src="//code.jquery.com/jquery.js"></script>
<script src="//netdna.bootstrapcdn.com/bootstrap/3.3.5/js/bootstrap.min.js"></script>
<script src="//d3js.org/d3.v3.min.js"></script>
<script src="//phylotree.hyphy.org/phylotree.js"></script>
<script src="//cdnjs.cloudflare.com/ajax/libs/underscore.js/1.8.3/underscore-min.js"></script>
<style>
.date-text {
padding: 0.25em;
text-align: center;
}
.dot {
stroke: #000;
opacity : 0.1;
}
.size-bubble, .size-bubble * {
fill: #CCC;
shape-rendering: crispEdges;
stroke: black;
stroke-width: 1px;
}
.size-label {
display: block;
text-align : center;
font-family: sans-serif;
font-size: 9pt;
}
.fit-label {
text-anchor : start;
font-family: sans-serif;
font-size: 7pt;
}
.axis text {
font: 14px sans-serif;
}
.axis path, .axis line {
fill: none;
stroke: #000;
shape-rendering: crispEdges;
}
.fit {
fill: none;
stroke-width: 2px;
stroke: red;
}
.fit-no-counts {
opacity: 0.1;
stroke: gray;
}
@media print
{
.no-print, .no-print *
{
display: none !important;
}
}
</style>
</head>
<body>
<nav class="navbar navbar-default" role="navigation">
<!-- Brand and toggle get grouped for better mobile display -->
<div class="navbar-header">
<button type="button" class="navbar-toggle" data-toggle="collapse" data-target=".navbar-ex1-collapse">
<span class="sr-only">Toggle navigation</span>
<span class="icon-bar"></span>
<span class="icon-bar"></span>
<span class="icon-bar"></span>
</button>
<a id = "menu-title" class="navbar-brand" href="#"></a>
</div>
<div class="collapse navbar-collapse">
<ul class="nav navbar-nav">
<li class="dropdown">
<a href="#" class="dropdown-toggle" data-toggle="dropdown">Tree to display <b class="caret"></b></a>
<ul class="dropdown-menu" id = 'tree_selector'>
<li><a href="#" data-tree = "js1774-gag.tre" data-tag = "JS1774 gag"></a></li>
<li><a href="#" data-tree = "sm1727-gag.tre" data-tag = "SM1727 gag"></a></li>
<li><a href="#" data-tree = "tm1679-gag.tre" data-tag = "TM1679 gag"></a></li>
<li><a href="#" data-tree = "js1774-pol-1.tre" data-tag = "JS1774 pol (PR)"></a></li>
<li><a href="#" data-tree = "js1774-pol-2.tre" data-tag = "JS1774 pol (RT)"></a></li>
<li><a href="#" data-tree = "sm1727-pol-1.tre" data-tag = "SM1727 pol (PR)"></a></li>
<li><a href="#" data-tree = "sm1727-pol-2.tre" data-tag = "SM1727 pol (RT)"></a></li>
<li><a href="#" data-tree = "tm1679-pol-1.tre" data-tag = "TM1679 pol (PR)"></a></li>
<li><a href="#" data-tree = "tm1679-pol-2.tre" data-tag = "TM1679 pol (RT)"></a></li>
</ul>
</li>
<li>
<a href="#" data-toggle="modal" data-target="#newick_modal">Input Newick</a>
</li>
</ul>
<form class="navbar-form navbar-left">
<div class="checkbox navbar-btn">
<label class="navbar-link">
<input type="checkbox" id="layout">
Radial layout
</label>
</div>
</form>
</div><!-- /.navbar-collapse -->
</nav>
<div class = "container">
<div class = "row" style = "margin-top: 1em">
<div class = "col-md-2" id="date-legend">
</div>
<div class = "col-md-2" id="compartment-legend">
</div>
<div class = "col-md-2" id="size-legend">
</div>
<div class = "col-md-6">
<svg id="rtt_display"></svg>
</div>
</div>
<div class = "row">
<div class = "col-md-12">
<svg id="tree_display"></svg>
</div>
</div>
</div>
<div class="modal" id = 'newick_modal'>
<div class="modal-dialog">
<div class="modal-content">
<div class="modal-header">
<button type="button" class="close" data-dismiss="modal" aria-hidden="true">&times;</button>
<h4 class="modal-title">Newick string to render</h4>
</div>
<div class="modal-body" id = 'newick_body'>
<textarea id = 'nwk_spec' autofocus = true placeholder = "" style = 'width: 100%; height: 100%' rows = 20 selectionStart = 1 selectionEnd = 1000>(a : 0.1, (b : 0.11, (c : 0.12, d : 0.13) : 0.14) : 0.15)</textarea>
</div>
<div class="modal-footer">
<button type="button" class="btn btn-primary" id = 'validate_newick'>Display this tree</button>
</div>
</div><!-- /.modal-content -->
</div><!-- /.modal-dialog -->
</div><!-- /.modal -->
<script>
// use these formats for parsing and displaying sampling dates
var date_in = d3.time.format("%Y%m%d"),
date_out = d3.time.format("%B %d, %Y");
// default scheme to color by date
var coloring_scheme = d3.scale.category10();
// different shapes for labeling compartments
var compartment_labels = d3.scale.ordinal().range(["circle", "square", "diamond", "cross", "triangle-up"]);
// global tree object
var tree;
// bubble size scale
var size_scale = d3.scale.pow ().exponent (0.4).range ([1,6]).clamp (false).domain ([0,1]);
var rescaled_size;
// determines the size of a node "bubble"
function bubbleSize(node) {
if (node && node.copy_number) {
return size_scale(node.copy_number);
}
return 1;
}
function nodeStyler(container, node) {
if (d3.layout.phylotree.is_leafnode(node)) {
var existing_circle = container.selectAll("circle");
if (existing_circle.size() == 1) {
existing_circle.remove();
}
if (node.copy_number) {
existing_circle = container.selectAll("path.node_shape").data([node.compartment]);
existing_circle.enter().append("path").classed("node_shape", true);
var bubble_size = tree.node_bubble_size(node);
var label = existing_circle.attr("d", function(d) {
return d3.svg.symbol().type(compartment_labels(d)).size(bubble_size * bubble_size)();
}).selectAll ("title").data ([node.copy_number]);
label.enter().append("title");
label.text ("" + node.copy_number + " copies");
existing_circle.style("stroke-width", "1px").style("stroke", "black");
if (node.text_angle) {
existing_circle.attr("transform", function(d) {
return "rotate(" + node.text_angle + ") translate(" + (node.text_align == "end" ? -1 : 1) * bubble_size / 2 + ",0)";
});
} else {
existing_circle.attr("transform", function(d) {
return "translate(" + bubble_size / 2 + ",0)";
});
}
}
}
if (node.date) {
var node_color = coloring_scheme (node.date);
container.selectAll("circle").style("fill", node_color);
container.selectAll("path").style("fill", node_color);
container.style("fill", node_color);
}
}
function drawATree(newick) {
tree = d3.layout.phylotree()
.svg(d3.select("#tree_display"))
.options({
'selectable': false,
// make nodes and branches not selectable
'collapsible': true,
// turn off the menu on internal nodes
'transitions': false,
// turn off d3 animations.
'draw-size-bubbles': true,
// draw node size bubbles
'annular-limit': 0.5,
// controls how much of the overall diameter can be taken by
// empty "annular" space
'max-radius': 512
// maximum radial layout radius
})
.node_span(bubbleSize)
.style_nodes(nodeStyler)
.node_circle_size (2) // do not draw clickable circles for internal nodes
.branch_name (function () {return ""}) // no leaf names
.layout_handler (function (tree) {
plotRTT (tree, "rtt_display", 350,350)
}
)
;
/* the next call creates the tree object, and tree nodes */
tree(d3.layout.newick_parser(newick));
/* parse tip names which are assumed to be like this:
6_PB_DNA_1_37-52
splitting on '_' yields
-- TimePoint; parsed and used to assign *colors* to edges/nodes and to draw root-to-tip distances
-- compartment; parsed and used to assign *shapes* to nodes in combination with the DNA/RNA source
-- source;
-- clone ID; ignored
-- another ID; ignored
-- copy number; parsed and used to scale "bubbles" on leaves
tips which do not conform to this naming convention are labeled as "reference",
and annotated using a special character.
*/
var unique_compartments = {},
unique_dates = {};
var copy_number_range = [1e100,0],
attributed_node = null;
_.each(tree.get_nodes(), function(value, key) {
var attributes = value.name.split('_');
if (attributes.length == 5) {
attributed_node = value;
value.compartment = attributes[1] + ' ' + attributes[2];
value.copy_number = parseFloat(attributes[4].split ('-')[1]);
copy_number_range[0] = Math.min (copy_number_range[0], value.copy_number);
copy_number_range[1] = Math.max (copy_number_range[1], value.copy_number);
value.numeric_date = parseInt (attributes[0]);
value.date = "Month " + attributes[0];
unique_compartments[value.compartment] = 1;
unique_dates[value.date] = 1;
} else {
value.is_reference = true;
}
});
size_scale.domain (copy_number_range);
coloring_scheme.domain (_.keys(unique_dates).sort());
compartment_labels.domain (_.keys(unique_compartments).sort());
if ($("#layout").prop("checked")) {
tree.radial (true);
}
tree.placenodes().layout();
var rescale = tree.node_bubble_size(attributed_node) / size_scale (attributed_node.copy_number) * 0.5;
rescaled_size = _.compose (function (d) {return d*rescale}, size_scale);
// create a legend for sizes
// power of 10 ticks
var size_bubbles = _.map (d3.range (Math.floor (Math.log(copy_number_range[0])/Math.log (10)), Math.ceil (Math.log(copy_number_range[1])/Math.log (10))),
function (v) {return Math.pow (10,v);});
d3.select ("#size-legend").selectAll (".row").remove();
var reference_sizes = d3.select ("#size-legend").selectAll (".row").data (_.map (size_bubbles, function (d) {return [d];}));
reference_sizes.enter().append ("div").classed ("row", true).selectAll (".col-md-2").data (function (d) {return [d];}).enter().append ("div").classed ("col-md-2", true);
var max_size = rescaled_size (size_bubbles[size_bubbles.length - 1]);
reference_sizes.selectAll (".col-md-2").each (function (d) {
var circle_size = rescaled_size (d);
d3.select (this).append ("svg").attr ("width", 2*max_size + 4).attr ("height", 2*circle_size + 4).append ("circle").classed("size-bubble", true).attr ("r", circle_size).attr ("cx", max_size+2).attr ("cy", circle_size+2);
d3.select (this).append ("span").classed ("size-label", true).style ("width", "" + (2*max_size+4) + "px").text (d);
});
// create a legend for compartments
d3.select ("#compartment-legend").selectAll (".row").remove();
var compartment_depictions = d3.select ("#compartment-legend").selectAll (".row").data (compartment_labels.domain());
compartment_depictions.enter().append ("div").classed ("row", true).selectAll (".col-md-2").data (function (d) {return [d];}).enter().append ("div").classed ("col-md-2", true);
compartment_depictions.exit().remove();
compartment_depictions.selectAll (".col-md-2").each (function (d) {
d3.select (this).append ("svg").attr ("width", 2*max_size + 4).attr ("height", 2*max_size + 4).append ("path").classed("size-bubble", true).attr("d", function(d) {
return d3.svg.symbol().type(compartment_labels(d)).size(max_size * max_size)();
}).attr ("transform", "translate( " + max_size + "," + max_size + ")");
d3.select (this).append ("span").classed ("size-label", true).style ("width", "" + (2*max_size+4) + "px").text (d);
});
// create a legend for dates
d3.select ("#date-legend").selectAll (".row").remove();
var date_colors = d3.select ("#date-legend").selectAll (".row").data (coloring_scheme.domain());
date_colors.enter().append ("div").classed ("row", true);
date_colors.exit().remove();
date_colors.each (function (d) {
d3.select(this).style ("background-color", coloring_scheme (d)).classed ("date-text", true).text (d).style ("color", d3.rgb(coloring_scheme (d)).hsl().l < 0.5 ? "white" : "black");
});
// UI handlers
$("#layout").on("click", function(e) {
tree.radial($(this).prop("checked")).placenodes().update();
});
var best_node = plotRTT (tree, "rtt_display", 350,350, true);
tree.reroot (best_node).update();
}
function computeRootToTip (node) {
if (node.parent) {
node.rtt = node.parent.rtt + tree.branch_length () (node);
}
else {
node.rtt = 0;
}
if (node.children) {
for (var c = 0; c < node.children.length; c++) {
computeRootToTip (node.children[c]);
}
}
}
function computeRootToTipOtherRoot (node, coming_from, shared_distance, distance_to_new_root) {
var my_bl = tree.branch_length () (node);
var go_up = false;
if (!coming_from) {
shared_distance = node.rtt;
distance_to_new_root = 0;
go_up = true;
}
if (node.children) {
for (var c = 0; c < node.children.length; c++) {
if (node.children[c] != coming_from) {
computeRootToTipOtherRoot (node.children[c], node, shared_distance, distance_to_new_root);
} else {
go_up = true;
}
}
}
node.rtta = node.rtt - shared_distance + distance_to_new_root;
if (go_up) {
shared_distance -= my_bl;
distance_to_new_root += my_bl;
}
if (node.parent && go_up) {
computeRootToTipOtherRoot (node.parent, node, shared_distance, distance_to_new_root);
}
}
function loadTreeFromURL(url) {
d3.text(url, function(error, newick) {
drawATree(newick);
});
}
function plotRTT (tree, tag, w, h, auto_root) {
computeRootToTip (tree.get_nodes()[0]);
if (auto_root) {
var max_r2 = 0;
var best_node = 0;
tree.traverse_and_compute (function (node) {
if (d3.layout.phylotree.is_leafnode (node)) {
if (node.numeric_date == 0) {
computeRootToTipOtherRoot (node, null, 0, 0);
linear_data = [];
tree.traverse_and_compute (function (node) {
if (d3.layout.phylotree.is_leafnode (node)) {
linear_data.push ([node.numeric_date, node.rtta, node.copy_number]);
}
});
var fit_with_counts = linear_fit (linear_data)["r2"];
if (fit_with_counts > max_r2) {
max_r2 = fit_with_counts;
best_node = node;
}
}
}
});
return best_node;
}
var margin = {top: 20, right: 80, bottom: 45, left: 65},
width = w - margin.left - margin.right,
height = h - margin.top - margin.bottom;
var x = d3.scale.linear()
.range([20, width-20]);
var y = d3.scale.linear()
.range([height, 0]);
var xAxis = d3.svg.axis()
.scale(x)
.orient("bottom");
var yAxis = d3.svg.axis()
.scale(y)
.orient("left");
var plot_here = d3.select ("#" + tag);
plot_here.attr ("width", w).attr("height", h);
plot_here.selectAll ("g").remove();
var svg = plot_here.append ("g")
.attr("transform", "translate(" + margin.left + "," + margin.top + ")");
var plot_data = [];
var linear_data = [];
tree.traverse_and_compute (function (node) {
if (d3.layout.phylotree.is_leafnode (node)) {
plot_data.push (node);
linear_data.push ([node.numeric_date, node.rtt, 1]);
}
});
xAxis.tickValues (_.keys(_.countBy (linear_data, function (d) {return d[0];})).sort());
x.domain(d3.extent(plot_data, function(d) { return d.numeric_date; })).nice();
y.domain(d3.extent(plot_data, function(d) { return d.rtt; })).nice();
svg.append("g")
.attr("class", "x axis")
.attr("transform", "translate(0," + height + ")")
.call(xAxis)
.append("text")
.attr("class", "label")
.attr("x", width)
.attr("y", 35)
.style("text-anchor", "end")
.text("Time");
svg.append("g")
.attr("class", "y axis")
.call(yAxis)
.append("text")
.attr("class", "label")
.attr("transform", "rotate(-90)")
.attr("y", -60)
.attr("dy", ".71em")
.style("text-anchor", "end")
.text("Root-to-tip") ;
svg.selectAll(".dot")
.data(plot_data)
.enter().append("circle")
.attr("class", "dot")
.attr("r", function (d) {return 2*rescaled_size(tree.node_bubble_size(d));})
.attr("cx", function(d) { return x(d.numeric_date); })
.attr("cy", function(d) {
return y(d.rtt); })
.style("fill", function(d,i) { return "#000000"; });
function linear_fit (data) {
var ss = linear_data.reduce (function (p,c) {return c[2] + p}, 0), // sample count
sx = linear_data.reduce (function (p,c) {return c[2]*c[0] + p}, 0), // sum X
sy = linear_data.reduce (function (p,c) {return c[2]*c[1] + p}, 0) // sum Y
sxoss = sx/ss,
syoss = sy/ss;
var fitB = 0,
st2 = 0,
vary = 0;
data.forEach (function (point) {
var t = point[0] - sxoss;
st2 += point[2] * t * t;
fitB += point[2] * t * point[1];
vary += point[2] * (point[1]-syoss)*(point[1]-syoss);
});
fitB /= st2;
a = (sy - sx * fitB) / ss;
var varres = 0;
data.forEach (function (point) {
var t = point[1] - a - fitB * point[0];
varres += point[2] * t * t;
});
var normvarres = Math.sqrt (varres / (ss - 2));
return {'intercept' : a,
'slope' : fitB,
'r2': 1 - varres / vary,
'var (intercept)' : Math.sqrt ((1+sx*sx/(ss*st2))/ss),
'var (slope)' : Math.sqrt (1/st2)
};
};
var fit_with_counts = linear_fit (linear_data);
//console.log (fit_with_counts["R"]);
var label_format = d3.format (".2r");
function format_line (i, s) {
return label_format (i) + (s >= 0 ? '+' : '') + label_format (s) + '*T';
}
var positive_slope = 0;
for (s = 0; s < 1000; s += 1) {
fit_without_counts = (linear_fit (linear_data.map (function (d) {return [d[0], d[1], Math.floor (d[2]*(0.05 + 9.95*Math.random (0,1)))];})));
if (fit_without_counts ["slope"] > 0) {
positive_slope ++;
}
svg.append ("line")
.attr ("x1", x (x.domain()[0]))
.attr ("y1", y (x.domain()[0]*fit_without_counts['slope'] + fit_without_counts ['intercept']))
.attr ("x2", x (x.domain()[1]))
.attr ("y2", y (x.domain()[1]*fit_without_counts['slope'] + fit_without_counts ['intercept']))
.attr ("class", "fit fit-no-counts");
}
svg.append ("text").attr ("x", x (x.domain()[0])).attr ( y (y.domain [1])).text (d3.format(".2p")(positive_slope/s) + " positive slope");
svg.append ("line")
.attr ("x1", x (x.domain()[0]))
.attr ("y1", y (x.domain()[0]*fit_with_counts['slope'] + fit_with_counts ['intercept']))
.attr ("x2", x (x.domain()[1]))
.attr ("y2", y (x.domain()[1]*fit_with_counts['slope'] + fit_with_counts ['intercept']))
.attr ("class", "fit fit-counts");
svg.append ("text").attr ("x", x (x.domain()[1]))
.attr ("y", y (x.domain()[1]*fit_with_counts['slope'] + fit_with_counts ['intercept']))
.text (format_line (fit_with_counts['intercept'],fit_with_counts['slope']))
.attr ("class", "fit-label")
.attr ("dx", "2em")
.attr ("dy", "-.25 em");
}
$("#validate_newick").on ("click", function (e) {
var res = d3.layout.newick_parser ( $('textarea[id$="nwk_spec"]').val(), true);
if (res["error"] || ! res["json"]) {
var warning_div = d3.select ("#newick_body").selectAll ("div .alert-danger").data ([res["error"]])
warning_div.enter ().append ("div");
warning_div.html (function (d) {return d;}).attr ("class", "alert-danger");
} else {
drawATree ($('textarea[id$="nwk_spec"]').val());
$('#newick_modal').modal('hide');
}
});
$("[data-tree]").each (function () {
$(this).text ($(this).data("tag"));
});
$("[data-tree]").on ("click", function (e) {
loadTreeFromURL ($(this).data("tree"));
$("#menu-title").text ($(this).data("tag"));
});
$("[data-tree='tm1679-gag.tre']").click();
</script>
</body>
</html>
((((0_LN_DNA_1_1-1688:0.00748,(((0_LN_DNA_1_1-609:0.0,0_PB_DNA_1_1-1687:0.0,0_PL_RNA_1_1-108:0.0):0.00057,((0_PB_DNA_1_1-205:0.00369,((0_LN_DNA_1_1-613:0.0,0_PB_DNA_1_1-1004:0.0,6_LN_DNA_1_1-979:0.0):0.00057,(0_LN_DNA_1_1-139:0.00373,0_LN_DNA_1_1-124:0.00372)0.751:0.00057)1.000:0.00056)0.144:0.00373,(0_PB_DNA_1_1-47:0.0,3_PB_DNA_1_1-24:0.0):0.00368)1.000:0.00057)0.553:0.00057,3_PB_RNA_1_1-547:0.00372)0.867:0.00373)1.000:0.00056,3_PB_DNA_1_1-2537:0.00371)0.354:0.00372,6_PB_DNA_1_1-1620:0.00749)0.915:0.00056,((((6_LN_DNA_1_1-2316:0.0,6_PB_DNA_1_1-2267:0.0):0.01500,3_LN_RNA_1_1-33:0.00057)0.224:0.00372,(3_PB_RNA_1_1-3607:0.00373,(0_PL_RNA_1_1-295:0.00057,(0_PL_RNA_1_1-18:0.01883,0_PL_RNA_1_1-7:0.00372)0.138:0.00057)0.798:0.00374)0.789:0.00374)0.979:0.00057,(0_LN_DNA_1_1-520:0.0,0_PB_DNA_1_1-2408:0.0,0_PL_RNA_1_1-17:0.0,3_LN_RNA_1_1-3415:0.0,6_LN_DNA_1_1-109:0.0,6_PB_DNA_1_1-109:0.0):0.00057)0.006:0.00057,(((6_LN_DNA_1_1-2063:0.0,6_PB_DNA_1_1-2063:0.0):0.00057,(6_LN_DNA_1_1-1185:0.0,6_PB_DNA_1_1-1185:0.0):0.00371)0.874:0.00372,3_LN_RNA_1_1-164:0.00372)0.951:0.00057);
((6_PB_DNA_1_1-347:0.00222,(6_PB_DNA_1_1-156:0.00053,(6_PB_DNA_1_1-142:0.00052,(6_PB_DNA_1_1-415:0.00211,(0_LN_DNA_1_1-98:0.00248,(0_PB_DNA_1_1-153:0.0,6_LN_DNA_1_1-172:0.0):0.00443)0.879:0.00697)1.000:0.00052)0.860:0.00229)1.000:0.00051)0.311:0.00230,((((6_PB_DNA_1_1-354:0.00229,6_PB_DNA_1_1-398:0.00053)0.916:0.00192,((6_PB_DNA_1_1-200:0.00229,6_PB_DNA_1_1-60:0.00053)0.916:0.00192,(6_PB_DNA_1_1-237:0.00053,(6_PB_DNA_1_1-1220:0.00052,6_PB_DNA_1_1-51:0.00230)0.858:0.00230)1.000:0.00051)0.162:0.00231)0.998:0.00052,(6_PB_DNA_1_1-140:0.00053,6_PB_DNA_1_1-106:0.00230)0.761:0.00052)0.602:0.00229,(6_PB_DNA_1_1-59:0.00229,6_PB_DNA_1_1-247:0.00053)1.000:0.00052)0.761:0.00210,(((((3_PB_RNA_1_1-37:0.01153,((0_LN_DNA_1_1-31:0.0,3_PB_DNA_1_1-172:0.0,3_PB_RNA_1_1-1084:0.0,6_LN_DNA_1_1-59:0.0):0.00052,((3_PB_RNA_1_1-47:0.00053,(3_PB_RNA_1_1-49:0.00220,(3_PB_RNA_1_1-188:0.00220,3_PB_RNA_1_1-72:0.00053)0.726:0.00457)1.000:0.00052)0.962:0.00692,((3_PB_RNA_1_1-67:0.00053,3_PB_RNA_1_1-23:0.00220)0.678:0.00458,(6_LN_DNA_1_1-57:0.00230,3_PB_RNA_1_1-97:0.00052)0.896:0.00220)1.000:0.00051)0.786:0.00051)0.918:0.00052)0.871:0.00216,(((0_LN_DNA_1_1-44:0.0,6_LN_DNA_1_1-415:0.0):0.00052,(0_PB_DNA_1_1-52:0.00229,3_PB_DNA_1_2-16:0.00690)0.327:0.00052)0.837:0.00052,((0_LN_DNA_1_1-135:0.00229,(0_LN_DNA_1_1-40:0.00229,(0_LN_DNA_1_1-33:0.0,3_PB_DNA_1_1-37:0.0,6_LN_DNA_1_1-208:0.0):0.00052)0.152:0.00052)0.921:0.00192,(((0_PL_RNA_1_1-93:0.0,3_PB_DNA_1_1-14:0.0):0.00052,(3_PB_RNA_1_1-36:0.00458,3_PB_DNA_1_2-15:0.00459)0.159:0.00052)0.921:0.00192,((((0_LN_DNA_1_2-93:0.00231,0_PB_DNA_1_1-65:0.00231)0.806:0.00229,(0_LN_DNA_1_1-35:0.00230,3_PB_DNA_1_1-21:0.00230)0.000:0.00052)0.000:0.00052,0_PL_RNA_1_1-72:0.00052)0.749:0.00052,0_PL_RNA_1_1-97:0.00461)1.000:0.00052)0.174:0.00230)1.000:0.00052)0.466:0.00230)1.000:0.00053,(6_PB_DNA_1_1-248:0.00052,(((((0_LN_DNA_1_1-34:0.0,0_PL_RNA_1_1-209:0.0,3_PB_DNA_1_1-43:0.0):0.00052,(3_PB_DNA_1_2-14:0.00464,(3_LN_RNA_1_1-36:0.00229,((3_LN_RNA_1_1-7481:0.00052,(3_LN_RNA_1_1-324:0.00229,3_LN_RNA_1_1-38:0.00461)0.542:0.00052)0.000:0.00052,3_LN_RNA_1_1-30:0.00052)0.806:0.00052)0.798:0.00233)0.804:0.00229)0.846:0.00230,(((0_LN_DNA_1_1-32:0.00229,(3_PB_DNA_1_1-23:0.00458,(0_PB_DNA_1_1-77:0.00920,((0_PB_DNA_1_1-219:0.00230,0_PB_DNA_1_1-85:0.00051)0.823:0.00227,((0_LN_DNA_1_1-93:0.0,0_PB_DNA_1_1-924:0.0):0.00052,(0_PB_DNA_1_1-551:0.00053,((0_PB_DNA_1_1-783:0.0,6_LN_DNA_1_1-53:0.0):0.00221,(0_PB_DNA_1_1-1719:0.00052,0_PB_DNA_1_1-72:0.00229)1.000:0.00051)0.121:0.00229)0.863:0.00222)0.983:0.00052)0.331:0.00228)0.855:0.00050)0.858:0.00228)0.464:0.00053,(3_PB_DNA_1_1-22:0.0,6_LN_DNA_1_1-38:0.0):0.00052)0.998:0.00052,((6_LN_DNA_1_1-42:0.01363,3_PB_DNA_1_1-25:0.00750)0.390:0.00461,(0_PB_DNA_1_1-116:0.00459,(((0_LN_DNA_1_1-144:0.00460,(6_LN_DNA_1_1-50:0.00052,(3_PB_DNA_1_1-15:0.00229,((6_LN_DNA_1_1-69:0.00448,(3_PB_DNA_1_1-20:0.00465,3_PB_DNA_1_1-17:0.00690)0.903:0.00052)0.204:0.00229,(0_LN_DNA_1_1-120:0.00229,3_PB_DNA_1_2-13:0.00052)0.922:0.00219)1.000:0.00052)0.388:0.00052)0.875:0.00051)0.845:0.00228,3_PB_DNA_1_1-31:0.00460)0.560:0.00053,(3_PB_DNA_1_1-13:0.0,6_LN_DNA_1_1-247:0.0):0.00052)0.943:0.00052)1.000:0.00053)0.978:0.00179)0.810:0.00228)0.890:0.00052,((3_PB_DNA_1_1-19:0.0,3_PB_RNA_1_1-235:0.0):0.00052,(3_PB_DNA_1_1-16:0.00229,3_PB_RNA_1_1-32:0.00230)0.911:0.00052)1.000:0.00052)0.870:0.00052,0_PL_RNA_1_1-69:0.00229)0.854:0.00230)0.859:0.00229)0.839:0.00052,(6_LN_DNA_1_1-75:0.00228,0_PL_RNA_1_2-69:0.00229)0.783:0.00052)0.000:0.00052,(0_LN_DNA_1_1-293:0.0,3_PB_DNA_1_1-64:0.0,6_LN_DNA_1_1-579:0.0,6_PB_DNA_1_1-1278:0.0):0.00052)1.000:0.00052);
(3_PB_DNA_1_1-22:0.00496,0_PL_RNA_1_1-209:0.00248,(((3_PB_RNA_1_1-598:0.00052,3_PB_RNA_1_1-21:0.00247)0.936:0.00230,(3_PB_RNA_1_2-21:0.00230,(0_PL_RNA_1_1-48:0.0,3_PB_RNA_1_1-209:0.0):0.00052)0.449:0.00249)1.000:0.00053,((0_LN_DNA_1_1-120:0.0,0_PB_DNA_1_1-65:0.0,0_PL_RNA_1_1-156:0.0,3_PB_DNA_1_1-33:0.0,3_PB_RNA_1_1-136:0.0,6_PB_DNA_1_1-194:0.0):0.00052,((((0_LN_DNA_1_1-122:0.0,3_PB_DNA_1_1-26:0.0):0.00052,3_PB_DNA_1_1-21:0.00748)0.838:0.00244,((3_PB_DNA_1_1-20:0.00245,(((((3_LN_RNA_1_1-47:0.00247,3_LN_RNA_1_1-29:0.00248)0.885:0.00052,((3_LN_RNA_1_1-67:0.00247,3_LN_RNA_1_1-129:0.00496)0.909:0.00052,3_LN_RNA_1_1-4859:0.00052)0.000:0.00052)0.000:0.00052,3_LN_RNA_1_1-44:0.00052)0.853:0.00248,(0_PB_DNA_1_1-384:0.00247,(0_PB_DNA_1_1-46:0.00052,((0_PB_DNA_1_1-140:0.00237,0_PB_DNA_1_1-54:0.00052)0.433:0.00251,(0_PB_DNA_1_1-42:0.00052,3_PB_DNA_1_1-52:0.01253)0.888:0.00052)0.892:0.00485)0.915:0.00501)0.795:0.00248)0.376:0.00052,6_LN_DNA_1_1-208:0.00745)0.953:0.00753)0.907:0.00495,3_PB_RNA_1_2-18:0.00051)0.352:0.00248)0.962:0.00051,(((0_PB_DNA_1_1-48:0.00051,(0_PB_DNA_1_1-87:0.00053,(0_PB_DNA_1_1-678:0.00052,0_PB_DNA_1_1-47:0.00052)0.853:0.00243)0.342:0.00245)0.914:0.00488,(((6_LN_DNA_1_1-84:0.0,6_PB_DNA_1_1-561:0.0):0.00052,0_PB_DNA_1_1-61:0.00248)0.861:0.00249,(0_PL_RNA_1_1-45:0.00052,(0_PL_RNA_1_1-98:0.00247,0_PL_RNA_1_1-64:0.00052)0.847:0.00246)0.854:0.00248)0.900:0.00049)0.601:0.00249,((3_PB_RNA_1_1-19:0.00052,(3_PB_RNA_1_1-48:0.00052,(3_PB_RNA_1_1-17:0.00051,(3_PB_RNA_1_1-18:0.00052,3_PB_RNA_1_1-267:0.00238)0.482:0.00246)0.915:0.00236)0.921:0.00216)0.855:0.00497,((0_PB_DNA_1_1-91:0.00051,0_PB_DNA_1_1-49:0.00489)0.460:0.00246,(0_LN_DNA_1_1-83:0.00991,((0_LN_DNA_1_1-69:0.00494,((0_LN_DNA_1_1-45:0.00246,6_LN_DNA_1_1-111:0.01234)0.801:0.00245,(0_PL_RNA_1_1-50:0.00052,(6_PB_DNA_1_1-526:0.00492,0_LN_DNA_1_1-48:0.00246)0.370:0.00052)0.955:0.00053)0.816:0.00246)0.586:0.00052,6_PB_DNA_1_1-133:0.00052)0.871:0.00052)1.000:0.00052)0.921:0.00214)1.000:0.00052)0.886:0.00052)0.921:0.00052)0.754:0.00052)0.536:0.00052);
((0_PB_DNA_1_21-1:0.0,6_LN_DNA_1_1-400:0.0):0.00052,((((0_PB_DNA_1_1-4:0.0,3_LN_RNA_1_1-332:0.0,6_LN_DNA_1_1-447:0.0):0.00052,(((0_PB_DNA_1_3-2:0.00920,0_PB_DNA_1_36-1:0.00052)0.466:0.00310,((0_PL_RNA_1_1-76:0.00052,(0_PB_DNA_1_2-1:0.00052,0_PB_DNA_1_22-1:0.00052)0.891:0.00613)0.415:0.00307,((0_PB_DNA_1_8-1:0.00306,0_PB_DNA_1_34-1:0.00306)0.418:0.00052,((0_PB_DNA_1_1-2:0.0,0_PL_RNA_1_1-253:0.0,3_LN_RNA_1_1-1745:0.0):0.00052,((3_LN_RNA_1_1-37:0.00052,3_LN_RNA_1_1-17:0.00306)0.857:0.00302,(3_LN_RNA_1_1-154:0.00051,(3_LN_RNA_1_1-69:0.00052,(3_LN_RNA_1_1-186:0.00052,3_LN_RNA_1_1-21:0.00310)0.843:0.00309)0.834:0.00303)0.108:0.00307)1.000:0.00051)0.818:0.00052)0.995:0.00051)0.778:0.00298)1.000:0.00052,0_PB_DNA_1_19-1:0.00936)0.724:0.00052)0.297:0.00052,(((3_LN_RNA_1_1-28:0.00052,(((0_PL_RNA_1_1-77:0.00052,0_PB_DNA_1_1-1:0.00916)0.786:0.00613,(0_PB_DNA_1_5-1:0.00306,0_PB_DNA_1_11-1:0.00923)0.987:0.00051)0.813:0.00305,(0_PB_DNA_1_4-1:0.00309,0_PB_DNA_1_29-1:0.00306)0.000:0.00052)0.925:0.00291)0.181:0.00309,(((0_PB_DNA_1_27-1:0.00306,6_LN_DNA_1_1-177:0.00052)0.000:0.00052,((3_PB_RNA_1_1-2510:0.00052,3_PB_RNA_1_1-243:0.00307)0.842:0.00305,(0_PB_DNA_1_20-1:0.00052,(0_PB_DNA_1_6-1:0.01042,(0_PB_DNA_1_12-1:0.00771,0_PB_DNA_1_26-1:0.00145)0.154:0.00308)0.732:0.00190)0.126:0.00051)0.782:0.00307)0.000:0.00052,(((0_PB_DNA_1_31-1:0.00308,((0_PB_DNA_1_2-2:0.0,0_PL_RNA_1_1-93:0.0):0.00052,0_PB_DNA_1_16-1:0.00308)0.646:0.00051)0.788:0.00306,(0_PL_RNA_1_1-147:0.00052,(0_PB_DNA_1_10-1:0.00621,0_PB_DNA_1_14-1:0.00306)0.746:0.00052)0.897:0.00051)0.139:0.00051,0_PB_DNA_1_35-1:0.00612)0.435:0.00306)0.864:0.00288)1.000:0.00052,0_LN_DNA_1_1-6706:0.00620)0.000:0.00052)0.831:0.00308,0_PB_DNA_1_32-1:0.00307)0.780:0.00051,((0_PB_DNA_1_30-1:0.00926,0_PB_DNA_1_38-1:0.00306)0.896:0.00052,(((0_PB_DNA_1_24-1:0.00914,(0_PB_DNA_1_3-1:0.00300,0_PB_DNA_1_28-1:0.00611)0.797:0.00326)0.207:0.00304,(((0_LN_DNA_1_1-43:0.0,6_LN_DNA_1_1-28:0.0):0.00052,(0_PB_DNA_1_9-1:0.00305,(0_PB_DNA_1_18-1:0.01226,(0_PL_RNA_1_1-410:0.00052,((0_PB_DNA_1_15-1:0.00303,0_PB_DNA_1_33-1:0.00052)0.862:0.00303,0_PB_DNA_1_37-1:0.00305)1.000:0.00052)1.000:0.00051)0.847:0.00303)0.861:0.00052)0.000:0.00052,(0_PB_DNA_1_23-1:0.00303,(0_PB_DNA_1_7-1:0.00052,(0_PB_DNA_1_13-1:0.00612,0_PB_DNA_1_17-1:0.00304)0.000:0.00052)0.852:0.00305)0.242:0.00052)1.000:0.00052)0.877:0.00294,0_PB_DNA_1_25-1:0.00305)0.922:0.00052)0.010:0.00052);
((((3_LN_RNA_1_1-177:0.00050,(3_LN_RNA_1_1-24:0.00050,(3_LN_RNA_1_1-276:0.00050,(3_LN_RNA_1_1-697:0.00052,3_LN_RNA_1_1-22:0.00271)0.436:0.00292)0.923:0.00266)0.448:0.00289)0.813:0.00285,0_PL_RNA_1_1-306:0.00291)0.874:0.00051,(((0_PL_RNA_1_1-157:0.00051,(0_PL_RNA_1_1-159:0.00051,0_PL_RNA_1_1-41:0.00289)0.847:0.00284)0.831:0.00581,(0_PL_RNA_1_1-244:0.00051,(0_PL_RNA_1_1-39:0.00289,0_PL_RNA_1_1-69:0.00288)0.397:0.00051)0.988:0.00050)0.848:0.00284,(((0_PL_RNA_1_1-100:0.00290,((0_LN_DNA_1_1-159:0.0,0_PL_RNA_1_1-1077:0.0):0.00051,(0_PL_RNA_1_2-40:0.00289,((0_LN_DNA_1_1-8809:0.0,6_PB_DNA_1_1-6204:0.0):0.00051,(0_LN_DNA_1_1-103:0.00288,(0_LN_DNA_1_1-1254:0.0,6_PB_DNA_1_1-623:0.0):0.00289)0.306:0.00051)0.863:0.00289)0.677:0.00051)0.836:0.00051)1.000:0.00051,0_PL_RNA_1_1-48:0.00270)0.104:0.00290,0_PL_RNA_1_1-219:0.00271)1.000:0.00050)0.898:0.00051)0.994:0.00051,(((0_PB_DNA_1_1-30:0.00584,(0_LN_DNA_1_1-182:0.0,0_PL_RNA_1_1-40:0.0):0.00581)0.045:0.00051,0_PL_RNA_1_3-40:0.00290)1.000:0.00051,(0_LN_DNA_1_1-1209:0.0,0_PB_DNA_1_1-749:0.0,0_PL_RNA_1_1-3677:0.0,3_LN_RNA_1_1-2619:0.0,6_LN_DNA_1_1-95:0.0):0.00051)0.000:0.00051,((0_PL_RNA_1_1-58:0.00290,0_PL_RNA_1_1-154:0.00581)0.915:0.00051,(((3_LN_RNA_1_1-129:0.00051,3_LN_RNA_1_1-62:0.00290)0.772:0.00289,((3_LN_RNA_1_1-20:0.00286,(3_LN_RNA_1_1-28:0.00051,((0_PL_RNA_1_1-89:0.00051,((((3_PB_RNA_1_1-23:0.00289,3_PB_RNA_1_1-39:0.00295)0.362:0.00051,(3_PB_RNA_1_1-159:0.00291,3_PB_RNA_1_1-18:0.00292)0.589:0.00051)0.000:0.00051,(((3_PB_RNA_1_1-442:0.00289,3_PB_RNA_1_1-15:0.00290)0.939:0.00051,3_PB_RNA_1_1-4606:0.00051)0.000:0.00051,(3_PB_RNA_1_1-17:0.00290,3_PB_RNA_1_1-16:0.00290)0.737:0.00051)0.000:0.00051)0.000:0.00051,(3_PB_RNA_1_1-518:0.00289,3_PB_RNA_1_1-21:0.00291)0.278:0.00051)0.835:0.00291)0.880:0.00292,3_LN_RNA_1_1-35:0.00291)0.867:0.00051)1.000:0.00051)0.259:0.00291,((3_LN_RNA_1_1-275:0.00273,(3_LN_RNA_1_1-45:0.00051,3_LN_RNA_1_1-18:0.00273)0.160:0.00292)1.000:0.00050,((0_LN_DNA_1_1-284:0.0,0_PL_RNA_1_1-763:0.0,3_LN_RNA_1_1-1669:0.0):0.00051,(3_LN_RNA_1_2-20:0.00292,3_LN_RNA_1_1-15:0.00290)0.344:0.00051)0.819:0.00051)0.774:0.00286)0.300:0.00051)0.820:0.00050,((3_LN_RNA_1_1-11:0.00286,(3_LN_RNA_1_1-12:0.00286,3_LN_RNA_1_1-25:0.00050)0.430:0.00294)1.000:0.00051,(((0_PL_RNA_1_1-57:0.0,3_LN_RNA_1_1-136:0.0):0.00051,0_LN_DNA_1_1-325:0.00291)0.844:0.00290,(((((0_PL_RNA_1_1-67:0.00271,((0_PL_RNA_1_1-274:0.00051,(0_PB_DNA_1_1-1186:0.0,0_PL_RNA_1_1-153:0.0):0.00290)1.000:0.00051,0_PL_RNA_1_1-168:0.00272)0.178:0.00291)1.000:0.00051,0_PL_RNA_1_1-97:0.00292)0.718:0.00051,0_PL_RNA_1_1-1398:0.00051)1.000:0.00051,0_PL_RNA_1_1-44:0.00275)0.183:0.00291,0_PL_RNA_1_1-635:0.00274)1.000:0.00051)0.445:0.00051)0.000:0.00051)0.000:0.00051)0.479:0.00051);
((((0_PB_DNA_1_1-2113:0.0,0_PL_RNA_1_1-212:0.0):0.00049,0_PB_DNA_1_1-213:0.00227)0.849:0.00228,((3_PB_RNA_1_1-20:0.00453,(3_PB_RNA_1_1-35:0.00049,3_PB_RNA_1_1-11:0.00224)0.879:0.00458)0.918:0.00050,(3_PB_RNA_1_1-3167:0.00049,(3_PB_RNA_1_1-43:0.00456,3_PB_RNA_1_1-70:0.00227)0.897:0.00049)0.584:0.00049)0.837:0.00227)0.893:0.00049,(0_LN_DNA_1_1-46:0.00914,0_PL_RNA_1_1-58:0.00459)0.541:0.00049,(((((0_LN_DNA_1_1-203:0.0,0_PL_RNA_1_1-205:0.0,3_LN_RNA_1_1-262:0.0):0.00049,((3_LN_RNA_1_1-192:0.00049,3_LN_RNA_1_2-21:0.00228)0.782:0.00229,(((0_LN_DNA_1_1-112:0.0,6_PB_DNA_1_1-61:0.0):0.00049,(((0_LN_DNA_1_1-89:0.0,6_PB_DNA_1_1-29:0.0):0.00460,(0_LN_DNA_1_1-191:0.0,6_PB_DNA_1_1-96:0.0):0.00926)0.341:0.00049,(6_PB_DNA_1_1-24:0.00049,(0_LN_DNA_1_1-2289:0.0,6_PB_DNA_1_1-2352:0.0):0.00049)0.000:0.00049)0.858:0.00229)0.900:0.00454,(3_LN_RNA_1_1-20:0.00225,(3_LN_RNA_1_1-44:0.00228,3_LN_RNA_1_2-44:0.00049)1.000:0.00050)0.427:0.00228)0.046:0.00049)0.739:0.00049)0.950:0.00048,(0_PL_RNA_1_1-42:0.00228,(0_PL_RNA_1_1-118:0.00049,0_PL_RNA_1_1-61:0.00228)0.168:0.00049)0.904:0.00456)0.427:0.00228,0_PL_RNA_1_1-67:0.00916)0.897:0.00048,(((3_LN_RNA_1_1-41:0.00049,(3_LN_RNA_1_1-263:0.00228,3_LN_RNA_1_1-21:0.00228)0.160:0.00049)0.854:0.00226,((0_PL_RNA_1_1-75:0.0,6_LN_DNA_1_1-2242:0.0):0.00049,6_LN_DNA_1_1-80:0.00227)0.853:0.00227)0.515:0.00049,(0_PB_DNA_1_1-1463:0.0,0_PL_RNA_1_1-94:0.0):0.00049)0.929:0.00049)0.000:0.00049);
(((0_PB_DNA_1_1-369:0.00666,(0_LN_DNA_1_1-124:0.00051,6_LN_DNA_1_1-77:0.02377)0.805:0.00332)0.808:0.00330,((0_PL_RNA_1_1-2637:0.00051,((0_PL_RNA_1_1-100:0.00050,(0_PL_RNA_1_1-623:0.00050,0_PL_RNA_1_1-24:0.00328)0.849:0.00331)0.850:0.00332,0_PL_RNA_1_1-360:0.00050)0.846:0.00331)0.173:0.00050,0_PB_DNA_1_1-105:0.00330)0.852:0.00329)0.444:0.00051,3_LN_RNA_1_1-13:0.00051,((6_PB_DNA_1_1-78:0.00328,0_LN_DNA_1_1-52:0.00328)0.922:0.00051,(((0_LN_DNA_1_1-3699:0.0,0_PB_DNA_1_1-136:0.0,3_LN_RNA_1_1-2407:0.0,6_PB_DNA_1_1-316:0.0):0.00051,3_LN_RNA_1_1-188:0.00331)0.000:0.00051,(0_LN_DNA_1_1-575:0.00329,0_LN_DNA_1_1-53:0.00329)0.854:0.00051)0.000:0.00051)0.000:0.00051);
(6_LN_DNA_1_1-4915:0.00228,6_PB_DNA_1_1-50:0.00227,(0_PB_DNA_1_1-53:0.00228,(((((3_LN_RNA_1_1-585:0.00052,3_LN_RNA_1_1-27:0.00229)0.824:0.00228,((3_LN_RNA_1_1-29:0.00220,(6_PB_DNA_1_1-288:0.00229,6_PB_DNA_1_1-137:0.00052)1.000:0.00052)0.745:0.00458,(6_LN_DNA_1_1-45:0.00457,((((((6_PB_DNA_1_1-39:0.00228,(0_PL_RNA_1_1-226:0.00228,0_PL_RNA_1_1-209:0.00052)0.812:0.00228)0.786:0.00229,(0_PL_RNA_1_1-240:0.00222,(0_PL_RNA_1_1-306:0.00052,0_PL_RNA_1_1-38:0.00221)0.402:0.00228)1.000:0.00051)0.748:0.00052,(0_LN_DNA_1_1-173:0.0,0_PB_DNA_1_1-49:0.0,0_PL_RNA_1_1-2144:0.0,3_LN_RNA_1_1-1230:0.0,6_LN_DNA_1_1-75:0.0,6_PB_DNA_1_1-1180:0.0):0.00052)0.994:0.00052,(6_LN_DNA_1_1-53:0.00229,0_PL_RNA_1_1-201:0.00228)0.734:0.00052)0.728:0.00052,(((6_PB_DNA_1_1-844:0.00210,(6_PB_DNA_1_1-376:0.00227,6_PB_DNA_1_1-230:0.00051)1.000:0.00050)0.457:0.00229,(6_PB_DNA_1_1-170:0.00227,6_PB_DNA_1_1-573:0.00052)0.926:0.00209)1.000:0.00052,(6_PB_DNA_1_1-104:0.00227,0_PL_RNA_1_1-198:0.00228)0.587:0.00052)0.000:0.00052)0.932:0.00053,0_PL_RNA_1_1-42:0.00228)0.477:0.00051)0.845:0.00219)1.000:0.00052)0.728:0.00052,((0_PL_RNA_1_1-635:0.00052,0_PL_RNA_1_1-77:0.00228)0.825:0.00228,(3_LN_RNA_1_1-83:0.00051,(3_LN_RNA_1_1-26:0.00221,(3_LN_RNA_1_2-73:0.00051,(3_LN_RNA_1_1-28:0.00229,(3_LN_RNA_1_1-1554:0.00052,3_LN_RNA_1_1-76:0.00231)0.265:0.00052)0.890:0.00221)0.178:0.00228)1.000:0.00052)0.943:0.00455)0.849:0.00051)0.000:0.00052,(0_LN_DNA_1_1-154:0.0,0_PB_DNA_1_1-190:0.0,0_PL_RNA_1_1-396:0.0,3_LN_RNA_1_1-73:0.0,6_PB_DNA_1_1-2420:0.0):0.00052)0.000:0.00052,(6_PB_DNA_1_1-174:0.00228,(0_LN_DNA_1_1-6843:0.00052,0_LN_DNA_1_1-37:0.00226)0.926:0.00454)0.877:0.00052)0.578:0.00052)0.525:0.00052);
((0_LN_DNA_1_1-164:0.0,0_PL_RNA_1_1-31:0.0,3_LN_RNA_1_1-952:0.0,6_LN_DNA_1_1-88:0.0):0.00054,((((((((0_LN_DNA_1_1-507:0.0,0_PL_RNA_1_1-139:0.0,6_PB_DNA_1_1-34:0.0):0.00054,6_LN_DNA_1_1-152:0.00248)0.852:0.00242,(6_PB_DNA_1_1-33:0.00054,(((0_LN_DNA_1_1-24:0.0,3_LN_RNA_1_1-32:0.0):0.00054,3_LN_RNA_1_1-31:0.00248)0.853:0.00248,6_PB_DNA_1_1-53:0.00249)0.308:0.00054)1.000:0.00053)0.849:0.00497,(6_LN_DNA_1_1-56:0.00249,((0_LN_DNA_1_1-577:0.0,0_PB_DNA_1_1-28:0.0,0_PL_RNA_1_1-238:0.0,6_LN_DNA_1_1-42:0.0,6_PB_DNA_1_1-146:0.0):0.00054,(6_PB_DNA_1_1-54:0.00250,(((0_LN_DNA_1_1-193:0.00242,(6_LN_DNA_1_1-149:0.00054,(0_LN_DNA_1_1-56:0.00248,6_LN_DNA_1_1-62:0.00054)0.893:0.00497)0.625:0.00054)0.941:0.00750,(6_LN_DNA_1_1-716:0.00054,(6_LN_DNA_1_1-86:0.00054,6_LN_DNA_1_1-36:0.00498)0.845:0.00247)1.000:0.00054)0.845:0.00240,0_PL_RNA_1_1-41:0.00249)0.390:0.00054)0.000:0.00054)0.914:0.00054)0.863:0.00242)1.000:0.00053,6_PB_DNA_1_1-41:0.00247)0.253:0.00054,(0_PB_DNA_1_1-45:0.0,6_LN_DNA_1_1-75:0.0,6_PB_DNA_1_1-152:0.0):0.00054)0.523:0.00051,((0_PB_DNA_1_1-65:0.0,0_PL_RNA_1_1-33:0.0,6_LN_DNA_1_1-24:0.0,6_PB_DNA_1_1-322:0.0):0.00054,(6_PB_DNA_1_1-139:0.00495,6_PB_DNA_1_1-92:0.00053)0.718:0.00495)0.902:0.00244)0.865:0.00246,0_PL_RNA_1_1-21:0.00246)0.231:0.00054,((0_PB_DNA_1_1-24:0.00248,6_LN_DNA_1_1-55:0.00497)0.900:0.00054,((3_LN_RNA_1_1-669:0.00054,3_LN_RNA_1_1-44:0.00247)0.859:0.00248,0_PL_RNA_1_1-25:0.00246)0.620:0.00054)1.000:0.00054);
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment