2f:T78c,summary_table = html`

${fel_results_json.input["number of sequences"]} sequences in the alignment	${fel_results_json.input["number of sites"]} codon sites in the alignment	${fel_results_json.input["partition count"]} partitions
${tested_branch_count} median branches/partition used for testing	${variable_site_count} non-invariant sites tested	${fel_results_json.simulated \|\| "N/A"} parametric bootstrap replicates
${_.filter (sites_table[1], (d)=>d.class == "Diversifying").length} Sites under diversifying positive selection at p≤${pv}	${_.filter (sites_table[1], (d)=>d.class == "Purifying").length} Sites under purifying selection at p≤${pv}

`30:T7d1,sites_table = { let results = []; const headers = _.clone(fel_results_json.MLE.headers); let format = {}; format [headers[0][0]] = (d)=>d.toFixed (3); format [headers[1][0]] = (d)=>d.toFixed (3); format [headers[2][0]] = (d)=>d.toFixed (3); format [headers[3][0]] = (d)=>d.toFixed (3); format [headers[4][0]] = (d)=>d <= pv ? html`${d.toFixed (4)}` : d.toFixed (4); if (has_pasmt) { format[headers[headers.length-1][0]] = format [headers[4][0]]; } format [headers[5][0]] = (d)=>d.toFixed (3); headers.push (["class","Site classification at p<=" + pv]); format["class"] = (d)=>html`${d}`; _.each (fel_results_json.MLE.content, (data, part)=> { const site_lookup = fel_results_json["data partitions"][part].coverage[0]; _.each (data, (row, i)=> { let row_object = { 'partition' : (+part) + 1, 'codon' : site_lookup [i] + 1 }; row_object[headers[0][0]] = +row[0]; row_object[headers[1][0]] = +row[1]; row_object[headers[2][0]] = +row[2]; row_object[headers[3][0]] = +row[3]; row_object[headers[4][0]] = +row[4]; if (has_T) { row_object[headers[5][0]] = +row[5]; } if (has_ci) { row_object[headers[6][0]] = row[6]; row_object[headers[7][0]] = row[7]; row_object[headers[8][0]] = row[8]; } if (has_pasmt) { row_object [headers[headers.length-2][0]] = row[ headers.length-2]; } row_object[headers[headers.length-1][0]] = row[4] <= pv ? (row[0] < row[1] ? "Diversifying" : "Purifying") : (row[0] + row[1] ? "Neutral" : "Invariable"); results.push (row_object); }); }); return [format, results, headers]; }31:T90c,dNdS_with_ci = (data, from, step)=> { let color_d = []; let color_r = []; _.each (table_colors, (v,c)=> {color_d.push (c); color_r.push (v);}); return { "width": {"step": 12}, "data" : {"values" : _.map ( _.filter (data, (d,i)=>i+1 >= from && i<= from + step), (d)=> { let dd = _.clone (d); _.each (["dN/dS LB","dN/dS UB", "dN/dS MLE"], (f)=> { dd[f] = Math.min (dyn_range_cap, dd[f]); }); return dd; })}, "transform" :[{"calculate": "1", "as": "neutral"}], "encoding": { "x": { "field": "codon", "type" : "nominal", "axis": {"grid" : false, "titleFontSize" : 14, "title" : "Codon"} } }, "layer": [ { "mark": {"opacity": 1., "type": "line", "color": "steelblue"}, "encoding": { "y": { "field": "dN/dS LB", "scale" : {"type" : "sqrt"}, "type" : "quantitative", "axis": {"titleColor": "black", "grid" : false, "titleFontSize" : 14, "title" : "dN/dS"} }, "y2": { "field": "dN/dS UB", "type" : "quantitative" } } }, { "mark": {"stroke": "black", "type": "point", "filled" : true, "size" : 64, "stroke" : null, "tooltip" : {"contents" : "data"}}, "encoding": { "y": { "field": "dN/dS MLE", "type" : "quantitative", }, "color" : {"field" : "class", "scale" : {"domain" : color_d, "range" : color_r}, "title" : "Selection class"} } }, /*{ "mark": {"stroke": "lightgrey", "type": "line", "size" : 1, "opacity" : 1.}, "encoding": { "y": { "field": "dN/dS MLE", "type" : "quantitative", }, } }*/ { "mark" : {"opacity": 0.5, "type": "line", "color": "gray"}, "encoding" : { "y": { "field": "neutral", "type" : "quantitative", }, "x": { "field": "codon", "type" : "nominal"}, "size": {"value": 2}, } } ] }; }32:Tc29,html` `33:[{"id":"f90985fa46d9ad5e9ea51b3c83628745ad9b6d5f3df4cd7e5c93f5a1ff0ee7fb765099785a5a127bfda55393d1200cc1d011b8acb463931ac69b6684c2cf57cb","href":"https://static.observableusercontent.com/files/f90985fa46d9ad5e9ea51b3c83628745ad9b6d5f3df4cd7e5c93f5a1ff0ee7fb765099785a5a127bfda55393d1200cc1d011b8acb463931ac69b6684c2cf57cb","name":"BCL6_orthologs.csv","type":"text/csv","lastModified":1713552411666.986,"size":75165},{"id":"b2eef377ec409c7797ab08e090dfeb69882061a2ffd07650ef5b09c973d7096d3c0bc038b85d6edd25ae9d7ab7d76397daa54096bf047879cf293c893d889cff","href":"https://static.observableusercontent.com/files/b2eef377ec409c7797ab08e090dfeb69882061a2ffd07650ef5b09c973d7096d3c0bc038b85d6edd25ae9d7ab7d76397daa54096bf047879cf293c893d889cff","name":"GnathoBCL6.RDP.fasta.dst","type":"application/octet-stream","lastModified":1713552491977.308,"size":16579128},{"id":"a6edc9bacb3664a28eb41fb1ad5b643e2b479181c4142ec63ddb2435e4a5cb65668046b7bca2d27c8861392e6fda2a7af7db7055938108a825e19ca24a4c1b21","href":"https://static.observableusercontent.com/files/a6edc9bacb3664a28eb41fb1ad5b643e2b479181c4142ec63ddb2435e4a5cb65668046b7bca2d27c8861392e6fda2a7af7db7055938108a825e19ca24a4c1b21","name":"GnathoBCL6.RDP.fasta","type":"application/octet-stream","lastModified":1713552592275.249,"size":2524087},{"id":"a0595fdf0efb56387d802e67359ebda22422797617f1b3bdf61ca0406ea7cfce5df791f4c1e368c30a1b65c8d84c32819b8aa9c53ee92ff25d9fef3a01dd8558","href":"https://static.observableusercontent.com/files/a0595fdf0efb56387d802e67359ebda22422797617f1b3bdf61ca0406ea7cfce5df791f4c1e368c30a1b65c8d84c32819b8aa9c53ee92ff25d9fef3a01dd8558","name":"GnathoBCL6.RDP.fasta.treefile","type":"application/octet-stream","lastModified":1713552744928.126,"size":44693},{"id":"2aa4959eeb69944ac9ba84220a36ef244a70cf7e07ee5401a247e7dbbfa15f4ac680e8c0c8e21707f648e2a92200e5bb6dbe6eab9fe4d8883c1f4d65b0f99850","href":"https://static.observableusercontent.com/files/2aa4959eeb69944ac9ba84220a36ef244a70cf7e07ee5401a247e7dbbfa15f4ac680e8c0c8e21707f648e2a92200e5bb6dbe6eab9fe4d8883c1f4d65b0f99850","name":"GnathoBCL6.RDP.fasta.FEL.json","type":"application/json","lastModified":1713552757113.175,"size":465315},{"id":"8bf539c1787d415d82b92105a32bcdd650288e0532f7e0a6fb6871af8419228acfa99aeaa2f80bce34a6d25d3b85e5457817187d9d34207d19bfb0ffd0ad0811","href":"https://static.observableusercontent.com/files/8bf539c1787d415d82b92105a32bcdd650288e0532f7e0a6fb6871af8419228acfa99aeaa2f80bce34a6d25d3b85e5457817187d9d34207d19bfb0ffd0ad0811","name":"humanBCL6_24983_unrelaxed_rank_001_alphafold2_ptm_model_3_seed_000.pdb","type":"application/octet-stream","lastModified":1714679733711.263,"size":446877}] 2d:{"title":"Selection analyses of GnathosomataBCL6","cells":[{"id":0,"value":"# Selection analyses of *${tag}${gene}*\n\nWe provide an evolutionary analysis for ${tag} ${gene}, which includes ${orthologs_table.length} species.\n\nData was retrieved from the NCBI Orthologs database: https://www.ncbi.nlm.nih.gov/gene/604/ortholog/?scope=7776&term=BCL6\n\nWorkflow: LIPSA","pinned":false,"hidden":false,"mode":"md"},{"id":421,"value":"## Table of species used (Orthologs in *${tag}*)","pinned":false,"hidden":false,"mode":"md"},{"id":9,"value":"//Inputs.table(btg1_orthologs_table)\n// Scientific name\tCommon name\tRefSeq Transcript accessions\tRefSeq Protein accessions\n\nInputs.table(orthologs_table, {\n columns: [\n \"Scientific name\",\n \"Common name\",\n \"RefSeq Transcript accessions\",\n \"RefSeq Protein accessions\"\n ],\n header: {}\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":36,"value":"## Genetic distances (nucleotide-level)","pinned":false,"hidden":false,"mode":"md"},{"id":762,"value":"This chart shows pairwise genetic distances calculated via the Tamura-Nei 1993 (TN93) method.","pinned":false,"hidden":false,"mode":"md"},{"id":42,"value":"Plot.plot({\n y: {grid: true},\n color: {legend: true},\n marks: [\n Plot.rectY(tn93_data, Plot.binX({y: \"count\"}, {x: \"Distance\"})),\n Plot.ruleY([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":165,"value":"## Evidence of selection on individual codon sites","pinned":false,"hidden":false,"mode":"md"},{"id":423,"value":"**High-level overview**","pinned":false,"hidden":false,"mode":"md"},{"id":196,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs"},{"id":425,"value":"**Site by site evolutionary rate**","pinned":false,"hidden":false,"mode":"md"},{"id":252,"value":"viewof table_filter = Inputs.checkbox([\"Diversifying\", \"Purifying\", \"Neutral\",\"Invariable\", \"HumanMapped\"], {value: [\"Diversifying\", \"Purifying\", \"Neutral\", \"Invariable\"], label: html`Show`, format: x => html`${x}`})","pinned":false,"hidden":false,"mode":"ojs"},{"id":237,"value":"// plot_options = [[\"Site-level dN/dS estimates\",(d)=>d[\"confidence interval\"]],[\"alpha/beta site-level estimates\", (d)=>1], [\"Bootstrap vs asymptotic p-value\", (d)=>has_pasmt], [\"Rate density plots\", (d)=>1], [\"Q-Q plots\", (d)=>has_pasmt], [\"Dense rate plot\", (d)=>1]]\n\nfigure1 = vegalite (plot_specs[\"Site-level dN/dS estimates\"])","pinned":false,"hidden":false,"mode":"ojs"},{"id":429,"value":"**Table of results for selection analysis**","pinned":false,"hidden":false,"mode":"md"},{"id":248,"value":"viewof table1 = Inputs.table (siteTableData, {\n rows : 40,\n format: sites_table[0]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":617,"value":"# Human Bcl6 Protein-Structure","pinned":false,"hidden":false,"mode":"md"},{"id":614,"value":"viewport = {\n let v = DOM.element(\"div\");\n v.id=\"viewport\";\n v.style.height = \"500px\"\n return v;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":438,"value":"# End of report ","pinned":false,"hidden":false,"mode":"md"},{"id":14,"value":"## Code behind the scenes -----------------------------------------","pinned":false,"hidden":false,"mode":"md"},{"id":494,"value":"##### File Imports","pinned":false,"hidden":false,"mode":"md"},{"id":7,"value":"orthologs_table = FileAttachment(\"BCL6_orthologs.csv\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs"},{"id":64,"value":"protein_file = FileAttachment(\"GnathoBCL6.RDP.fasta\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":127,"value":"nwk = FileAttachment(\"GnathoBCL6.RDP.fasta.treefile\").text()","pinned":false,"hidden":false,"mode":"ojs"},{"id":186,"value":"fel_results_json = FileAttachment(\"GnathoBCL6.RDP.fasta.FEL.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":324,"value":"PDB = FileAttachment(\"humanBCL6_24983_unrelaxed_rank_001_alphafold2_ptm_model_3_seed_000.pdb\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":297,"value":"##### TN93\n","pinned":false,"hidden":false,"mode":"md"},{"id":40,"value":"//tn93_data = FileAttachment(\"Primates_codons.SA.fasta.dst\").csv({typed: true})\n//Primates_codons.SA.TrimmedToHuman.fasta.dst\ntn93_data = FileAttachment(\"GnathoBCL6.RDP.fasta.dst\").csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs"},{"id":54,"value":"// Forked from: https://observablehq.com/@howardyao/sequence-alignment-dashboard\nprotein_data = {\n const url = '' \n const raw = (\n protein_file === null \n ? a.fastaParser(await d3.text(url)) \n : a.fastaParser(await protein_file.text())\n ).map(d => ({...d, \n 'seq_len': d.seq.length,\n 'seq_name': d.header // .split('|')[2].split(' ')[0] \n }));\n \n const max_len = d3.max(raw, d => d.seq_len);\n \n const raw1 = raw.map(d => ({...d, \n 'seq_full_str': d.seq + '-'.repeat(max_len - d.seq_len),\n 'seq_full': [...(d.seq + '-'.repeat(max_len - d.seq_len))].map((n,i) => ({'k': i, 'v': n, 'n': d.seq_name})) \n }));\n\n const derive = raw1.map(k => k.seq_full.map(j => ({'n': j.n, 'k': j.k + 1, 'v': j.v }))).flat()\n\n return {'raw': raw1, 'derive': derive, 'total_len': max_len}\n}\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":74,"value":"protein_data_viz = {\n const raw_len = protein_data.raw.length\n const countCol = d3.flatRollup(protein_data.derive, v => v.length, d => d.k, d => d.v)\n .map(([k, v, cnt]) => ({k, v, cnt}))\n\n const conversion = d3.groups(countCol, d=>d.k)\n .map(n => d3.greatest(n[1], a=>a.cnt))\n .map(d => ({...d, 'Conversion Rate': d.cnt / raw_len}))\n\n const derive2 = protein_data.derive.map(d => ({...d, 'mode': conversion.filter(j => j.k + 1 === d.k + 1)[0].v}))\n\n const data_viz = derive2.concat(conversion.map(d => ({'n': 'Consensus', 'k': d.k, 'v': d.v, 'mode': d.v})))\n\n return {'heatmap': data_viz, 'bar' : conversion}\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":612,"value":"height = 600;","pinned":false,"hidden":false,"mode":"ojs"},{"id":98,"value":"renderedTree = tree.render({\n height:height, \n width:width,\n 'left-right-spacing': 'fit-to-size', \n 'align-tips':true,\n 'node-styler': colorNodesByName,\n 'edge-styler': colorEdgesByTarget\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":102,"value":"colorEdgesByTarget = function (element, data) {\n edgeColorizer (element, data);\n var m = (data.target.data.name).split (\"_\");\n element.style(\"stroke-width\", 5);\n if (m.length > 10) {\n element.style (\"stroke\", color_scale(_.lowerCase(m[9])));\n }\n };","pinned":false,"hidden":false,"mode":"ojs"},{"id":105,"value":"colorNodesByName = function (element, data) {\n nodeColorizer (element, data);\n var m = (data.data.name).split (\"_\");\n if (m.length > 10) {\n element.style (\"stroke\", color_scale(_.lowerCase(m[9])));\n }\n };","pinned":false,"hidden":false,"mode":"ojs"},{"id":107,"value":"edgeColorizer = function(element, data) {\n //console.log (data[current_selection_name]);\ntry {\n var count_class = 0;\n\n selection_set.forEach (function (d,i) { if (data[d]) {count_class ++; element.style (\"stroke\", color_scale(i), 'important');}});\n\n if (count_class > 1) {\n element.classed (\"branch-multiple\", true);\n } else\n if (count_class == 0) {\n element.style (\"stroke\", null)\n .classed (\"branch-multiple\", false);\n }\n}\ncatch (e) {\n}\n\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":109,"value":"nodeColorizer = function (element, data) {\ntry{\n var count_class = 0;\n\n selection_set.forEach (function (d,i) { if (data[d]) {count_class ++; element.style (\"fill\", color_scale(i), 'important');}});\n\n if (count_class > 1) {\n\n } else {\n if (count_class == 0) {\n element.style (\"fill\", null);\n }\n }\n}\ncatch (e) {\n\n}\n\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":111,"value":"selection_set = \"\";","pinned":false,"hidden":false,"mode":"ojs"},{"id":113,"value":"color_scale = d3.scaleOrdinal(d3.schemeCategory10);","pinned":false,"hidden":false,"mode":"ojs"},{"id":125,"value":"tree = new phylotree.phylotree(nwk);","pinned":false,"hidden":false,"mode":"ojs"},{"id":292,"value":"##### Newick","pinned":false,"hidden":false,"mode":"md"},{"id":184,"value":"##### FEL - variables","pinned":false,"hidden":false,"mode":"md"},{"id":192,"value":"viewof pv = Inputs.text({label: html`p-value threshold`, value: \"0.1\", submit: \"Update\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":202,"value":"has_pasmt = fel_results_json.MLE[\"LRT\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":205,"value":"table_colors = ({\n 'Diversifying' : '#e3243b',\n 'Neutral' : '#444',\n 'Purifying' : 'green',\n 'Invariable' : '#CCC'\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":214,"value":"has_srv = _.chain (fel_results_json.MLE.content).some ((d)=>_.some (d,(dd)=>dd[0] > 0 && dd[0] != 1)).value()","pinned":false,"hidden":false,"mode":"ojs"},{"id":212,"value":"has_ci = fel_results_json [\"confidence interval\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":209,"value":"has_T = _.some (_.map (fel_results_json.MLE.content, (d)=>_.some(d, (dd)=>dd[5] > 0.)))","pinned":false,"hidden":false,"mode":"ojs"},{"id":198,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":227,"value":"tested_branch_count = d3.median (_.chain (fel_results_json.tested).map ().map ((d)=>_.map (d, (dd)=>_.filter (dd, ddd=>ddd == \"test\"))).map ((d)=>d.length).value())","pinned":false,"hidden":false,"mode":"ojs"},{"id":230,"value":"variable_site_count = d3.sum(_.chain (fel_results_json.MLE.content).map ((d)=>_.filter (d, (dd)=>dd[0]+dd[1] > 0)).map (d=>d.length).value())","pinned":false,"hidden":false,"mode":"ojs"},{"id":243,"value":"plot_specs = (\n{\n \"Site-level dN/dS estimates\" : {\n \"width\": 800, \"height\": 200, \n \"vconcat\" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {\n return dNdS_with_ci (fig1data, d, 70)\n })}\n}\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":246,"value":"fig1data = {\n let in_set = new Set (_.map (table1, (d)=>d.codon));\n return _.filter (siteTableData, (x)=>in_set.has (x.codon));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":250,"value":"siteTableData = {\n let data = \"\";\n data = _.filter (sites_table[1], (x)=>table_filter.indexOf (x.class)>=0);\n \n //let newEntry = selected_sites_human.codon;\n //let newData = data.map(entry => ({ ...entry, ...newEntry }));\n \n let newValues = newValuesMapped;\n let newData = data.map((entry, index) => ({ ...entry, HumanMapped: newValues[index] }));\n\n \n return newData;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":737,"value":"newValuesMapped = {\n let data = [];\n let site_table = sites_table;\n let counter = 1;\n function appendItem(item) {\n data.push(item);\n }\n \n site_table[1].forEach(site_datum => { \n if (protein_data[\"derive\"][site_datum.codon].v !== \"-\") {\n if (protein_data[\"derive\"][site_datum.codon].n === HumanRefSeq) {\n appendItem(counter);\n counter++\n }\n } else {\n appendItem(\"\");\n }\n });\n \n\n return data;\n};","pinned":true,"hidden":false,"mode":"ojs"},{"id":256,"value":"dyn_range_cap = 10","pinned":false,"hidden":false,"mode":"ojs"},{"id":254,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":47,"value":"### NOTE: Recombination detection (Found no breakpoints)","pinned":false,"hidden":false,"mode":"md"},{"id":323,"value":"##### Protein Structure","pinned":false,"hidden":false,"mode":"md"},{"id":328,"value":"color_hash = {\n return {\n blue: 0x00FFFF,\n red: 0xcf5b5b,\n yellow: 0xFFFF00,\n green: 0x5fcf5b,\n pink: 0xFFC0CB,\n gray: 0x808080\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":341,"value":"three_to_one = {\n return {\n \"ALA\": \"A\",\n \"CYS\": \"C\",\n \"ASP\": \"D\",\n \"GLU\": \"E\",\n \"PHE\": \"F\",\n \"GLY\": \"G\",\n \"HIS\": \"H\",\n \"ILE\": \"I\",\n \"LYS\": \"K\",\n \"LEU\": \"L\",\n \"MET\": \"M\",\n \"ASN\": \"N\",\n \"PRO\": \"P\",\n \"GLN\": \"Q\",\n \"ARG\": \"R\",\n \"SER\": \"S\",\n \"THR\": \"T\",\n \"VAL\": \"V\",\n \"TRP\": \"W\",\n \"TYR\": \"Y\"\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":344,"value":"stage.viewerControls.getOrientation()","pinned":false,"hidden":false,"mode":"ojs"},{"id":346,"value":"stage = {\n let s = new NGL.Stage( \"viewport\" ,{ backgroundColor: \"white\" })\n s.mouseControls.remove(\"hoverPick\");\n s.signals.hovered.add(function (pickingProxy) {\n if (pickingProxy && (pickingProxy.atom || pickingProxy.bond)){\n var atom = pickingProxy.atom || pickingProxy.closestBondAtom;\n var mp = pickingProxy.mouse.position;\n let residue = atom.residue;\n tooltip.innerText = residue.index + ' - index\\n' + \n residue.resno + ' - resno\\n' + \n residue.resname + ' - long resname\\n' +\n three_to_one[residue.resname] + ' - short resname\\n' +\n residue.chainname + ' - chain name\\n' +\n (selected_sites[residue.resno] ? selected_sites[residue.resno].annotation : '');\n tooltip.style.bottom = window.innerHeight - mp.y + 3 + \"px\";\n tooltip.style.left = mp.x + 3 + \"px\";\n tooltip.style.display = \"block\";\n } else {\n tooltip.style.display = \"none\";\n }\n })\n return s\n};","pinned":false,"hidden":false,"mode":"ojs"},{"id":360,"value":"pdbText = PDB.text()","pinned":false,"hidden":false,"mode":"ojs"},{"id":379,"value":"schemeId = {\n return NGL.ColormakerRegistry.addScheme(function (params) {\n this.atomColor = function (atom) {\n if(selected_sites[atom.residue.resno]) {\n \n console.log(\"# Site: \" + atom.residue.resno)\n \n let site_color = \"gray\"; // gray for Neutral\n \n console.log(\"Class: \" + selected_sites[atom.residue.resno].class)\n \n if (selected_sites[atom.residue.resno].class == \"Purifying\") {\n site_color = \"red\";\n }\n\n if (selected_sites[atom.residue.resno].class == \"Diversifying\") {\n site_color = \"blue\";\n }\n\n if (selected_sites[atom.residue.resno].class == \"Invariable\") {\n site_color = \"gray\";\n }\n\n return color_hash[site_color];\n //return color_hash[selected_sites[atom.residue.resno].color];\n }\n else {\n // gray\n //return 0x808080;\n return color_hash[\"gray\"]\n }\n }\n })\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":350,"value":"tooltip = {\n let tt = document.createElement(\"div\");\n document.body.appendChild(tt);\n Object.assign(tt.style, {\n display: \"none\",\n position: \"absolute\",\n zIndex: 10,\n pointerEvents: \"none\",\n backgroundColor: \"rgba(1, 1, 1, 1)\",\n color: \"lightgrey\",\n padding: \"0.5em\",\n fontFamily: \"sans-serif\"\n });\n return tt;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":382,"value":"component = {\n let str;\n stage.removeAllComponents();\n var stringBlob = new Blob( [ pdbText ], { type: 'text/plain'} );\n stage.loadFile(stringBlob, { ext: \"pdb\" }).then(d=>{\n d.addRepresentation(\"rope\", {color: schemeId})\n //d.addRepresentation(\"hyperball\", { sele: selection_string , color: schemeId});\n d.addRepresentation(\"licorice\", { sele: selection_string , color: schemeId});\n\n stage.viewerControls.orient([17.550527072492336, 77.73142323575034, -35.35459236339453, 0, 76.29412120673659, 1.939498913978275, 42.137723873208586, 0, 38.35789916992212, -39.42342635625532, -67.63584394065052, 0, 3.367788314819336, 0.33123016357421875, 3.3393964767456055, 1]);\n stage.setSpin(false);\n str = d;\n })\n return str;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":348,"value":"selected_sites = {\n let site_table = sites_table;\n let site_hash = {};\n \n site_table[1].forEach(site_datum => site_hash[+site_datum.codon] = site_datum);\n return site_hash;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":648,"value":"selected_sites_human = {\n let site_table = sites_table;\n let site_hash = {};\n let counter = 1;\n \n site_table[1].forEach(site_datum => {\n if (protein_data[\"derive\"][site_datum.codon].v !== \"-\") {\n site_hash[counter] = site_datum\n counter++\n }\n });\n return site_hash;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":386,"value":"selection_string = _.map (\n Object.keys(selected_sites_human), (i, s) => ['A'].map(chain => i + ':' + chain).join(\" OR \")\n).join (\" OR \")","pinned":false,"hidden":false,"mode":"ojs"},{"id":67,"value":"## Imports","pinned":false,"hidden":false,"mode":"md"},{"id":69,"value":"a = require('https://stephenshank.com/bundle.js')","pinned":false,"hidden":false,"mode":"ojs"},{"id":96,"value":"stylesheet = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":384,"value":"html`

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":100,"value":"phylotree = require('phylotree@1.0.0')","pinned":false,"hidden":false,"mode":"ojs"},{"id":115,"value":"d3 = require('d3')","pinned":false,"hidden":false,"mode":"ojs"},{"id":233,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":241,"value":"vegalite = require (\"vega-embed\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":330,"value":"NGL = require(\"ngl@2.0.0-dev.38\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":332,"value":"html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":334,"value":"_ = require (\"lodash\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":336,"value":"ReactDOM = r('react-dom')","pinned":false,"hidden":false,"mode":"ojs"},{"id":393,"value":"React = r('react')","pinned":false,"hidden":false,"mode":"ojs"},{"id":338,"value":"r = require.alias({\n react: \"react@16/umd/react.production.min.js\",\n \"react-dom\": \"react-dom@16/umd/react-dom.production.min.js\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":389,"value":"bioseq = require(\"bioseq\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":486,"value":"## Supplemental Figs\n","pinned":false,"hidden":false,"mode":"md"},{"id":487,"value":"Plot.plot({\n inset: 8,\n grid: true,\n color: {\n legend: true,\n },\n marks: [\n Plot.dot(protein_data[\"raw\"], {x: \"header\", y: \"seq_len\"})\n ]\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":161,"value":"## Gene-wide selection","pinned":false,"hidden":false,"mode":"md"},{"id":514,"value":"tag = \"Gnathosomata\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":523,"value":"gene = \"BCL6\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":626,"value":"### Human-protein mapping","pinned":false,"hidden":false,"mode":"md"},{"id":629,"value":"HumanRefSeq = \"NM_001706_5_Homo_sapiens_BCL6_transcription_repres\";","pinned":true,"hidden":false,"mode":"ojs"},{"id":632,"value":"CopyProteinData = protein_data[\"derive\"];","pinned":true,"hidden":false,"mode":"ojs"},{"id":634,"value":"HumanSites = CopyProteinData.filter(entry => entry.n === HumanRefSeq && entry.v !== \"-\");","pinned":true,"hidden":false,"mode":"ojs"},{"id":779,"value":"//console.log(HumanSites.length);","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}