{ "metadata": { "name": "", "signature": "sha256:813ccbb0d395bff71eecd0d0f92b28a4049bc7a014053bc72064b5243d792519" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Determing CpG ratio per Gene Function" ] }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Step 1: Obtain GO slim information (gene function)" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#from 00 notebook\n", "!head ../data/Piura_v1_GOslim.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Column1,GOSlim_bin\r", "\r\n", "PiuraChilensis_v1_contig_7798,protein metabolism\r", "\r\n", "PiuraChilensis_v1_contig_3751,other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_11851,other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_12118,protein metabolism\r", "\r\n", "PiuraChilensis_v1_contig_12118,transport\r", "\r\n", "PiuraChilensis_v1_contig_12118,cell organization and biogenesis\r", "\r\n", "PiuraChilensis_v1_contig_12118,signal transduction\r", "\r\n", "PiuraChilensis_v1_contig_12118,developmental processes\r", "\r\n", "PiuraChilensis_v1_contig_3855,developmental processes\r", "\r\n" ] } ], "prompt_number": 1 }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Step 2: Determine CpG ratio for each gene " ] }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Step 2a: Convert fasta to tab-delimited" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!perl -e '$count=0; $len=0; while(<>) {s/\\r?\\n//; s/\\t/ /g; if (s/^>//) { if ($. != 1) {print \"\\n\"} s/ |$/\\t/; $count++; $_ .= \"\\t\";} else {s/ //g; $len += length($_)} print $_;} print \"\\n\"; warn \"\\nConverted $count FASTA records in $. lines to tabular format\\nTotal sequence length: $len\\n\\n\";' \\\n", "../data/Piura_v1_contigs.fa > ../data/Piura_v1_contigs.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\r\n", "Converted 15022 FASTA records in 384556 lines to tabular format\r\n", "Total sequence length: 21729367\r\n", "\r\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "!head -1 ../data/Piura_v1_contigs.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "PiuraChilensis_v1_contig_1\t\tATTTACAATACGAAGTAAAATAGATAACGTGAAAATAATCTTGGTGCTGGATGATCGATCAAGTTCACCAATATTTTATTGTAAAAAATCATTCTAAACAGCATGAAATCGTGTACAATGTATAAACAAGCAAATATATAACACTAAAGCAAGAGGGCGTAAGTGGGGGGGTGGGTGAGAGTAAAAAATTCAAACATGTCAAATACCCCGGCGTTAGCCTTAAAAGCACCATGGACTTCTGCCTTCAATAAGCATAAAATTAAAACACCTAATACACAATGAATATACAGATAAAACAGATTTATGAATAGTTGGTGTTACATCTTTTACAGCCATAAGCCTTCATTTTGCTTCCAAACGTATAAAATCTGACTTGGAACAATATACAGCCATGAGATATGACACAGCGAGCACTACAATATATATTTATCTTGTACTATACAGCCTGTACAAGAAAATTCTGGAATTGTCTTCACAAGAGACAGAAAAATAGTTGCAATGTGAATGCTAGTCTACTATTTGATCACAATTGGATAGAAAAGTACAGCACATAAATGTTGGTGATACCTTAAAGAAAAGTGCAACAATATCAAAGGAATTAGTACCAGCATGCATTAGAAAAGTAAAAGTCTTGCTTATTACACAAAGCTGACTATATGATGTTCACCGCTTCTGGTGTGCAAAGAATTAAAAACAATGCAATTTCGGTCAGTTTTAACAAGGAATTAACAATTCATAGGAAAAATACAAGCATATGGTCTCAGGCCAATTGCTAGGACATAAAAAAAGCCTGCATATCACGAAAAGCCAAGTGCATGCATCGTATCCTGAAGACACCTTGATATTAACATGTAAGAAATTTAGCTTGCCACATTTCCATATTCCATAATTTCATTTTGAACACCGTGCCAGCAAATTCATCTGATATAAACACACAGGCAACTAATTTGGACACTTTCTAACTAGGTAGTTCAGAAAATACAGCTTTCAACAGGTACACATTTCTATAATAATAATAATAGCAAATGTCAGTGTGGCAGTTTTTGGACAAGTCCCTTTCAGAGGCCAAAATATCTATTTTGTATTAATTAATTATCCATTTTTTGGACTATACGGCTGTATCAAAACCATGGGTAACTGGGACTTGCCTAGCTTTGGGGGTGGAGTC\r\n" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "#temp replace name so c or g will not confound\n", "!sed 's/PiuraChilensis_v1_contig/999999/g' <../data/Piura_v1_contigs.tab> ../data/Piura_v1-99_contigs.tab\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 24 }, { "cell_type": "code", "collapsed": false, "input": [ "!head -1 ../data/Piura_v1-99_contigs.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "999999_1\t\tATTTACAATACGAAGTAAAATAGATAACGTGAAAATAATCTTGGTGCTGGATGATCGATCAAGTTCACCAATATTTTATTGTAAAAAATCATTCTAAACAGCATGAAATCGTGTACAATGTATAAACAAGCAAATATATAACACTAAAGCAAGAGGGCGTAAGTGGGGGGGTGGGTGAGAGTAAAAAATTCAAACATGTCAAATACCCCGGCGTTAGCCTTAAAAGCACCATGGACTTCTGCCTTCAATAAGCATAAAATTAAAACACCTAATACACAATGAATATACAGATAAAACAGATTTATGAATAGTTGGTGTTACATCTTTTACAGCCATAAGCCTTCATTTTGCTTCCAAACGTATAAAATCTGACTTGGAACAATATACAGCCATGAGATATGACACAGCGAGCACTACAATATATATTTATCTTGTACTATACAGCCTGTACAAGAAAATTCTGGAATTGTCTTCACAAGAGACAGAAAAATAGTTGCAATGTGAATGCTAGTCTACTATTTGATCACAATTGGATAGAAAAGTACAGCACATAAATGTTGGTGATACCTTAAAGAAAAGTGCAACAATATCAAAGGAATTAGTACCAGCATGCATTAGAAAAGTAAAAGTCTTGCTTATTACACAAAGCTGACTATATGATGTTCACCGCTTCTGGTGTGCAAAGAATTAAAAACAATGCAATTTCGGTCAGTTTTAACAAGGAATTAACAATTCATAGGAAAAATACAAGCATATGGTCTCAGGCCAATTGCTAGGACATAAAAAAAGCCTGCATATCACGAAAAGCCAAGTGCATGCATCGTATCCTGAAGACACCTTGATATTAACATGTAAGAAATTTAGCTTGCCACATTTCCATATTCCATAATTTCATTTTGAACACCGTGCCAGCAAATTCATCTGATATAAACACACAGGCAACTAATTTGGACACTTTCTAACTAGGTAGTTCAGAAAATACAGCTTTCAACAGGTACACATTTCTATAATAATAATAATAGCAAATGTCAGTGTGGCAGTTTTTGGACAAGTCCCTTTCAGAGGCCAAAATATCTATTTTGTATTAATTAATTATCCATTTTTTGGACTATACGGCTGTATCAAAACCATGGGTAACTGGGACTTGCCTAGCTTTGGGGGTGGAGTC\r\n" ] } ], "prompt_number": 25 }, { "cell_type": "code", "collapsed": false, "input": [ "#add column with length of sequence\n", "!perl -e '$col = 2;' -e 'while (<>) { s/\\r?\\n//; @F = split /\\t/, $_; $len = length($F[$col]); print \"$_\\t$len\\n\" } warn \"\\nAdded column with length of column $col for $. lines.\\n\\n\";' \\\n", "../data/Piura_v1-99_contigs.tab > ../data/Piura_v1-99-l_contigs.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\r\n", "Added column with length of column 2 for 15022 lines.\r\n", "\r\n" ] } ], "prompt_number": 26 }, { "cell_type": "code", "collapsed": false, "input": [ "!head -1 ../data/Piura_v1-99-l_contigs.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "999999_1\t\tATTTACAATACGAAGTAAAATAGATAACGTGAAAATAATCTTGGTGCTGGATGATCGATCAAGTTCACCAATATTTTATTGTAAAAAATCATTCTAAACAGCATGAAATCGTGTACAATGTATAAACAAGCAAATATATAACACTAAAGCAAGAGGGCGTAAGTGGGGGGGTGGGTGAGAGTAAAAAATTCAAACATGTCAAATACCCCGGCGTTAGCCTTAAAAGCACCATGGACTTCTGCCTTCAATAAGCATAAAATTAAAACACCTAATACACAATGAATATACAGATAAAACAGATTTATGAATAGTTGGTGTTACATCTTTTACAGCCATAAGCCTTCATTTTGCTTCCAAACGTATAAAATCTGACTTGGAACAATATACAGCCATGAGATATGACACAGCGAGCACTACAATATATATTTATCTTGTACTATACAGCCTGTACAAGAAAATTCTGGAATTGTCTTCACAAGAGACAGAAAAATAGTTGCAATGTGAATGCTAGTCTACTATTTGATCACAATTGGATAGAAAAGTACAGCACATAAATGTTGGTGATACCTTAAAGAAAAGTGCAACAATATCAAAGGAATTAGTACCAGCATGCATTAGAAAAGTAAAAGTCTTGCTTATTACACAAAGCTGACTATATGATGTTCACCGCTTCTGGTGTGCAAAGAATTAAAAACAATGCAATTTCGGTCAGTTTTAACAAGGAATTAACAATTCATAGGAAAAATACAAGCATATGGTCTCAGGCCAATTGCTAGGACATAAAAAAAGCCTGCATATCACGAAAAGCCAAGTGCATGCATCGTATCCTGAAGACACCTTGATATTAACATGTAAGAAATTTAGCTTGCCACATTTCCATATTCCATAATTTCATTTTGAACACCGTGCCAGCAAATTCATCTGATATAAACACACAGGCAACTAATTTGGACACTTTCTAACTAGGTAGTTCAGAAAATACAGCTTTCAACAGGTACACATTTCTATAATAATAATAATAGCAAATGTCAGTGTGGCAGTTTTTGGACAAGTCCCTTTCAGAGGCCAAAATATCTATTTTGTATTAATTAATTATCCATTTTTTGGACTATACGGCTGTATCAAAACCATGGGTAACTGGGACTTGCCTAGCTTTGGGGGTGGAGTC\t1168\r\n" ] } ], "prompt_number": 27 }, { "cell_type": "code", "collapsed": false, "input": [ "!awk -F\\CG '{print NF-1}' ../data/Piura_v1-99-l_contigs.tab > ../data/Piura_v1-99-l_contigs__CG.tab\n", "!awk -F\\C '{print NF-1}' ../data/Piura_v1-99-l_contigs.tab > ../data/Piura_v1-99-l_contigs__C.tab\n", "!awk -F\\G '{print NF-1}' ../data/Piura_v1-99-l_contigs.tab > ../data/Piura_v1-99-l_contigs__G.tab\n", "!paste ../data/Piura_v1-99-l_contigs.tab \\\n", "../data/Piura_v1-99-l_contigs__CG.tab \\\n", "../data/Piura_v1-99-l_contigs__C.tab \\\n", "../data/Piura_v1-99-l_contigs__G.tab \\\n", "> ../data/Piura_v1-99-l_contigs__C-G.tab\n", "!head -1 ../data/Piura_v1-99-l_contigs__C-G.tab\n" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "999999_1\t\tATTTACAATACGAAGTAAAATAGATAACGTGAAAATAATCTTGGTGCTGGATGATCGATCAAGTTCACCAATATTTTATTGTAAAAAATCATTCTAAACAGCATGAAATCGTGTACAATGTATAAACAAGCAAATATATAACACTAAAGCAAGAGGGCGTAAGTGGGGGGGTGGGTGAGAGTAAAAAATTCAAACATGTCAAATACCCCGGCGTTAGCCTTAAAAGCACCATGGACTTCTGCCTTCAATAAGCATAAAATTAAAACACCTAATACACAATGAATATACAGATAAAACAGATTTATGAATAGTTGGTGTTACATCTTTTACAGCCATAAGCCTTCATTTTGCTTCCAAACGTATAAAATCTGACTTGGAACAATATACAGCCATGAGATATGACACAGCGAGCACTACAATATATATTTATCTTGTACTATACAGCCTGTACAAGAAAATTCTGGAATTGTCTTCACAAGAGACAGAAAAATAGTTGCAATGTGAATGCTAGTCTACTATTTGATCACAATTGGATAGAAAAGTACAGCACATAAATGTTGGTGATACCTTAAAGAAAAGTGCAACAATATCAAAGGAATTAGTACCAGCATGCATTAGAAAAGTAAAAGTCTTGCTTATTACACAAAGCTGACTATATGATGTTCACCGCTTCTGGTGTGCAAAGAATTAAAAACAATGCAATTTCGGTCAGTTTTAACAAGGAATTAACAATTCATAGGAAAAATACAAGCATATGGTCTCAGGCCAATTGCTAGGACATAAAAAAAGCCTGCATATCACGAAAAGCCAAGTGCATGCATCGTATCCTGAAGACACCTTGATATTAACATGTAAGAAATTTAGCTTGCCACATTTCCATATTCCATAATTTCATTTTGAACACCGTGCCAGCAAATTCATCTGATATAAACACACAGGCAACTAATTTGGACACTTTCTAACTAGGTAGTTCAGAAAATACAGCTTTCAACAGGTACACATTTCTATAATAATAATAATAGCAAATGTCAGTGTGGCAGTTTTTGGACAAGTCCCTTTCAGAGGCCAAAATATCTATTTTGTATTAATTAATTATCCATTTTTTGGACTATACGGCTGTATCAAAACCATGGGTAACTGGGACTTGCCTAGCTTTGGGGGTGGAGTC\t1168\t15\t203\t202\r\n" ] } ], "prompt_number": 28 }, { "cell_type": "markdown", "metadata": {}, "source": [ "![](http://eagle.fish.washington.edu/cnidarian/skitch/BMC_Genomics___Full_text___DNA_methylation_patterns_provide_insight_into_epigenetic_regulation_in_the_Pacific_oyster__Crassostrea_gigas__1A0683A5.png)" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!awk '{print $1, \"\\t\", (($4)/($5*$6))*(($3**2)/($3-1))}' \\\n", "../data/Piura_v1-99-l_contigs__C-G.tab \\\n", "| sed 's/999999/PiuraChilensis_v1_contig/g' > ../data/Piura_v1_CpG.tab\n", "!head ../data/Piura_v1_CpG.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "PiuraChilensis_v1_contig_1 \t 0.427621\r\n", "PiuraChilensis_v1_contig_2 \t 0.600881\r\n", "PiuraChilensis_v1_contig_3 \t 0.750945\r\n", "PiuraChilensis_v1_contig_4 \t 1.22757\r\n", "PiuraChilensis_v1_contig_5 \t 0.848172\r\n", "PiuraChilensis_v1_contig_6 \t 0.86292\r\n", "PiuraChilensis_v1_contig_7 \t 0.718647\r\n", "PiuraChilensis_v1_contig_8 \t 1.03897\r\n", "PiuraChilensis_v1_contig_9 \t 0.97759\r\n", "PiuraChilensis_v1_contig_10 \t 0.817866\r\n" ] } ], "prompt_number": 29 }, { "cell_type": "code", "collapsed": false, "input": [ "!ls ../data" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Cgigas-HS-count.txt\r\n", "PiuraC_Coq_Trinity.fasta\r\n", "PiuraC_Val_Trinity.fasta\r\n", "PiuraC_Val_Trinity_2ndhalf.fasta\r\n", "PiuraC_Val_Trinity_uniprot_sprot.tab\r\n", "PiuraC_Val_Trinity_uniprot_sprot_2ndhalf.tab\r\n", "Piura_counts.txt\r\n", "Piura_v1-99-l_contigs.tab\r\n", "Piura_v1-99-l_contigs__C-G.tab\r\n", "Piura_v1-99-l_contigs__C.tab\r\n", "Piura_v1-99-l_contigs__CG.tab\r\n", "Piura_v1-99-l_contigs__G.tab\r\n", "Piura_v1-99_contigs.tab\r\n", "Piura_v1_CpG.tab\r\n", "Piura_v1_GOslim.csv\r\n", "Piura_v1_contigs.fa\r\n", "Piura_v1_contigs.tab\r\n", "Piura_v1_uniprot_sprot.tab\r\n", "Piura_v1_uniprot_sprot_sql.tab\r\n", "\u001b[34mwd\u001b[m\u001b[m\r\n" ] } ], "prompt_number": 30 }, { "cell_type": "code", "collapsed": false, "input": [ "!rm ../data/Piura_v1-99*" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 31 }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Join with GO Slim" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head ../data/Piura_v1_GOslim.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Column1,GOSlim_bin\r", "\r\n", "PiuraChilensis_v1_contig_7798,protein metabolism\r", "\r\n", "PiuraChilensis_v1_contig_3751,other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_11851,other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_12118,protein metabolism\r", "\r\n", "PiuraChilensis_v1_contig_12118,transport\r", "\r\n", "PiuraChilensis_v1_contig_12118,cell organization and biogenesis\r", "\r\n", "PiuraChilensis_v1_contig_12118,signal transduction\r", "\r\n", "PiuraChilensis_v1_contig_12118,developmental processes\r", "\r\n", "PiuraChilensis_v1_contig_3855,developmental processes\r", "\r\n" ] } ], "prompt_number": 32 }, { "cell_type": "code", "collapsed": false, "input": [ "!tr ',' \"\\t\" <../data/Piura_v1_GOslim.csv> ../data/Piura_v1_GOslim.tab" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 36 }, { "cell_type": "code", "collapsed": false, "input": [ "!sort ../data/Piura_v1_GOslim.tab | tail -n +2 > ../data/Piura_v1_GOslim.sorted\n", "!awk -F $'\\t' '{print $1, \"\\t\", $2}' ../data/Piura_v1_GOslim.sorted > ../data/Piura_v1_GOslim.sortedtab\n", "!head ../data/Piura_v1_GOslim.sortedtab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "PiuraChilensis_v1_contig_100 \t death\r", "\r\n", "PiuraChilensis_v1_contig_100 \t other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_100 \t signal transduction\r", "\r\n", "PiuraChilensis_v1_contig_100 \t stress response\r", "\r\n", "PiuraChilensis_v1_contig_100 \t transport\r", "\r\n", "PiuraChilensis_v1_contig_1000 \t other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_1000 \t other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_10003 \t RNA metabolism\r", "\r\n", "PiuraChilensis_v1_contig_10003 \t other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_10003 \t protein metabolism\r", "\r\n" ] } ], "prompt_number": 151 }, { "cell_type": "code", "collapsed": false, "input": [ "!sort ../data/Piura_v1_CpG.tab > ../data/Piura_v1_CpG.sorted\n", "!awk -F $'\\t' '{print $1, \"\\t\", $2}' ../data/Piura_v1_CpG.sorted > ../data/Piura_v1_CpG.sortedtab\n", "!head ../data/Piura_v1_CpG.sortedtab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "PiuraChilensis_v1_contig_1 \t 0.427621\r\n", "PiuraChilensis_v1_contig_10 \t 0.817866\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473\r\n", "PiuraChilensis_v1_contig_1000 \t 0.792597\r\n", "PiuraChilensis_v1_contig_10000 \t 0.947282\r\n", "PiuraChilensis_v1_contig_10001 \t 0.691634\r\n", "PiuraChilensis_v1_contig_10002 \t 0.936689\r\n", "PiuraChilensis_v1_contig_10003 \t 0.371111\r\n", "PiuraChilensis_v1_contig_10004 \t 0.89432\r\n", "PiuraChilensis_v1_contig_10005 \t 0.873687\r\n" ] } ], "prompt_number": 152 }, { "cell_type": "markdown", "metadata": {}, "source": [ "```\n", "SELECT * FROM [sr320@washington.edu].[Piura_v1_CpG.sorted]cpg\n", " left join\n", "[sr320@washington.edu].[Piura_v1_GOslim.sorted\u200b]go\n", " on\n", " cpg.Column1=go.Column1\n", "``` \n" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head ../data/Piura_v1_CpG-GOslim.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Column1,Column2,Column1,Column2\r", "\r\n", "PiuraChilensis_v1_contig_1,0.427621,,\r", "\r\n", "PiuraChilensis_v1_contig_10,0.817866,,\r", "\r\n", "PiuraChilensis_v1_contig_100,0.914473,PiuraChilensis_v1_contig_100,death\r", "\r\n", "PiuraChilensis_v1_contig_100,0.914473,PiuraChilensis_v1_contig_100,other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_100,0.914473,PiuraChilensis_v1_contig_100,signal transduction\r", "\r\n", "PiuraChilensis_v1_contig_100,0.914473,PiuraChilensis_v1_contig_100,stress response\r", "\r\n", "PiuraChilensis_v1_contig_100,0.914473,PiuraChilensis_v1_contig_100,transport\r", "\r\n", "PiuraChilensis_v1_contig_1000,0.792597,PiuraChilensis_v1_contig_1000,other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_1000,0.792597,PiuraChilensis_v1_contig_1000,other metabolic processes\r", "\r\n" ] } ], "prompt_number": 167 }, { "cell_type": "code", "collapsed": false, "input": [ "!tr ',' \"\\t\" <../data/Piura_v1_CpG-GOslim.csv> ../data/Piura_v1_CpG-GOslim.tab" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 168 }, { "cell_type": "code", "collapsed": false, "input": [ "!head ../data/Piura_v1_CpG-GOslim.tab" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Column1\tColumn2\tColumn1\tColumn2\r", "\r\n", "PiuraChilensis_v1_contig_1\t0.427621\t\t\r", "\r\n", "PiuraChilensis_v1_contig_10\t0.817866\t\t\r", "\r\n", "PiuraChilensis_v1_contig_100\t0.914473\tPiuraChilensis_v1_contig_100\tdeath\r", "\r\n", "PiuraChilensis_v1_contig_100\t0.914473\tPiuraChilensis_v1_contig_100\tother biological processes\r", "\r\n", "PiuraChilensis_v1_contig_100\t0.914473\tPiuraChilensis_v1_contig_100\tsignal transduction\r", "\r\n", "PiuraChilensis_v1_contig_100\t0.914473\tPiuraChilensis_v1_contig_100\tstress response\r", "\r\n", "PiuraChilensis_v1_contig_100\t0.914473\tPiuraChilensis_v1_contig_100\ttransport\r", "\r\n", "PiuraChilensis_v1_contig_1000\t0.792597\tPiuraChilensis_v1_contig_1000\tother biological processes\r", "\r\n", "PiuraChilensis_v1_contig_1000\t0.792597\tPiuraChilensis_v1_contig_1000\tother metabolic processes\r", "\r\n" ] } ], "prompt_number": 169 }, { "cell_type": "code", "collapsed": false, "input": [ "!awk -F $'\\t' '{print $1, \"\\t\", $2,\" \\t\", $4}' ../data/Piura_v1_CpG-GOslim.tab | tail -n +2 > ../data/Piura_v1_CpG-slim.tab\n", "!head ../data/Piura_v1_CpG-slim.tab\n" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "PiuraChilensis_v1_contig_1 \t 0.427621 \t \r", "\r\n", "PiuraChilensis_v1_contig_10 \t 0.817866 \t \r", "\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473 \t death\r", "\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473 \t other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473 \t signal transduction\r", "\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473 \t stress response\r", "\r\n", "PiuraChilensis_v1_contig_100 \t 0.914473 \t transport\r", "\r\n", "PiuraChilensis_v1_contig_1000 \t 0.792597 \t other biological processes\r", "\r\n", "PiuraChilensis_v1_contig_1000 \t 0.792597 \t other metabolic processes\r", "\r\n", "PiuraChilensis_v1_contig_10000 \t 0.947282 \t \r", "\r\n" ] } ], "prompt_number": 170 }, { "cell_type": "code", "collapsed": false, "input": [ "import pandas as pd" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 171 }, { "cell_type": "code", "collapsed": false, "input": [ "Piura = pd.read_table('../data/Piura_v1_CpG-slim.tab', header=None)\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 172 }, { "cell_type": "code", "collapsed": false, "input": [ "Piura" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "
\n",
        "<class 'pandas.core.frame.DataFrame'>\n",
        "Int64Index: 27267 entries, 0 to 27266\n",
        "Data columns (total 3 columns):\n",
        "0    27267  non-null values\n",
        "1    27267  non-null values\n",
        "2    27267  non-null values\n",
        "dtypes: float64(1), object(2)\n",
        "
" ], "metadata": {}, "output_type": "pyout", "prompt_number": 173, "text": [ "\n", "Int64Index: 27267 entries, 0 to 27266\n", "Data columns (total 3 columns):\n", "0 27267 non-null values\n", "1 27267 non-null values\n", "2 27267 non-null values\n", "dtypes: float64(1), object(2)" ] } ], "prompt_number": 173 }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 174 }, { "cell_type": "code", "collapsed": false, "input": [ "import matplotlib.pyplot as plt " ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 177 }, { "cell_type": "code", "collapsed": false, "input": [ "Piura.groupby(2)[1].mean().plot(kind='barh', color=list('myb'))\n", "plt.axis([0.7, 0.9, 0, 15])" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 186, "text": [ "[0.7, 0.9, 0, 15]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAhwAAAD5CAYAAABs6OxkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdYVMfaAPB3QRClCBqN5UZRUBTYCiwddlVQA0RExaBB\nsMYKsQavDVtCoqjYO6ixAZYEsaDIKlV6VxAU42dFQMoCUna+P/Dcu9kL2BZ3gfeX5zwye87Oec+4\nkdk5c+alEUIAIYQQQqgtKcg6AIQQQgh1fNjhQAghhFCbww4HQgghhNocdjgQQggh1Oaww4EQQgih\nNocdDoQQQgi1uS6yDgChL4lGo+Fz4Agh9AkIIbTPeT+OcKAP9ssvv/xb1jF8iJ07d/5UU1PTreUj\nCG5S29bLQQwdZZN2WwIQQjrttn79epnH0JE2acAOB/pgv/7666rmXieE0Mhn9nylpbGxUTEgIMC7\nurq6u6xj6RyKZB1AB1Ik6wA6lKKiIlmHgCRghwN9EB8fH7+amppubDY7zd3d/eTjx48H6enp5Xl4\neByn0+lZT548+WbBggX7TExMkgwNDbN9fX19qfdqa2sX+fr6+hoZGaUwGIzMvLw8PQCA27dv27LZ\n7DQ2m53G4XBSq6qq1AQCAc/GxuaOo6Pj5eHDh9+fP3/+fqozc+bMGTcGg5FJp9OzfHx8/Kj61dTU\nqpYvX76NxWKl//LLL/9+9uxZfz6fHzVq1KjIL95QCCGEmifrYRrc2s+mpqZWSf386NEjbQUFhca7\nd+9yqddKS0u1CCHQ0NCgyOPxorKysgwJIaCtrf1oz549CwkhsG/fvvmzZ88+TAgBJyenv+Li4swJ\nISAUCrs3NDQoRkVF8VRUVGoePXqk3djYqGBnZxcRGho68enTp/0HDhz4+PXr170aGhoUR44cGXnp\n0qXxhBCg0WiikJCQSVQc2traj0pKSno2dw0AQAAIblLbouQgho6ySbstgXRmUVFRsg6hQ3n3eYLP\n2XDSKPpkgwYNeszlchOp8rlz56YcPnx4TkNDQ5fnz5/3y83N1Tc0NMwGAHBxcbkAAMDhcFIvXLjg\nAgBgaWkZu2TJkh3Tpk075eLicmHAgAFPAQC4XG6itrZ2EQCAm5vbmZiYGCslJaV6Ho8n6NWrVwkA\nwLRp007duXPHZvz48X8qKio2Tpw48fyHR+4JANrvftYEABYA8N6VBe/+xPKHlanX5CWe9lzmSb1+\ngaCpzON1vjKPx5OreNpbWSAQQFBQEAAAaGtrg1R8bo8Ft86zSY5wGBoaZlHlhw8fDtbV1X3w5s2b\nHoQQ8PT0DDx+/Ph0Qv454pCUlGTM4/GiqPdlZ2cb/PbbbysHDRpUdP/+fb2oqCiera2tgNp/9OjR\nmUuWLNn+559/fjd9+vTj1OtHjhyZtWzZsm2ScUmeT3IDwBEO+f5W3pk3abclkM4MRzik693nCT5n\nwzkc6IMpKSnVNzQ0NDsqVlFRoaGqqirU0NCoePny5ddXr14d9776CgsLdQwMDHJWrlz5u4mJSRI1\ntyMxMZFbVFSkLRKJFIKDg12tra2juVxu4u3bt21LSkp6NTY2Kp49e/Z7W1vb283Vq66uXllRUaHx\neVeLEEJImrDDgT7Y3LlzDzEYjEx3d/eTNBqNiK9pwWQyM9hsdtrw4cPvT5s27ZSVlVVMS/VQ7wsI\nCPCm0+lZTCYzQ1lZuW7cuHFXAQBMTEySFi1atEdfXz93yJAhDydMmHCxb9++L/z8/Hz4fH4Ui8VK\nNzY2TnZycgoTr088zrFjx17DSaNfAk/WAXQgPFkH0KFQtwmQ/KA1jZQgJB8EAgHP399/WVhYmFNb\n1N/UOcHPPOoMaID/viNpodFoQHDhL9SRSI6cIHknkHUAHYhA1gF0KNQESCQ/cIQDdSrYmUGdhbq6\nFlRUlMo6DJkRCAR4W0WKpDHCgR0O1KnQaDSCn3mEEPo4eEsFIYQQQu0CLvyFOh0aTS7SviD0D2pq\n3aCyslrWYXQYeEtF/mCHo507fvy4h729fUS/fv2eyzqW9iIqStYRdBzp6QAslqyj6Bj4/BpZh4BQ\nm8JbKu1cUFCQ57Nnz/o3t08kErXZ329LC4ChzgU7G0he4eiG/MEORzvR2Nio6OnpGUSn07MYDEbm\nzp07fzp//vzE5ORk42nTpp3icDiptbW1Ktra2kU+Pj5+RkZGKSEhIZMjIiLsLSws4oyMjFJcXV2D\nhUKhKkBT9lcDA4McJpOZsXLlyt8BAEJCQibT6fQsFouV3twqngKBgGdtbR09fvz4Pw0NDbNFIpHC\nihUrtnK53EQmk5lx6NChuQAAz58/72djY3OHzWan0en0rNjYWEuApqyuS5cu3W5oaJg9evTom69f\nv/4KACA9PZ1lZmaWwGQyM1xcXC68efNGEwCAx+MJfHx8/ExNTe/q6enlxcTEWAEA5OTkGJiamt5l\ns9lpTCYzo7CwUAcA4I8//viBen3evHkH2rLDhRBC6OPgP8jtRHp6OuvZs2f9s7Ky6JmZmYyZM2ce\nmzhx4nljY+Pk06dPT01NTeWoqKjU0mg08tVXX71OSUkxGjVqVOSWLVtWR0ZGjkpJSTEyMjJK2b59\n+9LS0tKely5dcs7JyTHIyMhgrlmzZjMAwKZNm9ZGRETYp6ens1paeCstLY29a9cur/v37w8/cuTI\nbE1NzTeJiYncxMRE7uHDh+cUFRVpnzlzxm3s2LHX0tLS2BkZGUwmk5kBAFBdXd3dxMQkKTs729DW\n1vb2hg0b1gMATJ8+/cTWrVtXZGRkMOl0ehb1Oo1GI42NjYp379413blz50/U6wcOHJjn7e0dkJaW\nxk5JSTEaMGDA03v37o0IDg52jYuLs0hLS2MrKCiITp06Ne3L/O10Xunpso4AoebhOhzyB4fF2wkd\nHZ3Chw8fDvHy8trl4OAQbm9vH0Htk3xUacqUKecAABISEsxyc3P1LSws4gAA6urqlC0sLOJ69OhR\nrqKiUjtr1qyjjo6Olx0dHS8DNGVv9fDwOO7q6hpMZXeVxOVyEwcNGvQYACAiIsI+KyuLHhoaOgmg\nKZ9KQUGBromJSdLMmTOP1dfXKzk7O1+iOhwKCgoiKrYffvjhDxcXlwsVFRUa5eXlPaytraMBADw8\nPI5Pnjw5hDqfeJbZoqIibQAACwuLuC1btqz+v//7v3+5uLhc0NXVLaA6VcbGxskAADU1Nd369u37\nQgpNjxBCSAqww9FOaGpqvsnMzGRcu3Zt7IEDB+YFBwe7Hj16dBbA/y5mpaqqKqR+trOzu3H69Omp\nkvUlJiZyIyMjR4WGhk7as2fPosjIyFH79++fn5iYyA0PD3cwMjJKSUlJMerZs2dpS3UDAOzZs2eR\nnZ3dDcn6o6OjrS9fvuzo6ekZtHTp0u3u7u4nxfcTQmjNLcIl2Xnq2rXrWwAARUXFRmreiJub2xkz\nM7OEy5cvO3777bdXDh48+CNAU2fll19++XdLbUjx8wPo27fpZzU1AF3d/85FoL6xY/nDytRr8hJP\ney/LU3ry9l7G9PTyl54eF/5qJ0pKSnopKSnVa2hoVGRnZxtOnz79RGpqKue77777a+nSpdt5PJ4A\nAGDw4MGPqI5CcXFxb2Nj4+Rbt26N1NHRKRQKharPnj3r379//2dCoVC1T58+r8rLy3vo6OgUvn79\n+qvCwkIdHR2dQoCmkYwjR47MZjAYmVQMknlODh8+POfKlSvfhoSETO7SpUtDfn7+sH/961//9/r1\n668GDBjwVFFRsXHPnj2LHj58OGT79u1LFRQURGfOnHGbMmXKuc2bN68pLi7uHRAQ4M1isdL37Nmz\nyMrKKsbX19e3srJS3d/ffxmfz4/y9/dfxuFwUl+/fv2ViYlJ0qNHjwY/fPhwyJAhQx4CAKxYsWLr\nN99888TOzu7G+PHj/4yNjbXs3bt3cWlpac+qqiq1gQMH/i3ejjQajeBTKkge8fmAuU+Q3JLGwl84\nwtFOPH36dMCMGTMCqYmQfn5+PgAAnp6eQfPmzTvQvXv36ri4OAvx9/Tu3bs4KCjI083N7czbt2+7\nAgBs2bJltbq6euX48eP/rK2tVSGE0Hbs2LEEAGDlypW/P3jwYCghhDZ69Oib4p0NgP/NczJ79uwj\nRUVF2hwOJ5UQQuvTp8+rixcvThAIBLytW7euUFJSqldXV688ceLEdICm0ZHExETu5s2b13z99dcv\nz507NwWg6dHeefPmHaiuru6uo6NTGBgYOKO5NqDOHRwc7PrHH3/8oKSkVN+vX7/nq1ev3qKpqflm\n8+bNa+zt7SNEIpGCkpJS/b59+xZIdjiQdOFjsUhe4Toc8gdHONAXo66uXllZWakuyxhwhEO6sMMh\nPTjCIV3Y4ZAuXNoctSuYOK3jwc4GklfY2ZA/2OFAX0xFRYWGrGNACCEkG3hLBXUqOMqC5BXmUpEu\nvKUiXThpFKFPgJ1s6cF/1KUHF6pCHR2OcKBOhUajEfzMI4TQx8ERDoQ+AaanR0h+aaprQllFmazD\nQG0AOxyo04kCfC5WWtIhHViAj6pIA7ZlE34lXyr14O0++YNPqcjQ+vXrN0RGRo6Sdr3a2tpFpaWl\nPSVf/5Blv6XN09Mz6Pz58xM/5b2S8VpaWsZKJyqEEEJfGs7h6IDElzcXf72lhbeo+3Jt8QTHjBkz\nAp2cnMJaSgbXmrZYKIxGoxEc4UBIfvGBjxO75RAu/NVONDY2Knp6egbR6fQsBoORGRAQ4A3wz2//\nV65c+XbEiBH3jI2Nk728vHY5OTmFAQD4+vr6zpw58xifz4/S0dEp3L1792Kq3gkTJlw0NjZONjQ0\nzD58+PCc1mLw8fHxq6mp6cZms9Pc3d1PPn78eJCenl6eh4fHcTqdnvXkyZNvFixYsM/ExCTJ0NAw\n29fX15d6r7a2dpGvr6+vkZFRCoPByMzLy9MDALh9+7Ytm81OY7PZaRwOJ1UoFKoSQmiLFi3aM3z4\n8Pt2dnY3Xr161Ue8HmrkJTk52ZjP50cBAFRVVanNmDEjkMFgZDKZzIwLFy64rFq16lfxeAEA1NTU\nqgCaOkgrVqzYSrVncHCwK0BTrhcejyeYPHlyyIgRI+798MMPf0jhrw8hhJAU4ByOLyA9PZ317Nmz\n/llZWXSA/y6AReUmqa2tVZk3b96B6Oho60GDBj2eOnXqafHRhvz8/GFRUVH8iooKDT09vbwFCxbs\nU1RUbDx27NhMLS2tspqamm5cLjdx0qRJoVpaWs3OtvLz8/PZu3fvwrS0NDYAQFFRkXZBQYHuyZMn\n3blcbiJAU54VLS2tssbGRsXRo0ffzM7ONjQ0NMym0Wikd+/exSkpKUb79++fv23btuWHDx+e4+/v\nv2zfvn0LzM3N46urq7t37dr17cWLFyfk5+cPu3fv3ogXL1701dfXz501a9ZR6nqbi23Tpk1rtbS0\nyjIzMxkAAG/evNF0cXG5sGfPnkVUvOLvv3DhgktGRgYzMzOTUVxc3NvExCTJxsbmDtXWubm5+v36\n9XtuaWkZGxsba4m3YtoOzjuQHmxL6cI5HPIHRzi+AB0dncKHDx8O8fLy2nX9+vUx6urqldQ+Qgjt\n/v37w4cMGfJw0KBBjwGa0q+L3+ZwcHAIV1JSqu/Vq1dJnz59Xr18+fJrAAAq06q5uXn8kydPvnnw\n4MHQj4lr0KBBj6nOBgDAuXPnphgZGaVwOJzUnJwcg9zcXH1qH3VLhMPhpBYVFWkDNM2pWLJkyY7d\nu3cvLisr01JUVGyMjo62pjpM/fr1ez5y5Mhb74sjMjJy1MKFC/dSZU1NzTetHR8TE2NFnaNPnz6v\nbG1tbyclJZnQaDTC5XIT+/fv/4xGoxEWi5VOxSrOD/wg6N1/oRAK6ZD+n33p7/7D8oeVC6BAruLB\ncvsvixMIBP9YnwTLX64sEAjA09MTPD09QWzA+7PgCMcXoKmp+SYzM5Nx7dq1sQcOHJgXHBzsevTo\n0VnUfslv/pL3yZSVleuonxUVFRsbGhq6CAQCXmRk5KiEhAQzFRWVWj6fH1VbW6vyMXGpqqoKqZ8f\nPXo02N/ff1lycrJxjx49ymfMmBEoXl/Xrl3fip8fAODnn3/+zdHR8XJ4eLiDpaVl7PXr18c0Fz+l\nS5cuDVS2W8lYP+be4Lu1NGiSr4nHKRmrOB/wabFuyW+YWG69PAkmyVU87bnc3OiGPMUni7LkCMXH\nlJsb3fic+jpbmcfj/aO8YcMG+Fw4wvEFlJSU9GpoaOji4uJyYdOmTWslbxPo6enlPXz4cMjjx48H\nATSNNFC/QJv7RUwIoVVUVGhoaWmVqaio1N6/f394QkKC2fviUFJSqm/uFzBA020eVVVVoYaGRsXL\nly+/vnr16rj31VdYWKhjYGCQs3Llyt9NTEyS7t+/P9zGxubOuXPnpohEIoXnz5/3i4qK+s8zbtra\n2kXJycnGAADiT67Y2dnd2Lt370Kq/ObNG83W4rW2to6mzlFcXNz7zp07NlwuN/FzJzQhhBBqO9jh\n+AKePn06gM/nR1ETIH/99ddV4vtVVFRq9+3bt2Ds2LHXjI2NkzU0NCp69OhRDvDfeR7ix9NoNDJ2\n7NhrDQ0NXfT19XNXrVr1q7m5efz74pg7d+4hBoOR6e7uflKyXiaTmcFms9OGDx9+f9q0aaesrKxi\nWqqHel9AQIA3nU7PYjKZGcrKynXffvvtlQkTJlwcOnToA319/VwPD4/jFhYWcdT71q9fv8Hb2zvA\nxMQkqUuXLg1UPWvWrNlcVlamRafTs1gsVrpAIOBJxit+3gkTJlykJpiOGjUqcuvWrSv69OnzqqW2\nel+7oE8nOQSOPh22pXThUvHyBx+LlRNCoVCVusWxcOHCvcOGDcv39vYOkHVcHQ0+Fitd6YATHaUF\n27KJtB6LxUmj0oWPxXYghw8fnsNms9MMDAxyKioqNH788ceDso4JoffBX5DSg20pXdjZkD84woE6\nFbzFgpB8w1wq8gmTtyH0CbCTLT04bC092JbShe0pf/CWCkIIIYTaHN5SQZ0K3lJBHZG6uhZUVJS+\n/0CEPhFOGu0kdu7c+VNNTU239x03Z86cw/fu3RvxJWICADh+/LjH8+fP+73vOB6PJ0hJSTH60HqD\ngoI8Fy9evBsA4ODBgz+ePHnS/XPi/F8EN9w61FZZiXMekPzDDoecaGxsVGxpX0BAgHd1dXX399Vx\n+PDhOSNGjLgn3chaFhQU5Pns2bP+7zvuY0cVxI//8ccfD1LrcCB5JJB1AB2IQNYBdCi4Dof8wQ7H\nF6Cmpla1dOnS7YaGhtmjR4+++fr1668Amr75L1myZIeJiUlSQECAd2Rk5CgOh5PKYDAyZ82adbSu\nrk55165dXs+ePevP5/OjRo0aFQkAEBERYW9hYRFnZGSU4urqGiwUClWp+lJTUznUOdesWbOZyrUi\nnrWV4uvr6+vh4XHcxsbmjra2dtGFCxdcli9fvo3BYGSOGzfuKrXKZ0pKihGPxxMYGxsnjx079tqL\nFy/6hoaGTkpOTjaeNm3aKQ6Hk1pbW6uycePGdVwuN5FOp2dJPtZ78uRJdzabnUan07OSkpJMAABK\nS0t7Ojs7X2IymRnm5ubxVHI7yRj9/f2XAQDs2rXLy8DAIIfJZGZMnTr19MdcA0IIIRkjhODWxhuN\nRhOdPn3ajRACGzduXLto0aLdhBDg8XhRCxcu3EMIgZqaGpVvvvnm7wcPHugSQmD69OnHd+7c6U0I\nAW1t7UclJSU9CSFQXFz8lY2Nze3q6upuhBDw8/P7eePGjWup+lJSUjjUOS9fvuxACIGVK1f+tnnz\n5tWSca1fv97X2tr6TkNDg2JGRgajW7du1deuXRtDCIEJEyZcuHTp0vi6ujolc3PzuNevX/cihMDZ\ns2enzJw586jk+QghUFpaqkX97O7ufiIsLMyREAK2traCuXPnHiSEwJ07d6wNDQ2zCCGwaNGi3VTs\nt27d4rNYrDRCCAQGBnpSbeTr67ve399/KSEE+vfv/7Surk6JEALl5eUaH3oN4tcMAASA4IZbB9uA\nINSW3n3G4HM2HOH4AhQUFERTpkw5BwDwww8//BETE2NF7aNez8vL0xs8ePAjXV3dAgAADw+P43fu\n3LGRrCshIcEsNzdX38LCIo7NZqedOHFi+t9//z1Q8jhlZeU6BweHcAAAIyOjlOayptJoNDJu3Lir\nioqKjYaGhtkikUhhzJgx1wEA6HR6VlFRkXZ+fv6wnJwcg9GjR99ks9lpW7ZsWf306dMBVB1EbBLR\nrVu3RpqZmSUwGIzMW7dujaSyzdJoNOLm5nYGoCkPSkVFhUZ5eXmP2NhYS+p2CZ/PjyopKelVWVmp\n3lI7MhiMzKlTp54+derUNEVFxcYPvYaW/2YQQgh9KTjc/IURQmjicxTEM7a2dpw4Ozu7G6dPn57a\n2nmUlJTqqZ8VFBRELd1aoDLRKigoiJp7DyGEZmBgkBMXF2fR3PupGGtra1UWLly4NyUlxWjAgAFP\nN2zYsL617LUtJadr7pqpY8LDwx3u3LljExYW5rRly5bV1C2Y913D/57dEwC03/2sCQAsAOC9Kwve\n/YnlDyvvBGw/aZWpnz/l/e9K7+YtUOtPdOay+BwOeYinvZUFAgEEBQUBAIC2tjZIxecOkeD2YbdU\nzp49O4UQAps2bVrj5eUVQEjTLYnk5GQjQppuqQwcOPBxQUGBDiEEPDw8gnbt2rWYEAJ0Oj3z0aNH\n2oQQePXqVW/x46qqqlTz8/OHUvVRtzjU1NQqqfOHhIRM8vT0DJSMy9fXd/22bduWUWXx91D76urq\nlHR1dR/Ex8ebEUKgrq5OKScnR58QAk5OTn9FRUXxCCFQVlam+fXXX7+oqalRqaysVDMwMMjesGHD\nOkKabqnMmzdvPyEEoqOjrRgMRgYhBLy8vAI2bdq0hhACUVFRPA6Hk0LIP2+prF+/3nfbtm3LRCIR\njWqDuro6pf79+z998+ZND2p/a9cgfs0AeEtFuluUHMTQUbbPaUsg6J+ioqJkHUKH8u4zBp+z4QjH\nF6CqqipMTEzkbt68ec3XX3/98ty5c1OofdQ3ehUVldrAwMAZkydPDmloaOjC5XIT582bdwCgKWvq\n2LFjrw0YMOBpZGTkqKCgIE83N7czb9++7QoAsGXLltVDhw59IH5O8ZGC5rKotnSc5D4lJaX60NDQ\nSV5eXrvKy8t7NDQ0dFmyZMkOfX39XE9Pz6B58+Yd6N69e3VcXJzFnDlzDhsaGmb37dv3hamp6V3x\nelRUVGo5HE5qQ0NDl2PHjs0EaJrwOXPmzGNMJjNDVVVVePz4cQ/JeKmfGxsbFd3d3U+Wl5f3IITQ\nvL29A3r06FEueW2YLfZL48k6gA6EJ+sAOhRcZVT+4MJfX4C6unpla3MT0JfT1AHBzzzqaGiA/5aj\ntoQLf7UT+C0bdVwCWQfQgQhkHUCHgutwyB+8pfIFVFRUaMg6BiTuszrpCMkddXUtWYeA0HvhLRXU\nqdBoNIKfeYQQ+jh4SwUhhBBC7QJ2OBBCnwzvk0sPtqV0YXvKH5zDgTodGg3ncKD2Q1NTHcrKKmQd\nBkKfrdPN4fjll1/+/e9///sXAICioiJtJyensOaShsnCn3/+OX7YsGH578v46unpGeTk5BQ2ceLE\n8x9Sr0Ag4Pn7+y8LCwtzCgsLc8rNzdX/+eeff5NO1O0LjUYjUVGyjgKhD8fnAz7yimQO53B8gl9/\n/XWVtOpqLaX8p7h48eIEKv9Iaz7nMVsnJ6cwaXc2pN0OCCGEOp4O2+E4c+aMG4PByKTT6Vk+Pj5+\nAAA+Pj5+NTU13dhsdpq7u/tJagXLuXPnHjI0NMweM2bMdSr/R2Fhoc64ceOuGhsbJ9vY2NzJy8vT\nA2gaXZg3b94BMzOzBMlf3EFBQZ7Ozs6X7O3tIwYPHvxoz549i7Zt27acw+Gkmpubx5eVlWm1VHdc\nXJxFWFiY04oVK7ZyOJzUhw8fDjl8+PAcLpebyGKx0idNmhRaU1PTjTrXzZs3R5uYmCTp6enlhYeH\nOwA05TOZMWNGIIPByORwOKkCgYAn2S5BQUGeixcv3g0A8PLly68nTJhwkcVipbNYrPT4+HhzyePV\n1NSqli5dut3Q0DB79OjRN1+/fv0VAACPxxMsWbJkh4mJSVJAQIB3ZGTkKA6Hk8pgMDJnzZp1tK6u\nThkAICkpycTS0jKWxWKlm5qa3hUKhaqNjY2KK1as2MrlchOZTGbGoUOH5gIAPH/+vJ+Njc0dKo19\nbGyspUgkUvD09Ayi0+lZDAYjc+fOnT+19vcTEhIymU6nZ7FYrHRbW9vbn/kxQu+Rni7rCDoObEvp\nwjkc8qdDzuF49uxZfx8fH7/U1FSOpqbmG3t7+4g///xzvJ+fn8/evXsXpqWlsQGabqk8ePBg6Nmz\nZ78/dOjQ3ClTppw7f/78xGnTpp2aO3fuoYMHD/6oq6tbcPfuXdMFCxbsi4yMHEXVHx8fb97cSENO\nTo5Beno6q6amppuOjk7h1q1bV6SmpnKWLl26/cSJE9O9vb0DWqr7u++++8vJySnMxcXlAgCApqbm\nmzlz5hwGAFi7du2mo0ePzlq0aNEeQgjt8ePHg5KSkkwKCgp0+Xx+VEFBge7evXsXKioqNmZmZjLy\n8vL07O3tI/Lz84e11E5eXl67+Hx+1MWLFyeIRCKFqqoqNcljqquru5uYmCRt37596aZNm9Zu2LBh\n/e7duxfTaDRSX1+vlJSUZFJbW6sybNiw/Fu3bo3U1dUt8PDwOL5///758+fP3//999+fDQ4OdjUy\nMkqpqqpSU1FRqT169OgsTU3NN4mJidy3b992tbKyirG3t4+4cOGCy9ixY6/9+9///oUQQhMKhapp\naWnsZ8+e9adue1FrmrTUhps2bVobERFh369fv+e4/glCCMmPDtnhSEpKMuHz+VG9evUqAQCYNm3a\nqTt37tiMHz/+T8ljBw8e/IjBYGQC/DeNu1AoVI2Li7OYPHlyCHUc9Y2dRqORyZMnhzTX2aDRaITP\n50epqqoKVVVVhZqamm+cnJzCAJpSpWdmZjJaqxsA/nGPLCsri75mzZrN5eXlPaqqqtTGjh17jTqP\nq6trMACwjR86AAAgAElEQVSArq5uwZAhQx7ev39/eGxsrKWXl9cuAAA9Pb28QYMGPW6twxEVFcX/\n448/fgBoyqyqoaHxPzPTFBQURFOmTDkHAPDDDz/8QXWGAACo1/Py8vQGDx78SFdXtwAAwMPD4/je\nvXsXjho1KrJfv37PjYyMUgCaRksAACIiIuyzsrLooaGhkwCaOhEFBQW6JiYmSTNnzjxWX1+v5Ozs\nfInJZGbo6OgUPnz4cIiXl9cuBweHcHt7+4iqqiq1+Ph48+ba0NLSMtbDw+O4q6trsHisqG2wWLKO\noOPAtpQuzKUifzpkh+Pd4k7/+cVNWkn13rVr17fUz4qKio21tbUqIpFIQUtLq4waCZHUvXv36pbO\nLV6fgoKCiCpTqdLfV7d4nJ6enkF//fXXd3Q6Pev48eMezd0ikXyf5KSe9833+JhJQJLtqKqqKvzU\nOvfs2bPIzs7uhuTr0dHR1pcvX3b09PQMWrp06XZ3d/eTGRkZzOvXr485cODAvODgYNedO3f+pKmp\n+aa5Nty/f//8xMREbnh4uIORkVFKSkqKUc+ePUvFj/HzA+jbt+lnNTUAXd3//mNPDWtjGcvyVKbI\nU/pyLHfsclukp++QT6k8f/68n7m5eXxKSoqRpqbmm7Fjx17z8vLa5eTkFNazZ8/SV69e9enSpUuD\n5FMq/v7+y6qqqtTWr1+/wdLSMnbJkiU7Jk2aFEoIoWVlZdEZDEbmjBkzAh0dHS8394RIUFCQZ0pK\nitHu3bsXAzSNnlC/8MT3tVS3l5fXLg6Hk+rp6RkEANC7d+/i3NxcfU1NzTfffvvtlW+++ebJsWPH\nZnp6egYVFxf3vnz5suPDhw+H8Hg8QWFhoc7evXsX5uTkGBw5cmR2fn7+MHt7+4gHDx4MjY2NtaSe\nUhGPw83N7YyZmVmCt7d3QGNjo6JQKFSVHOVQUFAQnTlzxm3KlCnnNm/evKa4uLh3QECAN5/Pj9q2\nbdtyIyOjlNraWhU9Pb28W7dujdTR0Sn09PQMMjIySpk3b96B4cOH3z937twUY2Pj5MrKSvXu3btX\nHzt2bOaVK1e+DQkJmdylS5eG/Pz8Yf/617/+7/Xr118NGDDgqaKiYuPevXsXFhQU6K5Zs2azkpJS\nvYaGRkV2drbh9OnTT6SmpnJaasPCwkIdHR2dQgAALpebeOTIkdnUCBYAPqUibenp+M1cWlpqS3xK\n5dMIBAIc5ZAifEqlBf369Xvu5+fnw+fzo1gsVrqxsXEydWtj7ty5hxgMRiY1abSldOanTp2advTo\n0VksFivd0NAw+6+//vpO8hhJraVKF9/XUt3ff//92a1bt64wMjJKefjw4ZBNmzatNTU1vWtlZRUj\n/qgsjUYjAwcO/JvL5SZ+++23Vw4ePPijsrJy3YIFC/aJRCIFBoOR+f333589fvy4h5KSUn1z6d4B\nAAICAryjoqL4DAYj09jYOPnevXsjJK9JVVVVmJiYyKXT6VkCgYC3bt26jZLXp6KiUhsYGDhj8uTJ\nIQwGI7NLly4N8+bNO6CkpFR/7ty5KYsXL97NYrHSx4wZc/3t27ddZ8+efURfXz+Xw+Gk0un0rPnz\n5+9vaGjoIhAIeCwWK53D4aQGBwe7/vTTTzufPn06gM/nR1ETfamnjFpqw5UrV/5OTRa2tLSMFe9s\nIIQQkp0OOcKBpEddXb2ysrJSXdZxSAuOcKD2Bkc4kDzAEQ7U5j5nzQ+EEEKIgiMcqFPBDhRqb3Bp\n80+DczikSxojHB3yKRWEWoOdbOnBf9SlB9sSdXQ4woE6lXePTMs6DIQQaldwDgdCCCGE2gW8pYI6\nHUxPj+SdpromlFWUyTqMdg1vUcmfTtnhaMsU9ZKLf4mztLSMjY2NtWzpvZ8bi4ODQ/iZM2fcmlui\n/FNj7oiiAJ+LlZZ0SAcW4Mpf0iDelvxKvoyjQUj6OuUtlbZMUd/aUxCtdTakITw83OFjOxsAn/fk\nRkNDQ6fstKIm2NmQHmxL6cLRDfnToTscskhRDwDw5MmTb/h8ftSwYcPyN27cuI56nUpeRgihrVix\nYiuVcj04ONhVso6WUs1XV1d3d3V1DTYwMMhxcXG5YGZmlpCamsoBANDW1i4qLS3tCQBw4sSJ6Uwm\nM4PFYqV7eHgcBwAICwtzMjMzS+BwOKl2dnY3Xr161ae19vP19fV1d3c/aWFhETds2LD8I0eOzAYA\nEAgEPGtr6+jx48f/aWhomP327duuzcXa2NiouHz58m10Oj2LyWRm7NmzZxEAQEpKihGPxxMYGxsn\njx079tqLFy/6AgDs2rXLy8DAIIfJZGa4ubmdAQC4ffu2LZvNTmOz2WkcDidVKBSqAgBs3bp1BZXe\n3tfX1xcAQCgUqjo4OISzWKx0Op2e1Vy7IoQQko0O++1UVinqCSG0xMREbk5OjkG3bt1qTExMkhwd\nHS9zOJxU6tgLFy64ZGRkMDMzMxnFxcW9TUxMkmxtbW+L19NSqvl9+/Yt6NWrV0lOTo5BTk6OAYvF\n+k96J6r+nJwcgy1btqyOj48379mzZ2lZWZkWAIC1tXV0QkKCGQDAkSNHZv/+++8rt23btry1mcfZ\n2dmGCQkJZlVVVWpsNjvNwcEhHAAgLS2NnZOTYzBo0KDH/v7+y5qL9dixYzP//vvvgRkZGUwFBQVR\nWVmZVn19vdLixYt3h4WFOfXq1avk3LlzU1avXr3l6NGjs3777befi4qKtJWUlOqp1PL+/v7L9u3b\nt8Dc3Dy+urq6e9euXd9GRETYFxQU6CYmJnJFIpHC+PHj/4yOjrYuLi7uPWDAgKfh4eEOAP9NZY/a\nDt5SkR5sS+nCORzyp8N2OGSVoh4AwN7ePkJLS6sMAMDFxeVCdHS0NYfDSaX2x8TEWE2dOvU0jUYj\nffr0eWVra3ubyldCHdNSqvnY2FjLn376aScAgIGBQY5krhBCCO3WrVsjXV1dg6ksqVQsT548+cbV\n1TX4xYsXfevq6pSHDBnysLU2pNFoZPz48X927dr1bdeuXd/y+fyoxMRErqam5hsul5s4aNCgx63F\nGhkZOWr+/Pn7FRQURFQc2dnZhjk5OQajR4++CdA0CtK/f/9nAAAMBiNz6tSpp52dnS85OztfAmia\n97JkyZId06ZNO+Xi4nJhwIABTyMiIuwjIiLs2Wx2GkDTyEZBQYGulZVVzLJly/x9fHz8HB0dL1tZ\nWcW0dn0IIYS+nA7b4ZBVivrmRjyoX7gtxdbc+6j3NneO9z0L3Vz9AACLFy/evXz58m2Ojo6Xb9++\nbUvdivgY1LVIpqb/0FgJITQDA4OcuLg4C8ljw8PDHe7cuWMTFhbmtGXLltXZ2dmGP//882+Ojo6X\nw8PDHSwtLWOvX78+BgBg1apVv86dO/eQZB1paWns8PBwhzVr1mweNWpU5Nq1azdJHuMHftAXmvLT\nq4Ea6ILuf75ZpkPTgBGWP6xMvSYv8bTnMgtY/ylT5CldeXsrUynW5SWe9lbG9PQfQZYp6levXr0l\nOzvbUEVFpdbMzCwhMDBwBofDSaUSoV28eHHCwYMHf7xy5cq3JSUlvUxMTJISExO51dXV3alYduzY\nsUQy1Xx+fv6wXbt2eT18+HDIvn37FuTm5uqzWKz0hIQEMw6Hkzp48OBHKSkpRi9evOg7YcKEi+K3\nVLS0tMo4HE7qkSNHZnM4nNQZM2YEFhUVaUdFRfFbekrF19fX988//xxP3VLhcDipd+/eNb1///5w\nKt09AEBzsT548GDosWPHZt68eXP02bNnv1dUVGwsKyvTUlNTq9LX1889efKku5mZWUJ9fb3SgwcP\nho4YMeLe48ePB2lraxfV19craWtrF927d29EcXFxbyrd/OTJk0Pc3d1PduvWrWbt2rWbIiMjR6mq\nqgqfPn06QFlZua6hoaGLlpZWmYqKSu3ly5cdjx49OuvixYsTxK+JRqMRfEoFyTs+8HFFXCRXcOGv\nVsgyRT2Xy02cOHHieSaTmTFp0qRQ6nYK9Z4JEyZcZDAYmUwmM2PUqFGRW7duXdGnT59X4sc0l2qe\nSkFfXFzc28DAIGft2rWbDAwMcnr06FEuHoO+vn7u6tWrt9ja2t5msVjpy5Yt8wdo6kBMnjw5xNjY\nOLl3797FzaWsl7wWBoORyefzo8zNzePXrVu3sW/fvi8kj28uViUlpfrZs2cfGThw4N8MBiOTxWKl\nnzlzxk1JSak+NDR00s8///wbi8VKZ7PZafHx8eaNjY2K7u7uJ6mJp97e3gEaGhoVAQEB3tSkU2Vl\n5bpx48ZdtbOzuzF16tTT5ubm8QwGI9PV1TW4srJSPSsri25qanqXzWanbdq0aW1zoxtIuiS/kaNP\nh20pXdS3dSQ/OuwIR0clEokU6uvrlbp27fq2sLBQx87O7kZ+fv6wLl26NEj7XBs2bFivpqZWRXVY\nOgIc4ZAunOgoPf9YhwNHOD4bThqVLkze1gkJhULVkSNH3qqvr1cihND2798/vy06GxTMropag50N\n6cG2lC7sbMgfHOFAnQp2oFB7gEubI3mDczgQ+gSEENyktEVFRck8ho6yibcldjY+H87hkD/Y4UAI\nIYRQm8NbKqhTebdGiazDQAihdgUnjSL0CTA9PUJtT11dCyoqSmUdBpIjnfKWSlBQkOfixYt3y3ud\nn+vx48eDzpw54/a+44qKirTFl1Xv+AhuUtui5CCGjrJ1rLasrJTtPBScwyF/OmWHoy3I49MPjx49\nGnz69OmpbX0ekUiEnyOEEEKt6jS/KAIDA2fo6enlmZqa3hXP41FcXNx70qRJoVwuN5HL5SbGxcVZ\niEQihcGDBz8qLy/vQR03dOjQB8XFxb2bO17yXEVFRdojR468xWQyM0aPHn3zyZMn3wD8N629iYlJ\nkp6eXh6V1TQoKMjT2dn5kr29fcTgwYMf7dmzZ9G2bduWczicVHNz83gq22thYaHOuHHjrhobGyfb\n2NjcycvL06Pq9fb2DrC0tIzV0dEpPH/+/EQAAB8fH7/o6GhrNpudFhAQ4P348eNBNjY2d4yMjFKM\njIxS4uPjzVtrM4FAwLOxsbnj6Oh4efjw4ffnz5+/n7qHp6amVrV8+fJtLBYrPT4+3nz79u1L6XR6\nFp1OzwoICPCm6jhx4sR0JpOZwWKx0qdPn36ipTYHaD4V/fPnz/vZ2NjcYbPZaXQ6PSsmJsYKACAi\nIsLewsIizsjIKMXV1TWYSlvv4+PjR6W4X7FixdZP+aygj8GTdQAdCE/WAXQouA6HHJL1o2BfYnv2\n7Fm/gQMHPn79+nWvuro6JUtLy5jFixfvIoSAm5vb6ZiYGEtCCDx+/HjgiBEjcgkh4O3tvTMwMNCT\nEAIJCQmmdnZ2Ea0dHxgY6Llo0aLdhBBwdHQMO3HihDshBI4dOzbD2dn5IiEEPDw8gsaNG3eFEAIP\nHjzQ/de//vWktra2a2BgoKeuru6Dqqoq1eLi4q80NDTKDx48OJcQAkuWLNm+c+dOb0IIjBw5MvLB\ngwe6VEwjR46MpOp1dXU9RwiB3NzcEbq6ug8IISAQCGwdHR3DqHaorq7uVltb25UQAvn5+UONjY2T\nCCHw6NEjbUNDw6xmHtPjqaio1Dx69Ei7sbFRwc7OLiI0NHQiIQRoNJooJCRkEiEEkpOTjeh0emZ1\ndXW3qqoqVQMDg+y0tDRWdna2wbBhw/JKSkp6EkKgrKxMs7U2dHJy+isuLs6cEAJCobB7Q0ODor+/\n/9ItW7b8mxACjY2NCpWVlWrFxcVf2djY3K6uru5GCAE/P7+fN27cuLakpKSnnp7efSr+8vJyDclr\nAgACQHDDDbc234CgjuPd3yd8ztYpJo3evXvXVDxV/ZQpU87l5+cPAwC4efPm6Hv37o2gjq2srFSv\nrq7uPmXKlHMbN25c5+npGXT27Nnvp0yZcq6l46lv15SEhASzS5cuOQMA/PDDD3+sXLnyd4Cm2y6u\nrq7BAAC6uroFQ4YMeXj//v3hNBqN8Pn8KFVVVaGqqqpQU1PzDZX3hU6nZ2VmZjKEQqFqXFycxeTJ\nk0Oo89TV1SlT9VLp3EeMGHHv5cuXXwMAEIkZxXV1dcqLFi3ak5GRwVRUVGyk2qA1XC43UVtbuwgA\nwM3N7UxMTIzVxIkTzysqKjZSyetiYmKsXFxcLnTr1q0GAMDFxeVCdHS0NXW9PXv2LAUA0NTUfNNa\nGzaXit7ExCRp5syZx+rr65WcnZ0vMZnMDIFAwMvNzdW3sLCIo67LwsIirkePHuUqKiq1s2bNOuro\n6HjZ0dHxcvNX5QkA2u9+1gQAFvz326Xg3Z9Y/rDyTsD2k1aZ+lle4vn8sqyznVLkIftqeyu3RbbY\nTtHhaC1VPSGEdvfuXVNlZeU68feYmZklFBQU6L5+/fqrP//8c/y6des2tnZ8c2npPzQ2AICuXbu+\npV5TUFAQUWUFBQVRQ0NDF5FIpKClpVWWlpbGbq4e8XhaOveOHTuW9OvX7/nJkyfdGxsbFVVUVGo/\nND6qXio9vYqKSq148jfJ9m0tlpbasLlU9NbW1tHR0dHWly9fdvT09AxaunTpdi0trTI7O7sbzc1P\nSUxM5EZGRo4KDQ2dtGfPnkWRkZGj/veqglq5Yh6WP6rMknhN1vFgWZ7Kkrc1sNx+yjwe7x/lDRs2\nwOdS+Owa2gEul5t4+/Zt29LS0p719fVKISEhk6l99vb2Ebt27fKiyunp6SyApl+iEyZMuLhkyZId\n+vr6uVpaWmWtHS/+i9XCwiLu7Nmz3wM0ZZy1sbG5Qx0TEhIymRBCKyws1Hn48OGQ4cOH32+tc0Lt\nU1dXrxw8ePCj0NDQSdTrmZmZjNauW11dvbKyslKdKldUVGj07dv3BUDT3IrGxkbF97VdYmIit6io\nSFskEimcO3duipWVVYzkMdbW1tGXLl1yrqmp6SYUClUvXbrkbGNjc2fkyJG3QkJCJpeWlvYEAKDm\norTUhoWFhToGBgY5K1eu/N3ExCQpLy9P7++//x7Yu3fv4tmzZx+ZPXv2kbS0NLaZmVlCbGysZWFh\noQ5AU36ZBw8eDBUKhapv3rzRHDdu3NXt27cvzcjIYL7v+tDn4sk6gA6EJ+sAOhScwyF/OkWHo1+/\nfs99fX19zc3N462srGIMDAxyqH27du3ySk5ONmYymRkGBgY5hw4dmkvtmzJlyrlTp05No26ntHa8\neMr23bt3Lw4MDJzBZDIzTp06NY2aREmj0cjAgQP/5nK5id9+++2VgwcP/qisrFwnme5d8meqfOrU\nqWlHjx6dxWKx0g0NDbP/+uuv71p6DwAAk8nMUFRUbGSxWOkBAQHeCxYs2Hf8+HEPFouVnpeXp6em\nplbV3PvFXzMxMUlatGjRHn19/VwdHZ3CCRMmXJQ8ns1mp3l6egZxudxEMzOzhDlz5hxmMpkZ+vr6\nuatXr95ia2t7m8VipVNZZ1tqQ8lU9GPHjr0mEAh4LBYrncPhpAYHB7t6e3sHfPXVV6+DgoI83dzc\nzjCZzAwLC4u4vLw8vcrKSnUnJ6cwJpOZYW1tHb1jx44lH/dJQQgh1FZwpdEvaMaMGYFOTk5hLi4u\nF2Qdy4cQCAQ8f3//ZWFhYU6yjkVamjpK+JmXHgHgN3NpEUDHaksayPL3C6anly5M3obalOTIC0II\nIfSpcIQDdSrYgULoy8ClzTsWzKWC0CfATjZCCH15eEsFIfTJMF+F9GBbShe2p/zBDgdCCCGE2hzO\n4UCdCs7hQB2BpqY6lJVVyDoM1IngHA6EPkFUlKwjQOjz8PmVsg4BoY+Gt1SQTAUFBXkuXrx496e8\nNyMjg3n16tVxVNnX19fX399/mfSiQ++Tni7rCDoObEvpwjkc8gc7HKjdSktLY1+5cuVbqoy3SxBC\nSH5hhwN9cYGBgTP09PTyTE1N78bFxVlQrxcXF/eeNGlSKJfLTeRyuYnUvsTERK6FhUUch8NJtbS0\njM3Pzx9WV1envG7duo3nzp2bwmaz04KDg10BAHJzc/X5fH6Ujo5O4e7duxfL6ho7CxZL1hF0HNiW\n0oWrjMqfVieN3rt3b8SzZ8/6m5qa3hXPu3Ht2rWxY8eOvfZFIkQdyvPnz/uZmZklpKamcjQ0NCr4\nfH4Uh8NJ3bVrl9fUqVNPL1y4cK+lpWXs33//PXDs2LHXcnNz9SsrK9W7d+9eraio2Hjz5s3RBw4c\nmBcaGjrp+PHjHikpKUZUIjhfX1/fGzdu2AkEAl5FRYWGnp5e3suXL79WVFRspM5Po9HImDEAffs2\nldXUAHR1//uPPTWsjWUsy3N5yZKm9WTkKZ05ljtWWTI9/YYNGz570miLHY5du3Z57d27d+GIESPu\npaWlsQMCArydnZ0vATQl62opTTpCrbl06ZLzpUuXnIOCgjwBmhLd5efnD9u9e/fiPn36vBowYMBT\n6tjXr19/lZeXp1dSUtLLy8trV0FBgS6NRiMNDQ1dcnNz9YOCgjxTUlKMqJGMDRs2rFdWVq5btWrV\nrwAA+vr6uTdv3hzdv3//Z1SdNBqN4KRR6UlPx2/m0vIxbcnn4wJ274O5VKSrTZ9SOXTo0NyUlBQj\nNTW1qqKiIu1JkyaFFhUVaf/00087P+eEqHOj0WhE/ENLCKFRcy8IIbS7d++aKisr14m/Z8GCBftG\njRoVefHixQmPHz8exOPxBC3VL/5eRUXFxoaGBnwSCyGE5ECLczgIITTqNoq2tnaRQCDgXb16ddyS\nJUt2fG4vB3VeXC438fbt27alpaU96+vrlUJCQiZT++zt7SOo2yMATU+hAABUVFRoUKMUgYGBM6j9\nGhoaFZWVlepfMn70Tzi6IT3YltKFoxvyp8UOR58+fV6lp6f/538BNTW1qsuXLzuWlJT0yszMZHyZ\n8FBH069fv+e+vr6+5ubm8VZWVjEGBgY51L5du3Z5JScnGzOZzAwDA4OcgwcP/ggAsHLlyt9XrVr1\nK4fDSW1sbFSkRkT4fH5Ubm6uvvikUXxSBSGE5FOLcziePHnyjZKSUn3fvn1fiL9OCKHFxsZaWllZ\nxXyRCBGSIpzDIV04h0N6cA6HdOEcDumSxhwOXNocdSo4AoI6Alza/P2wwyFd2OFA6CO9m7Qq6zAQ\nQqhdkUaHAxf+QgghhFCbww4HQuiTYb4K6cG2lC5sT/mDaxSgTodGw6e6EUJtR1NdE8oqymQdhtzB\nDkc7JRAIeP7+/svCwsKcJFfc/NJ1Pnv2rL+3t3eA+Joa0iAew8GDB3/s3r17tbu7+8nPrTcK8DEV\nhFDb4VfyZR2CXMIORwfQFk9efEyd/fv3fybtzoZkDD/++ONBadePEELoy8E5HHLk2rVrY42MjFJY\nLFb66NGjbwIACIVC1ZkzZx4zNTW9y+FwUv/666/vJN/X2sxhadd5+/ZtWzabncZms9M4HE6qUChU\nLSoq0qbT6VkAANXV1d1dXV2DDQwMclxcXC5QidoAmhaPW7NmzWYWi5Vubm4e/+rVqz4AAGFhYU5m\nZmYJHA4n1c7O7gb1ujhfX19ff3//ZQAAPB5P4OPj42dqanpXT08vLyYmxqqlc6ekpBi9v+XRp0qH\ndFmH0GFgW0oXtqf8wREOOVFcXNx77ty5h6Kjo60HDRr0+M2bN5oAAFu2bFk9atSoyGPHjs188+aN\npqmp6V2q4yCLOv39/Zft27dvgbm5eXx1dXX3rl27vhXfv2/fvgW9evUqycnJMcjJyTFgsVj/+b++\nurq6u7m5efzmzZvX/Pzzz78dPnx4zurVq7dYW1tHJyQkmAEAHDlyZPbvv/++ctu2bcvFOz00Go1Q\nIx40Go00NjYq3r171/Tq1avjNmzYsP7GjRt2zZ0b191ACCH5gB0OOZGQkGBma2t7e9CgQY8BADQ1\nNd8AAERERNiHhYU5bdu2bTkAwNu3b7s+efLkG1nVaWlpGbtkyZId06ZNO+Xi4nJBPLsrAEBsbKwl\nleDPwMAgh8FgZFL7lJWV6xwcHMIBAIyMjFJu3LhhB9C0qq2rq2vwixcv+tbV1SkPGTLk4fvicHFx\nuQAAwOFwUouKirTfd25xfuAHfaEpP70aqIEu6AILmpZ4pL4VYfnDytRr8hJPey6zgCVX8bT3sizb\nkyJP6eY/tiyZnl4asMMhJySzqIq7cOGCy9ChQx+Iv/b8+fN+kseJRCIFDoeTSqPRyHffffeXiYlJ\n0ufWKennn3/+zdHR8XJ4eLiDpaVl7PXr18dIjnK0dE4lJaV66mcFBQURlcl18eLFu5cvX77N0dHx\n8u3bt219fX193xcHdU7JjLAfsjCND/i0uE/8FymWsYxlLH9OWXKl0/ZU5vF4/yhv2LABPhfO4ZAT\npqamd+/cuWNDfVsvLS3tCQAwZsyY6+IZVNPS0tgt1aGgoCBKT09npaWlsTds2LBeGnVKKiws1DEw\nMMhZuXLl7yYmJkl5eXl64vstLS1jqURqubm5+llZWfT31SmeDTYoKMizuWMIIbT3dSY+5dzo8+B9\ncunBtpQubE/5gx0OOdG7d+/iQ4cOzXVxcbnAYrHS3dzczgAArF27dlN9fb0Sg8HINDQ0zF6/fv1/\nupnicxqam6vQFnUGBAR40+n0LCaTmaGsrFw3bty4q+LvW7Bgwb7i4uLeBgYGOWvXrt1kYGCQ06NH\nj3LxYyTr9/X19Z08eXKIsbFxcu/evYubi6GleD703AghhGQLc6kgqRKJRAr19fVKXbt2fVtYWKhj\nZ2d3Iz8/f1iXLl0a5OHcNBqN4DocCKG2xAd+h8vmK41cKjiHA0mVUChUHTly5K36+nolQght//79\n879EZ0PW50YIIdQ6HOFAnQo+JosQamsdcWlzHOFA6BNgJ1t6BALB/8x0R58G21K6sD3lD45woE7l\n3ePHsg4DIYTaFWmMcOBTKgghhBBqc3hLBXU6mJ4etXfq6lpQUVEq6zDkGt5SkT9yNcIhEAh4Tk5O\nYQBNC0AtXrx4t6xjAgA4ePDgjydPnnT/2Pfdvn3bNj4+3vxz65E28WRr4sTbX5KDg0N4RUWFRttH\n97Sqkq8AACAASURBVPHmzJlz+N69eyM+/B0EN6ltUXIQQ0fZPrwtKys71oRE1DnI7QjH5zxN0NDQ\n0EWaj0N+amr0qKgovrq6eqW5uXn859QjD8LDwx1kHUNLDh8+PEfWMXRePFkH0IHwZB1Ah4KjG/Kn\nzUc4pJ0evbS0tKezs/MlJpOZYW5uHk8tX+3r6+vr7u5+0srKKsbDw+P469evv7Kzs7thaGiYPWfO\nnMPa2tpF1NLeEyZMuGhsbJxsaGiYLf7LqqX06VRq9OfPn/ejUrOz2ey0Ll26NDx58uSb5tKrFxUV\naR88ePDHHTt2LGGz2WkxMTFW4inW09PTWWZmZglMJjPDxcXlApXJtaXU6+KEQqHq6NGjbxoZGaUw\nGIxMqv2Kioq0R4wYcW/u3LmHDA0Ns8eMGXO9trZWBQAgJSXFiMlkZrBYrPR9+/YtaK5taTQaqaio\n0HB0dLw8fPjw+/Pnz99P/T2It9/27duX0un0LDqdnhUQEOBNvX/Tpk1rhw8fft/a2jp66tSpp6lr\nLSws1Bk3btxVY2PjZBsbmzvUcuienp5B3t7eAZaWlrE6OjqF58+fn0jVtXXr1hVcLjeRyWRmULlV\nhEKhqoODQziLxUqn0+lZISEhk6k2S01N5YhEIgVPT88gOp2exWAwMnfu3PlTc9eJEEJIBgghbba9\nevWq9zfffPN3UVHRIEIIlJWVaRJCYNWqVb/88ccf06jXhg0blicUCrtHRUXxHB0dwwghEBgY6Llo\n0aLdknUuWrRo98aNG9cSQuDWrVt8FouVRgiB9evX+xobGyfV1tZ2JYTAwoUL9/j5+f1MCIFr166N\nodFoopKSkp6EECgtLdUihEB1dXU3Q0PDLKpMo9FEly9fdiCEwMqVK3/bvHnzakII+Pr6rt+2bdsy\n8Tj27NmzcMqUKWfFr4sQAocPH569bNmybdT7/P39l1L7xMt0Oj3zzp071oQQWLdu3YaffvppByEE\neDxe1PLly7cSQuDKlSvjRo8efUOyDRoaGhQrKirUCSFQXFz8la6u7gNCCDx69Ei7S5cu9RkZGQxC\nCLi6up6j2plOp2dGR0dbEUJgxYoVvxsaGmZJ1hsVFcVTUVGpefTokXZjY6OCnZ1dRGho6ERCCGhr\naz8qKSnpmZycbESn0zOrq6u7VVVVqRoYGGSnpaWxEhMTTVgsVtrbt2+VKysr1YYOHZpPXevIkSMj\nHzx4oEsIgYSEBNORI0dGEkLAw8MjyNXV9RwhBHJzc0dQ13H9+nX7uXPnHiSEQGNjo4Kjo2PYnTt3\nrM+fP+8yZ86cQ1S85eXlGlSbpaSkcJKTk43s7OwiqP1v3rzpIXmNAEAACG5S26LkIIaOsn1MWwJB\nrYuKipJ1CB3Ku88cfM7WprdU2iI9emxsrOWFCxdcAAD4fH5USUlJr8rKSnUqQyqVRTQ2Ntby0qVL\nzgBNycq0tLT+c9MzICDAm9r35MmTbx48eDCUy+UmtpQ+vbkYjhw5Mjs2NtaSqqOl9OqkmZGaiooK\njfLy8h7W1tbRAAAeHh7HJ0+eHELtby71ujiRSKSwatWqX6Ojo60VFBREz54960+NxgwePPgRlZbd\nyMgopaioSLu8vLxHeXl5DysrqxgAAHd395NXr14d19y1cbncRG1t7SIAADc3tzMxMTFWEydOPE9d\nS0xMjJWLi8uFbt261VCxRkdHW4tEIgVnZ+dLysrKdcrKynXUXBChUKgaFxdnIX59dXV1ygBNIyrO\nzs6XAABGjBhx7+XLl18DNH0+IiIi7NlsdhpVR0FBga6VlVXMsmXL/H18fPwcHR0vU9dD0dHRKXz4\n8OEQLy+vXQ4ODuH29vYRzV0jgCcAUM2qCQAs+O9wtuDdn1j+sHL6e/ZjuW3K70pylM4cyx2r3O7S\n00sj5XpzWqqze/fu1e87TiAQ8CIjI0clJCSYqaio1PL5/CjqtkNL6dMlY5w9e/aRsLAwJ+p8n5Je\nvbU4W0q9Tjl16tS0169ff5WamspRVFRsHDx48CPqGsRTxSsqKjbW1NR0e9/5xInPnSGE0BQUFESS\n+8Xf/76fRSKRgpaWVllLGWmVlZXrmnv/qlWrfp07d+4hyePT0tLY4eHhDmvWrNk8atSoyLVr126i\n9mlqar7JyMhgXr9+fcyBAwfmBQcHux49enTW/541qKXLh/+9j47l1suSd61kHU97Lkvue//75Smd\nubyVm5vDIU/xyXuZx+P9oyz36enbIj26tbV19KlTp6YBNHUeevfuXayurl4p+UtUPFV5RESEfVlZ\nmRZA0+iClpZWmYqKSu39+/eHJyQkmL3vnFTdDQ0NXSZPnhzy+++/r9TV1S2g9reUXl1dXb2ysrJS\nXbIuDQ2NCi0trTJqfsbJkyfdeTye4EPboKKiQqNPnz6vFBUVG6OioviPHz8e1NrxPXr0KNfU1HxD\njchQ7decxMREblFRkbZIJFI4d+7cFPFRBBqNRqytraMvXbrkXFNT000oFKpeunTJ2cbG5o6lpWVs\nWFiY09u3b7tWVVWpUZNM1dXVKwcPHvwoNDR0EnX9mZmZjNbiHTNmzPVjx47NFAqFqgAAT58+HVBc\nXNz7+fPn/VRUVGqnTZt2avny5dvEPzeEEFpJSUmvxsZGRRcXlwubNm1am5qayvmQ9kQIIdT22nSE\nQzw9ukgkUvj6669fXr9+fczatWs3/fTTTzsZDEamSCRSGDJkyENq4uP70pH7+vr6zpw58xiTycxQ\nVVUVHj9+3KO549evX7/Bzc3tzMmTJ93Nzc3j+/bt+0JdXb1y7Nix1w4cODBPX18/V09PL496gkT8\n3JL1UT/HxcVZpKSkGK1bt27junXrNgIAXLly5VsqvbqWllbZyJEjb1EdACcnp7BJkyaF/vXXX99R\nHSyqzuPHj3vMmzfvQHV1dXcdHZ3CwMDAGc21YXNtMG3atFNOTk5hDAYj09jYOHnEiBH3WjqeKgcG\nBs6YOXPmMRqNRuzt7SOaq5dGoxETE5OkRYsW7SkoKNAdOXLkrQkTJlwUr4fNZqd5enoGcbncRICm\nR1KZTGYGAMB33333F4PByPz6669f0un0LCo1/KlTp6bNnz9//+bNm9fU19crubm5naFu+0i2OQCA\nnZ3djXv37o2g/m7U1dUrT5486V5QUKC7YsWKrQoKCiIlJaX6AwcOzBN/79OnTwfMmDEjUCQSKQAA\n+Pn5+TTXpkiaBND8N3P08QSAbSk9uA6H/OmwS5vX1dUpKyoqNioqKjbGx8ebL1y4cC9+421bQqFQ\nVVVVVVhdXd3d1tb29uHDh+ewWKz097/zy2nq1HTMz7xsCAB/SUqLAD68LWnQUf/tlhbscEgXJm9r\nxd9//z3Q1dU1WCQSKSgrK9fhWg1tb+7cuYdyc3P/v707j2ri+h4AfsOmFlBAsbiVKEgEQhISIIIg\nRJTqqbhERahVQCtVqliXilYpUMWlghZcv1AF3BeK2nRxo0YRBBRIWEUFoqhdIpKyqAgyvz/o/BpT\n1ABBAtzPOXPkJTOTN9cxPt57867V8+fPe/v5+SWoW2MDdQS3zq5AN+LW2RXoVrCxoX66bQ8HQi3B\n9PSoO8ClzdG7hsnbEGqD9j5Ljtsra7d0eh26y9aaWGJj4+3IRzyR+sAGB0IIIYQ6HA6poB4Fh1TQ\nu2RgoA9VVdWdXQ2E2g0njSLUBpcvd3YNUE/B49V0dhUQUhs9fkhFPiV7QkKC39KlS3d2xOfIJ257\nFzrqWuTjJRAIPLdu3RoMACCVSo25XG4mh8PJJhcYa4+zZ89OlU85HxoaGp6SkuLe3vMi1RLhc0gq\ng3MOVAvjqX6wh0NOR3a3d6Wu/MbGRi0tLa3Gt+3n6ekpIBsfKSkp7gwGI681jx83NTVpKC6dTjp9\n+vR0T09PAbmoWXh4eKiy50UIIaR+umUPx7lz5yZyOJxsFoslGj9+/CWA5kWp5s+ff4DL5Way2ewc\ncmVTea8bn6qtrdXz9/ePZzAYeUwmU5ycnMyPj4/3X758+Q5yn7i4uIUrVqzYDgBw8ODBeWQqeF9f\n30TF870uXbu8rKwsBycnp3Q2m50zZsyYtNu3b1sANPdc8Pn85EmTJv1qYWFxOzg4eCt5THx8vD+N\nRivhcrmZ6enpTi1dS1hYWNjcuXMPOTk5pVtYWNz+/vvvPwVo7rlwcXFJnTp16lk6nV5QX1/fi7xm\nNpudIxQK3RTPRfaiiMViZnBw8NazZ89OZbPZOc+fP+994cIFDycnp3QOh5Pt5eV1klymnEqlStas\nWbOFw+Fknzp1atb333//qYODQxaLxRLNnDkz6dmzZ33S09OdBAKB55dffrmNzWbnlJWVjfDz80sg\n09enpKS4s9nsHAaDkbdgwYL9ZDI4KpUqCQsLC+NwONkMBiOvpbgi1WKxOrsG3QeuG6FaGE/10+16\nOKRSqXFAQEBsamqqi6mp6T2ZTGYAABAREbHO3d095cCBA/NlMpkBl8vNJBsjb7Nhw4YQQ0PDKjIH\niEwmM9DW1m6IiIhYFxkZuUpTU/NlQkKCX2xsbEBhYaF1RETEuuvXrzsaGRk9IT8f4N9ejoCAgNj/\n/e9/n5mbm9/NzMzkBgYG7lEcLrC0tCxOTU110dTUfHnp0qXxX3311SYyH4lYLGaKRCKWjo7OCxqN\nVhIUFBSjoaHRFBYWFpaTk8Pu27dvNY/Hu8xms3Naup6CggJ6RkbG6NraWj1bW9tcMkNubm6ubWFh\nobWpqem9qKiolZqami/z8vIYJSUlNA8Pjwtko0cRk8kUf/PNN19nZ2dzYmJigh4/fjwgIiJiXUpK\ninufPn2ebd26NXj79u0rQkJCNlAoFGLAgAGPs7OzOQDN+XU+/fTT7wEAQkJCNuzfv3/BkiVLdk2Z\nMuVHT09PAZk5l1xe/vnz5739/f3jf/vtt3Hm5uZ3fX19E/fu3bt42bJl0RQKhTA2NpZmZ2dz9u7d\nuzgyMnIVLviGEELqods1ODIyMka7urpeMTU1vQfQnEEUoDmBm0Ag8IyMjFwFAFBfX9+roqJimDLn\nTElJcT9x4sRsskyec9y4cb8JBALPUaNG3WpoaNC2trYu3Llz51IvL6+TRkZGT+T3Jb0pXbs8mUxm\nMG/evIN37941p1AohHzWWHd39xR9ff0aAAArK6siiURClUqlxm5ubsL+/ftXAgDMnj37REsNBAqF\nQkydOvVsr1696nv16lXP4/EuZ2VlORgYGMgcHByyyLilpaWNCQoKigEAoNFoJaampvde1+AAaO4d\nInuIMjIyRhcVFVk5OTmlk9dH/kzWjfw5Pz/fZv369Rv//vvvfrW1tXoTJ048J39Oxc8oKSmhDR8+\nvJxMnufr65u4e/fuz5ctWxYNAEA2UNhsdk5ycjK/pbpu2QJgYtL8s54egLn5v7+pk3MSsKxcOSkJ\n4/emMsCrS2y/LR04SR3Sk3f1Msazh6Wn7wyK6dPlJScn80eOHHlH/rXff/99kDLnbemcn3766fcR\nERHrLC0ti+fPn3/gbZ8P8PZ07aSQkJAN7u7uKadPn55+7949U/lssoop6BsbG7UU54i05vElch6F\nrq5u3ZvO0Zp5KBMmTLh49OjRj1t6T/5z/Pz8En788ccpNjY2+YmJib7yQzevSzCnWEf518jYkHFp\n6fPXvCGlm+IQAZbfXJZvbKhDfdStDPDmlOlYxrK6lt26Wnr6zsDlcjOvXr06ViKRUAGau+wBmlOe\nkxlbAZqHD5Q954QJEy7u3r37c7JMDpM4ODhkPXjwYOjRo0c/9vHxOQbQ3Otx6tSpWeTnVlVVGZLH\nEQRBUTZdu3zK+9dlkiVRKBSCy+VmXrlyxfXJkydGDQ0N2qdOnZrV0r4EQVDOnj07tb6+vldlZWV/\noVDoZm9vf0OxceHi4pJKprG/ffu2xf379z+g0WglysSLy+VmpqWljSktLTUDaO7VuXPnzsiW9q2t\nrdUzMTH5o6GhQfvw4cOfkI0HfX39murq6r6K10mj0UokEgmVPPehQ4fmurq6XlGmXkj1cA6H6ih+\n+aP2wXiqn27X4DA2NpbGxsYG8Pn8ZBaLJSIbAiEhIRsaGhq0GQxGHp1OLwgNDf3/5ppiGnrFc65f\nv35jVVWVoY2NTT6LxRLJ/xbu5eV10tnZ+RqZit3Kyqpo3bp1Ea6urldYLJZo5cqVUYqfc+TIkTn7\n9+9fwGKxRHQ6vaClCayrV6/+du3atZvZbHbOy5cvNd9WRxMTkz/CwsLCHB0drzs7O1+ztrYufF0P\nAYPByOPxeJcdHR2vf/3119+YmJj8oXjewMDAPU1NTRoMBiPP29v7eGJioq+2tnaD/H6v+9nY2Fia\nkJDg5+Pjc4zJZIqdnJzSXzeBc8OGDSFcLjfT2dn5GvlECgCAt7f38W3btn3J4XCyy8rKRpCv9+rV\nqz4+Pt5/1qxZpxgMRp6WllbjokWL9snH901xQggh1DlwpdF28vT0FKxYsWI7j8frEstJhYeHh+rp\n6dXKN4R6EgqFQuDCX6ojEmEvx5vweM25e5SB6dRVC+OpWpi8rRPJZDIDGo1W8t577z3tKo0NEv7m\njxBC6F3DHg7Uo2BjC71LmEsFdReYSwWhNsBGNkIIvXs4pIIQajPMV6E6GEvVwniqH+zhQD0OhdKu\nXkGEkBow0DeAquqqzq4GagWcw4F6FAqFQlyGLjXHFyHUAh7wcHj0HcKnVJBKyaeeb296ezc3N2FO\nTg4bAEBPT6+2vXV79OjRYPnl4BFCCHUt2OBALWrv0xyKi3C1tz6DBw9+9LrVU1HnEYGos6vQbWAs\nVQvncKgfbHB0c+fOnZvI4XCyWSyWiMyOW1dXpzt//vwDXC43k81m57S00unrus6ysrIcnJyc0tls\nds6YMWPSyIRuz5496+Pt7X3cysqqiM/nJz979qyP/HHr16/fyGKxRI6Ojtf/+uuvgQDNmX1nzpyZ\n5ODgkOXg4JCVnp7uBABw5coVV1tb21xbW9tcNpudU1dXpyuRSKg2Njb5AABkxlgGg5HHZrNzyJVf\nExIS/Ph8fvKkSZN+tbCwuB0cHLxVhaFECCHUDjhptBuTSqXGAQEBsampqS6mpqb3yBwwERER69zd\n3VMOHDgwXyaTGXC53EyyMfI2lpaWxampqS6ampovL126NP6rr77alJSUNHPv3r2L9fT0aouKiqzy\n8/Nt2Gx2DnlMXV2drqOj4/WNGzeuDw4O3hoXF7dw3bp1EcuWLYtevnz5jjFjxqTdv3//g4kTJ54r\nKiqyioqKWrlnz55AR0fH60+fPn1PPlkdAMDu3bs/19TUfJmXl8coKSmheXh4XCAbPmKxmCkSiVg6\nOjovaDRaSVBQUMyQIUMeqjKu6F8swGVGVQVjqVq4yqj6wQZHN5aRkTHa1dX1Cply3sDAQAYAcOHC\nBQ+BQOAZGRm5CgCgvr6+V0VFxTBlzimTyQzmzZt38O7du+YUCoUgM7Kmpqa6kCnibWxs8hkMRh55\njI6OzouPPvroZwAADoeTffHixQkAAJcuXRpfXFxsSe5XU1OjX1dXpztmzJi05cuX75gzZ84RPp+f\nrNhgSEtLGxMUFBQDAECj0UpMTU3v3b5924JCoRDu7u4p+vr6NQDNeW0kEglV8fgtsAVMoDk/vR7o\ngTmY//+XPdmtjWUsY1n9y+qUzr27lTsiPT0+pdKN/fTTT5OPHz/uffjw4U/kX7ezs7t57Ngxn5Ej\nR96Rf10oFLpFRUWtFAgEngkJCX7Z2dmcnTt3LpXfx8/PL8HOzu7mkiVLdkkkEiqPx7tcXl4+fPr0\n6aeDgoJiyGXeORxOdlxc3EI2m52jr69fU1NTow8AkJSUNPPnn3/+KD4+3t/Y2Fj68OHDITo6Oi8U\n615YWGj9888/f7Rnz57A8+fPf9irV696T09PQX5+vg2fz09eunTpTvKzxo4de3X37t2f5+TksG/e\nvGlH1tnT01Pw5Zdfbhs7duxV8rz4lIpqiUCEv5mrCMaydd72lArmUlEtfEoFvRGXy828evXqWIlE\nQgUAePLkiREAwIcffng+JiYmiNwvNzfXVtlzVldX9x08ePAjgOY5E+TrY8eOvXr06NGPAQAKCgro\neXl5jLedy8PD44J8PUQiEQsAoLS01Mza2rpw9erV39rb299QzDTr4uKSeuTIkTkAALdv37a4f//+\nB6NGjbrV0j+G9v4DQQghpBrY4OjGjI2NpbGxsQF8Pj+ZxWKJfHx8jgEAhISEbGhoaNBmMBh5dDq9\nIDQ0NJw8pqV08/JWr1797dq1azez2eycly9fapL7LF68eG9tba2elZVVUWhoaLidnd1NxXMqnjcm\nJibo5s2bdkwmU2xtbV0YGxsbAAAQHR29zMbGJp/JZIp1dHReTJo06Vf58wQGBu5pamrSYDAYed7e\n3scTExN9tbW1G1qqM+ZO6Vj4G7nqYCxVC3s31A8OqaAeBYdUEOoecOGvd0sVQyrY4EA9CvZ4INQ9\nvG1pc5zDoVqYLRahNsBGturgl7rqYCxRd4c9HKhHoVAoBN7zCCHUOviUCkIIIYS6BGxwoB6HQqHg\n1sW3vn2NOvs2UjnM/aFaGE/1g3M4UA+EQyqqIwQAt3f+qTU1uLwKQl0N9nCoMSqVKmEwGHlk4rVH\njx4NJt/T0NBoWrVqVSRZjoyMXBUeHh4qf7z82hsd4bvvvvtCMUlbS6hUqoRcdEwZYWFhYVFRUSsB\nAEJDQ8NTUlLc21NP1JHcOrsC3QZOGFUtjKf6wQaHGqNQKIRQKHQTiUQsZ2fna5s3b15Lvqejo/Pi\n9OnT0ysrK/uT+8ofW1xcbNm7d+/nmZmZ3KdPn77XEfWLjo5epsy5/5moqfSvpPLXEh4eHuru7p7S\n1joihBBSD9jg6CJGjx6dUVpaakaWtbW1GwICAmJ37NixvKX9jx075uPj43PMw8PjwtmzZ6e2tI+f\nn19CYGDgHkdHx+tmZmalQqHQzdfXN9HKyqrI398/ntzvwoULHk5OTukcDifby8vrZF1dnW5MTEzQ\no0ePBvN4vMtkg2Dx4sV77e3tb9Dp9IKwsLAw+c/69ttvVzMYjDwul5tJXodEIqGOGzfuNyaTKR4/\nfvyllhLI+fn5Jfzwww8zAADWrFmzxdraupDJZIpXr179bWuuAXUUYWdXoNvAOQeqhfFUP9jgUHNk\nz8C5c+cm0un0Avn3AgMD9xw5cmROdXV1X8XjTp486eXl5XXSy8vr5LFjx3xaOjeFQiFkMpnB9evX\nHXfs2LF8ypQpP65evfrbwsJC6/z8fBuxWMx8/PjxgIiIiHUpKSnu2dnZHA6Hk719+/YVQUFBMYMH\nD34kFArdyCGPTZs2fXXjxg17sVjMvHLlimtBQQGd/CwDAwNZXl4eY8mSJbu++OKL7wAAli5dutPf\n3z9eLBYz58yZc4TMAKtYRwqFQlRWVvY/c+bMtMLCQmuxWMxcv379RmWvoT3xRwghpBo4aVTN8Xi8\ny0+ePDHS0tJqlP8PHABAX1+/Zt68eQdjYmKC+vTp84x8/ebNm3bGxsbSQYMG/T5w4MC//Pz8Eqqq\nqgwNDQ3/syyfp6enAACATqcXmJiY/GFtbV0IAGBtbV0okUioFRUVw4qKiqycnJzSAQBevHihQ/6s\n6MSJE7Pj4uIWNjY2av3++++DioqKrMhGEjmXxNvb+/jy5ct3AABkZGSMPnPmzDQAgE8++eQw2WvR\nEgMDA1nv3r2fL1iwYP/kyZN/mjx58k/KXgOTyRS/ejY/AKCSZwYAFvw7F0H4z59YVq5MvvauP/+f\nkhql825vmUwJri716epljKf6pacHgiBwU9ONSqWWV1ZWGjU2NmrOnj37+Pbt25eT7+np6dUQBAFP\nnjwxpFKp5eHh4V+HhYWFEgQBK1asiOrfv/9jKpVaTqVSy/v16yeLi4v7VPH8fn5+8UlJSTMIgoDy\n8nIqnU7PV3xPIBBM9vHxOfqm+hEEAWVlZcPNzc3vyGSyfuTxiYmJ88j9ysvLqQRBwIsXL7QHDBgg\nJQgCBgwYIG1oaNBSfD0sLCw0KipqhWId6+vrdX755ZdJ8+fP3z9u3LgUZa9Bvs4AQAAQuHX5DQiE\n0Lvzz7+5//w/0JoNh1S6AE1NzZfffffdF1FRUStra2v15N8zNDSs8vLyOrl///4FFAqFaGpq0jh1\n6tSsgoICenl5+fDy8vLhZ86cmfa6YZU3oVAoxOjRozPS0tLGkPMu6urqdO/cuTMSoLmHhRzOqa6u\n7qurq1vXt2/f6j///PP9X3/9dRJ5HoIgKCdOnJgN0NwLQvaQODk5pR8/ftwbAODIkSNzxo4de5Xc\nn1CYZFpXV6crk8kMJk2a9Ov27dtX4FCJuhB2dgW6DZxzoFoYT/WDQypqTP5pDRMTkz/4fH7y7t27\nPw8ODt4q/97KlSujdu3atQQA4Nq1a85Dhw59YGJi8gf5vouLS2pRUZHVn3/++f7777//5+s+Q/FJ\nFwCAAQMGPE5ISPDz8fE5Vl9f3wsAICIiYt3IkSPvBAQExE6cOPHckCFDHqakpLjb2trmjho16taw\nYcMqnJ2dr8mft6qqypDJZIp79+79nGz87Ny5c6m/v3/8tm3bvhw4cOBf8fHx/uT+ivWqqanRnzp1\n6tnnz5/3JgiCIj9Z9k3X0NI1IYQQevcwlwrqUZobIHjPd30UwO8uhN4dCgWzxSLUBrhKZVenr2/Y\n2VVACLUSzuFAPU57Jz7h9u92+fLlTvnc6uonnX0bqRzOOVAtjKf6wQYHQgghhDoczuFAPco/y6x3\ndjUQQqhLwTkcCLUBhYJzONrKwEAfqqqqO7saCKEuCBscqMe5fLmza9B18Xg1r5SFQiFm5VQRjKVq\nYTzVD87hUFNkanoGg5FnbW1dGBISsoFcB0MikVA1NDSayLU3AACWLFmyKzEx0ZcsNzY2ahkbG0vX\nrl27uaPquGnTpq+U2U9PT6+2NeeVT9i2cOHCuOLiYsu21A8hhJD6wAaHmiJT0+fl5TGysrIcF6ox\nlgAAC7VJREFUysrKRnz22Wf/I98fOHDgXzExMUENDQ3a5P7yx1+8eHECh8PJJv/j7gibN29eq8x+\nrV18S37hr7i4uIWWlpbFbakf6nj4G6TqYCxVC+OpfrDB0QXo6urW7du3b9GZM2emyWQyAwAAY2Nj\nqbu7e4p8r4a848ePey9evHjviBEjyq5fv+7Y0j5ubm7CFStWbLe3t79haWlZfOPGDfvp06eftrCw\nuB0SErKB3O/w4cOfcLncTFtb29xFixbta2pq0lizZs2WZ8+e9bG1tc2dO3fuIQCAadOmnbGzs7tJ\np9ML4uLiFsp/1ooVK7bT6fSC8ePHX3r8+PEAAACRSMQaPXp0BpPJFPP5/GTy2hTrmJOTw25qatLw\n8/NLsLGxyWcwGHnR0dHLWnMNCCGEOhc2OLoIfX39muHDh5eTeUwAAFavXv1tZGTkqqamplf+Hp8/\nf977t99+Gzdp0qRf35aevlevXvU3btywX7x48d6pU6ee3bdv36KCggJ6QkKCX1VVlWFxcbHlyZMn\nvdLT051yc3NtNTQ0mo4cOTJny5Yta/r06fMsNzfX9tChQ3MBAOLj4/1v3rxpd+PGDfuYmJigqqoq\nQ4DmPCj29vY3CgoK6K6urlfCw8NDAQDmzZt3cNu2bV+KxWKmjY1NPvm6Yh0BAHJzc20fPXo0OD8/\n3yYvL4/h7+8fr+w1qOrvAP0XrnWgOhhL1cJ4qh+cNNqFKD6SNHz48HIul5t59OjRj+Vf/+mnnya7\nubkJdXR0XkybNu1MWFhYWHR09LKWhjamTJnyI0Bzanc6nV5A5loZMWJE2f379z9ITU11yc7O5tjZ\n2d0EAHj27Fkf+Twt8qKjo5eR6eYrKiqG3blzZ6SDg0OWhoZG0+zZs08ANKeh5/P5ydXV1X3//vvv\nfi4uLqkAAL6+vomzZs069bprNzMzKy0rKxsRFBQU89FHH/3s4eFxQdlrMDQ0rJI/15YtACYmzT/r\n6QGYmwOwWM1lkaj5Tyy3XAZ4dTKe6J8d1CGdNpaxjGXVlYUdkJ4e1+FQU8OHDy/Pzs7mGBkZPQEA\nqKmp0R82bFjFvXv3TKuqqgw9PT0F+fn5NiUlJbSZM2cmubq6XrG3t7/h6+ubOGPGjB/S0tLG9OnT\n5xkAgFQqNT5z5sy08ePHX5L/DB6PdzkqKmolm83OEQqFblFRUSsFAoEn+V5kZOSq69evOz569Ghw\nSxNE9fX1a2pqavQBAIRCoVtISMiGixcvTujdu/dzHo93OTw8PHTs2LFXtbS0Gl+8eKGjoaHRVFZW\nNmLmzJlJQqHQzcbGJv/evXumAAClpaVmXl5eJ7Ozszn+/v7xnp6eAj6fnyxfx7q6Ot3z589/eOjQ\noblGRkZP9u/fv+Bt10C+R9aZQqEQ+JRK2/F4gDlMEOqBVLEOBw6pqDHyL7e2tlYvMDBwz/Tp00/3\n69fvb/l9aDRaiZWVVZFAIPCkUChEdXV132vXrjlXVFQMI9PT79q1a0lb09O7u7unJCUlzZRKpcYA\nAE+ePDG6f//+BwAA2traDY2NjVoAzenpDQ0Nq3r37v381q1bozIyMkaT52lqatI4derULACAo0eP\nfuzi4pLat2/fakNDw6pr1645AwAcOnRorpubm/B1caisrOz/8uVLTT6fn7xhw4aQ3Nxc29ZeD0II\noc6DQypqjMfjXSYIgtLU1KTB5/OT5SdByg+PrFu3LoL8Lf7MmTPT3N3dU7S1tRvI96dMmfJjcHDw\n1oaGBm351+UppoQnWVpaFm/cuHG9h4fHhaamJg1tbe2GPXv2BH7wwQf3AwICYhkMRh6Hw8nev3//\ngn379i2ysrIqotFoJY6OjtfJc+jq6tZlZWU5bNy4cf3777//54kTJ2YDACQmJvouWrRo39OnT98z\nMzMrJdPTt1S3hw8fDvH3948n56ts2bJljbLXgDqOENc6UBmMpWphPNUPDqmgHgWHVNpHcUgFv9RV\nB2OpWhhP1VLFkAo2OFCPgj0g7YNLmyPUM2EuFYTaABvZCCH07uGkUYRQm5GP0aH2w1iqFsZT/WCD\nAyGEEEIdDudwoB4F53CgrsxA3wCqqqveviNCKoZzOBBqg8uAj6mgrolXw+vsKiDUZjikghBqMxGI\nOrsK3QbGUrVwDof6wQYHQgghhDoczuFAPQqFQiFwSAV1VTzg4WPdqFNgLhWEEEIIdQk4aRT1OFtg\nC5hAc356PdADczAHFjTnXyfH0bGsXDkJkjB+KirLz+F43f4Ary7ZrU7pzNWtLD+HQx3q09XKmJ4e\noXbCIRXVEoHolf8MUdspE0scUlEe5lJRLcylglArYYMDdWXY4ECdBedwIIQQQqhLwAYHQqjNcO0I\n1cFYqhauw6F+cNIo6nF4gKs1oq7JQN+gs6vQZYhEIpzDoWawhwP1OARB4KaiLTQ0tNPr0F02ZWKJ\neVSUJ5PJOrsKSAE2OBBCCCHU4bDBgRBqM4lE0tlV6DYwlqqF8VQ/+Fgs6lEwPT1CCLUNgetwIIQQ\nQkjd4ZAKQgghhDocNjgQQggh1OGwwYG6jXPnzk0cNWrUrZEjR97ZunVrsOL7kZGRq2xtbXNtbW1z\nbWxs8rW0tBplMpmBMsf2NO2JJZVKlTAYjDxbW9tcBweHrHdfe/Xztng+fvx4wMSJE8+xWCwRnU4v\nSEhI8FP22J6oPfHE+/NVb4tlVVWV4fTp008zmUwxl8vNLCwstFb22P/o7GfPccNNFVtjY6OmmZnZ\n3fLycuqLFy+0mUymqKioyPJ1+wsEgsnu7u6X2nJsd9/aE0uCIIBKpZZXVlYadfZ1qMumTDxDQ0PD\n1qxZs5kgCJBKpQOMjIwqGxoatPDeVG08CQLvz9bGctWqVdu++eabEIIg4NatW7T2fG9iDwfqFrKy\nshzMzc3vUqlUiba2doO3t/fxs2fPTn3d/kePHv3Yx8fnWFuO7e7aE0sS0c7Z7N2JMvEcNGjQ79XV\n1X0BAKqrq/v279+/UktLqxHvzf9qTzzJ9/H+bKZMLIuLiy15PN5lAAAajVYikUiof/3118C23JvY\n4EDdwsOHD4cMGzasgiwPHTr0wcOHD4e0tO/Tp0/fO3/+/IczZsz4obXH9gTtiSVA86PH48ePv2Rn\nZ3czLi5u4buoszpTJp4LFy6MKywstB48ePAjJpMpjo6OXqbssT1Ne+IJgPenPGViyWQyxcnJyXyA\n5gbKvXv3TB88eDC0Lfcm5lJB3UJr1tcQCASezs7O1wwMDGStPbYnaE8sAQDS0tLGDBo06HepVGo8\nYcKEi6NGjbrl4uKS2jG1VX/KxHPTpk1fsVgskVAodCstLTWbMGHCRbFYzHwX9etq2hNPfX39Grw/\n/6VMLNesWbNl2bJl0eR8LVtb21xNTc2XbfnexB4O1C0MGTLkYUVFxTCyXFFRMWzo0KEPWtr3+PHj\n3vJDAK05tidoTywBmruzAQCMjY2l06dPP52VleXQsTVWb8rEMz093WnWrFmnAADMzMxKhw8fXl5S\nUkIbOnToA7w3X9WeeALg/SlPmVjq6+vXHDhwYH5ubq7twYMH50mlUmMzM7PSNn1vdvakFdxwU8XW\n0NCgNWLEiNLy8nJqfX29zusmMMlksn5GRkaVT58+7dPaY3vK1p5Y1tXVvVddXa1PEATU1tbqOjk5\npZ0/f96js69J3eO5fPny7WFhYaEEQcAff/zx/pAhQx5UVlYa4b2p2nji/dn6WMpksn719fU6BEFA\nbGzsQl9f3wRlj1XcOv2CccNNVdsvv/wyycLCosTMzOzupk2b1hIEAfv27fts3759n5H7JCQk+Pr4\n+BxV5tievLU1lmVlZcOZTKaIyWSKrK2tCzCWysVTKpUOmDx5soDBYIjpdHr+kSNHPn7TsT19a2s8\nS0tLR+D92bpYpqenO1pYWJTQaLRbM2bMSJLJZP3edOybNlzaHCGEEEIdDudwIIQQQqjDYYMDIYQQ\nQh0OGxwIIYQQ6nDY4EAIIYRQh8MGB0IIIYQ6HDY4EEIIIdTh/g/7bDzG3ZNWIwAAAABJRU5ErkJg\ngg==\n", "text": [ "" ] } ], "prompt_number": 186 }, { "cell_type": "code", "collapsed": false, "input": [ "# pandas density plot\n", "Piura[1].plot(kind='kde', linewidth=3);\n", "plt.axis([0.2, 1.5, 0, 1.9])" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 191, "text": [ "[0.2, 1.5, 0, 1.9]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAD5CAYAAAAgGF4oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XdYFFf3B/DvIogliiKWCChG7NLsJSq+xqjY2xuNDQtR\nY4lGI7ZETSxYY4tKNPb6xmhsaKLGVWMjgmDBgglEQEVFEBWVNr8/rvsbQMqyOzt3dvZ8nmefzF2W\n9ZwMHGbP3LmjEQQBhBBCLIcV7wAIIYTIiwo/IYRYGCr8hBBiYajwE0KIhaHCTwghFoYKPyGEWBhr\n3gHoQ6PR0JxTQggxgCAImpzPmc0RvyAIqnjMmjWLewyUj/pzUVs+aspFznzyYjaFXy2io6N5hyAp\nNeWjplwAdeWjplwA/vlQ4SeEEAtDhV9mvr6+vEOQlJryUVMugLryUVMuAP98NPn1gZRCo9EI5hAn\nIYQoiUajgWDOJ3fVQqvV8g5BUmrKR025AOrKR025APzzocJPCCEWhlo9hBCiUtTqIYQQAoAKv+x4\n9/akpqZ81JQLoK581JQLwD8fKvzEogkC8Po1+y8hloJ6/MTi3LsHbNsG7N0L3L0LvHgBvP8+0KIF\nMHYs4O3NO0JCpJFXj58KP7EYDx4A334LrF8PZGTk/bqePYF164AKFeSLjRBToJO7CsG7tyc1c8gn\nIwNYuRKoVYsV9JxFX/P/vxZaAMD+/UDLlsA//8gZpfTMYd/oS025APzzocJPVC00FGjWDPjiC+D5\nc/H51q2BnTuBhAQgLQ0ICwPatxe/fvcua/38/bf8MRNiatTqIar04gXwzTfAihVAZqb4fO3awPLl\nwMcfZz3SF+3fD/TvD7x5I77+wgWgTBl54iZEStTjJxbj4EF2kjYmRnzO1haYMQOYMoVt5+f0aaBD\nB7H4f/QRcOwYUKSI6WImxBSox68QvHt7UlNSPrGxQK9eQPfu2Yt+u3bAtWvA11/nX/R1ubRpA2ze\nLD5/4gSwZIlJQjYpJe0bY6kpF4B/PlT4idnTnbytU4e1anQcHNi0zePHgRo1Cvee/fqxTwg6X3/N\nzhcQogbU6iFm7fp1YOhQ4PLl7M8PHw4sXAiUK2f4e6elAR9+CAQHs3GtWqz4lyhh+HsSIidq9RBV\nEQRg1SqgUaPsRb92bdaj37DBuKIPADY2wI4dQMmSbHz7NjB5snHvSYgSUOGXGe/entR45BMfD3Tu\nDIwfL56AtbUFvvuOTcts3dqw980tF1dXNjNIZ+1a4PBhw95fbmr6WVNTLgD/fKjwE7Ny9izg4QEc\nPSo+5+HBWjAzZxY8Y8cQw4axq3l1hg4FHj6U/t8hRC7U4ydmQdfamTQJSE8Xn580CZg3zzQFP6uE\nBMDdHbh/n407dACCggArOnQiCkY9fmK2UlKAQYPY1be6ou/gAPz2G5tmaeqiD7DzBVu3ihd9/fYb\nm0lEiDmiwi8z3r09qZk6n/v3Wc9+xw7xucaNWWvn44+l/bcKyqVdu+wnd/392TkFpVLTz5qacgH4\n50OFnyhWaCgr8iEh4nMjRgBnzgDOznximjsXaNCAbaemsuUdkpP5xEKIoajHTxRp/35g4EDW5gHY\ncgmrVgGjR/ONC2DTOhs0EGPz8WHLRNCSDkRpqMdPzMaGDUDv3mJhtbNjs3iUUPQBdiHXjz+K46Ag\nYOJEuosXMR9U+GXGu7cnNanzWb0a8PMTi2j16sDFi9mXTDaVwuQyYAAwdao4XrWKLfGgpOKvpp81\nNeUC8M+HCj9RjB07gHHjxHHDhqzo167NL6b8zJsH9O0rjhcsAL76Kvsy0IQoEfX4iSKcPAl07ChO\n12zRgrVQ7Oz4xlWQ1FSgTx/g0CHxuf/+F9i4UVzqgRBeaD1+olgPHrCrbx8/ZuN69dgVumXL8o1L\nX69fs9k9v/4qPlerFrBrF+DlxS8uQujkrkLw7u1Jzdh8MjPZxVm6ol+pErvpCY+ib2guxYoBe/cC\nY8aIz92+zW75uHw5v76/mn7W1JQLwD8fKvyEqx9+YG0egF0Vu3074OTENyZD6Kabbt4stnhSU9ls\nn48+AiIjuYZHSDbU6iHcPHoE1KwJPHvGxjNmsAukzN2dO6z1k/XGLba27GYuX30FFC3KLzZiWajH\nTxRnxAjgp5/Ydo0a7PaIcqy7I4c3b1ihX7o0+yyfevXYNQAtWvCLjVgO6vErBO/entQMzefqVbHo\nA2zNe95FX8p9Y2sLLFrE7t6lW+IBAG7cAFq2BEaNAhITJfvncqWmnzU15QLwz4cKP+Fi3jxxu0sX\noFMnfrGYUsOGwKVLwLJl2ad3BgaymT9btyrroi9iGajVQ2R38yZreeh2aUhI9qNitbp3Dxg7Nvuc\nfwDw9gbWrGE3iydEStTqIYoxf75Y9Dt3toyiDwBVqgAHDrAF6LLOXNJq2XUM06eL6xMRYkpU+GXG\nu7cntcLmc/8+sHu3OJ45U9p4jCHHvtFogB492KeeyZPFFT3T0tiSD3XrZr+tpDHU9LOmplwA/vmY\ntPAPGzZsY8WKFePd3Nyu5fZ1rVbrbWdn98zLy+uKl5fXlblz5yqoDBBTWL9eXJahVSt2kZMleu89\nYPFi4MoVdrJX599/2TLPkyax6wAIMQWT9vjPnj3b6r333nsxePDgrdeuXXPL+XWtVuu9bNmyLw8e\nPNgt3yCpx68KaWmAi4t439rdu4FPPuEakiJkZrILv6ZMYff21WnThrWGlL5eEVEuLj3+Vq1anS1b\ntmy+k9ZyC4qo04EDYtGvVAno2ZNvPEphZQUMGwbcusVmOOmcPs1uO/nwIb/YiDpx7fFrNBrh/Pnz\nLTw8PMJ9fHyCIiIi6vKMRw68e3tSK0w+gYHitp+f8q5g5b1vHBzYnbyyTnW9epVNdTXk9o6885GS\nmnIB+OdjzfMfb9CgQWhMTIxziRIlUo4ePdqpR48ev965c6dmbq/19fWFi4sLAKBMmTLw9PSEt7c3\nAPF/Io2VO378GDh5ko0BLerVAwDlxJcVz3g0GqBFCy2mTgUWL/ZGRgYQFqbFf/4DXLjgDRsb88pH\nqnFYWJii4lFqPlqtFps3bwaA/6+XuTH5PP7o6GiXrl27Hsqtx59TtWrVokJCQhra29s/zfo89fjN\nX0AAMG0a227XDjhxgm885mDjRmD4cHHs78/+PxKiL0XO44+Pj6+oCyo4OLiJIAianEWfmD9BYFeo\n6gwZwi8WczJsGDB7tjhetEhcyZQQY5i08Pfv339XixYtzt++fbuWs7NzzMaNG4cFBgaODAwMHAkA\ne/fu7ePm5nbN09MzbMKECct3797dz5TxKEHOj+HmTp98Ll9m89YBtmyBUk/qKnHffP018PHHbFsQ\n2L0LdKuZFkSJ+RhKTbkA/PMxaY9/165d/fP7+pgxY34YM2bMD6aMgfC3Y4e43acPm8NO9GNlBWzZ\nwq7sffSI3a1s1ix2gxdCDEVr9RCTysgAnJ1ZwQKA338H2rfnG5M5+t//xGserKzYWv8eHnxjIsqn\nyB4/Ub8//xSLfvnyQNu2fOMxV337spPiALvga9w4WtWTGI4Kv8x49/akVlA+e/aI2336ANZcJxDn\nT8n7RqMBVq8GbGzY+OxZ9ukpP0rOp7DUlAvAPx8q/MRk0tPZTch1+qn+1L1p1a7NLnzTmTmTjvqJ\nYajHT0zm+HFxRkrlykBMDOtPE8Pdvw9Urw68fs3G+/Ypd5YU4Y96/ER2WZdf/u9/qehLoXJl4PPP\nxfG8eXTUTwqPfhVlxru3J7W88klNZUejOuawCqe57Bt/f6BYMbYdEsJu5JIbc8lHH2rKBeCfDxV+\nYhK//w4kJbHtqlWBpk35xqMmFSoAvr7ieNEibqEQM0U9fmISgwYB27ez7SlTgIUL+cajNnfvspu1\nZ2aycXg44O7ONyaiPNTjJ7JJSQF+/VUc02we6bm6Ar17i+PVq/nFQswPFX6Z8e7tSS23fA4dAl68\nYNs1awKenvLGZChz2zfjx4vbO3a8u4aPueWTHzXlAvDPhwo/kVzWtXkGDGAXHxHptWwJ1K/PtlNS\ngG3b+MZDzAf1+ImkEhLYbRV1N1SPjGRtCWIaa9YAY8aw7bp1gevX6Q8tEVGPn8ji55/Fot+0KRV9\nUxs4UFztNCKCLeVASEGo8MuMd29Pajnz0c3kAVibx5yY474pXZoVf501a8Rtc8wnL2rKBeCfDxV+\nIpnoaODcObZdpIh5XLSlBqNHi9v79gHx8fxiIeaBevxEMvPnAzNmsO1OnYCgIL7xWJKWLYHz59n2\nvHnA9Ol84yHKQD1+YlKC8O5sHiKfrOv3/PijeGEXIbmhwi8z3r09qenyCQ9nJxcBoEQJoHt3fjEZ\nypz3TZ8+gL092/73X+DUKfPOJyc15QLwz4cKP5HE5s3ido8edF9dudnaZv+UtWkTv1iI8lGPnxjt\n1SvA0RFITGTj334T1+En8rlyBWjQgG0XKwY8fAjY2fGNifBFPX5iMr/8Ihb9atWAjz7iG4+l8vIS\nb8D++nX2214SkhUVfpnx7u1JTavVYv16ceznZ743XFHDvhk6VNxesULLLQ6pqWHfZMU7HzP9FSVK\nER0NnDnDtq2tsxceIr8BA8QbskdEADdv8o2HKBP1+IlRRo0CAgPZdq9erO1D+OrdW7z7Gd0LwbLl\n1eOnwk8MlpAAODuzk7sAcPo00Lo135gIcPgw0LUr265Uid3k3tqab0yEDzq5qxC8e3tSWr8eePVK\nC4CdWGzVim88xlLLvunYEahYEQC0ePiQzbIyd2rZNzq886HCTwySmpr9rk8TJtBywEphbc1ufalD\nc/pJTtTqIQbZuBEYPpxtV6zIrha1teUbExFFRAD16rFtGxvg/n3AwYFvTER+1OohksnIAAICxPHE\niVT0laZuXaBJE7adlgbs3Mk3HqIsVPhlxru3J4V9+9idtQCgZElttmWBzZka9k1WLVtq/3/b3Ns9\nats3vPOhwk8KJTMT+O47cdyzJ7sZCFGe//yHLd0AAGFhbEkHQgDq8ZNC2r0b6N+fbZcsCURFAeXL\n842J5G3AALHNM24csHIl33iIvKjHT4yWng7MmiWOJ0ygoq90w4aJ29u3szV8CKHCLzPevT1jbN8O\n3LnDtu3sgEmTzDufnNSUC8DyadsWcHFh48RE4MABriEZTI37hicq/EQvqanAnDniePJkoGxZfvEQ\n/VhZZV8/aeNGfrEQ5aAeP9HL2rXi7f3KlWO9/VKl+MZE9HPvHjvqFwR2kV1UFFC1Ku+oiByox08M\n9uoVMHeuOJ46lYq+OalSBWjfnm0LArBlC994CH9U+GXGu7dniLVr2ZWfAPD++9lv7G2O+eRFTbkA\n2fPJepJ30ybzuxm7mvcND1T4Sb5evAAWLBDHM2awm6kT89K9u3hOJjqa3YydWC7q8ZN8zZ/Pij3A\nWgZ37tDyDOZq/Hhg1Sq2/emnwI4dfOMhpkfr8ZNCS0pi99BNSmLjDRvEhdmI+QkLY8tnA+yP94MH\nNDNL7Qw+udurV699R44c6ZyZmUltIQnw7u0VxtKlYtF3dQUGD373NeaUT0HUlAvwbj6enmLhf/MG\n2LVL/pgMpfZ9I7cCi/no0aPX7tixY4Crq+vdqVOnBty+fbuWHIERvh4/BpYvF8dz5oj3ciXmK+sn\nNprTb7n0bvUkJSWV2b17d7+5c+fOrFKlyj0/P7/1AwcO3G5jY5Nm4hip1cPBlCnA4sVsu149IDwc\nKFKEb0zEeImJbGbWmzdsHBYGeHjwjYmYjlHz+BMSEspt3rzZd8OGDSMaNGgQOn78+JUhISEN27dv\nf1z6UAlvCQnAmjXieM4cKvpqUbYs0KuXOA4M5BcL4afAwt+zZ8/9H3744Z8pKSklDh061PXgwYPd\n+vXrt3v16tVjnz9/TpfxFBLv3p4+VqwAXr5k225ubOnlvJhDPvpSUy5A3vl89pm4vW0bkJwsTzzG\nsJR9Ixfrgl7g5+e33sfHJyjrc2/evLG1tbV9ExIS0tB0oREekpPFKX8AMH06W++FqEebNuwOXRER\n7DqNbduAMWN4R0XkVGCP38vL68qVK1e8sj7XoEGD0NDQ0AYmjSwL6vHLJyAAmDaNbbu6ArduUZtH\njX74ARg7lm3XqQPcuMHW8SHqklePP88j/gcPHrx///79yq9evSoeGhraQBAEjUajEZKTk0unpKTQ\ntZsqlJICLFsmjqdOpaKvVoMGsf374gVw8yZw+jTg7c07KiKXPD/E//bbbx0mT568JC4uznHSpElL\nJ0+evGTSpElLly1b9uX8+fOn6/Pmw4YN21ixYsV4Nze3a3m9Zvz48Str1KgR6eHhEZ7zk4Ua8e7t\n5WfDBjaNEwCcnVlxKIiS8yksNeUC5J9P6dLZ92/Wk/lKZEn7RhaCIOT72Lt3b++CXpPX48yZM61C\nQ0O96tevfy23rx85csSnU6dOQYIg4OLFi02bNm16MbfXsTDV4dSpU7xDyNWbN4Lg5CQIbP1GQVi5\nUr/vU2o+hlBTLoJQcD7Xron729paEOLi5InLEJa2b6Tytna+U1Pz7PFv27Zt0KBBg7YtXbp0kkaj\nEbL8odBoNBrhyy+/XJbrN+YQHR3t0rVr10PXrl1zy/m1UaNGrWvbtu2pTz75ZA8A1K5d+9bp06fb\nVKxYMT7r66jHb3obNgB+fmy7QgW2kFfx4lxDIjJo0wY4c4Ztz56d/daaxPwVusev6+M/f/68VG6F\nX4qg4uLiHJ2dnWN0Yycnp9jY2FinnIUfAHx9feHy9h5yZcqUgaenJ7zfNiV1H5tobNj45Ent2194\nNu7RQ4tLl5QTH41NN/78c+DMGTZet84bU6cCFy4oJz4aF26s1WqxefNmAPj/epmr3D4GSPmIiopy\nyavV06VLl0N//vlnS924Xbt2J0JCQhrkfB2o1WNSO3eKH/nLlBGEZ8/0/14l5mMoNeUiCPrl8+aN\nILz/vrj/t2wxfVyGsMR9IwXk0eopcIb2lClTFiUnJ5dOS0uzadeu3UkHB4cn27Zt0+O0X8EcHR3j\nYmJinHXj2NhYJ0dHxzgp3pvoJzOTLb2sM348O/FHLEPRosC4ceJ42TL2J4CoXG5/DbI+3N3dwwVB\nwL59+3oOGzbsp6SkJDs3N7erBX2f7pHfEX/Wk7sXLlxoZgknd5Xm11/Fo72SJQXhyRPeERG5JSQI\nQokS4s/ByZO8IyJSQR5H/AVeuZuenm4NAIcPH+7Sp0+fvXZ2ds/07fH3799/1+nTp9s8efLEwdnZ\nOWbOnDmz0tLSbABg5MiRgT4+PkFBQUE+rq6ud0uWLPly06ZNQ434G0YKSRCAefPE8ejR7EbqxLLY\n2wO+vuKUzmXLgP/8h2tIxNRy+2uQ9eHv7x9Qq1atWx4eHmFv3rwpGh8fX6FJkyaXCvo+KR9Q0RG/\nknqVv/8uHuXZ2grC/fuFfw8l5WMsNeUiCIXL5/ZtQdBoxJ+HiAjTxWUIS943xoChPf6AgICp586d\naxkSEtKwaNGiqSVLlnx54MCB7ib/i0RMLuvR/vDhbLleYplq1gS6dhXHWe/FQNRHr/X4z5071/Lf\nf/+tqmvTaDQaYfDgwVtNHt1bNI9feufOAR9+yLatrYG7d4GqVfnGRPg6c4bN6weAYsWAe/eA8uX5\nxkSMU+h5/DoDBw7c/s8//3zg6ekZVqRIkQzd83IWfiK9rEf7AwdS0SdAq1ZAw4ZASAjw+jVbyG32\nbN5REZPIrf+T9VG7du2bmZmZmoJeZ8oHqMcvqZAQsZer0QjCrVuGv5cS8pGKmnIRBMPyyXpNR7ly\ngvDihfRxGYL2jWFgaI+/fv361x88eEDdXxXJOm+/Tx+gFt1FmbzVty+gu+AzIYHuy6tWBfb4vb29\ntWFhYZ5NmjQJtrW1fQOwnvvBgwe7yRIhqMcvpZs32T10df87r1wBPD35xkSUZfVq8aKuqlWByEjA\nxoZvTMQwefX4Cyz8Wq3W++0bCLo30Gg0Qps2bU6bItDcUOGXzpAhwNa3Z2c6dwYOH+YbD1GelBRW\n8J88YePt24EBA/jGRAxj8M3Wvb29tS4uLtFpaWk23t7e2iZNmgR7eXldMU2Y6qdbUImHqChgxw5x\nPGOG8e/JMx+pqSkXwPB8SpTIvozDokX8l3GgfSOtAgv/jz/++Fnfvn1/HjlyZCDA1tPp2bPnftOH\nRqS2aBGQ8XZeVtu2QPPmfOMhyjVmDPsDAABXrwLHjvGNh0irwFaPh4dHeHBwcJNmzZpd1N0hy83N\n7Vpu6+ubCrV6jHf/PlCtGpCaysYnTgDt2vGNiSjbhAnAihVs29sbOHWKazjEAAa3emxtbd/oTuoC\nbO0eqdbjJ/JZulQs+k2b0lospGBffinec1mrBYKDuYZDJFRg4W/Tps3pefPmzUhJSSlx/Pjx9n37\n9v25a9euh+QITo149PaePAHWrRPHM2YAmneOAQzDu1cpJTXlAhifT5UqQP/+4njhQuPiMQbtG2np\ntVZP+fLlH7u5uV0LDAwc6ePjEzR37tyZcgRHpLFiBZupAQDu7kCXLnzjIeZjyhRxe/9+4PZtfrEQ\n6ei1Vs+jR48qAECFChUemTyiXFCP33DPnrGpec+esfHu3cAnn/CNiZiXzp2BoCC2PWIEsH4933iI\n/grd4xcEQTN79uzZDg4OT2rVqnW7Vq1atx0cHJ7MmTNnVm5vRJRpzRqx6Nesya7UJaQw/P3F7a1b\ngQcP+MVCpJFn4f/+++8nnjt3ruVff/3VODExsWxiYmLZ4ODgJufOnWv5/fffT5QzSDWRs7f34gW7\nqYbO1KniyTqp8O5VSklNuQDS5dOqFdCsGdtOTeWzZDPtG2nlWfi3bt06eOfOnZ9Wq1YtSvfcBx98\n8M+OHTsGbN26dbA84RFjrFkjXn1ZtSpdfUkMo9FkP+pft078FEnMU549/vr161+/fv16/cJ+zRSo\nx194L1+yefuPH7NxYCDw2Wd8YyLmKzMTqFtXPLm7cGH2E79EmQrd47exsUkz5GtEGdauFYt+lSrs\nnqqEGMrKCvjqK3G8fDnw5k3eryfKlmfhv3r1qnupUqWe5/aQ86pdtZGjt5eSAixeLI6nTQOKFjXN\nv8W7VyklNeUCSJ/PwIFA5cps+8EDYNs2Sd8+X7RvpJVn4c/IyCjy/PnzUrk90tPTC7xzF+Fn3Trg\n0duJt87OwNChfOMh6mBry5Zx0Fm8WFz7iZgXvebx80Y9fv2lpAAffADEx7PxmjXA6NF8YyLqkZzM\nWoe6k7s//0xThJXM4LV6iHlZuVIs+k5OwLBhfOMh6lK6NFu5U2fuXP5LNpPCo8IvM1P29p4+BQIC\nxPGMGezjuSnx7lVKSU25AKbLZ8IEccnm8HDgyBGT/DPZ0L6RFhV+FVmwQPwIXqMGMHw433iIOpUv\nD4waJY7pqN/8UI9fJWJiWLHXTbH73//YjbMJMYWc93c4fhz46CO+MZF3UY9f5WbNEot+48Z0wo2Y\nVuXK2T9RzpvHLxZSeFT4ZWaK3t6NG8CWLeI4IEC69fYLwrtXKSU15QKYPp8pUwBra92/Bfz5p+n+\nLdo30qLCrwL+/uySegDo0IHurkXk4eICDBokjufM4RYKKSTq8Zu5oCC2XjrAjvJDQgAvL74xEcsR\nGQnUri0eeGi1QJs2XEMiWVCPX4VSU7NfSTlsGBV9Iq8aNYAhQ8TxzJk0w8ccUOGXmZS9vRUr2BEX\nwC6smT9fsrfWG+9epZTUlAsgXz7ffAPY2LDtP/8Efv9d+n+D9o20qPCbqQcPgG+/Fcdz5gAVKvCL\nh1guFxfAz08c01G/8lGP30wNGcJugwewddLDwsSjLkLkdv8+UL068Po1G+/fD/TowTcmQj1+Vblw\nQSz6AGv5UNEnPFWunH0Nn6+/Fk/4EuWhwi8zY3t7mZnAuHHiuFcvvldM8u5VSklNuQDy5+PvD7z3\nHtu+fh3YsUO696Z9Iy0q/GZm0yY2ZRMAihUDli7lGw8hOuXLAxMniuOZM8XWD1EW6vGbkaQkoGZN\n8ZaK33xDF80QZXn+nPX6dT+jS5YAkybxjcmSUY9fBebMEX+hnJ3ZR2tClKRUKXZAojNvHpCYyC8e\nkjsq/DIztLcXEQGsWiWOly4V10TniXevUkpqygXgl89nn7GjfoAV/QULjH9P2jfSosJvJiZOFO9v\n6u1Nq28S5SpaNHuxX7kSuHePXzzkXdTjNwO//QZ07Mi2rayAK1cAd3e+MRGSH0EAmjUDgoPZeMgQ\nYPNmriFZJOrxm6mMDGDyZHE8YgQVfaJ8Gg2waJE43roVuHqVXzwkOyr8Mitsb2/TJjYnGgBKllTe\nLB7evUopqSkXgH8+bdqIK8cKAjB1quHvxTsXqfHOhwq/gr14wa6A1PH3BypV4hcPIYUVEMDakwBw\n9Cjwxx984yEM9fgVbPZs8Qi/cmXgzh121E+IORk+HNi4kW03bMj6/lZ0yCmLvHr8VPgV6sEDwNUV\nSElh402bAF9friERYpDYWLZuv+4q3l27gH79+MZkKejkrkLo29tbvFgs+h4e2W9xpyS8e5VSUlMu\ngHLycXLKfsOgGTPYTYQKQym5SIV3PlT4FSg+Hli3Thx/9x1QpAi/eAgxlr8/YG/Ptv/5J/vPN5Ef\ntXoUaMoUdsQPsFsphoSw6XGEmLPvvwe+/JJtlysH/P03YGfHNya149LqOXbsWMfatWvfqlGjRuTC\nhQvfWVlGq9V629nZPfPy8rri5eV1Ze7cuTNNGY85ePwY+OEHcfzNN1T0iTp8/jm7WxcAJCRkn+dP\nZCYIgkke6enpRapXr343KirKJTU11cbDwyMsIiKiTtbXnDp1yrtr164HC3ovFqY6nDp1Kt+vT5sm\nCGzWsyC4uwtCRoY8cRmqoHzMiZpyEQRl5rNjh/jzXby4IMTG6vd9SszFGHLl87Z2vlNTTXbEHxwc\n3MTV1fWui4tLtI2NTVq/fv12HzhwoHsuf3joePatp0+zL8T29dc07Y2oS79+QIMGbPvVKzZlmcjP\n2lRvHBcX5+js7ByjGzs5OcVeunSpadbXaDQa4fz58y08PDzCHR0d45YsWTK5bt26Ebm9n6+vL1ze\nfk4sU6ZDG7EpAAASEElEQVQMPD094e3tDUA8Q24OY29v7zy//scf3njxAgC0qFoV6NWLf7zG5ENj\nGuccW1kBn36qRWgoAHhj40agZUstXFwK/n4dJeVjzNgU+Wi1Wmx+uyiSrl7mKrePAVI89u7d23vE\niBHrdeNt27YNHDt27Kqsr0lOTi718uXLEoIgICgoqFONGjXu5PZeUFGrJy+JiYJQurT4MXjXLt4R\nEWI6H38s/qx37co7GvWC3K0eR0fHuJiYGGfdOCYmxtnJySk262tKlSr1vESJEikA0KlTp6NpaWk2\nT58+tTdVTEqQ86+9zsqVQHIy265VC+jbV76YjJFXPuZITbkAys5n4UJx0sKhQ8DZs/m/Xsm5GIJ3\nPiYr/I0aNbocGRlZIzo62iU1NbXonj17PunWrdvBrK+Jj4+vKLzt8QcHBzcRBEFjb2//1FQxKVVy\nMpvqpjNzJs3bJ+rm6QkMHCiOp0xhx/9EJrl9DJDqERQU1KlmzZq3q1evfnf+/PnTBEHAunXrRq5b\nt26kIAhYvXr1mHr16l338PAIa968+fkLFy40y+19oPJWz7x54sdeV1dBSEvjHREhphcdLQhFi4o/\n+3v38o5IfZBHq4cu4OLs+XM2t/np2885tCYPsSSTJ7PbiAJsPZ8bNwAbG74xqQmt1aMQOXt7a9eK\nRb9aNWDAAPljMgbvXqWU1JQLYB75TJ8OlCnDtiMjgfXrc3+dOeRSGLzzocLP0cuXwJIl4nj6dDra\nIZbF3h6YNk0cz5kjTnIgpkOtHo6WLQMmTWLbVaqwI56iRfnGRIjcXr1iM9li3l71M3myuFYVMQ61\nehTm1avsa5VMm0ZFn1im4sXZnbp0li8HInK9jJNIhQq/zHS9vQ0b2PLLAODoCAwdyi8mY/DuVUpJ\nTbkA5pVP//5A69ZsOz0dGD8++/ROc8pFH7zzocLPwatX2Y9wpk4FbG35xUMIbxoNW6dKd/3KyZPA\nL7/wjUnNqMfPQdbefqVK7MYUxYvzjYkQJfjiC3YVOwA4OwM3b9J9po1BPX6FeP4cWLBAHM+cSUWf\nEJ05c4AKFdh2TAy7+xyRHhV+mX3xhRZPnrDtqlUBPz++8RiLd69SSmrKBTDPfMqUYev46CxeDPz1\nl3nmkh/e+VDhl9HTp8CePeJ49myayUNIToMHA23bsu3MTHYle2Fvzk7yRz1+GU2dKh7N1KoFXL8O\nWJvsjgiEmK+oKMDNjV3kCLDe//LlfGMyR9Tj5+zBA/GkFcB6l1T0CcldtWrZL+JasQI4fJhfPGpD\nhV8ms2ezaZyAFp6eQO/enAOSCO9epZTUlAtg/vmMGgV07aobaeHrC0RH84tHSrz3DRV+GYSFZV98\nav58upcuIQXRaNhqtY6ObJyQwP4QPH/ONy41oB6/iQkCO1F1+jQbd+wIBAWJdx8ihOTv/Hn2O6Q7\nwdulC/Drr3SzIn1Qj5+TPXvEom9tzS7eoqJPiP5atMj+ifnwYTZRghiOCr8JJSay2Qg6Y8cC8fFa\nbvGYAu9epZTUlAugrnyqVNFmK/ZLlgA//sgvHmPx3jdU+E3I3x949IhtOzqyqxIJIYaZNw/o3l0c\njxoF/O9//OIxZ9TjN5Hffwc6dBDH+/YBPXvyi4cQNXjxAmjTBggNZWNra+DAAcDHh29cSpVXj58K\nvwkkJgL16wP377Nxjx7A/v18YyJELR4/Zks437rFxsWKAceOsT8IJDs6uSsTQQBGjhSLfvnyQGCg\n+HXevT2pqSkfNeUCqCufrLmULw+cOAG4uLDx69dsps+FC1xCMwjvfUOFX2IrVwI//yyO168XVxsk\nhEjD0ZEV/0qV2PjFC9ZavXiRb1zmglo9EvrzTzbfOD2djUePBtas4RsTIWp24wb7nXv8mI1LlwaO\nHweaNOEbl1JQj9/EIiOBZs3YCpwA+8E7c4burEWIqV2/zoq/brlzOztW/Bs35huXElCP34Ti49ms\nAl3Rr1CBTTPLrejz7u1JTU35qCkXQF355JdL/frAH38A5cqx8bNnQPv2wOXL8sRmCN77hgq/kRIS\n2A/Z3btsXKwYcPAgu8kKIUQebm7sPr05i/9ff/GNS6mo1WOE+Hi29k5YGBtbWQF799J8fUJ4CQsD\n2rUTP32XKsWWeGjdmm9cvFCrR2KRkWwNEV3R12iALVuo6BPCk6cnm+2jO/J//pwdnB07xjcupaHC\nb4DLl4GWLYF//mFjKytg40Zg4MCCv5d3b09qaspHTbkA6sqnMLl4ebGFEXVTPV+9Arp1Y5/GlYL3\nvqHCX0i//w54e4vTx4oVY0vE+vryjIoQklW9esDZs+K5trQ04L//BRYtYhdZWjrq8RfCjh2swOvm\n6Zcty/qHLVpwDYsQkoeYGOCjj4A7d8Tnhg4F1q0DihblF5dcqMdvpKVLWStHV/SdnYFz56joE6Jk\nut/TVq3E5zZtAj7+GHj4kF9cvFHhL4AgAF99BUyeLD5Xrx67K1CdOoV/P969PampKR815QKoKx9j\ncnFwYBd0DRkiPnf6NJsCeuCA8bEZgve+ocKfj8xMtub3kiXicx9+yHqHTk784iKEFI6tLTvSDwgQ\nn3vyhK2c6+fH1vqxJNTjz0N6OusFbt8uPtejB7BzJ1C8uKyhEEIkdPIkO/qPixOfc3ICFiwAPv2U\nzdJTC+rxF0JqKvDJJ9mL/sCBbNVNKvqEmLd27YCrV9ksH53YWGDQIKB5c9bGVTsq/Dm8esWO7Pft\nE58bOZJdnGVtbfz78+7tSU1N+agpF0Bd+Uidi709sHs3O7irWFF8PjiYXaPTpYtpl3vgvW+o8Gfx\n/DnQuTNw9Kj43MSJwNq16vr4RwhhV9sPGMCuwp82LfuiikeOsBV2Tf0HgBfq8b+VkAB06pR9J3/9\nNbtBuuadDhkhRG2io4Hp09kngZzlpnNnYNYs81vqmdbjz0dcHJvXGxEhPhcQAPj7m+yfJIQo1M2b\nwHffqeMPAJ3czUNkJOvp6Yq+RsOu6jNV0efd25OamvJRUy6AuvKRM5c6ddjsvRs3gP79s3/i17WA\nunY1br1/3vvGogv/8eNA06bAv/+ysY0NsGsXO5lLCLFsuj8A168D/fpl/wNw+DA76jf2DwAvFtnq\nSU8H5s0Dvv2WXaQFsGma+/axJVwJISSniAjWAtqz590WUJcuwDffKK8FRD3+ty5fBsaMYdO2dCpX\nZitsKm2nEUKUJ78/AC1aAJ9/DvTpo4z7bVt8j//aNTZ1q3Hj7EW/VSv2x0Cuos+7tyc1NeWjplwA\ndeWjpFzq1mUt4WvX2IWeWVtA58+ziz2dnYGxY9nU8Fev3n0P3vmouvDHxwMbNrDbrrm7s36dTtGi\n7BLtU6eA99/nFyMhxDzVq8dm/ly7xpZ6sLERv/b4MfDDD4CPD7sbWKdOwMyZrJ0cFSW2mHlRTasn\nMxO4dQu4dIk9Ll4EwsNzf23XrsCyZYCrqwmCJYRYpIcPgZ9+YrMCY2Pzf22xYkD16kCNGuzh6ipu\nV64s3QWjquvxP3okFvlLl1j7Jjk57/coUgTo3p1N02zSxMQBE0IsVno6cOYMm/p55Ahw+3bhvr94\ncfZHoVq1dx8uLkDp0vq/l9kX/iNHBISGAiEhQGgocO9ewd9XpAhbRrlHDzYfN+uaHLxotVp4e3vz\nDkMyaspHTbkA6srHnHP5+292cHrlCqtd4eFAQoIWgLdB71e2LGtPV6rEHhUritsVKrCvly0LlCkD\nVKyYe+GXYNmxvB07dqzjhAkTlmdkZBQZMWLEBn9//4U5XzN+/PiVR48e7VSiRImUzZs3+3p5eV3J\n7b06dy7436tQgc3L1z0aNwbs7IzPQ0phYWFm+wOcGzXlo6ZcAHXlY865VK/OHp9+Kj43f34Y2rf3\nRmQkcPcuu5BUt52QkP/7JSayR9aVBgrLZIU/IyOjyNixY1efOHHiI0dHx7jGjRv/1a1bt4N16tS5\nqXtNUFCQz927d10jIyNrXLp0qeno0aPXXrx4sZk+729rC3h5Ac2aiYXexUX56+okJSXxDkFSaspH\nTbkA6spHTbkAQGpqEho3zn024dOn7FNCdDQ7EZz1ER3Nlo03lskKf3BwcBNXV9e7Li4u0QDQr1+/\n3QcOHOietfAfPHiw25AhQ7YAQNOmTS8lJSWViY+Pr1ixYsX4nO/XogXQoIH4qFs3+1l0QghRA3t7\n9sjtj0JmJju/GR/PTibn/O/jx+zTQFIS+29enx5MVvjj4uIcnZ2dY3RjJyen2EuXLjUt6DWxsbFO\nuRX+c+dMFam8oqOjeYcgKTXlo6ZcAHXlo6ZcAMPzsbIS+/keHgW/Pq8OiMkKv0aj0euscc4TD3l9\nn0bpPZxC2LJlC+8QJKWmfNSUC6CufNSUC8A3H5MVfkdHx7iYmBhn3TgmJsbZyckpNr/XxMbGOjk6\nOsYhh9zOShNCCDGMya7cbdSo0eXIyMga0dHRLqmpqUX37NnzSbdu3Q5mfU23bt0Obt26dTAAXLx4\nsVmZMmWScmvzEEIIkY7Jjvitra3TV69ePbZDhw6/ZWRkFBk+fPhPderUuRkYGDgSAEaOHBno4+MT\nFBQU5OPq6nq3ZMmSLzdt2jTUVPEQQgh5SxAExTyOHj3asVatWrdcXV0jAwIC/HN+ffv27QPc3d3D\n3dzcrrZo0eJceHi4O++YDc1F9wgODm5cpEiR9F9++aUX75iNyeXUqVPenp6eV+rVq3e9TZs2Wt4x\nG5PP48ePHTp06HDMw8MjrF69etc3bdrkyzvmvB5Dhw7dWKFChfj69etfy+s148aNW+nq6hrp7u4e\nHhoa6sU7ZkNzMafff333jSDwqQHc/+foHunp6UWqV69+NyoqyiU1NdXGw8MjLCIiok7W15w/f755\nUlKSnSCwX96mTZte5B23obnoXte2bds/OnfufHjv3r29ecdtaC6JiYll6tateyMmJsZJEFjh5B23\nMfnMmjVr9tSpUxfocrG3t09IS0uz5h17bo8zZ860Cg0N9cqruBw5csSnU6dOQYIg4OLFi02V+juj\nTy7m8vuvbz6CwK8GKGZ1zqzz/m1sbNJ08/6zvqZ58+YX7OzsngFs3n9sbKwTn2jzp08uALBq1apx\nffr02Vu+fPnHPOLUhz657Ny589PevXv/ojt57+Dg8IRPtAXTJ5/333//QXJycmkASE5OLl2uXLkE\na2vrdD4R569Vq1Zny5Ytm5jX1/O6Vka+CPVXUC7m8vuvU1A+AL8aoJjCn9uc/ri4OMe8Xv/TTz8N\n9/HxCZInusLRJ5e4uDjHAwcOdB89evRaQP/pr3LTJ5fIyMgaT58+tW/btu2pRo0aXd62bdsg+SPV\njz75+Pn5rb9x40a9ypUr3/fw8AhfsWLFF/JHKo28rpXhGZMUlPz7ry+eNcCka/UURmGSPnXqVNuN\nGzcOO3fuXEtTxmQofXKZMGHC8oCAgKlvVx7VCAqdsqpPLmlpaTahoaENTp482S4lJaVE8+bNLzRr\n1uxijRo1IuWIsTD0yWf+/PnTPT09w7Rarffff/9dvX379sfDw8M9SpUq9VyOGKWW82dLqQcZ+lL6\n77++eNYAxRR+feb9A8DVq1fd/fz81h87dqxjQR+jeNEnl5CQkIb9+vXbDQBPnjxxOHr0aCcbG5u0\nnFNeedMnF2dn5xgHB4cnxYsXf1W8ePFXrVu3PhMeHu6hxMKvTz7nz59vMWPGjHkAUL169b+rVasW\ndfv27VqNGjUyu9tq63utjLkwh99/fXGtAbxPgOgeaWlp1h988MHfUVFRLm/evCma20m3f//9t0r1\n6tXvXrhwoRnveI3NJevD19d3k1Jn9eiTy82bN2u3a9fuRHp6epGXL1+WqF+//rUbN27U5R27oflM\nnDhx2ezZs2cJgoCHDx9WdHR0jE1ISLDnHXtej6ioKBd9Tu5euHChmdJPiOaXi7n8/uubT9aH3DVA\nMUf8+sz7//bbb79JTEwsq+uJ2djYpAUHByvutir65MI7Rn3pk0vt2rVvdezY8Zi7u/tVKyurTD8/\nv/V169Y1YtFY09Enn+nTp88fOnToJg8Pj/DMzEyrRYsWTbG3t3/KO/bc9O/ff9fp06fbPHnyxMHZ\n2Tlmzpw5s9LS0mwA87tWpqBczOX3X6egfHjGZhY3YiGEECIdxczqIYQQIg8q/IQQYmGo8BNCiIWh\nwk8IIRaGCj8hhFgYKvyEEGJh/g81XEuNi6738AAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 191 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }