{ "metadata": { "name": "", "signature": "sha256:583ab56294995ef0cf4b6dd3627bbede6d2ff040c58107eedd471d6f1bc3c36f" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "#RNA-Seq data corresponding to BS-Seq Oyster Sperm sample \n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Shortcuts \n", "\n", "- [RNAseq on genes](#RNA-seq-on-Genes-Only)\n", "- [March2014](#FastQC-in-iPlant)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "Data\n", "\n", "RNA-seq data from both the sperm and pooled gill samples were provided by the core facility. \n", "Originally BAM files were provided, later followed up with fastq files.\n", "\n", "\"72975_GACTAAGA%20et%20al%20and%20(274)%20Discovery%20Environment\"\n", "\n", "\n", "\n", "\n", "Sperm data = 72976/GTGTCTAC\n", "\n", "\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "##CLC \n", "fastq import (s_1.GTGTCTAC_1 (paired)-1) \n", "51,084,360 sequences\n", "\n", "`\n", " Discard read names = Yes\n", " Discard quality scores = No\n", " Paired orientation = Paired reads (forward-reverse)\n", " Minimum distance = 180\n", " Maximum distance = 250\n", " Original sequence resource = /Volumes/NGS Drive/NGS Raw Data/RNA_complement_BS/C25GNACXX.oyster.fastqs/s_1.GTGTCTAC_1.fastq\n", " Original sequence resource = /Volumes/NGS Drive/NGS Raw Data/RNA_complement_BS/C25GNACXX.oyster.fastqs/s_1.GTGTCTAC_2.fastq\n", " Quality score = NCBI/Sanger or Illumina Pipeline 1.8 and later\n", " MiSeq de-multiplexing = No\n", "`\n", "\n", "####QC Report\n", "\n", "\n", "####Supp QC Report\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "Trimming \n", "\n", "\"CLC%20Genomics%20Workbench%206.0.2\"\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "---\n", "\n", "###RNA-seq \n", "\n", "Input file\n", "\n", "\"CLC%20Genomics%20Workbench%206.0.2\" \n", "\n", "\n", "\n", "*Need to modify GFF so that CLC recognizes to annotate*\n", "\n", "" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/oyster.v9.gene_mRNA.gff" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "C16582\tGLEAN\tgene\t35\t385\t0.555898\t-\t.\tID=CGI_10000001;\r\n", "C16582\tGLEAN\tmrna\t35\t385\t.\t-\t0\tParent=CGI_10000001;\r\n", "C17212\tGLEAN\tgene\t31\t363\t0.999572\t+\t.\tID=CGI_10000002;\r\n", "C17212\tGLEAN\tmrna\t31\t363\t.\t+\t0\tParent=CGI_10000002;\r\n", "C17316\tGLEAN\tgene\t30\t257\t0.555898\t+\t.\tID=CGI_10000003;\r\n", "C17316\tGLEAN\tmrna\t30\t257\t.\t+\t0\tParent=CGI_10000003;\r\n", "C17476\tGLEAN\tgene\t34\t257\t0.998947\t-\t.\tID=CGI_10000004;\r\n", "C17476\tGLEAN\tmrna\t104\t257\t.\t-\t0\tParent=CGI_10000004;\r\n", "C17476\tGLEAN\tmrna\t34\t74\t.\t-\t2\tParent=CGI_10000004;\r\n", "C17998\tGLEAN\tgene\t196\t387\t1\t-\t.\tID=CGI_10000005;\r\n" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/oyster.v9.gene_mRNA.gff" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 224718 2022462 14376214 /Volumes/web/cnidarian/oyster.v9.gene_mRNA.gff\r\n" ] } ], "prompt_number": 4 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Annotated version 9 of oyster genome then ran RNA-seq with Exon Discovery and Expression value as transcripts. Only using paired data.\n", "\n", "\n", "\"RNA-Seq%20Analysis\"\n", "\n", "\n", "##_Concerned about the number of mapped reads_ \n", "\n", "\n", "\"CLC_Genomics_Workbench_6.0.2_179F262E.png\"/\n", "\n", "Full Report \n", "\n", "\n" ] }, { "cell_type": "heading", "level": 2, "metadata": {}, "source": [ "Redux on RNA-seq on annotated genome" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#September 12, 2013" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"RNA-Seq_Analysis_17E25F01.png\"/" ] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"CLC_Genomics_Workbench_6.0.2_17AB0078.png\"/" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#export of transcript-based table\n", "!wc /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 196692 2753710 15000644 /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt\r\n" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\"Feature ID\"\t\"Expression values\"\t\"Gene name\"\t\"Transcripts annotated\"\t\"Transcript length\"\t\"Transcript ID\"\t\"Unique transcript reads\"\t\"Total transcript reads\"\t\"Ratio of unique to total (transcript reads)\"\t\"Exons\"\t\"RPKM\"\t\"Relative RPKM\"\t\"Chromosome\"\t\"Chromosome region start\"\t\"Chromosome region end\"\r\n", "CGI_10000780.1\t0\tCGI_10000780\t2\t1248\t\t0\t0\t?\t1\t0\t1\tscaffold350\t1959\t3206\r\n", "CGI_10000780.2\t0\tCGI_10000780\t2\t102\t\t0\t0\t?\t1\t0\t1\tscaffold350\t1105\t1206\r\n", "CGI_10000456.1\t0\tCGI_10000456\t4\t97\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t349\t445\r\n", "CGI_10000456.2\t0\tCGI_10000456\t4\t151\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t943\t1093\r\n", "CGI_10000456.3\t0\tCGI_10000456\t4\t124\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t1481\t1604\r\n", "CGI_10000456.4\t0\tCGI_10000456\t4\t66\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t1973\t2038\r\n", "CGI_10000457.1\t0\tCGI_10000457\t7\t94\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t5599\t5692\r\n", "CGI_10000457.2\t0\tCGI_10000457\t7\t138\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t5175\t5312\r\n", "CGI_10000457.3\t0\tCGI_10000457\t7\t50\t\t0\t0\t?\t1\t0\t1\tscaffold28720\t4978\t5027\r\n" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "from pandas import *\n", "\n", "# read data from data file into a pandas DataFrame \n", "BiGoRNAseq_exon = read_table(\"/Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt\", # name of the data file\n", " #sep=\",\", # what character separates each column?\n", " na_values=[\"\", \" \"]) # what values should be considered \"blank\" values?" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "IOError", "evalue": "File /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt does not exist", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mIOError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[1;32m 4\u001b[0m BiGoRNAseq_exon = read_table(\"/Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt\", # name of the data file\n\u001b[1;32m 5\u001b[0m \u001b[0;31m#sep=\",\", # what character separates each column?\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 6\u001b[0;31m na_values=[\"\", \" \"]) # what values should be considered \"blank\" values?\n\u001b[0m", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/io/parsers.pyc\u001b[0m in \u001b[0;36mparser_f\u001b[0;34m(filepath_or_buffer, sep, dialect, compression, doublequote, escapechar, quotechar, quoting, skipinitialspace, lineterminator, header, index_col, names, prefix, skiprows, skipfooter, skip_footer, na_values, true_values, false_values, delimiter, converters, dtype, usecols, engine, delim_whitespace, as_recarray, na_filter, compact_ints, use_unsigned, low_memory, buffer_lines, warn_bad_lines, error_bad_lines, keep_default_na, thousands, comment, decimal, parse_dates, keep_date_col, dayfirst, date_parser, memory_map, nrows, iterator, chunksize, verbose, encoding, squeeze)\u001b[0m\n\u001b[1;32m 399\u001b[0m buffer_lines=buffer_lines)\n\u001b[1;32m 400\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 401\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0m_read\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilepath_or_buffer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 402\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 403\u001b[0m \u001b[0mparser_f\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__name__\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/io/parsers.pyc\u001b[0m in \u001b[0;36m_read\u001b[0;34m(filepath_or_buffer, kwds)\u001b[0m\n\u001b[1;32m 207\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 208\u001b[0m \u001b[0;31m# Create the parser.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 209\u001b[0;31m \u001b[0mparser\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mTextFileReader\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilepath_or_buffer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 210\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 211\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mnrows\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/io/parsers.pyc\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, f, engine, **kwds)\u001b[0m\n\u001b[1;32m 507\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moptions\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'has_index_names'\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mkwds\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'has_index_names'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 508\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 509\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_make_engine\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mengine\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 510\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 511\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_get_options_with_defaults\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mengine\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/io/parsers.pyc\u001b[0m in \u001b[0;36m_make_engine\u001b[0;34m(self, engine)\u001b[0m\n\u001b[1;32m 609\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m_make_engine\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mengine\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'c'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 610\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mengine\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'c'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 611\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_engine\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mCParserWrapper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moptions\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 612\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 613\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mengine\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'python'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/io/parsers.pyc\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, src, **kwds)\u001b[0m\n\u001b[1;32m 891\u001b[0m \u001b[0;31m# #2442\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 892\u001b[0m \u001b[0mkwds\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'allow_leading_cols'\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mindex_col\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 893\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_reader\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0m_parser\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mTextReader\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0msrc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwds\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 894\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 895\u001b[0m \u001b[0;31m# XXX\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/_parser.so\u001b[0m in \u001b[0;36mpandas._parser.TextReader.__cinit__ (pandas/src/parser.c:2771)\u001b[0;34m()\u001b[0m\n", "\u001b[0;32m/Users/sr320/anaconda/lib/python2.7/site-packages/pandas/_parser.so\u001b[0m in \u001b[0;36mpandas._parser.TextReader._setup_parser_source (pandas/src/parser.c:4803)\u001b[0;34m()\u001b[0m\n", "\u001b[0;31mIOError\u001b[0m: File /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_1.txt does not exist" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "BiGoRNAseq_exon.dtypes" ], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "code", "collapsed": false, "input": [ "print BiGoRNAseq_exon" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'BiGoRNAseq_exon' is not defined", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;32mprint\u001b[0m \u001b[0mBiGoRNAseq_exon\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'BiGoRNAseq_exon' is not defined" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "BiGoRNAseq_exon['RPKM'].hist(bins=100);\n", "#Axis limits are changed using the axis([xmin, xmax, ymin, ymax]) function.\n", "axis([0, 1000, 0, 1000000])\n" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'BiGoRNAseq_exon' is not defined", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mBiGoRNAseq_exon\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'RPKM'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mhist\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbins\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m100\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m;\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2\u001b[0m \u001b[0;31m#Axis limits are changed using the axis([xmin, xmax, ymin, ymax]) function.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 3\u001b[0m \u001b[0maxis\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m1000\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m1000000\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mNameError\u001b[0m: name 'BiGoRNAseq_exon' is not defined" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "BiGoRNAseq_exon(x='Transcripts annotated', y='Transcripts length', style='o');" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "TypeError", "evalue": "'DataFrame' object is not callable", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mTypeError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mBiGoRNAseq_exon\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mx\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Transcripts annotated'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0my\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'Transcripts length'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstyle\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'o'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m;\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mTypeError\u001b[0m: 'DataFrame' object is not callable" ] } ], "prompt_number": 44 }, { "cell_type": "markdown", "metadata": {}, "source": [ "--- \n", "###Upload to SQLShare\n", "\n", "Query to get Genome track\n", "\n", "```\n", "SELECT\n", " [\"Chromosome\"] as seqname,\n", " 'CLC_RNAseq' as source,\n", " 'TranscriptExp' as feature, \n", " [\"Chromosome region start\"] as start,\n", " [\"Chromosome region start\"] as [end],\n", " [\"RPKM\"] as score,\n", " '.' as strand,\n", " '.' as frame,\n", " [\"Gene name\"] as attribute\n", " FROM [sr320@washington.edu].[table_BiGoRNAseq_exon_exp_1.txt]\u200b\n", "``` \n", "\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "seqname,source,feature,start,end,score,strand,frame,attribute\r", "\r\n", "C16582,CLC_RNAseq,TranscriptExp,35,35,31.6,.,.,CGI_10000001\r", "\r\n", "C17212,CLC_RNAseq,TranscriptExp,31,31,2.776,.,.,CGI_10000002\r", "\r\n", "C17316,CLC_RNAseq,TranscriptExp,30,30,0,.,.,CGI_10000003\r", "\r\n", "C17476,CLC_RNAseq,TranscriptExp,217,217,0,.,.,CGI_10000004\r", "\r\n", "C17476,CLC_RNAseq,TranscriptExp,34,34,0,.,.,CGI_10000004\r", "\r\n", "C17998,CLC_RNAseq,TranscriptExp,196,196,0,.,.,CGI_10000005\r", "\r\n", "C18346,CLC_RNAseq,TranscriptExp,174,174,14.671,.,.,CGI_10000009\r", "\r\n", "C18428,CLC_RNAseq,TranscriptExp,286,286,240.811,.,.,CGI_10000010\r", "\r\n", "C18964,CLC_RNAseq,TranscriptExp,203,203,0,.,.,CGI_10000011\r", "\r\n" ] } ], "prompt_number": 49 }, { "cell_type": "code", "collapsed": false, "input": [ "!tr ',' \"\\t\" /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track.gff" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 52 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track.gff" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "seqname\tsource\tfeature\tstart\tend\tscore\tstrand\tframe\tattribute\r", "\r\n", "C16582\tCLC_RNAseq\tTranscriptExp\t35\t35\t31.6\t.\t.\tCGI_10000001\r", "\r\n", "C17212\tCLC_RNAseq\tTranscriptExp\t31\t31\t2.776\t.\t.\tCGI_10000002\r", "\r\n", "C17316\tCLC_RNAseq\tTranscriptExp\t30\t30\t0\t.\t.\tCGI_10000003\r", "\r\n", "C17476\tCLC_RNAseq\tTranscriptExp\t217\t217\t0\t.\t.\tCGI_10000004\r", "\r\n", "C17476\tCLC_RNAseq\tTranscriptExp\t34\t34\t0\t.\t.\tCGI_10000004\r", "\r\n", "C17998\tCLC_RNAseq\tTranscriptExp\t196\t196\t0\t.\t.\tCGI_10000005\r", "\r\n", "C18346\tCLC_RNAseq\tTranscriptExp\t174\t174\t14.671\t.\t.\tCGI_10000009\r", "\r\n", "C18428\tCLC_RNAseq\tTranscriptExp\t286\t286\t240.811\t.\t.\tCGI_10000010\r", "\r\n", "C18964\tCLC_RNAseq\tTranscriptExp\t203\t203\t0\t.\t.\tCGI_10000011\r", "\r\n" ] } ], "prompt_number": 53 }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track.gff" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 196692 1770228 14001279 /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track.gff\r\n" ] } ], "prompt_number": 54 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track2.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "seqname,start,end,feature,score\r", "\r\n", "C16582,35,385,tr_expression,31.6\r", "\r\n", "C17212,31,363,tr_expression,2.776\r", "\r\n", "C17316,30,257,tr_expression,0\r", "\r\n", "C17476,217,257,tr_expression,0\r", "\r\n", "C17476,34,187,tr_expression,0\r", "\r\n", "C17998,196,387,tr_expression,0\r", "\r\n", "C18346,174,551,tr_expression,14.671\r", "\r\n", "C18428,286,546,tr_expression,240.811\r", "\r\n", "C18964,203,658,tr_expression,0\r", "\r\n" ] } ], "prompt_number": 60 }, { "cell_type": "code", "collapsed": false, "input": [ "!tr ',' \"\\t\" /Volumes/web/cnidarian/BiGoRNAseq_exon_exp_track2.igv" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 61 }, { "cell_type": "markdown", "metadata": {}, "source": [ "IGV file based on exons / transcripts\n", "\n", "" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#WHY not just run bedtools on Accepted hit bam?\n", "\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 62 }, { "cell_type": "heading", "level": 2, "metadata": {}, "source": [ "RNA-seq on Genes Only" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"CLC_Genomics_Workbench_6.5_17E2669D.png\"/" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGo_RNAseq_genes" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\"Feature ID\"\t\"Expression values\"\t\"Gene length\"\t\"Unique gene reads\"\t\"Total gene reads\"\t\"RPKM\"\r\n", "CGI_10000780\t0\t1350\t0\t0\t0\r\n", "CGI_10000456\t7.892\t438\t8\t9\t7.892\r\n", "CGI_10000457\t7.643\t603\t6\t12\t7.643\r\n", "CGI_10000774\t0\t375\t0\t0\t0\r\n", "CGI_10000917\t0\t426\t0\t0\t0\r\n", "CGI_10000861\t0\t2004\t0\t0\t0\r\n", "CGI_10000994\t16.913\t1635\t64\t72\t16.913\r\n", "CGI_10000643\t0.696\t552\t1\t1\t0.696\r\n", "CGI_10000763\t0\t249\t0\t0\t0\r\n" ] } ], "prompt_number": 6 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "\n", "---\n", "\n", "\n", "\n", "##RefMap (to see what proportion map)\n", "\n", "\"greenbird_17A05407.png\"/\n", "\n", "Report \n", "\n", "\n", " \n", "Only about 50% of reads map?\n", " \n", "Will run again an try to assemble unmapped reads to determine what they are..\n", " \n", "Running- reduced stringency to default, ignore vs random, kept unmapped reads, and output set as track \n", "\n", "\"CLC_Genomics_Workbench_6.0.2_17A1B50E.png\"/ \n", " \n", "\"CLC_Genomics_Workbench_6.0.2_17A1B4E8.png\"/ \n", " \n", " \n", "---\n", "\n", " \n", " \n" ] }, { "cell_type": "heading", "level": 2, "metadata": {}, "source": [ "REDUX on Reference assembly" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#September 11 2013" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"CLC_Genomics_Workbench_6.5_17E0D7AD.png\"/" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Big Difference! Seems to be the fact that mapped non-specific reads randomly...." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"CLC_Genomics_Workbench_6.5_17E0D854.png\"/" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "--- \n", "\n", "##De novo assembly of reads that did not map\n", " \n", "Fast de novo on 21M unmapped \n", " \n", " \n", "\"CLC_Genomics_Workbench_6.5_Beta_3_17A6B13D.png\"/\n", " " ] }, { "cell_type": "code", "collapsed": false, "input": [ "!tail /Volumes/web/cnidarian/BiGo_RNAseq_unmapped.fa" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ ">s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_3093\r\n", "GTGGTGCTCGGTCTGTATAACATCTGGACAGACTCCTGTTTTGTCAATCAGTTCTTGTAT\r\n", "GTGGTGTTGAAGAAAGTTAGAAACTCCAATATTCTTTACTTTACCTTCCTTCTGTAACTT\r\n", "GATCAGATCCCGCCAGCTCCCCACCCTCAGCTCCCCCTGATCTGGGTCATCTGGTTTCTT\r\n", "GCCCTGGGTCCCTGGCCAGTGAATCAGGTACAGGTCCAGGTACTGGCAGTCTAAGTTAG\r\n", ">s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_3094\r\n", "GGGCTGCACAGATCTCTCAACAATCTCTTTCGTGAGAGTTAAGCGACTTTTAGGTTTGTT\r\n", "GTACATAACTGACAACAAACAGACAAAAGAACTTGGATTTTGGCACGGGGTTTCACTGTG\r\n", "TTACAAAGATCATGATACATTTTTGAGCAAATATTTTCACAACTGCATGAAAAATACATG\r\n", "TTTTTAACCGTTGGTGCCCTTGGTTGCATAGGTTGCCGCCATTTTGAAACATAGGAAGC\r\n" ] } ], "prompt_number": 67 }, { "cell_type": "code", "collapsed": false, "input": [ "!blastx -query /Volumes/web/cnidarian/BiGo_RNAseq_unmapped.fa -db /Volumes/web/whale/fish546/blast/db/swissprot -out /Volumes/web/cnidarian/BiGoRNAseq_unmapped_swissprot_blastout -outfmt 6 -evalue 1E-10 -max_target_seqs 1 -num_threads 2" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 690 replaced by X\r\n", "Selenocysteine (U) at position 690 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 690 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 140 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 196 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 204 replaced by X\r\n", "Selenocysteine (U) at position 196 replaced by X\r\n", "Selenocysteine (U) at position 193 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 493 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 24 replaced by X\r\n", "Selenocysteine (U) at position 63 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 25 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 63 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 60 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 73 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "Selenocysteine (U) at position 43 replaced by X\r\n", "Selenocysteine (U) at position 73 replaced by X\r\n" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_unmapped_swissprot_blastout" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_1\tgi|3121895|sp|Q37705.1|COX1_ARTSF\t55.23\t507\t224\t2\t4105\t2594\t5\t511\t2e-140\t 457\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_4\tgi|3914777|sp|O61463.1|RLA2_CRYST\t56.25\t112\t47\t1\t392\t57\t1\t110\t9e-21\t86.7\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_5\tgi|54041577|sp|P13418.2|POLS_CRPVC\t22.81\t754\t440\t25\t2769\t787\t43\t747\t1e-25\t 118\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_6\tgi|135489|sp|P11833.1|TBB_PARLI\t99.29\t420\t3\t0\t1425\t166\t11\t430\t0.0\t 839\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_7\tgi|341941223|sp|P29341.2|PABP1_MOUSE\t85.23\t88\t13\t0\t1063\t800\t536\t623\t2e-38\t 150\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_8\tgi|2500286|sp|Q26481.3|RL5_STYCL\t75.93\t295\t70\t1\t1168\t284\t1\t294\t1e-134\t 395\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_9\tgi|109940067|sp|P00428.2|COX5B_BOVIN\t52.94\t51\t23\t1\t430\t281\t72\t122\t7e-13\t65.9\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_10\tgi|71151870|sp|Q8TC29.1|ENKUR_HUMAN\t49.79\t235\t112\t2\t891\t187\t26\t254\t1e-67\t 218\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_13\tgi|115305840|sp|Q2NKZ1.1|TCPH_BOVIN\t82.76\t116\t20\t0\t1089\t742\t299\t414\t3e-108\t 212\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_14\tgi|6094094|sp|O57592.3|RL7A_FUGRU\t77.08\t240\t53\t1\t770\t51\t29\t266\t1e-125\t 366\r\n" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/BiGoRNAseq_unmapped_swissprot_blastout" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 1600 19200 271972 /Volumes/web/cnidarian/BiGoRNAseq_unmapped_swissprot_blastout\r\n" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "##Blast againa and get taxonomy info\n", "# ./blastn -query /Volumes/web/cnidarian/BiGo_RNAseq_unmapped.fa -db /Volumes/CLC_blastdatabases/nt -out /Volumes/web/cnidarian/BiGo_RNAseq_unmapped_nt_blastout_taxa2 -outfmt \"6 std stitle staxids sscinames scomnames sblastnames\" -evalue 1E-20 -max_target_seqs 1 -task blastn -num_threads 6" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "^C\r\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGo_RNAseq_unmapped_nt_blastout_taxa2" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_1\tgi|6636083|gb|AF177226.1|AF177226\t99.54\t11122\t47\t2\t1476\t12593\t18224\t7103\t0.0\t1.982e+04\tCrassostrea gigas mitochondrial DNA, complete genome\t29159\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_1\tgi|6636083|gb|AF177226.1|AF177226\t99.86\t1474\t2\t0\t2\t1475\t1474\t1\t0.0\t2650\tCrassostrea gigas mitochondrial DNA, complete genome\t29159\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_2\tgi|187762791|gb|EU672831.1|\t99.56\t901\t4\t0\t4\t904\t8445\t7545\t0.0\t1608\tCrassostrea gigas isolate ORCg-4 mitochondrion, complete genome\t29159\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_3\tgi|375073746|gb|JF744732.1|\t72.50\t240\t49\t5\t36\t264\t326\t93\t1e-26\t 129\tOstrea edulis clone M13-F.OEI1 CD650310-like protein mRNA, partial sequence\t37623\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_4\tgi|212374119|emb|FM165434.1|\t74.27\t342\t82\t2\t54\t392\t339\t1\t5e-53\t 217\tMytilus galloprovincialis mRNA for ribosomal protein (p2 gene)\t29158\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_6\tgi|194068374|dbj|AB374929.1|\t95.00\t1260\t63\t0\t166\t1425\t1356\t97\t0.0\t1988\tSaccostrea kegaki mRNA for beta-tubulin, complete cds\t182713\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_6\tgi|194068374|dbj|AB374929.1|\t74.74\t190\t12\t4\t5\t158\t1613\t1424\t1e-29\t 141\tSaccostrea kegaki mRNA for beta-tubulin, complete cds\t182713\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_7\tgi|224081794|ref|XM_002196958.1|\t77.39\t283\t61\t1\t800\t1079\t1857\t1575\t1e-53\t 221\tPREDICTED: Taeniopygia guttata similar to poly A binding protein, cytoplasmic 4 (LOC100228528), mRNA\t59729\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_8\tgi|40642989|emb|AJ563456.1|\t99.76\t829\t2\t0\t367\t1195\t829\t1\t0.0\t1487\tCrassostrea gigas partial mRNA for ribosomal protein L5 (rpl5 gene)\t29159\tN/A\tN/A\tN/A\r\n", "s_1.GTGTCTAC_1_(paired)-1_trimmed_(paired)_un-mapped_reads_[no_read_group]_(single)_contig_10\tgi|259013499|ref|NM_001165022.1|\t69.58\t503\t147\t3\t224\t723\t863\t364\t2e-50\t 210\tSaccoglossus kowalevskii enkurin, TRPC channel interacting protein (enkur), mRNA\t10224\tN/A\tN/A\tN/A\r\n" ] } ], "prompt_number": 69 }, { "cell_type": "code", "collapsed": false, "input": [ "###Get Taxonomic Distribution\n", "#SQLSshare\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 70 }, { "cell_type": "markdown", "metadata": {}, "source": [ "```\n", "SELECT [Column14], count([Column1])\n", " From\n", "[sr320@washington.edu].[BiGo_RNAseq_unmapped_nt_blastout_taxa2]\n", "Group by Column14\n", "``\n", " \n", " " ] }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNAseq_taxid_unmapped.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Column14,\r", "\r\n", "7091,3\r", "\r\n", "345175,1\r", "\r\n", "198214,2\r", "\r\n", "9601,1\r", "\r\n", "6945,7\r", "\r\n", "358574,2\r", "\r\n", "100272,1\r", "\r\n", "45954,2\r", "\r\n", "1132475,1\r", "\r\n" ] } ], "prompt_number": 72 }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/BiGoRNAseq_taxid_unmapped.csv" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 142 142 1308 /Volumes/web/cnidarian/BiGoRNAseq_taxid_unmapped.csv\r\n" ] } ], "prompt_number": 73 }, { "cell_type": "code", "collapsed": false, "input": [ "from pandas import *\n", "\n", "# read data from data file into a pandas DataFrame \n", "BiGoRNAsequm = read_csv(\"http://eagle.fish.washington.edu/cnidarian/BiGoRNAseq_taxid_unmapped.csv\", # name of the data file\n", " sep=\",\", # what character separates each column?\n", " na_values=[\"\", \" \"]) # what values should be considered \"blank\" values?" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 74 }, { "cell_type": "code", "collapsed": false, "input": [ "print BiGoRNAsequm" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\n", "Int64Index: 141 entries, 0 to 140\n", "Data columns:\n", "Column14 141 non-null values\n", "Unnamed: 1 141 non-null values\n", "dtypes: int64(2)" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "\n" ] } ], "prompt_number": 76 }, { "cell_type": "code", "collapsed": false, "input": [ "BiGoRNAsequm['Unnamed: 1'].hist(bins=50);\n", "#Axis limits are changed using the axis([xmin, xmax, ymin, ymax]) function.\n", "#plt.axis([0, 1, 0, 40])" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 83, "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAD9CAYAAABdoNd6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGutJREFUeJzt3X9M1Pf9B/DnGWg2t2phkQ8dx3aNgniIQKtCm5nZssPY\nFUbVkaKxTLzuD1Mzu6aZa7JYTCrXGrPRVv8xrFzafotmC8IapWrctX5rq+uQ1okB1x3JoXCJPa6l\n2g7B9/cPw/VLsFSONz/er8/zkZD087kfn/ezd31y97rPUYdSSoGIiESbNd0LICKiyceyJyKyAZY9\nEZENsOyJiGyAZU9EZAMseyIiGxiz7KuqqmBZFnJyckZdtmfPHsyaNQuRSCS2r6amBhkZGcjKysLR\no0f1r5aIiOIyZtlv2rQJLS0to/aHQiEcO3YMP/7xj2P72tvbceDAAbS3t6OlpQVbtmzBjRs39K+Y\niIjGbcyyX7FiBZKSkkbt/+1vf4sXX3xxxL6mpiZUVFQgMTERLpcLCxYswJkzZ/SuloiI4pIw3hs0\nNTXB6XRiyZIlI/ZfvnwZhYWFsW2n04lLly6NuI7D4YhzmURE9jbRP3Ywrg9or127hl27dqG6uvq2\nFnCrcldKif3ZsWPHtK+B+ZjPjvkkZ1NKz1+0Gdcr+08++QRdXV3Izc0FAHR3d+O+++7D6dOnkZaW\nhlAoFLtud3c30tLStCzSFF1dXdO9hEnFfGaTnE9yNl3G9co+JycH4XAYwWAQwWAQTqcTra2tsCwL\npaWlaGhowMDAAILBIC5evIjly5dP1rqJiGgcxiz7iooKPPDAA+js7ER6ejpeffXVEZf//zGN2+1G\neXk53G43Vq9ejX379tluRv+rX/1qupcwqZjPbJLzSc6mi0PpGgjdzsEcDm3zJyIiu9DRnfwGrUaB\nQGC6lzCpmM9skvNJzqYLy56IyAY4xiEimuE4xiEiotvCstdI+tyQ+cwmOZ/kbLqw7ImIbIAzeyKi\nGY4zeyIiui0se42kzw2Zz2yS80nOpgvLnojIBjizJyKa4TizJyKi28Ky10j63JD5zCY5n+RsurDs\niYhsgDN7IqIZjjN7IiK6LeP6f9DqcPr0aezf/9qo/Q4H8PTTTyIrK2uql6RNIBDAypUrp3sZk4b5\nzCY5n+Rsukx52b/99tuoq/s3gEdG7L/jjjdQUPC/Rpc9EdFMNeVlf9NyAE+OXEhC2/QsRSPpryyY\nz2yS80nOpgtn9kRENsCy10j6ub7MZzbJ+SRn04VlT0RkAyx7jaTPDZnPbJLzSc6my5hlX1VVBcuy\nkJOTE9v3zDPPYNGiRcjNzcWaNWvw2WefxS6rqalBRkYGsrKycPTo0clbNRERjcuYZb9p0ya0tLSM\n2FdcXIzz58/jo48+QmZmJmpqagAA7e3tOHDgANrb29HS0oItW7bgxo0bk7fyGUj63JD5zCY5n+Rs\nuox56uWKFSvQ1dU1Yp/H44n9c0FBAf76178CAJqamlBRUYHExES4XC4sWLAAZ86cQWFh4YjbNzY2\nApgD4DkAdwHIA7ASANDR0THiyxHDD6Ap221tbTNqPczHfHbKJ2k7EAigvr4eAOByuaDDt/5tnK6u\nLpSUlODcuXOjLispKUFFRQXWr1+PrVu3orCwEBs2bAAAeL1erF69GmvXrv36YA4HqqursWPHIICd\nI+5r9mwvamsL4fV6NcQiIpJjWv82zvPPP4877rgD69ev/8brOByOeO+eiIg0iqvs6+vrcfjwYbzx\nxhuxfWlpaQiFQrHt7u5upKWlTXyFBhl+GyYV85lNcj7J2XQZd9m3tLRg9+7daGpqwne+853Y/tLS\nUjQ0NGBgYADBYBAXL17E8uXLtS6WiIjiM+YHtBUVFXjnnXdw5coVpKeno7q6GjU1NRgYGIh9UHv/\n/fdj3759cLvdKC8vh9vtRkJCAvbt22e7Mc7wBy1SMZ/ZJOeTnE2XMcv+zTffHLWvqqrqG6//7LPP\n4tlnn534qoiISCt+g1Yj6XND5jOb5HySs+nCsicisgGWvUbS54bMZzbJ+SRn04VlT0RkAyx7jaTP\nDZnPbJLzSc6mC8ueiMgGWPYaSZ8bMp/ZJOeTnE0Xlj0RkQ2w7DWSPjdkPrNJzic5my4seyIiG2DZ\nayR9bsh8ZpOcT3I2XVj2REQ2wLLXSPrckPnMJjmf5Gy6sOyJiGyAZa+R9Lkh85lNcj7J2XRh2RMR\n2QDLXiPpc0PmM5vkfJKz6cKyJyKyAZa9RtLnhsxnNsn5JGfThWVPRGQDLHuNpM8Nmc9skvNJzqYL\ny56IyAbGLPuqqipYloWcnJzYvkgkAo/Hg8zMTBQXFyMajcYuq6mpQUZGBrKysnD06NHJW/UMJX1u\nyHxmk5xPcjZdxiz7TZs2oaWlZcQ+n88Hj8eDzs5OFBUVwefzAQDa29tx4MABtLe3o6WlBVu2bMGN\nGzcmb+VERHTbxiz7FStWICkpacS+5uZmVFZWAgAqKytx6NAhAEBTUxMqKiqQmJgIl8uFBQsW4MyZ\nM5O07JlJ+tyQ+cwmOZ/kbLokjPcG4XAYlmUBACzLQjgcBgBcvnwZhYWFses5nU5cunRp1O0bGxsB\nzAHwHIC7AOQBWAkA6OjoQCAQiL0lG34ATdlua2ubUethPuazUz5J24FAAPX19QAAl8sFHRxKKTXW\nFbq6ulBSUoJz584BAJKSktDX1xe7PDk5GZFIBFu3bkVhYSE2bNgAAPB6vXj44YexZs2arw/mcKC6\nuho7dgwC2DniOLNne1FbWwiv16slGBGRFA6HA99S1d9q3GfjWJaF3t5eAEBPTw9SUlIAAGlpaQiF\nQrHrdXd3Iy0tbUKLIyIiPcZd9qWlpfD7/QAAv9+PsrKy2P6GhgYMDAwgGAzi4sWLWL58ud7VznDD\nb8OkYj6zSc4nOZsuY87sKyoq8M477+DKlStIT0/Hzp07sX37dpSXl6Ourg4ulwsHDx4EALjdbpSX\nl8PtdiMhIQH79u2Dw+GYkhBERDS2b53Zaz0YZ/ZEROM2LTN7IiIyD8teI+lzQ+Yzm+R8krPpwrIn\nIrIBlr1Gw1+OkIr5zCY5n+RsurDsiYhsgGWvkfS5IfOZTXI+ydl0YdkTEdkAy14j6XND5jOb5HyS\ns+nCsicisgGWvUbS54bMZzbJ+SRn04VlT0RkAyx7jaTPDZnPbJLzSc6mC8ueiMgGWPYaSZ8bMp/Z\nJOeTnE0Xlj0RkQ2w7DWSPjdkPrNJzic5my4seyIiG2DZayR9bsh8ZpOcT3I2XVj2REQ2wLLXSPrc\nkPnMJjmf5Gy6sOyJiGyAZa+R9Lkh85lNcj7J2XRh2RMR2UDcZV9TU4Ps7Gzk5ORg/fr1+O9//4tI\nJAKPx4PMzEwUFxcjGo3qXOuMJ31uyHxmk5xPcjZd4ir7rq4u7N+/H62trTh37hyGhobQ0NAAn88H\nj8eDzs5OFBUVwefz6V4vERHFIa6ynzNnDhITE3Ht2jUMDg7i2rVr+OEPf4jm5mZUVlYCACorK3Ho\n0CGti53ppM8Nmc9skvNJzqZLQjw3Sk5OxtNPP40f/ehH+O53v4tVq1bB4/EgHA7DsiwAgGVZCIfD\no27b2NgIYA6A5wDcBSAPwEoAQEdHBwKBQOwt2fADaMp2W1vbjFoP8zGfnfJJ2g4EAqivrwcAuFwu\n6OBQSqnx3uiTTz5BSUkJTp48iblz5+KXv/wl1q5di61bt6Kvry92veTkZEQika8P5nCguroaO3YM\nAtg54j5nz/aitrYQXq83/jRERAI5HA7EUdUjxDXG+fDDD/HAAw/gBz/4ARISErBmzRq8//77SE1N\nRW9vLwCgp6cHKSkpE1ocERHpEVfZZ2Vl4YMPPsCXX34JpRSOHz8Ot9uNkpIS+P1+AIDf70dZWZnW\nxc50w2/DpGI+s0nOJzmbLnHN7HNzc/H4449j6dKlmDVrFu699178+te/Rn9/P8rLy1FXVweXy4WD\nBw/qXi8REcUhrpl93AfjzJ6IaNymbWZPRERmYdlrJH1uyHxmk5xPcjZdWPZERDbAstdo+MsRUjGf\n2STnk5xNF5Y9EZENsOw1kj43ZD6zSc4nOZsuLHsiIhtg2WskfW7IfGaTnE9yNl1Y9kRENsCy10j6\n3JD5zCY5n+RsurDsiYhsgGWvkfS5IfOZTXI+ydl0YdkTEdkAy14j6XND5jOb5HySs+nCsicisgGW\nvUbS54bMZzbJ+SRn04VlT0RkAyx7jaTPDZnPbJLzSc6mC8ueiMgGWPYaSZ8bMp/ZJOeTnE0Xlj0R\nkQ2w7DWSPjdkPrNJzic5my4seyIiG4i77KPRKNatW4dFixbB7Xbj9OnTiEQi8Hg8yMzMRHFxMaLR\nqM61znjS54bMZzbJ+SRn0yXusv/Nb36Dhx9+GBcuXMDHH3+MrKws+Hw+eDwedHZ2oqioCD6fT+da\niYgoTnGV/WeffYaTJ0+iqqoKAJCQkIC5c+eiubkZlZWVAIDKykocOnRI30oNIH1uyHxmk5xPcjZd\nEuK5UTAYxLx587Bp0yZ89NFHuO+++/CnP/0J4XAYlmUBACzLQjgcHnXbxsZGAHMAPAfgLgB5AFYC\nADo6OhAIBGJvyYYfQFO229raZtR6mI/57JRP0nYgEEB9fT0AwOVyQQeHUkqN90Yffvgh7r//fpw6\ndQrLli3Dtm3bcOedd+KVV15BX19f7HrJycmIRCJfH8zhQHV1NXbsGASwc8R9zp7tRW1tIbxeb/xp\niIgEcjgciKOqR4hrjON0OuF0OrFs2TIAwLp169Da2orU1FT09vYCAHp6epCSkjKhxRERkR5xlX1q\nairS09PR2dkJADh+/Diys7NRUlICv98PAPD7/SgrK9O3UgMMvw2TivnMJjmf5Gy6xDWzB4CXX34Z\nGzZswMDAAObPn49XX30VQ0NDKC8vR11dHVwuFw4ePKhzrUREFKe4yz43Nxf/+Mc/Ru0/fvz4hBZk\nsuEPWqRiPrNJzic5my78Bi0RkQ2w7DWSPjdkPrNJzic5my4seyIiG2DZayR9bsh8ZpOcT3I2XVj2\nREQ2wLLXSPrckPnMJjmf5Gy6sOyJiGyAZa+R9Lkh85lNcj7J2XRh2RMR2QDLXiPpc0PmM5vkfJKz\n6cKyJyKyAZa9RtLnhsxnNsn5JGfThWVPRGQDLHuNpM8Nmc9skvNJzqYLy56IyAZY9hpJnxsyn9kk\n55OcTReWPRGRDbDsNZI+N2Q+s0nOJzmbLix7IiIbYNlrJH1uyHxmk5xPcjZdWPZERDbAstdI+tyQ\n+cwmOZ/kbLqw7ImIbCDush8aGkJ+fj5KSkoAAJFIBB6PB5mZmSguLkY0GtW2SFNInxsyn9kk55Oc\nTZe4y762thZutxsOhwMA4PP54PF40NnZiaKiIvh8Pm2LJCKiiYmr7Lu7u3H48GF4vV4opQAAzc3N\nqKysBABUVlbi0KFD+lZpCOlzQ+Yzm+R8krPpkhDPjZ566ins3r0bn3/+eWxfOByGZVkAAMuyEA6H\nb3nbxsZGAHMAPAfgLgB5AFYCADo6OhAIBGJvyYYfQFO229raZtR6mI/57JRP0nYgEEB9fT0AwOVy\nQQeHGn5pfpveeustHDlyBHv37kUgEMCePXvwt7/9DUlJSejr64tdLzk5GZFIZOTBHA5UV1djx45B\nADtHXDZ7the1tYXwer3xpyEiEsjhcGCcVT3KuF/Znzp1Cs3NzTh8+DC++uorfP7559i4cSMsy0Jv\nby9SU1PR09ODlJSUCS2MiIj0GffMfteuXQiFQggGg2hoaMBDDz2E1157DaWlpfD7/QAAv9+PsrIy\n7Yud6YbfhknFfGaTnE9yNl0mfJ798Nk427dvx7Fjx5CZmYkTJ05g+/btE14cERHpMe6Z/YQOxpk9\nEdG46ZjZ8xu0REQ2wLLXSPrckPnMJjmf5Gy6sOyJiGyAZa/R8JcjpGI+s0nOJzmbLix7IiIbYNlr\nJH1uyHxmk5xPcjZdWPZERDbAstdI+tyQ+cwmOZ/kbLqw7ImIbIBlr5H0uSHzmU1yPsnZdGHZExHZ\nAMteI+lzQ+Yzm+R8krPpwrInIrIBlr1G0ueGzGc2yfkkZ9OFZU9EZAMse42kzw2Zz2yS80nOpgvL\nnojIBlj2GkmfGzKf2STnk5xNF5Y9EZENsOw1kj43ZD6zSc4nOZsuLHsiIhtg2WskfW7IfGaTnE9y\nNl1Y9kRENhBX2YdCITz44IPIzs7G4sWL8dJLLwEAIpEIPB4PMjMzUVxcjGg0qnWxM530uSHzmU1y\nPsnZdImr7BMTE/HHP/4R58+fxwcffIC9e/fiwoUL8Pl88Hg86OzsRFFREXw+n+71EhFRHOIq+9TU\nVOTl5QEAvv/972PRokW4dOkSmpubUVlZCQCorKzEoUOH9K3UANLnhsxnNsn5JGfTJWGid9DV1YWz\nZ8+ioKAA4XAYlmUBACzLQjgcHnX9xsZGAHMAPAfgLgB5AFYCADo6OhAIBGJvyYYfQFO229raZtR6\nmI/57JRP0nYgEEB9fT0AwOVyQQeHUkrFe+MvvvgCP/3pT/GHP/wBZWVlSEpKQl9fX+zy5ORkRCKR\nrw/mcKC6uho7dgwC2DnivmbP9qK2thBerzfe5RARieRwODCBqgYwgbNxrl+/jrVr12Ljxo0oKysD\ncPPVfG9vLwCgp6cHKSkpE1ocERHpEVfZK6WwefNmuN1ubNu2Lba/tLQUfr8fAOD3+2O/BOxi+G2Y\nVMxnNsn5JGfTJa6Z/XvvvYfXX38dS5YsQX5+PgCgpqYG27dvR3l5Oerq6uByuXDw4EGtiyUiovjE\nVfY/+clPcOPGjVtedvz48QktyGTDH7RIxXxmk5xPcjZd+A1aIiIbYNlrJH1uyHxmk5xPcjZdWPZE\nRDbAstdI+tyQ+cwmOZ/kbLqw7ImIbIBlr5H0uSHzmU1yPsnZdGHZExHZAMteI+lzQ+Yzm+R8krPp\nwrInIrIBlr1G0ueGzGc2yfkkZ9OFZU9EZAMse42kzw2Zz2yS80nOpgvLnojIBlj2GkmfGzKf2STn\nk5xNF5Y9EZENsOw1kj43ZD6zSc4nOZsuLHsiIhtg2WskfW7IfGaTnE9yNl1Y9kRENsCy10j63JD5\nzCY5n+RsurDsiYhsYMaU/bVr/4MnnngCDodjxM+cOcnTvbTbJn1uyHxmk5xPcjZdZkzZA18CUKN+\n+vv7pnVV49HW1jbdS5hUzGc2yfkkZ9NFe9m3tLQgKysLGRkZeOGFFzTcY8KoV/s3f+645f7pfCcQ\njUan7dhTgfnMJjmf5Gy6aC37oaEhPPnkk2hpaUF7ezvefPNNXLhwYYL3OohbveIHrt9yf39//23/\nYpiJvzDsbM6cZC2Ph677Gc9987lkhsl8bsx0CTrv7MyZM1iwYAFcLhcA4LHHHkNTUxMWLVqk8zDf\nYviXw//nuMW+b97f3++I68hdXV1x3c4Uk53v5shu4o9HvPdzO/m+6b51P5cmg+Tn5+1m0/UcM5FD\nKXWrZ25c/vKXv+Dtt9/G/v37AQCvv/46Tp8+jZdffvnmwRzy/4USEU2GiVa11lf231bmGn+vEBHR\nOGid2aelpSEUCsW2Q6EQnE6nzkMQEVEctJb90qVLcfHiRXR1dWFgYAAHDhxAaWmpzkMQEVEctI5x\nEhIS8Morr2DVqlUYGhrC5s2bp/jDWSIiuhXt59mvXr0aHR0d+Pe//43f//73sf36z7+felVVVbAs\nCzk5ObF9kUgEHo8HmZmZKC4uHnG+b01NDTIyMpCVlYWjR49Ox5LHJRQK4cEHH0R2djYWL16Ml156\nCYCcjF999RUKCgqQl5cHt9sde35KyQfcPP05Pz8fJSUlAGRlc7lcWLJkCfLz87F8+XIAsvJFo1Gs\nW7cOixYtgtvtxunTp/XmU1NgcHBQzZ8/XwWDQTUwMKByc3NVe3v7VBxaq3fffVe1traqxYsXx/Y9\n88wz6oUXXlBKKeXz+dTvfvc7pZRS58+fV7m5uWpgYEAFg0E1f/58NTQ0NC3rvl09PT3q7NmzSiml\n+vv7VWZmpmpvbxeV8erVq0oppa5fv64KCgrUyZMnReXbs2ePWr9+vSopKVFKyXp+ulwu9emnn47Y\nJynf448/rurq6pRSN5+f0WhUa74pKftTp06pVatWxbZrampUTU3NVBxau2AwOKLsFy5cqHp7e5VS\nN8ty4cKFSimldu3apXw+X+x6q1atUu+///7ULnaCfvGLX6hjx46JzHj16lW1dOlS9a9//UtMvlAo\npIqKitSJEyfUI488opSS9fx0uVzqypUrI/ZJyReNRtU999wzar/OfFPyt3EuXbqE9PT02LbT6cSl\nS5em4tCTLhwOw7IsAIBlWQiHwwCAy5cvjzgTybTMXV1dOHv2LAoKCkRlvHHjBvLy8mBZVmxkJSXf\nU089hd27d2PWrK//s5aSDbh5avfPfvYzLF26NPZdHin5gsEg5s2bh02bNuHee+/FE088gatXr2rN\nNyVlb5cvUw1/9Xqsy03wxRdfYO3ataitrcWdd9454jLTM86aNQttbW3o7u7Gu+++i7///e8jLjc1\n31tvvYWUlBTk5+d/4/dZTM027L333sPZs2dx5MgR7N27FydPnhxxucn5BgcH0draii1btqC1tRXf\n+9734PP5RlxnovmmpOwln39vWRZ6e3sBAD09PUhJSQEwOnN3dzfS0tKmZY3jcf36daxduxYbN25E\nWVkZAHkZAWDu3Ln4+c9/jn/+858i8p06dQrNzc245557UFFRgRMnTmDjxo0isg27++67AQDz5s3D\no48+ijNnzojJ53Q64XQ6sWzZMgDAunXr0NraitTUVG35pqTsJZ9/X1paCr/fDwDw+/2xgiwtLUVD\nQwMGBgYQDAZx8eLF2BkEM5VSCps3b4bb7ca2bdti+6VkvHLlSuxshi+//BLHjh1Dfn6+iHy7du1C\nKBRCMBhEQ0MDHnroIbz22msisgHAtWvX0N/fDwC4evUqjh49ipycHDH5UlNTkZ6ejs7OTgDA8ePH\nkZ2djZKSEn35tH3C8C0OHz6sMjMz1fz589WuXbum6rBaPfbYY+ruu+9WiYmJyul0qj//+c/q008/\nVUVFRSojI0N5PB7V19cXu/7zzz+v5s+frxYuXKhaWlqmceW35+TJk8rhcKjc3FyVl5en8vLy1JEj\nR8Rk/Pjjj1V+fr7Kzc1VOTk56sUXX1RKKTH5hgUCgdjZOFKy/ec//1G5ubkqNzdXZWdnxzpESj6l\nlGpra1NLly5VS5YsUY8++qiKRqNa82n9Q2hERDQzzaD/UxUREU0Wlj0RkQ2w7ImIbIBlT0RkAyx7\nIiIbYNkTEdnA/wEfVfa46/PXTQAAAABJRU5ErkJggg==\n" } ], "prompt_number": 83 }, { "cell_type": "raw", "metadata": {}, "source": [ "SELECT Column1,\n", " Column13\n", "FROM [sr320@washington.edu].[BiGo_RNAseq_unmapped_nt_blastout_taxa2]\u200b" ] }, { "cell_type": "code", "collapsed": false, "input": [ "@fu" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "SyntaxError", "evalue": "unexpected EOF while parsing (, line 1)", "output_type": "pyerr", "traceback": [ "\u001b[0;36m File \u001b[0;32m\"\"\u001b[0;36m, line \u001b[0;32m1\u001b[0m\n\u001b[0;31m @fu\u001b[0m\n\u001b[0m ^\u001b[0m\n\u001b[0;31mSyntaxError\u001b[0m\u001b[0;31m:\u001b[0m unexpected EOF while parsing\n" ] } ], "prompt_number": 10 }, { "cell_type": "markdown", "metadata": {}, "source": [ "--- \n", "\n", "\n", "##Tophat Analysis \n", "\n", "\n", "--\n", "\n", "\"(274)%20Discovery%20Environment\"\n", "\n", "\n", "TopHat Output\n", "\n", "The tophat script produces a number of files in the directory in which it was invoked. Most of these files are internal, intermediate files that are generated for use within the pipeline. The output files you will likely want to look at are:\n", "\n", "accepted_hits.bam. A list of read alignments in SAM format. SAM is a compact short read alignment format that is increasingly being adopted. The formal specification is here. \n", "\n", "junctions.bed. A UCSC BED track of junctions reported by TopHat. Each junction consists of two connected BED blocks, where each block is as long as the maximal overhang of any read spanning the junction. The score is the number of alignments spanning the junction. \n", "\n", "insertions.bed and deletions.bed. UCSC BED tracks of insertions and deletions reported by TopHat. \n", "\n", "Insertions.bed - chromLeft refers to the last genomic base before the insertion. \n", "\n", "Deletions.bed - chromLeft refers to the first genomic base of the deletion. \n", "\n", "---\n", "\n", "Data @ \n", "\n", "\n" ] }, { "cell_type": "code", "collapsed": false, "input": [ "ls /Volumes/web/cnidarian/tophat_071313" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\u001b[31maccepted_hits.bam\u001b[m\u001b[m* \u001b[31mprep_reads.info\u001b[m\u001b[m*\r\n", "\u001b[31mdeletions.bed\u001b[m\u001b[m* \u001b[31ms_1.bam\u001b[m\u001b[m*\r\n", "\u001b[31minsertions.bed\u001b[m\u001b[m* \u001b[31ms_1.bam.bai\u001b[m\u001b[m*\r\n", "\u001b[31mjunctions.bed\u001b[m\u001b[m* \u001b[31munmapped.bam\u001b[m\u001b[m*\r\n", "\u001b[31mlookup_roberts_grc_oyster.xls\u001b[m\u001b[m*\r\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/tophat_071313/junctions.bed" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "track name=junctions description=\"TopHat junctions\"\r\n", "C14162\t6\t239\tJUNC00000001\t12\t-\t6\t239\t255,0,0\t2\t36,49\t0,184\r\n", "C14832\t78\t229\tJUNC00000002\t4\t-\t78\t229\t255,0,0\t2\t43,27\t0,124\r\n", "C17940\t206\t349\tJUNC00000003\t1\t-\t206\t349\t255,0,0\t2\t31,19\t0,124\r\n", "C17998\t32\t262\tJUNC00000004\t734\t-\t32\t262\t255,0,0\t2\t49,49\t0,181\r\n", "C18186\t378\t530\tJUNC00000005\t1\t+\t378\t530\t255,0,0\t2\t29,21\t0,131\r\n", "C18450\t359\t554\tJUNC00000006\t6\t+\t359\t554\t255,0,0\t2\t38,49\t0,146\r\n", "C18694\t57\t245\tJUNC00000007\t1\t-\t57\t245\t255,0,0\t2\t30,20\t0,168\r\n", "C18754\t51\t472\tJUNC00000008\t6\t+\t51\t472\t255,0,0\t2\t47,32\t0,389\r\n", "C19028\t89\t279\tJUNC00000009\t76\t+\t89\t279\t255,0,0\t2\t40,49\t0,141\r\n" ] } ], "prompt_number": 2 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"IGV%20-%20Session:%20http://eagle.fish.washington.edu/cnidarian/oyster_v9_igv_session.xml%20and%20Home%20%7C%20TWiT.TV\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "URL to load IGV Session corresponding to screenshot above. \n", "\n", "`http://eagle.fish.washington.edu/cnidarian/BiGo_RNAseq_072213_igv_session.xml`" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "## Using Samtools to interogate TopHat Bam file\n", "\n", "reference: " ] }, { "cell_type": "code", "collapsed": false, "input": [ "cd /Volumes/Bay3/Software/BSMAP/bsmap-2.74/samtools" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "/Volumes/Bay3/Software/BSMAP/bsmap-2.74/samtools" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "\n" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "ls" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "AUTHORS bam_index.o bedidx.o kstring.c\r\n", "COPYING bam_lpileup.c bgzf.c kstring.h\r\n", "ChangeLog bam_lpileup.o bgzf.h kstring.o\r\n", "INSTALL bam_mate.c bgzf.o libbam.a\r\n", "Makefile bam_mate.o bgzip.c \u001b[34mmisc\u001b[m\u001b[m/\r\n", "Makefile.mingw bam_md.c cut_target.c phase.c\r\n", "NEWS bam_md.o cut_target.o phase.o\r\n", "bam.c bam_pileup.c errmod.c razf.c\r\n", "bam.h bam_pileup.o errmod.h razf.h\r\n", "bam.o bam_plcmd.c errmod.o razf.o\r\n", "bam2bcf.c bam_plcmd.o \u001b[34mexamples\u001b[m\u001b[m/ razip.c\r\n", "bam2bcf.h bam_reheader.c faidx.c sam.c\r\n", "bam2bcf.o bam_reheader.o faidx.h sam.h\r\n", "bam2bcf_indel.c bam_rmdup.c faidx.o sam.o\r\n", "bam2bcf_indel.o bam_rmdup.o kaln.c sam_header.c\r\n", "bam2depth.c bam_rmdupse.c kaln.h sam_header.h\r\n", "bam2depth.o bam_rmdupse.o kaln.o sam_header.o\r\n", "bam_aux.c bam_sort.c khash.h sam_view.c\r\n", "bam_aux.o bam_sort.o klist.h sam_view.o\r\n", "bam_cat.c bam_stat.c knetfile.c sample.c\r\n", "bam_cat.o bam_stat.o knetfile.h sample.h\r\n", "bam_color.c bam_tview.c knetfile.o sample.o\r\n", "bam_color.o bam_tview.o kprobaln.c \u001b[31msamtools\u001b[m\u001b[m*\r\n", "bam_endian.h bamtk.c kprobaln.h samtools.1\r\n", "bam_import.c bamtk.o kprobaln.o \u001b[34mwin32\u001b[m\u001b[m/\r\n", "bam_import.o \u001b[34mbcftools\u001b[m\u001b[m/ kseq.h\r\n", "bam_index.c bedidx.c ksort.h\r\n" ] } ], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "/bin/sh: samtools: command not found\r\n" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "cp samtools /usr/local/bin" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\r\n", "Program: samtools (Tools for alignments in the SAM format)\r\n", "Version: 0.1.18 (r982:295)\r\n", "\r\n", "Usage: samtools [options]\r\n", "\r\n", "Command: view SAM<->BAM conversion\r\n", " sort sort alignment file\r\n", " mpileup multi-way pileup\r\n", " depth compute the depth\r\n", " faidx index/extract FASTA\r\n", " index index alignment\r\n", " idxstats BAM index stats (r595 or later)\r\n", " fixmate fix mate information\r\n", " flagstat simple stats\r\n", " calmd recalculate MD/NM tags and '=' bases\r\n", " merge merge sorted alignments\r\n", " rmdup remove PCR duplicates\r\n", " reheader replace BAM header\r\n", " cat concatenate BAMs\r\n", " targetcut cut fosmid regions (for fosmid pool only)\r\n", " phase phase heterozygotes\r\n", "\r\n" ] } ], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools view -c /Volumes/web/cnidarian/tophat_071313/s_1.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "164718919\r\n" ] } ], "prompt_number": 11 }, { "cell_type": "code", "collapsed": false, "input": [ "#only mapped reads\n", "!samtools view -c -F 4 /Volumes/web/cnidarian/tophat_071313/s_1.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "164718919\r\n" ] } ], "prompt_number": 12 }, { "cell_type": "code", "collapsed": false, "input": [ "#unmapped reads\n", "!samtools view -c -f 4 /Volumes/web/cnidarian/tophat_071313/s_1.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "0\r\n" ] } ], "prompt_number": 13 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools view -c /Volumes/web/cnidarian/tophat_071313/accepted_hits.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "164718919\r\n" ] } ], "prompt_number": 19 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools view -c /Volumes/web/cnidarian/tophat_071313/unmapped.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "7263587\r\n" ] } ], "prompt_number": 22 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools flagstat /Volumes/web/cnidarian/tophat_071313/s_1.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "164718919 + 0 in total (QC-passed reads + QC-failed reads)\r\n", "0 + 0 duplicates\r\n", "164718919 + 0 mapped (100.00%:nan%)\r\n", "164718919 + 0 paired in sequencing\r\n", "82181831 + 0 read1\r\n", "82537088 + 0 read2\r\n", "42418114 + 0 properly paired (25.75%:nan%)\r\n", "161646608 + 0 with itself and mate mapped\r\n", "3072311 + 0 singletons (1.87%:nan%)\r\n", "1264082 + 0 with mate mapped to a different chr\r\n", "395312 + 0 with mate mapped to a different chr (mapQ>=5)\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "164718919 + 0 in total (QC-passed reads + QC-failed reads)\r\n", "0 + 0 duplicates\r\n", "164718919 + 0 mapped (100.00%:nan%)\r\n", "164718919 + 0 paired in sequencing\r\n", "82181831 + 0 read1\r\n", "82537088 + 0 read2\r\n", "42418114 + 0 properly paired (25.75%:nan%)\r\n", "161646608 + 0 with itself and mate mapped\r\n", "3072311 + 0 singletons (1.87%:nan%)\r\n", "1264082 + 0 with mate mapped to a different chr\r\n", "395312 + 0 with mate mapped to a different chr (mapQ>=5)\r\n" ] } ], "prompt_number": 18 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools flagstat /Volumes/web/cnidarian/tophat_071313/accepted_hits.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "164718919 + 0 in total (QC-passed reads + QC-failed reads)\r\n", "0 + 0 duplicates\r\n", "164718919 + 0 mapped (100.00%:nan%)\r\n", "164718919 + 0 paired in sequencing\r\n", "82181831 + 0 read1\r\n", "82537088 + 0 read2\r\n", "42418114 + 0 properly paired (25.75%:nan%)\r\n", "161646608 + 0 with itself and mate mapped\r\n", "3072311 + 0 singletons (1.87%:nan%)\r\n", "1264082 + 0 with mate mapped to a different chr\r\n", "395312 + 0 with mate mapped to a different chr (mapQ>=5)\r\n" ] } ], "prompt_number": 15 }, { "cell_type": "code", "collapsed": false, "input": [ "!samtools flagstat /Volumes/web/cnidarian/tophat_071313/unmapped.bam" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "7161175 + 102412 in total (QC-passed reads + QC-failed reads)\r\n", "0 + 0 duplicates\r\n", "0 + 0 mapped (0.00%:0.00%)\r\n", "7161175 + 102412 paired in sequencing\r\n", "3670234 + 76916 read1\r\n", "3490941 + 25496 read2\r\n", "0 + 0 properly paired (0.00%:0.00%)\r\n", "0 + 0 with itself and mate mapped\r\n", "0 + 0 singletons (0.00%:0.00%)\r\n", "0 + 0 with mate mapped to a different chr\r\n", "0 + 0 with mate mapped to a different chr (mapQ>=5)\r\n" ] } ], "prompt_number": 16 }, { "cell_type": "markdown", "metadata": {}, "source": [ "##Tophat Data Into Cufflinks (IPlant)" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#cufflinks failed" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "##Quality Trimming in iPlant\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"_264__Discovery_Environment_17A1B77F.png\"/ \n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Pre-Process Complete, generating one file \n", "\n", "Downloading and running through Fastx quality Trimmer" ] }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/BiGo_RNAseq_fastx_qual.fastq" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 190998292 190998292 6469564027 /Volumes/web/cnidarian/BiGo_RNAseq_fastx_qual.fastq\r\n" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "!wc /Volumes/web/cnidarian/BiGo_RNAseq_PreProcess.fastq" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ " 191927856 191927856 6501051882 /Volumes/web/cnidarian/BiGo_RNAseq_PreProcess.fastq\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ " 191927856 191927856 6501051882 /Volumes/web/cnidarian/BiGo_RNAseq_PreProcess.fastq\r\n" ] } ], "prompt_number": 6 }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "FastQC in iPlant" ] }, { "cell_type": "code", "collapsed": false, "input": [ "cd /Volumes/web/cnidarian/" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "/Volumes/web/cnidarian\n" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "!iget -r /iplant/home/sr320/analyses/FastQC_0.10.1__BiGoRNA-2014-03-24-10-33-44.916" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"BiGoRNA_GTGTCTAC_1_fastq_FastQC_Report_18E0A92A.png\"/" ] }, { "cell_type": "code", "collapsed": false, "input": [ "Run on filtered file.." ], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"fastq_quality_filter_out_fastq_FastQC_Report_18E0DD4F.png\"/" ] }, { "cell_type": "raw", "metadata": {}, "source": [ "Still not great...." ] }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "In CLC" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "\"CLC_Genomics_Workbench_7_0_18E1E1D9.png\"/" ] }, { "cell_type": "code", "collapsed": false, "input": [ "new file" ], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "code", "collapsed": false, "input": [ "!head /Volumes/web/cnidarian/BiGoRNA_genetable_clc" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "\"Name\"\t\"Chromosome\"\t\"Region\"\t\"Expression value\"\t\"Gene length\"\t\"RPKM\"\t\"Unique gene reads\"\t\"Total gene reads\"\t\"Transcripts annotated\"\t\"Detected transcripts\"\t\"Exon length\"\t\"Exons\"\t\"Unique exon reads\"\t\"Total exon reads\"\t\"Ratio of unique to total (exon reads)\"\t\"Unique exon-exon reads\"\t\"Total exon-exon reads\"\t\"Unique intron reads\"\t\"Total intron reads\"\t\"Ratio of intron to total gene reads\"\r\n", "CGI_10000780\tscaffold350\tcomplement(1105..3206)\t0\t2102\t0\t0\t0\t2\t0\t1350\t2\t0\t0\t?\t0\t0\t0\t0\t?\r\n", "CGI_10000456\tscaffold28720\t349..2038\t0\t1690\t0.323\t2\t2\t4\t0\t438\t4\t0\t0\t?\t0\t0\t2\t2\t1\r\n", "CGI_10000457\tscaffold28720\tcomplement(3692..5692)\t1\t2001\t0.469\t0\t4\t7\t0\t603\t7\t0\t1\t0\t0\t0\t0\t3\t0.75\r\n", "CGI_10000774\tscaffold31684\tcomplement(6549..7563)\t0\t1015\t0.189\t1\t1\t2\t0\t375\t2\t0\t0\t?\t0\t0\t1\t1\t1\r\n", "CGI_10000917\tscaffold32586\t1736..2161\t0\t426\t0\t0\t0\t1\t0\t426\t1\t0\t0\t?\t0\t0\t0\t0\t?\r\n", "CGI_10000861\tscaffold900\tcomplement(6516..8519)\t1\t2004\t0.035\t1\t1\t1\t1\t2004\t1\t1\t1\t1\t0\t0\t0\t0\t0\r\n", "CGI_10000994\tscaffold33132\tcomplement(5089..12821)\t16\t7733\t1.385\t27\t32\t11\t2\t1635\t11\t14\t16\t0.875\t0\t0\t13\t16\t" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "0.5\r\n", "CGI_10000643\tscaffold740\t109..2125\t0\t2017\t0\t0\t0\t4\t0\t552\t4\t0\t0\t?\t0\t0\t0\t0\t?\r\n", "CGI_10000763\tscaffold31534\tcomplement(1194..3724)\t0\t2531\t0\t0\t0\t3\t0\t249\t3\t0\t0\t?\t0\t0\t0\t0\t?\r\n" ] } ], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "#Location of SQLShare python tools: you can empty (\"\") if tools are in PATH\n", "spd=\"/Users/sr320/sqlshare-pythonclient/tools/\"" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 10 }, { "cell_type": "code", "collapsed": false, "input": [ "!python {spd}singleupload.py -d BiGoRNA_genetable_clc /Volumes/web/cnidarian/BiGoRNA_genetable_clc\n" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "processing chunk line 0 to 28028 (0.0644998550415 s elapsed)\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "pushing /Volumes/web/cnidarian/BiGoRNA_genetable_clc...\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "parsing FA491187...\r\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "finished BiGoRNA_genetable_clc\r\n" ] } ], "prompt_number": 11 }, { "cell_type": "code", "collapsed": false, "input": [ "%pylab inline" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "from pandas import *\n", "\n", "# read data from data file into a pandas DataFrame \n", "BiGoRNA = read_table(\"/Volumes/web/cnidarian/BiGoRNA_genetable_clc\", # name of the data file\n", " #sep=\"\\t\", # what character separates each column?\n", " #na_values=[\"\", \" \"], # what values should be considered \"blank\" values?\n", " #header=None\n", " )" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "ImportError", "evalue": "No module named pandas", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mImportError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;32mfrom\u001b[0m \u001b[0mpandas\u001b[0m \u001b[0;32mimport\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 3\u001b[0m \u001b[0;31m# read data from data file into a pandas DataFrame\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 4\u001b[0m BiGoRNA = read_table(\"/Volumes/web/cnidarian/BiGoRNA_genetable_clc\", # name of the data file\n\u001b[1;32m 5\u001b[0m \u001b[0;31m#sep=\"\\t\", # what character separates each column?\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mImportError\u001b[0m: No module named pandas" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "BiGoRNA\n" ], "language": "python", "metadata": {}, "outputs": [] }, { "cell_type": "code", "collapsed": false, "input": [ "plot.BiGoRNA[\"Expression value\"]\n", "#Axis limits are changed using the axis([xmin, xmax, ymin, ymax]) function.\n", "#plt.axis([0, 200, 0, 30000]);\n", "#plt.title('Expression value');" ], "language": "python", "metadata": {}, "outputs": [ { "ename": "AttributeError", "evalue": "'function' object has no attribute 'BiGoRNA'", "output_type": "pyerr", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mAttributeError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mplot\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mBiGoRNA\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"Expression value\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2\u001b[0m \u001b[0;31m#Axis limits are changed using the axis([xmin, xmax, ymin, ymax]) function.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 3\u001b[0m \u001b[0;31m#plt.axis([0, 200, 0, 30000]);\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 4\u001b[0m \u001b[0;31m#plt.title('Expression value');\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", "\u001b[0;31mAttributeError\u001b[0m: 'function' object has no attribute 'BiGoRNA'" ] } ], "prompt_number": 58 }, { "cell_type": "code", "collapsed": false, "input": [ "\n", "plt.hist(BiGoRNA[\"Total intron reads\"], 20)\n", "plt.axis([0, 100, 0, 2000])\n", "plt.show()\n" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAD9CAYAAABazssqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEhRJREFUeJzt3FFsU+fdx/Hf2WKpUlfSbWqczGYLSpwYk5RkoimaVKld\n5rBWwQujQmQTZMCqiWgVnXq1qxcuIK6mqqNbkaaJShmTluxmhFbDyqq+Wau1GEFAmzBavMmRHCex\nNkIQdC2h8LwXpe6bEUxJ0jjt//uRHqk+9jl+zqPir7F98JxzTgAAkz5X6gkAAEqHCACAYUQAAAwj\nAgBgGBEAAMOIAAAYVjQC2WxWjz32mNasWaOGhga9+OKLkqSpqSlFo1HV1dWpra1N09PThX16enoU\nCoUUDoc1ODhY2H769Gk1NjYqFAppz549n9DpAADuRtEI+Hw+vfDCCzp37pxOnDihl156SefPn1c8\nHlc0GtXIyIhaW1sVj8clSalUSv39/UqlUkokEuru7taHlyHs3r1bhw8fVjqdVjqdViKR+OTPDgBQ\nVNEIVFZWqqmpSZL0hS98QatXr1Yul9OxY8fU1dUlSerq6tLRo0clSQMDA+rs7JTP51N1dbVqa2uV\nTCY1MTGhy5cvq6WlRZK0ffv2wj4AgNIp+7gPHB0d1ZkzZ/Twww8rn8/L7/dLkvx+v/L5vCRpfHxc\n69evL+wTDAaVy+Xk8/kUDAYL2wOBgHK53Kzje563oBMBAKsW8g8/fKwvhq9cuaLNmzfr4MGDuu++\n+2bd53neor2Al5e3S3IM/c8ymMNyGawFa8FaFB8Lc8cIXLt2TZs3b9a2bdvU0dEh6YN3/5OTk5Kk\niYkJVVRUSPrgHX42my3sOzY2pmAwqEAgoLGxsVnbA4HAgicPAFiYohFwzmnXrl2KRCJ65plnCttj\nsZh6e3slSb29vYU4xGIx9fX1aWZmRplMRul0Wi0tLaqsrNSKFSuUTCblnNORI0cK+wAASsgV8eab\nbzrP89zatWtdU1OTa2pqcsePH3cXLlxwra2tLhQKuWg06i5evFjYZ//+/a6mpsbV19e7RCJR2H7q\n1CnX0NDgampq3NNPP33Lc0ly5eXtTnIM/e8ymMNyGawFa8FaFB8q9jJ+R95y+aekPc9TeXm7Ll16\npdRTAYBPEU+f+BfDAIDPJiIAAIYRAQAwjAgAgGFEAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEB\nADCMCACAYUQAAAwjAgBgGBEAAMOIAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEQAAw4gA\nABhGBADAMCIAAIYRAQAwjAgAgGFEAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEBADCMCACAYUQA\nAAwjAgBgGBEAAMOIAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEQAAw4pGYOfOnfL7/Wps\nbCxs27t3r4LBoJqbm9Xc3Kzjx48X7uvp6VEoFFI4HNbg4GBh++nTp9XY2KhQKKQ9e/Z8AqcBAJgX\nV8Qbb7zhhoeHXUNDQ2Hb3r173fPPP3/LY8+dO+fWrl3rZmZmXCaTcTU1Ne7GjRvOOeceeughl0wm\nnXPOPf744+748eO37C/JlZe3O8kxGAwG42MPFXsZv6OyYoF45JFHNDo6Olc4btk2MDCgzs5O+Xw+\nVVdXq7a2VslkUl/72td0+fJltbS0SJK2b9+uo0eP6tvf/vYtx3jvvRFJe2/eevTmAAB8ZOjmWBxF\nI3A7v/jFL/Sb3/xG69at0/PPP6/7779f4+PjWr9+feExwWBQuVxOPp9PwWCwsD0QCCiXy8153Hvu\nqdPVq3vnMyUAMOJRzX6DvG9BR7vrL4Z3796tTCajs2fPqqqqSs8+++yCJgAAKJ27jkBFRYU8z5Pn\nefrhD3+okydPSvrgHX42my08bmxsTMFgUIFAQGNjY7O2BwKBRZg6AGCh7joCExMThf/+wx/+UPjl\nUCwWU19fn2ZmZpTJZJROp9XS0qLKykqtWLFCyWRSzjkdOXJEHR0di3cGAIB5K/qdQGdnp/785z/r\n3//+t1auXKl9+/ZpaGhIZ8+eled5WrVqlX71q19JkiKRiLZs2aJIJKKysjIdOnRInudJkg4dOqQf\n/OAHevfdd/XEE0/M+aUwAGDpeW6un/qUgOd5Ki9v16VLr5R6KgDwKeLN+YvNj4srhgHAMCIAAIYR\nAQAwjAgAgGFEAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEBADCMCACAYUQAAAwjAgBgGBEAAMOI\nAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEQAAw4gAABhGBADAMCIAAIYRAQAwjAgAgGFE\nAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEBADCMCACAYUQAAAwjAgBgGBEAAMOIAAAYRgQAwDAi\nAACGEQEAMIwIAIBhRAAADCMCAGBY0Qjs3LlTfr9fjY2NhW1TU1OKRqOqq6tTW1ubpqenC/f19PQo\nFAopHA5rcHCwsP306dNqbGxUKBTSnj17PoHTAADMR9EI7NixQ4lEYta2eDyuaDSqkZERtba2Kh6P\nS5JSqZT6+/uVSqWUSCTU3d0t55wkaffu3Tp8+LDS6bTS6fQtxwQAlEbRCDzyyCP64he/OGvbsWPH\n1NXVJUnq6urS0aNHJUkDAwPq7OyUz+dTdXW1amtrlUwmNTExocuXL6ulpUWStH379sI+AIDSKrvb\nHfL5vPx+vyTJ7/crn89LksbHx7V+/frC44LBoHK5nHw+n4LBYGF7IBBQLpeb89jvvTciae/NW4/e\nHACAjwzdHIvjriPw/3meJ8/zFmsuuueeOl29unfRjgcAnz2PavYb5H0LOtpd/zrI7/drcnJSkjQx\nMaGKigpJH7zDz2azhceNjY0pGAwqEAhobGxs1vZAILCgSQMAFsddRyAWi6m3t1eS1Nvbq46OjsL2\nvr4+zczMKJPJKJ1Oq6WlRZWVlVqxYoWSyaScczpy5EhhHwBAibkitm7d6qqqqpzP53PBYNC9/PLL\n7sKFC661tdWFQiEXjUbdxYsXC4/fv3+/q6mpcfX19S6RSBS2nzp1yjU0NLiamhr39NNPz/lcklx5\nebuTHIPBYDA+9lCxl/E78tyHv+MsMc/zVF7erkuXXin1VADgU8TTQl7GuWIYAAwjAgBgGBEAAMOI\nAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEQAAw4gAABhGBADAMCIAAIYRAQAwjAgAgGFE\nAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEBADCMCACAYUQAAAwjAgBgGBEAAMOIAAAYRgQAwDAi\nAACGEQEAMIwIAIBhRAAADCMCAGAYEQAAw4gAABhGBADAMCIAAIYRAQAwjAgAgGFEAAAMIwIAYBgR\nAADDiAAAGEYEAMAwIgAAhhEBADBs3hGorq7Wgw8+qObmZrW0tEiSpqamFI1GVVdXp7a2Nk1PTxce\n39PTo1AopHA4rMHBwYXPHACwYPOOgOd5Ghoa0pkzZ3Ty5ElJUjweVzQa1cjIiFpbWxWPxyVJqVRK\n/f39SqVSSiQS6u7u1o0bNxbnDAAA87agj4Occ7NuHzt2TF1dXZKkrq4uHT16VJI0MDCgzs5O+Xw+\nVVdXq7a2thAOAEDplM13R8/z9K1vfUuf//zn9aMf/UhPPfWU8vm8/H6/JMnv9yufz0uSxsfHtX79\n+sK+wWBQuVzulmO+996IpL03bz16cwAAPjJ0cyyOeUfgL3/5i6qqqvSvf/1L0WhU4XB41v2e58nz\nvNvuP9d999xTp6tX9853SgBgwKOa/QZ534KONu+Pg6qqqiRJDzzwgDZt2qSTJ0/K7/drcnJSkjQx\nMaGKigpJUiAQUDabLew7NjamQCCwkHkDABbBvCLwn//8R5cvX5YkvfPOOxocHFRjY6NisZh6e3sl\nSb29vero6JAkxWIx9fX1aWZmRplMRul0uvCLIgBA6czr46B8Pq9NmzZJkt5//319//vfV1tbm9at\nW6ctW7bo8OHDqq6u1u9//3tJUiQS0ZYtWxSJRFRWVqZDhw4V/agIALA0PPffP/EpEc/zVF7erkuX\nXin1VADgU8S75Zead4MrhgHAMCIAAIYRAQAwjAgAgGFEAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAA\nhhEBADCMCACAYUQAAAwjAgBgGBEAAMOIAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEQAA\nw4gAABhGBADAMCIAAIYRAQAwjAgAgGFEAAAMIwIAYBgRAADDiAAAGEYEAMAwIgAAhhEBADCMCACA\nYUQAAAwjAgBgGBEAAMOIAAAYRgQAwDAiAACGEQEAMIwIAIBhRAAADCMCAGAYEViWhko9gWVkqNQT\nWEaGSj2BZWSo1BP4zFjSCCQSCYXDYYVCIT333HNL+dSfMkOlnsAyMlTqCSwjQ6WewDIyVOoJfGYs\nWQSuX7+uH//4x0okEkqlUvrd736n8+fPL9XTAwDmsGQROHnypGpra1VdXS2fz6etW7dqYGBgqZ4e\nADCHsqV6olwup5UrVxZuB4NBJZPJWY+5dOlVSd5STWmZ21fqCSwjrMVHWIuPsBaLYcki4HnFX9yd\nc0s0EwDAh5bs46BAIKBsNlu4nc1mFQwGl+rpAQBzWLIIrFu3Tul0WqOjo5qZmVF/f79isdhSPT0A\nYA5L9nFQWVmZfvnLX2rDhg26fv26du3apdWrVy/V0wMA5rCk1wk8/vjj+vvf/65//OMf+ulPf1rY\nbvn6gWw2q8cee0xr1qxRQ0ODXnzxRUnS1NSUotGo6urq1NbWpunp6RLPdOlcv35dzc3N2rhxoyS7\nazE9Pa0nn3xSq1evViQSUTKZNLsWPT09WrNmjRobG/W9731PV69eNbMWO3fulN/vV2NjY2FbsXPv\n6elRKBRSOBzW4ODgHY9f8iuGrV8/4PP59MILL+jcuXM6ceKEXnrpJZ0/f17xeFzRaFQjIyNqbW1V\nPB4v9VSXzMGDBxWJRAo/JrC6Fnv27NETTzyh8+fP669//avC4bDJtRgdHdWvf/1rDQ8P629/+5uu\nX7+uvr4+M2uxY8cOJRKJWdtud+6pVEr9/f1KpVJKJBLq7u7WjRs3ij+BK7G33nrLbdiwoXC7p6fH\n9fT0lHBGpfWd73zH/elPf3L19fVucnLSOefcxMSEq6+vL/HMlkY2m3Wtra3u9ddfd+3t7c45Z3It\npqen3apVq27ZbnEtLly44Orq6tzU1JS7du2aa29vd4ODg6bWIpPJuIaGhsLt2537gQMHXDweLzxu\nw4YN7u233y567JL/TWCu6wdyuVwJZ1Q6o6OjOnPmjB5++GHl83n5/X5Jkt/vVz6fL/HslsZPfvIT\n/exnP9PnPvfR/5oW1yKTyeiBBx7Qjh079PWvf11PPfWU3nnnHZNr8aUvfUnPPvusvvrVr+orX/mK\n7r//fkWjUZNr8aHbnfv4+PisX11+nNfTkkfgTtcPWHHlyhVt3rxZBw8e1H333TfrPs/zTKzTq6++\nqoqKCjU3N9/2uhEra/H+++9reHhY3d3dGh4e1r333nvLxx1W1uKf//ynfv7zn2t0dFTj4+O6cuWK\nfvvb3856jJW1mMudzv1O61LyCHD9gHTt2jVt3rxZ27ZtU0dHh6QP6j45OSlJmpiYUEVFRSmnuCTe\neustHTt2TKtWrVJnZ6def/11bdu2zeRaBINBBYNBPfTQQ5KkJ598UsPDw6qsrDS3FqdOndI3vvEN\nffnLX1ZZWZm++93v6u233za5Fh+63Z+J/349HRsbUyAQKHqskkfA+vUDzjnt2rVLkUhEzzzzTGF7\nLBZTb2+vJKm3t7cQh8+yAwcOKJvNKpPJqK+vT9/85jd15MgRk2tRWVmplStXamRkRJL02muvac2a\nNdq4caO5tQiHwzpx4oTeffddOef02muvKRKJmFyLD93uz0QsFlNfX59mZmaUyWSUTqfV0tJS/GCL\n/QXGfPzxj390dXV1rqamxh04cKDU01lSb775pvM8z61du9Y1NTW5pqYmd/z4cXfhwgXX2trqQqGQ\ni0aj7uLFi6We6pIaGhpyGzdudM45s2tx9uxZt27dOvfggw+6TZs2uenpabNr8dxzz7lIJOIaGhrc\n9u3b3czMjJm12Lp1q6uqqnI+n88Fg0H38ssvFz33/fv3u5qaGldfX+8SicQdj+85xz/aAwBWlfzj\nIABA6RABADCMCACAYUQAAAwjAgBgGBEAAMP+D3dxABU5HEp1AAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 94 }, { "cell_type": "code", "collapsed": false, "input": [ "#generate some data\n", "x = BiGoRNA[\"Expression value\"]\n", "y = BiGoRNA[\"Total intron reads\"]\n", "\n", "#plot the data\n", "plt.plot(x, y, 'bo')\n", "plt.show()" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAD9CAYAAACBdWEIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X9QVHee7/9nG5hYMxklsiM4tDNo0x1FCJIV5G5Vdkkx\ngJoJmpBIzAzgxOSbJd8xJpuyrCRjxJ01mJqyatXRu7lztYJmN2TK2qg1DghlhpjKTLD8tbeu7s2w\nbKvQIDtB8JJEQOR9/yD0iOCPqKdp8PWosqr59Pmc8z4knFd/zvmc0y4zM0RERG6zcSNdgIiIjE0K\nGBERcYQCRkREHKGAERERRyhgRETEEQoYERFxxDUDpquri7lz5zJ79mwSExN55ZVXACgtLcXtdpOa\nmkpqaiqVlZXBPmVlZXi9XmbMmEF1dXWw/ciRIyQnJ+P1elmxYkWwvbu7m4KCArxeLxkZGZw+fTr4\nXnl5OT6fD5/Px44dO27bTouISAjYdXzxxRdmZnbx4kWbO3euffTRR1ZaWmobNmwYsuyJEycsJSXF\nenp6zO/3m8fjsb6+PjMzS0tLs7q6OjMzmz9/vlVWVpqZ2ZYtW6ykpMTMzCoqKqygoMDMzNra2mz6\n9OnW3t5u7e3twdciIjI6XPcU2Te/+U0Aenp6uHTpEvfee+9AMA1Zds+ePSxZsoTIyEji4+NJSEig\nrq6OlpYWOjs7SU9PB6CoqIjdu3cDsHfvXoqLiwHIz8/nwIEDAOzfv5+cnByioqKIiooiOzubqqqq\n2xCpIiISChHXW6Cvr48HHniAhoYGSkpKmDVrFrt27WLz5s3s2LGDOXPmsGHDBqKiomhubiYjIyPY\n1+12EwgEiIyMxO12B9vj4uIIBAIABAIBpk6d2l9MRAQTJ06kra2N5ubmQX0G1nU5l8t1a3svInKH\nGm6QcLtddwQzbtw4jh8/TlNTEwcPHqS2tpaSkhL8fj/Hjx9nypQpvPzyy44XejVmFvb/1qxZM+I1\nqE7VqDpV58C/ULnhWWQTJ07k4Ycf5vDhw0yePBmXy4XL5eKZZ57h0KFDQP/IpLGxMdinqakJt9tN\nXFwcTU1NQ9oH+pw5cwaA3t5ezp8/T3R09JB1NTY2DhrRiIhIeLtmwHz22Wd0dHQAcOHCBWpqakhN\nTeXs2bPBZd5//32Sk5MByMvLo6Kigp6eHvx+P/X19aSnpxMbG8uECROoq6vDzNi5cycLFy4M9ikv\nLwdg165dZGVlAZCTk0N1dTUdHR20t7dTU1NDbm7u7f8NiIiII655DaalpYXi4mL6+vro6+ujsLCQ\nrKwsioqKOH78OC6Xi2nTpvHWW28BkJiYyOLFi0lMTCQiIoKtW7cGr5Ns3bqVpUuXcuHCBRYsWMC8\nefMAWLZsGYWFhXi9XqKjo6moqABg0qRJrF69mrS0NADWrFlDVFSUY78IJ2VmZo50CTdEdd4+o6FG\nUJ2322ipM1RcFsoTcreZy+UK6flEEZGxIFTHTt3JLyIijlDAiIiIIxQwIiLiCAWMiIg4QgEjIiKO\nUMCIiIgjFDAiIuIIBYyIiDhCASMiIo5QwIiIiCMUMCIi4ggFjIiIOEIBIyIijlDAiIiIIxQwIiLi\niFEfMJmZpeTm/ox9+w6OdCkiInKZa36j5Wjw4YelADQ0vAbAww//9QhWIyIiA0b9CGZAQ8M6Nm+u\nGekyRETkK2MmYAC6uu4a6RJEROQrYypgxo+/NNIliIjIV8ZMwHg8r7J8efZIlyEiIl8Z9Rf5/+Zv\nShk//hLLl8/TBX4RkTDiMjMb6SJulsvlYhSXLyIyIkJ17LzmKbKuri7mzp3L7NmzSUxM5JVXXgHg\n3LlzZGdn4/P5yMnJoaOjI9inrKwMr9fLjBkzqK6uDrYfOXKE5ORkvF4vK1asCLZ3d3dTUFCA1+sl\nIyOD06dPB98rLy/H5/Ph8/nYsWPHbdtpEREJAbuOL774wszMLl68aHPnzrWPPvrIVq5caW+++aaZ\nma1fv95WrVplZmYnTpywlJQU6+npMb/fbx6Px/r6+szMLC0tzerq6szMbP78+VZZWWlmZlu2bLGS\nkhIzM6uoqLCCggIzM2tra7Pp06dbe3u7tbe3B19f7gbKFxGRK4Tq2HndazDf/OY3Aejp6eHSpUvc\ne++97N27lw8//BCA4uJiMjMzWb9+PXv27GHJkiVERkYSHx9PQkICdXV1fP/736ezs5P09HQAioqK\n2L17N/PmzWPv3r2sXbsWgPz8fH76058CsH//fnJycoiKigIgOzubqqoqnnzyyUH1lZaWBl9nZmaS\nmZl5K3krIjLm1NbWUltbG/LtXjdg+vr6eOCBB2hoaKCkpIRZs2bR2tpKTEwMADExMbS2tgLQ3NxM\nRkZGsK/b7SYQCBAZGYnb7Q62x8XFEQgEAAgEAkydOrW/mIgIJk6cSFtbG83NzYP6DKzrSpcHjIiI\nDHXlh++BD/VOu27AjBs3juPHj3P+/Hlyc3P53e9+N+h9l8uFy+VyrEARERmdbvg+mIkTJ/Lwww9z\n5MgRYmJiOHv2LAAtLS1MnjwZ6B+ZNDY2Bvs0NTXhdruJi4ujqalpSPtAnzNnzgDQ29vL+fPniY6O\nHrKuxsbGQSMaEREJb9cMmM8++yw4Q+zChQvU1NSQmppKXl4e5eXlQP9Mr0WLFgGQl5dHRUUFPT09\n+P1+6uvrSU9PJzY2lgkTJlBXV4eZsXPnThYuXBjsM7CuXbt2kZWVBUBOTg7V1dV0dHTQ3t5OTU0N\nubm5zvwWRETktrvmKbKWlhaKi4vp6+ujr6+PwsJCsrKySE1NZfHixWzbto34+Hh+/etfA5CYmMji\nxYtJTEwkIiKCrVu3Bk+fbd26laVLl3LhwgUWLFjAvHnzAFi2bBmFhYV4vV6io6OpqKgAYNKkSaxe\nvZq0tDQA1qxZE7zgLyIi4U83WoqI3GHC4kZLERGRm6WAERERRyhgRETEEQoYERFxhAJGREQcoYAR\nERFHKGBERMQRChgREXGEAkZERByhgBEREUcoYERExBEKGBERcYQCRkREHKGAERERRyhgRETEEQoY\nERFxhAJGREQcoYARERFHKGBERMQRChgREXGEAkZERByhgBEREUdcM2AaGxt56KGHmDVrFklJSWza\ntAmA0tJS3G43qamppKamUllZGexTVlaG1+tlxowZVFdXB9uPHDlCcnIyXq+XFStWBNu7u7spKCjA\n6/WSkZHB6dOng++Vl5fj8/nw+Xzs2LHjtu20iIiEgF1DS0uLHTt2zMzMOjs7zefz2cmTJ620tNQ2\nbNgwZPkTJ05YSkqK9fT0mN/vN4/HY319fWZmlpaWZnV1dWZmNn/+fKusrDQzsy1btlhJSYmZmVVU\nVFhBQYGZmbW1tdn06dOtvb3d2tvbg68vd53yRURkGKE6dl5zBBMbG8vs2bMBuOeee5g5cyaBQGAg\nmIYsv2fPHpYsWUJkZCTx8fEkJCRQV1dHS0sLnZ2dpKenA1BUVMTu3bsB2Lt3L8XFxQDk5+dz4MAB\nAPbv309OTg5RUVFERUWRnZ1NVVXVbQlVERFxXsSNLnjq1CmOHTtGRkYGH3/8MZs3b2bHjh3MmTOH\nDRs2EBUVRXNzMxkZGcE+brebQCBAZGQkbrc72B4XFxcMqkAgwNSpU/uLiYhg4sSJtLW10dzcPKjP\nwLquVFpaGnydmZlJZmbmDe+8iMidoLa2ltra2pBv94YC5vPPP+fxxx9n48aN3HPPPZSUlPD6668D\nsHr1al5++WW2bdvmaKFXc3nAiIjIUFd++F67dm1ItnvdWWQXL14kPz+fH//4xyxatAiAyZMn43K5\ncLlcPPPMMxw6dAjoH5k0NjYG+zY1NeF2u4mLi6OpqWlI+0CfM2fOANDb28v58+eJjo4esq7GxsZB\nIxoREQlv1wwYM2PZsmUkJiby4osvBttbWlqCr99//32Sk5MByMvLo6Kigp6eHvx+P/X19aSnpxMb\nG8uECROoq6vDzNi5cycLFy4M9ikvLwdg165dZGVlAZCTk0N1dTUdHR20t7dTU1NDbm7u7d17ERFx\nzDVPkX388ce888473H///aSmpgLwxhtv8O6773L8+HFcLhfTpk3jrbfeAiAxMZHFixeTmJhIREQE\nW7duxeVyAbB161aWLl3KhQsXWLBgAfPmzQNg2bJlFBYW4vV6iY6OpqKiAoBJkyaxevVq0tLSAFiz\nZg1RUVHO/BZEROS2c9lw08FGCZfLNexsNhERubpQHTt1J7+IiDhCASMiIo5QwIiIiCMUMCIi4ggF\njIiIOEIBIyIijlDAiIiIIxQwIiLiCAWMiIg4QgEjIiKOUMCIiIgjFDAiIuIIBYyIiDhCASMiIo5Q\nwIiIiCMUMCIi4ggFjIiIOEIBIyIijlDAiIiIIxQwIiLiCAWMiIg4QgEjIiKOUMCIiIgjrhkwjY2N\nPPTQQ8yaNYukpCQ2bdoEwLlz58jOzsbn85GTk0NHR0ewT1lZGV6vlxkzZlBdXR1sP3LkCMnJyXi9\nXlasWBFs7+7upqCgAK/XS0ZGBqdPnw6+V15ejs/nw+fzsWPHjtu20yIiEgJ2DS0tLXbs2DEzM+vs\n7DSfz2cnT560lStX2ptvvmlmZuvXr7dVq1aZmdmJEycsJSXFenp6zO/3m8fjsb6+PjMzS0tLs7q6\nOjMzmz9/vlVWVpqZ2ZYtW6ykpMTMzCoqKqygoMDMzNra2mz69OnW3t5u7e3twdeXu075IiIyjFAd\nO685gomNjWX27NkA3HPPPcycOZNAIMDevXspLi4GoLi4mN27dwOwZ88elixZQmRkJPHx8SQkJFBX\nV0dLSwudnZ2kp6cDUFRUFOxz+bry8/M5cOAAAPv37ycnJ4eoqCiioqLIzs6mqqrqtgesiIg4I+JG\nFzx16hTHjh1j7ty5tLa2EhMTA0BMTAytra0ANDc3k5GREezjdrsJBAJERkbidruD7XFxcQQCAQAC\ngQBTp07tLyYigokTJ9LW1kZzc/OgPgPrulJpaWnwdWZmJpmZmTe6SyIid4Ta2lpqa2tDvt0bCpjP\nP/+c/Px8Nm7cyLe//e1B77lcLlwulyPF3YjLA0ZERIa68sP32rVrQ7Ld684iu3jxIvn5+RQWFrJo\n0SKgf9Ry9uxZAFpaWpg8eTLQPzJpbGwM9m1qasLtdhMXF0dTU9OQ9oE+Z86cAaC3t5fz588THR09\nZF2NjY2DRjQiIhLerhkwZsayZctITEzkxRdfDLbn5eVRXl4O9M/0GgievLw8Kioq6Onpwe/3U19f\nT3p6OrGxsUyYMIG6ujrMjJ07d7Jw4cIh69q1axdZWVkA5OTkUF1dTUdHB+3t7dTU1JCbm3v7fwMi\nIuKMa80A+Oijj8zlcllKSorNnj3bZs+ebZWVldbW1mZZWVnm9XotOzt70OyudevWmcfjsfvuu8+q\nqqqC7YcPH7akpCTzeDy2fPnyYHtXV5c98cQTlpCQYHPnzjW/3x98b/v27ZaQkGAJCQn29ttvD6nv\nOuWLiMgwQnXsdH21sVHJ5XIxissXERkRoTp26k5+ERFxhAJGREQcoYARERFHKGBERMQRChgREXGE\nAkZERByhgBEREUcoYERExBEKGBERcYQCRkREHKGAERERRyhgRETEEQoYERFxhAJGREQcoYARERFH\nKGBERMQRChgREXGEAkZERByhgBEREUcoYERExBEKGBERcYQCRkREHKGAERERR1wzYJ5++mliYmJI\nTk4OtpWWluJ2u0lNTSU1NZXKysrge2VlZXi9XmbMmEF1dXWw/ciRIyQnJ+P1elmxYkWwvbu7m4KC\nArxeLxkZGZw+fTr4Xnl5OT6fD5/Px44dO27LzoqISAjZNRw8eNCOHj1qSUlJwbbS0lLbsGHDkGVP\nnDhhKSkp1tPTY36/3zwej/X19ZmZWVpamtXV1ZmZ2fz5862ystLMzLZs2WIlJSVmZlZRUWEFBQVm\nZtbW1mbTp0+39vZ2a29vD76+0nXKFxGRYYTq2HnNEcyDDz7IvffeO1woDWnbs2cPS5YsITIykvj4\neBISEqirq6OlpYXOzk7S09MBKCoqYvfu3QDs3buX4uJiAPLz8zlw4AAA+/fvJycnh6ioKKKiosjO\nzqaqqupWclREREIs4mY6bd68mR07djBnzhw2bNhAVFQUzc3NZGRkBJdxu90EAgEiIyNxu93B9ri4\nOAKBAACBQICpU6f2FxIRwcSJE2lra6O5uXlQn4F1Dae0tDT4OjMzk8zMzJvZJRGRMau2tpba2tqQ\nb/drB0xJSQmvv/46AKtXr+bll19m27Ztt72wG3V5wIiIyFBXfvheu3ZtSLb7tWeRTZ48GZfLhcvl\n4plnnuHQoUNA/8iksbExuFxTUxNut5u4uDiampqGtA/0OXPmDAC9vb2cP3+e6OjoIetqbGwcNKIR\nEZHw97UDpqWlJfj6/fffD84wy8vLo6Kigp6eHvx+P/X19aSnpxMbG8uECROoq6vDzNi5cycLFy4M\n9ikvLwdg165dZGVlAZCTk0N1dTUdHR20t7dTU1NDbm7uLe+siIiEzjVPkS1ZsoQPP/yQzz77jKlT\np7J27Vpqa2s5fvw4LpeLadOm8dZbbwGQmJjI4sWLSUxMJCIigq1bt+JyuQDYunUrS5cu5cKFCyxY\nsIB58+YBsGzZMgoLC/F6vURHR1NRUQHApEmTWL16NWlpaQCsWbOGqKgox34JIiJy+7lsuClho4TL\n5Rp2RpuIiFxdqI6dupNfREQcoYARERFHKGBERMQRChgREXGEAkZERByhgBEREUcoYERExBEKGBER\ncYQCRkREHKGAERERRyhgRETEEQoYERFxhAJGREQccVNfmSyj2759B9m0qZru7gjuvruXF17I4eGH\n/3qkyxKRMUYBc4fZt+8gK1bsp6FhXbCtoeE1AIWMiNxWOkV2h9m0qXpQuAA0NKxj8+aaEapIRMYq\nBcwdprt7+EFrV9ddIa5ERMY6Bcwd5u67e4dtHz/+UogrEZGxTgFzh3nhhRw8ntcGtXk8r7J8efYI\nVSQiY5XLRvGX2ofqe6XHmn37DrJ5cw1dXXcxfvwlli/P1gV+kTtIqI6dChgRkTtMqI6dOkUmIiKO\nUMCIiIgjrhkwTz/9NDExMSQnJwfbzp07R3Z2Nj6fj5ycHDo6OoLvlZWV4fV6mTFjBtXV1cH2I0eO\nkJycjNfrZcWKFcH27u5uCgoK8Hq9ZGRkcPr06eB75eXl+Hw+fD4fO3bsuC07KyIioXPNgPnJT35C\nVVXVoLb169eTnZ3NH//4R7Kysli/fj0AJ0+e5L333uPkyZNUVVXx/PPPB8/xlZSUsG3bNurr66mv\nrw+uc9u2bURHR1NfX89LL73EqlWrgP4Q+/u//3sOHTrEoUOHWLt27aAgExGR8HfNgHnwwQe59957\nB7Xt3buX4uJiAIqLi9m9ezcAe/bsYcmSJURGRhIfH09CQgJ1dXW0tLTQ2dlJeno6AEVFRcE+l68r\nPz+fAwcOALB//35ycnKIiooiKiqK7OzsIUEnIiLh7Ws/i6y1tZWYmBgAYmJiaG1tBaC5uZmMjIzg\ncm63m0AgQGRkJG63O9geFxdHIBAAIBAIMHXq1P5CIiKYOHEibW1tNDc3D+ozsK7hlJaWBl9nZmaS\nmZn5dXdJRGRMq62tpba2NuTbvaWHXbpcLlwu1+2q5aZcHjAiIjLUlR++165dG5Ltfu1ZZDExMZw9\nexaAlpYWJk+eDPSPTBobG4PLNTU14Xa7iYuLo6mpaUj7QJ8zZ84A0Nvby/nz54mOjh6yrsbGxkEj\nGhERCX9fO2Dy8vIoLy8H+md6LVq0KNheUVFBT08Pfr+f+vp60tPTiY2NZcKECdTV1WFm7Ny5k4UL\nFw5Z165du8jKygIgJyeH6upqOjo6aG9vp6amhtzc3NuywyIiEiJ2DU8++aRNmTLFIiMjze122/bt\n262trc2ysrLM6/Vadna2tbe3B5dft26deTweu++++6yqqirYfvjwYUtKSjKPx2PLly8Ptnd1ddkT\nTzxhCQkJNnfuXPP7/cH3tm/fbgkJCZaQkGBvv/32sPVdp3wRERlGqI6delSMiMgdRo+KERGRUU0B\nIyIijlDAiIiII27pPhgZGfv2HWTTpmq6uyO4++5eXnghR9/nIiJhRwEzyuzbd5AVK/bT0LAu2NbQ\n0P8NlQoZEQknOkU2ymzaVD0oXAAaGtaxeXPNCFUkIjI8Bcwo0909/KCzq+uuEFciInJtCphR5u67\ne4dtHz/+UogrERG5NgXMKPPCCzl4PK8NavN4XmX58uwRqkhEZHi6k38U2rfvIJs319DVdRfjx19i\n+fJsXeAXkRsWqmOnAkZE5A6jR8WIiMiopoARERFHKGBERMQRChgREXGEAkZERByhZ5GJiIxxVz4g\nN1Q0TVlEZAwb7gG5oGnKIiJyi4Z7QG6oKGBERMawqz0gNxQUMCIiY1gor7lcSQEjIjKGDfeA3FDR\nLDIRkTFs4EG4mzevDj4gd//+0Gz7pmeRxcfHM2HCBO666y4iIyM5dOgQ586do6CggNOnTxMfH8+v\nf/1roqKiACgrK2P79u3cddddbNq0iZycHACOHDnC0qVL6erqYsGCBWzcuBGA7u5uioqKOHr0KNHR\n0bz33nt8//vfH1y8ZpGJiHxtYf+wS5fLRW1tLceOHePQoUMArF+/nuzsbP74xz+SlZXF+vXrATh5\n8iTvvfceJ0+epKqqiueffz64cyUlJWzbto36+nrq6+upqqoCYNu2bURHR1NfX89LL73EqlWrbnVf\nRUQkhG7pGsyVCbh3716Ki4sBKC4uZvfu3QDs2bOHJUuWEBkZSXx8PAkJCdTV1dHS0kJnZyfp6ekA\nFBUVBftcvq78/HwOHDhwK6WKiEiI3fQ1GJfLxQ9+8APuuusunnvuOZ599llaW1uJiYkBICYmhtbW\nVgCam5vJyMgI9nW73QQCASIjI3G73cH2uLg4AoEAAIFAgKlTp/YXGRHBxIkTOXfuHJMmTRpUR2lp\nafB1ZmYmmZmZN7tLIiKOuvKO+hdeyAnJlwXW1tZSW1vr+HaudNMB8/HHHzNlyhT+9Kc/kZ2dzYwZ\nMwa973K5cLlct1zg9VweMCIi4Wq4O+obGvpndzkdMld++F67dq2j2xtw06fIpkyZAsB3vvMdHn30\nUQ4dOkRMTAxnz54FoKWlhcmTJwP9I5PGxsZg36amJtxuN3FxcTQ1NQ1pH+hz5swZAHp7ezl//vyQ\n0YuIyGgx3B31DQ3r2Ly5ZoQqct5NBcyXX35JZ2cnAF988QXV1dUkJyeTl5dHeXk5AOXl5SxatAiA\nvLw8Kioq6Onpwe/3U19fT3p6OrGxsUyYMIG6ujrMjJ07d7Jw4cJgn4F17dq1i6ysrFveWRGRkXK1\nO+q7uu4KcSWhc1OnyFpbW3n00UeB/tHFj370I3JycpgzZw6LFy9m27ZtwWnKAImJiSxevJjExEQi\nIiLYunVr8PTZ1q1bWbp0KRcuXGDBggXMmzcPgGXLllFYWIjX6yU6OpqKiorbsb8iIiPianfUjx9/\nKcSVhI6epiwiEgLDXYPxeF5l48Z5IbnQf7lQHTsVMCIiIbJv30E2b64J3lG/fHl2yMMFFDA3xOVy\nkZPzWsim+kn4G6lpoCKjSagCZtQ/i6y6+h9CNtVPwttITgMVkaHGxNOUx/pUP7kxd+I0UJFwNiYC\nBsb2VD+5MXfiNFCRcDYGAuZnwEH+9//+d/btOzjSxcgIuhOngYqEszEQMP8A/AttbX/DihX7FTJ3\nsOG+WMnjeZXly7NHqCKRO9uon0UGA+U/C8Rwzz0N3H13H7GxUcTFfUeziO4w4TINVCScaZryDfhz\nwBwE3gX++2Xvvgbk4vHsZ+PGXB1kRES+EvZfOBZeqhkcLgDrgBrNIhIRGSFjJGAar9LeP3tIs4hE\nREJvDATMauCeq7zXP3tIs4hEREJvDATMz4ECxo1bdkX7q0C2ZhGJiIyQMXKR/1W+9716Zs6cQXPz\n57S0tBATMxG3e/Jtm0WkZ1yJyFihZ5HdsEeARGbOnEFV1c8d2YKecSUi8vWNkRHMs8ya9RlxcbMc\nGWHk5v6M6up/GKZ9tWOhdi0aTYnIrdAI5mv5FSdO5HLixHeB5wH4X/9rGd/85v+ktfU83d19jB9/\nD17vJH7+84KrHoyvduAOp2dcaTQlIqPFGAkYgEn032z5W+BbnD3bBeQCAWAdn38Ox47BY4/9f3i9\n/zzkLv9rHbhv9BlXoRhZXP2JwasVMCISVsZAwJQCvcCXwHeB/5/+Gy8jgEpg5aCle3r+BydOrObE\niZ8P+uR/rQP3Cy/k0NDw2pCvOl2+fF7w51CNLMJpNCUici1jYJry/6H/RstewA9soT9ceukPl/30\nP0rmcv0H48vv8r/egXvChA7uvbeYe+99kgceeHbI92iH6rtI9MRgERktxsAI5nn6RyzTgGPAn4Dx\n9N98+UsgE6gBLh9F/PlgfOjQGfbtO3jVA/eRI4d45pn/y9mzW4Jtd9/9d6xevYNf/OKD4KmwqwXU\nJ580kpv7s2FPl93MKbUbGU2JiISDMTCL7FX6nzs24G+Bp+gPlL+jf4QTCbz81ftb6D+V9i0gB6jB\n4+njxz+O4513AleMQp4GPgd+PczWV9N/kyd4PK8xYUIHx45tuepyHs9rgx66OdwptdjYv2PKlP/L\nhAnuawaOnhgsIrciVLPIsFEMMLBh/j1j8JrBGoNHDH5gsMig+IrlfvLVe/MM5hrkfrVcwVfr+NDg\niats48VBPz/wwDPm8bx6xTKvfLWO3xmY5eb+LFh7Ts5rV1nvz4KvPZ5X7Te/+TBkv8/f/e53IdvW\nrRgNdY6GGs1U5+02WuoM1aE/rK/BVFVVMWPGDLxeL2+++eZVlpoPJF3280H6z/z9A/0TAPYC3wfu\nBt6+ou92IIP+yQDJ9I+G3gcqgMlfLeO6ynZbBv307W/HsXFjLhERjwBL6R+5zKN/JFULDL4Qf7VT\nagPXhyD03ydfW1sbsm3ditFQ52ioEVTn7TZa6gyVsA2YS5cu8dOf/pSqqipOnjzJu+++y7//+78P\ns2Ql8N/UVBmQAAAIk0lEQVT4c8gM9+j+/wlcuMqWBg7ov6L/Ws2AdV/9HEX/d8tc7lVg4qCW8eMv\n8fDDf823vhUNuOk/ffbXQ5YZcLVrPpdfHwLNDhOR0StsA+bQoUMkJCQQHx9PZGQkTz75JHv27LnK\n0r8Cpn71uuEqy9x9lfbLD+hXHszvAr5D//00q+kfEQ2MTCYHl7r8gZoRERfov7YzOJTGjXt60EM3\nh/t634EHdF5Os8NEZLQK24v8u3btYv/+/fzqV78C4J133qGuro7NmzcHl+m/yC8iIl9XKA79YTtN\n+UbCI0yzUURECONTZHFxcTQ2/vmbKhsbG3G73SNYkYiIfB1hGzBz5syhvr6eU6dO0dPTw3vvvUde\nXt5IlyUiIjcobE+RRURE8Mtf/pLc3FwuXbrEsmXLmDlz5kiXJSIiNyhsRzAA8+fP59NPP+U//uM/\neOWVVwa9d2P3yNw+jY2NPPTQQ8yaNYukpCQ2bdoEwLlz58jOzsbn85GTk0NHR0ewT1lZGV6vlxkz\nZlBdXR1sP3LkCMnJyXi9XlasWBFs7+7upqCgAK/XS0ZGBqdPn77pei9dukRqaiqPPPJI2NbZ0dHB\n448/zsyZM0lMTKSuri4s6ywrK2PWrFkkJyfz1FNP0d3dHRZ1Pv3008TExJCcnBxsC1Vd5eXl+Hw+\nfD4fO3bs+Fo1rly5kpkzZ5KSksJjjz3G+fPnR7TGq9U5YMOGDYwbN45z586FbZ2bN29m5syZJCUl\nsWrVqhGvMygkt3PeZr29vebxeMzv91tPT4+lpKTYyZMnHd1mS0uLHTt2zMzMOjs7zefz2cmTJ23l\nypX25ptvmpnZ+vXrbdWqVWZmduLECUtJSbGenh7z+/3m8Xisr6/PzMzS0tKsrq7OzMzmz59vlZWV\nZma2ZcsWKykpMTOziooKKygouOl6N2zYYE899ZQ98sgjZmZhWWdRUZFt27bNzMwuXrxoHR0dYVen\n3++3adOmWVdXl5mZLV682N5+++2wqPPgwYN29OhRS0pKCraFoq62tjabPn26tbe3W3t7e/D1jdZY\nXV1tly5dMjOzVatWjXiNV6vTzOzMmTOWm5tr8fHx1tbWFpZ1fvDBB/aDH/zAenp6zMzsv/7rv0a8\nzgGjMmB+//vfW25ubvDnsrIyKysrC2kNCxcutJqaGrvvvvvs7NmzZtYfQvfdd5+Zmb3xxhu2fv36\n4PK5ubn2hz/8wZqbm23GjBnB9nfffdeee+654DKffPKJmfUfcP/iL/7ipmprbGy0rKws++CDD+yH\nP/yhmVnY1dnR0WHTpk0b0h5udba1tZnP57Nz587ZxYsX7Yc//KFVV1eHTZ1+v3/QwSYUdf3Lv/yL\n/e3f/m2wz3PPPWfvvvvuDdd4uX/913+1H/3oRyNe49XqfPzxx+3f/u3fBgVMuNX5xBNP2IEDB4Ys\nN9J1moX5o2KuJhAIMHXq1ODPbrebQCAQsu2fOnWKY8eOMXfuXFpbW4mJiQEgJiaG1tZWAJqbmwfN\nehuo8cr2uLi4YO2X71dERAQTJ04cNCy/US+99BK/+MUvGDfuz/95w61Ov9/Pd77zHX7yk5/wwAMP\n8Oyzz/LFF1+EXZ2TJk3i5Zdf5nvf+x7f/e53iYqKIjs7O+zqHOB0XW1tbVdd183Yvn07CxYsCMsa\n9+zZg9vt5v777x/UHm511tfXc/DgQTIyMsjMzOTw4cNhU+eoDJiRvMHy888/Jz8/n40bN/Ltb397\n0Hsul2vEb/78zW9+w+TJk0lNTb3qfULhUGdvby9Hjx7l+eef5+jRo3zrW99i/fr1g5YJhzobGhr4\nx3/8R06dOkVzczOff/4577zzzqBlwqHO4YRrXQPWrVvHN77xDZ566qmRLmWIL7/8kjfeeIO1a9cG\n26729zTSent7aW9v55NPPuEXv/gFixcvHumSgkZlwIzUPTIXL14kPz+fwsJCFi1aBPR/Sjx79iwA\nLS0tTJ48edgam5qacLvdxMXF0dTUNKR9oM+ZM2eA/v9pzp8/z6RJk75Wjb///e/Zu3cv06ZNY8mS\nJXzwwQcUFhaGXZ1utxu3201aWhoAjz/+OEePHiU2Njas6jx8+DB/9Vd/RXR0NBERETz22GP84Q9/\nCLs6Bzj93zk6Ovq2/P29/fbb/Pa3v+Wf//mfg23hVGNDQwOnTp0iJSWFadOm0dTUxF/+5V/S2toa\nVnVC/9/SY489BkBaWhrjxo3js88+C486r3sSLQxdvHjRpk+fbn6/37q7u0Nykb+vr88KCwvtxRdf\nHNS+cuXK4HnOsrKyIRcsu7u77T//8z9t+vTpwQts6enp9sknn1hfX9+QC2wD5znffffdW7rIb2ZW\nW1sbvAYTjnU++OCD9umnn5qZ2Zo1a2zlypVhV+fx48dt1qxZ9uWXX1pfX58VFRXZL3/5y7Cp88rz\n8aGoq62tzaZNm2bt7e127ty54OsbrbGystISExPtT3/606DlRrLG4eq83HAX+cOlzn/6p3+y119/\n3czMPv30U5s6dWpY1Gk2Si/ym5n99re/NZ/PZx6Px9544w3Ht/fRRx+Zy+WylJQUmz17ts2ePdsq\nKyutra3NsrKyzOv1WnZ29qBf+rp168zj8dh9991nVVVVwfbDhw9bUlKSeTweW758ebC9q6vLnnji\nCUtISLC5c+ea3++/pZpra2uDs8jCsc7jx4/bnDlz7P7777dHH33UOjo6wrLON9980xITEy0pKcmK\nioqsp6cnLOp88sknbcqUKRYZGWlut9u2b98esrq2b99uCQkJlpCQYG+//fYN17ht2zZLSEiw733v\ne8G/o4FZSyNV4+V1fuMb3wj+Li83bdq0YMCEW509PT324x//2JKSkuyBBx4Y9J00I1XngLB92KWI\niIxuo/IajIiIhD8FjIiIOEIBIyIijlDAiIiIIxQwIiLiCAWMiIg44v8BUdu3mbOzYzkAAAAASUVO\nRK5CYII=\n", "text": [ "" ] } ], "prompt_number": 77 }, { "cell_type": "code", "collapsed": false, "input": [ "fig = plt.loglog(x, y, 'rs')\n" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEHCAYAAABCwJb2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvX9sHNd57/1dUj/WCh0yurKkWGTuuJFbSlYQ3hcy+sKx\nQrpwqyiBZbmvRFP5wx5KboIUUuI/ijpNWvDYcFM7CPDeNrq1EaQp0waRlVzkSk2VMG+UmEwEBxIu\ngk1t0wvLtY9NypFkyZeUKWYpcTnvH2fOzpmz58zO7A/ucvl8AGJmzpxfs9x9zpnnPOd5Up7neSAI\ngiCWBS317gBBEASxeJDQJwiCWEaQ0CcIglhGkNAnCIJYRpDQJwiCWEaQ0CcIglhGkNAnCIJYRpDQ\nJwiCWEasqHaFnufhr//6r/Hee+9h+/bteOihh6rdBEEQBFEmVZ/pHz9+HOfPn8eqVavQ2dlZ7eoJ\ngiCICogl9A8cOIANGzbgIx/5SCh9ZGQE3d3duP322/H0008DAF599VV87GMfw9e+9jU888wz1e8x\nQRAEUTaxhP7g4CBGRkZCafl8HocOHcLIyAjGx8dx9OhRvPLKK+js7ERHR4eovIWWDAiCIBqJWDr9\nHTt2gHMeSjt79iw2b94Mx3EAAAMDAzhx4gS+8IUv4PDhw/jlL3+Jvr4+Y32pVKqSPhMEQSxbKvWR\nWfZU/Pz58+jq6ipcd3Z24vz587jpppvwzW9+E//wD/+Az33uc9byQ0NDeP755+F5Xll/Q0NDFecz\n3dPToq7jnNfy+ZI+W6M8X63+d9V4vsX635XzfKXS4jxrMz5fPX571Xi+JN/N559/HkNDQ+WK6xCt\njDEWJ+PU1BSOHj2KP//zPwcAvPLKK3jttdewe/duAMB//Md/4Pz58/jkJz9Zsq7HH38co6OjhbeE\ncolbPiqf6Z6eFnVtOh8dHbW+5SQhzvMlfTZTej2er1b/O1N6kudbzP9dqXxxvpt6Wqlnbdbnq8dv\nr1S+asqWvr4+9PX14fHHH0dMkW3Hi8kbb7zhbdu2rXD9q1/9ytu5c2fh+itf+Yr31FNPxaoLgDc0\nNOQ9//zzcZtfUgwNDdW7CzWlmZ+vmZ/N8+j5lirPP/+8NzQ05CUQ2VbKVu9s374d586dA+cc169f\nx7Fjxwqz/jgwxqoyIjcizfpckmZ+vmZ+NoCeb6nS19dX+QzfJ+V5XslVgf3792NsbAxXrlzB+vXr\n8cQTT2BwcBA//vGP8eijjyKfz+PgwYP4q7/6q3iNplKI0SxBEAShUA3ZGUvoV5tUKoWhoaGCnoog\nCIKwMzo6itHRUTz++ONLV+jTTJ8gCCIZ1ZCddds9xRjD6OhovZonCIJYMoyOji6uTr/a0EyfIAgi\nOUt6pk8QBEEsPqTeIQiCaHBIvUMQBLEMIfUOQRAEkQgS+gRBEMsI0ukTBEE0OKTTJ5oe5rqAFsMB\nAOA4YMPDi9wbgmgMqiE7qx4YnSCqAudgY2NFyWzxe0IQTQXp9AmCIJYRpNMnCIJocEinTzQ9rK/P\nrN7p7QWjyQKxTCE7fYIgCCIRtJBLNCaOY160rTCuMkEsd0i9QxAEsUQg9Q5BEASRCBL6BEEQywgy\n2SQIgmhwyGSTIAhiGUI6fYIgCCIRJPQJgiCWEST0CYIglhEk9AmCIJYRtCOXIIiaQ/ERGoeqC/3R\n0VH8zd/8DbZt24aBgQH09vZWuwmCIJYaFB+hYai6eqelpQU333wz5ubm0NnZWe3qCYIgiAqIJfQP\nHDiADRs24CMf+UgofWRkBN3d3bj99tvx9NNPAwB27NiBH/3oR3jqqacwNDRU/R4TBEEQZRNL6A8O\nDmJkZCSUls/ncejQIYyMjGB8fBxHjx7FK6+8glQqBQDo6OjA3Nxc9XtMEARBlE0snf6OHTvAtUWY\ns2fPYvPmzXB8V7cDAwM4ceIEstksfvKTn2BqagqHDx+21qluKe7r60NfX1/SvhME0WDYFmwz2eyi\n96UZGB0drbq7mrIXcs+fP4+urq7CdWdnJ86cOYMvfvGLeOCBB0qWr5YfCYIgGgjLgq27YQOYyaiD\n4iNEok+IH3/88YrrLFvoSzVOuTDGaIZPEMsEp7ubwlxWQDVn/GUL/U2bNmFiYqJwPTExkchah2b6\nBEEQ8ZAT5GrM9Ms22dy+fTvOnTsHzjmuX7+OY8eOYffu3bHLk2tlgiCIeFTTtXKsmf7+/fsxNjaG\nK1euoKurC0888QQGBwdx5MgR7Ny5E/l8HgcPHsSWLVtiN0wzfYIgFpulujO4mjP9WEL/6NGjxvRd\nu3Zh165dZTVMOn2CaEIaPaD9Et0ZXE2dPgVRIQhi2cD6+sxCv7d3SSw0L+kgKqTTJwiCiAeFSyQI\ngigDmumTa2WCIJYBcgGXZzL17krdqZvQp4VcgiAWDX8BlyG8aMvb2+H09DTOQrMFWsglCIJIwFJX\n60iW9EIuQRAEsfiQTp8oyVLd0EIsXeg7VztIp0+UZoluaCGWMPSdC9EQDtcqhdwwEASxaDT6TuES\nLLobBoIgiKUMqYQCSOgTyx7SHxPLCdLpE0SV9cdxBhEaaIgkkE6fWFyWuD500YkziNBCZTT0nQvR\nFDp9Jmf4NLNpeOj/Qyw29J2rHfUT+v4shyUpQ6/EBEEQFbG0FnLplZggCKIilpbQJ4haQPpjYhlB\nQp9Y9lRdNRhnEKGBhqgT9dPpA+irV+MEUUPiDCK0BkUkoZomm/Xzstnbi9HeXprZEARBlKCvr69q\nZu71m+mXM2pV4ZWYLIAIIgz9JpYXS8pOvypfQLIAIogw9JtYViwpO32CIAiiMsh6h2gqSFVBENGQ\n0CcS09CCtcFUFY30Wdn6wrPZmrfREN8NAkCNhP61a9cKq82f+tSnatEEsQjYfsCZbBbHL14szl/z\nHtWGmgqqRhqELH1x29tr3garXgtEhdRE6H/1q1/Fgw8+WIuqK4c2xcRnMYREAmomnBtMUKnPybNZ\nIJcDALw+N4ffW71aZEqn4XR3i/MKn38qnQbr6Sm+Qb+JpiSW0D9w4ABOnjyJ9evX48UXXyykj4yM\n4NFHH0U+n8cjjzyCxx57DD/96U+xdetW5Pwvqg3W2ytOFvmLRa+YS5gGE85xYYY0nsmAua75+2gb\nbAEMy9/V9DTgv22Z6k9CT3d3eSbUxJIkltAfHBzE4cOH8dBDDxXS8vk8Dh06hFOnTmHTpk248847\nsXv3boyNjeHatWsYHx/HTTfdhE9+8pNIpVJFdTbrl6zUbHQ56DyXwzNWijM9DX78eEHwh2b3mUyi\nungmE5hAFxqo/LOm/2NzEkvo79ixA1z75589exabN2+G48/UBwYGcOLECTz55JMAgG9/+9u45ZZb\njAIfCAdRaaoIWqVmo0t0tpqICp4xStDEIqb6LqmQlflVdQsgVCM93d2xBGFRv6anweSzKp+Zsf8R\nONPTRZ930jqMLIfvaoNTTfcLkrJ1+ufPn0dXV1fhurOzE2fOnClcP/zww5HlKXLWEqZMHXD29Gm4\nHR3R+ugKBU3sGWhSIWvr1/Q02MWL9jocBzyTgTM9ba7Wn6XzTAasRF+iFaYxWIz1LFozqyr6hLiu\nkbNsM/i4UIzcJYDlB+yU+Xrfnc+DTU+LGa6vj+aZDNyNGwuz58nZ2fL7q5BUNcH8I29vhyMHtCoI\nKjY8LFQvhgEDCM/SmTFHQLoKfak1pPapDQ0RI3fTpk2YmJgoXE9MTKCzszN2eZrpNz61/gEzQCxI\nKrNgN6pAxCxSF/Lq7FotYyofut/T0xDrTcw/vrRyJVYsLKA7n8cUwn3lvhWV7S2CaB4aIkbu9u3b\nce7cOXDOceutt+LYsWM4evRo7PJ9jgOnowNOTw/NDggAQqDlYJ/xRn1PWF+fUa/Ntfqk1UxkPwxv\nCXHUL1Z8FQ+mp4vKvwQx0E0B6EDxAMXuugs8my28Gam4vprMpHYimotFn+nv378fY2NjuHLlCrq6\nuvDEE09gcHAQR44cwc6dO5HP53Hw4EFs2bIldsOjb74JvPkmWEdH2Z1vSErpNJeDztN/RnW2zQE4\nMYo+Z0hjFXRlWE+YnoZ7/Dim0hHKkgTrCtxPjzLBlNY5jHPwTAbDhpk50+p329uxJ51GRzaL3OXL\nRW0zALnLl5HJZsW+CXWdBKjO92k5fFeXCIs+07fN4Hft2oVdu3aV1TBDcwZRKWnBUYO3mkpN6yop\nH1XWAUL6alayJ2Z4e7tYR6jUskdmn54GT6cLe0Uy2Sw6pEWOb50j+8qUcry9HW46jdzly+jO50Vd\nMo9qiWNAfo5R+v1QH/11hajBJ53P4zlpq18DtRS9gTcODaHTrxTmH0fr1YEq0DB2zJWa1lVSPqqs\nMlPMZLNwpWBNp8XCbVxddDotzCUNVjAsbj91ZF8cBz1QhKtSv16v4wtWkyppMckC2FO31ol60BA6\n/UphaIKZfo3smBtmMKmQSB28r+4AAlcDuZkZwJ9Bh8jlqv45S6sZN5MRg0qDwSP61W1MjYf+3Sq8\n5RjUQ0vpu9bsLNuZ/p7Nm9Fx+XJR+tS6dTj+2mtV6FWDUOdNMbKdop2eVRQERt13TNWHDs9kMJVO\nh30CxXibYAAwPS0GmxpgWxCOgzM9DV4yVxlo3y0GxYpKWSxmtWibKJvmmOmX4Xun4/Jl4yKYW50u\nNTw8my3ebo/quMZlan0QumpnehoYGzPqt2tBZnLS6Mxtcm4uspwzPS3s/xVYTw/gOHCPHzeaNHL/\nOAyAmd4uIHT4oU1oSRcwDYM3A3DvypXoXLMGAJCbncXCjRtYgFDZGLa8matGvIXxalFkEuu/nRV2\nJEvoDaHhqZ8/fX/kWuzNWYuhOgnNkKv5I7CoOarl9ZLFTFPRBxyZP/vCC4nfEno6O8H+8z+L0t32\n9pDrA50sgoF/CkJwZrJZHB8dDdwnaOsBMl8Ujm1xNKFVi573A/PzsSx49qTTcGGxw29vB3p6gvxV\ntqjR3/J4NothkzttZaMdQG8ItaI51Dv12py1CKqTJLssy0Gvc2Z2Fm57e8EMUbVEiTX4KHbkSSna\nMATfTPLGjZCqxo3yKlkmUtjnEOi5uX+Un4FqNaOrNVTU68KuXJsQL+MZ1PqZ58UqU3gGU32GAama\nExrdn0+93GkTguZQ75QxE56xbNG3pcfuC8Izm9ivrsqML8q/SrXIzcwUNhw5UITBjRvCZFA384vp\nfreUq4Aoerq74830NK+SSZj0B7UCvr4+Nz0d367fceBmMpiankYHwgOFnt/V66vEpBWIrZvnCM/2\nTTuKIwckcpBGxKB+Qr+MmfCV+Xlj/ivz87HrsOm/h339tUqpV1f1B19VMz7LYNKdzxv162pe273F\nglvSHfifZ4Qtu427b9wI6ewLevYSn7e+8J+bncUqlFbtyFkukwk2//YxF7qdyF6G88l1BPl/Z1oe\nfYavewtV8+tlC31ULgumtL71zmJMXoj6sqRMNjtSKeOr8QNJnL+VCO7SCJQzmDCgSD3DqtCX0GJm\njXdi2gZkHqMsU86zABauXsWq6eniHbkIz+SZdo8jgZDWVCD3vvAC3OPHC9czs7NoSaWwNaY6B4Dw\nYJpkL0CFs/sezduprd24Fk6FN2Za0K0qzaHT94+jCcq033yzUe/cfvPNVehRc1My0EaVPWrG7odS\nv74+IJmKWXeorOfFFnwcgaB3lHpUt8dx2l1340bsWTKDZYCp9aSkzEEip7xlAqLvudZWXG5pAbtx\no5Be9JZEVIXm0On7x0wVzA2TMJVOF5n38WpUXIafkjiBOapFqUAbZQt27bknX3gBe2/cwAqYNxHp\n/VDL6usDxn5CCOPJuTmsVe5x5VzNb0Ov01QmiQAz5dHbAAJBL496ObdEnZFhFquB4Xt8+oUX0KkI\ndkD0PbNuHe4mh29LjroLfXeR1S0mwcLMWRNRsBLRF/04t7/ulhuYA2G1i83yRtcPV4uohc27AfCx\nMaNAK1WnrpO2Mezruk1545SfRSBcZwBMAvgNgNVKnosQXi9dRMNQetLAtHP1T2f66lW4HR2YnJ3F\nWgDfK8owDXdkJDA6sLyFJFXLRZkZux0d5v0xS0BVShRTPzt9SZIt8Om02aywCtvoXQSv2tw/Tq5c\nCXbXXUGmODPvaltQxFC72CxvpJ25u3Fj8eYlGGzpZXulZpGWZ9zjb+l3lDZCxWDXl/OREfvGuxgD\nXBK2an0bQLF3T6bkUc+zEMFMHCWvTS2VgVgwHli5Eml/M9asxbfPOIAFAO95Hj46PV14G2CGvLnL\nl0P/b1Pb1j0GFmptZkw0DnVfyA35+yiB84lPGK0/nCRqEIMQ1S0lCvf9H2o90QVwYZYt3yDgW2Bs\n2FD8WcrPxTIje06zpQcSzs616/Tly0A+j9MATlny62Vk5KxZg1oHELPw+KF5inEhBHQWQqhuLaMO\nZqhPwlHsB1/NywC4a9YUBrSPQwwyOpMAboX9c1PpVnYQD5va9b2SEs3Dsl3IrYYe0+bzxajv1sw4\nWcWtVwHOiwR1D/woSnWwmGDqhS+M+hOUly4UXMv9tQib0zJLPkkWQqjKdz+5CNwNX+VheEuQLhA6\nEH7bk23pbTItTQ/8Is9nAbgbNog++O3+kaVOINnnFoUtMBFz3cD/PoRFTlpxEW0l6g27jLUsIjlN\nsZBLVAbTEyJs4JMsXlfDyZrt/Uj3ZSMXCBlKRMyKuKfTDaFWMb0/2jbxydm/qQ1buxzmWbaa/zqA\nNbkc0jMzhTeEDMTgYqKltdXoB4gj2WdghXMcN6xnceUcKF4sdrq7Q87YJE53N5llLkHqJvQHWlsB\nAG//+td1n0Hrjr5sLn4rEYgmZ2k8m8WeDRuEENCsd1iVrHeY66JDCQoi2QOLeWSEdU2idv1jFkDa\n30Wqm3+qC4RuRHt6GgPw//l16+RgV7f0axYo5TCzciXcNWtCunkbPYBxLaUUpjxuzLJFdWmbt0wM\n6wn6BIJm801F3YT+c75QdVtaYpeplbM03dEXs+SLMjcM9UW5lOaYk+++C0dTywwDYL29tQ3EzbnR\ncicNc/9N+uYi/GcstajK/Pr0RUX5f1Q3/Ewh2SYsDvtM21bPdeW+E6MtQAwsDGKhdWtvL7b537WB\nVauE+4sqwQAsLCwUxfSF5bowSSkx+PCRkYKjNL2e2H2j2XxTsbTUO0vAt0jRwqtlhyOrZpv+0fQm\nInG1MjZjuzmE+xYKJJ7AhQL361kA8PMzZwDFHFD+H9V2elCd/RIc9mfrgd1cUr+WyLWBG6mU2Cnr\numB9fViwCHyO6IEl29qK+ZYW44DxPcPOXdkvxz934S/U7tmDzOnTxrWQqcnJ4ILMKgmNulvvJCGT\nzZpN42qwwYshcG6l2rmr7UfZNqsufauKxSMmA6wLzwzFws72ibVDE4D+q34mm0WPptedArAvlcId\nvrDi8P3wa3W4uVyhX6p/IPWTYbAL3iyAvQC2+dfckg+ojY95DmCF54H19eH0mTM4ZVCXqTDYn6c7\nn8fPfRURIPwBYWEBCxaf/joOIFwqS5cJ/huq2hYHCt/LyQqdERKNQVNY7wDCcifJl9KkmwaKN4lU\nHChcHg0xUZlyH1FWJYbZbDVgw8NFoQZzvqlkURcsAw6DfdORaUGVZzJYMTtrVge1tBQWH5mhrNom\ngNBgZeuDTjeKrWkqHeb7EV5wngTwAICP+tcchgFsbKywFgUUC1oH8VxG/N7q1UV7Epg5axHZ1lZ0\n21w+yxPle9mP4HPOQQye8wDSra3ovvvuqux7IGpPU1jvMP84sLAQu8ykRfAUDRwNpAbKWNqV6SV9\n4hjQHbLZvE2WuwuX6QnT0yFhlxSbgmFGaYsDeB1Cd65b/0wBWOHPjuUzldqOxxAMDNKSx/a/mAJw\nt1JOHk15F5QBTm8PCBaRMxDCdhZCxdXm33e0/BJuaEumq3nTbW2JdOz6RjTZNm9rC9RVlgkS0ZzU\nXac/v3Jl7LzrFhbMs80EA4cpzu7k7Cz2ptPY9od/GM4sv/iaq+PIxUtNrWOzJBlobQVvazO7dI73\nKCXhSlsXFxbgel5B0E4imAE6St6Lhvaj+hNXLWGiTbv+PYTVM0w9+jujTW9cLysqJrWci/ACsTpQ\nXEdgqsmVPLJeDjOlvmnDhjSGcH9zMzOxF5MdaJ9/FXafcwQhKGmRdvlR95n+fJKIPC0tRjUGElgA\nmeLsMgDZBAMHs5wDiK3WSa9bJ+yfa+isyoHfP99DoqlPrn+chcHPi6WMirSIke4JkjALoa54zvA/\ndZW2OYQjt861a0OmtdIx3fuyWaOjNiBYvI2qX8WUplLOD4Zr9ZqeV5rQcoQH4Sn4C7f+/oakO23V\nttU+tFEkrGVL3YX+wLvvxi6TXrPGOMtOa+4SrH7ZIxZ8n8vnrTNu1ezNdD8Kx5aewP1EHJh2nYV5\ng1JRPxCt4y9gGXBVixhjP/y3DBgsU9YA5kEcxTPcvTduwNH+BxwAHAeO42DgO98JuSfgEG8zjrF2\nO+MobbbKEP/z1cm2tkZ+jg/IBEXQ29R96v4S294SB8WqKgcAqvz9I5YOdVfvIMEMO67DtZxlIMm9\n+27RABELi9lbtrUV7G5fG2yYgTFE73yVefQySTGVuVe7x2EWzBxhtwU2cmvXig1jCtnTpwuCVjoX\n07nc0oLWlSvBcrmi3ahTAFZp+U39BYArMDynb1nERkdF8BLlu+GUeB4baxDsnC1qD8FC9wOpFLIt\nLVjI59ECsTiqq6vUvsi63La2SPXgSgDo7Y0Vx0DdX8IQflMgCBtVF/rZbBZ///d/jytXrmDnzp04\nePBgZP5cAp1w1HZwlcsw/2Avw+y8i1vai3ozAPxFNYsZFdOOOsZQeJGt2SoK1htOnzmDdb799we0\noBdRbVxGsJCpwyHMMmfefVcMVH5YvUw2i5Q/YDPYN3u5N24Ix3W5HDpQrPOO28d+S/r4L36BgRUr\nsJDPhwQes+SPA4d9LUamfdTzQi4TXFRZ2HKeyG8+U47McF8GkQdoUFjuVF3od3d345lnnsHCwgIG\nBgZKCn19phdJzO3gnWvWGB1ruUln+YYZvtp+bmYm8HQ5OYmezs4in/A/hxAi11IptPprDwv5vNHb\n4zgA1tsrnjNGgBU2PGy05Mkgnn5d9rMDxTNrply3eh7u9oOvq8HW2xEMmFHva1PpNLKzs+jWNiQx\n2AdcHZsVylZN+ELLZ6vfZijsIHrAyCIsYDMQVkhrIf5//QA8APJJo2LxqnAIvX4KCO1pkLEYRAIP\nl/G/a7a+qnSj+P9LLE9iCf0DBw7g5MmTWL9+PV588cVC+sjICB599FHk83k88sgjeOyxxwAAP/zh\nD/GP//iP+LM/+7OSdSeIHorM6dNFljeAtgMR0aadbR/6UJH+WvejIstKoZ6bmTH/YJR1gD0tLcB/\n/mfRLKo1lcKw5wmddgl79pbW1sKbg3Unb4wAK1GCRvVCmUOxH3kgWOTMRtzXZ+2upb0ZCJfLNisf\nB+HPYgzmz8a2xY0b8ss0+dcPoAVhHbzqkkH2Q6cfYjDbhkCHr+5jYAg+a1OfXaU/LoQ6a/rqVbjt\n7Zi9ejUUO1e2H6rHt+xyfVWgKd6AqV0TurM7MslcvsQS+oODgzh8+DAeeuihQlo+n8ehQ4dw6tQp\nbNq0CXfeeSd2796NLVu24L777sN9992H+++/H3/6p39qrJP5x3cSdHbmrbfQY9i+flqz048y7Xzu\ntdeK0vds3Biy/uDwBZq2uOvCbnO+ytamzeLIwEI+D9bXh8zkpHFwqwaq4FNnrZIoL5CSSQDrDGWZ\nIe0ygL58HlmEBbRsx9HyR/mmN2HzzsmUo5xBm/Ldi2KVn8ynvl24SK4yyiH8ppEGAM+zRv6y1Z10\nv4Vq7VOoow5ut4nGJJbQ37FjB7j2ann27Fls3ry5YEI2MDCAEydO4NKlS/jBD36AXC6He+65x1rn\nqH/kEFuM+/QNSgbi2unPWwTtvMW0Uw+heC/MP0CeSmFjTCEuyy9YdOumvACAsTGkLRYe1vLK3gCn\nRF6OYOCyzeKHTf1SWGspa8q/3pKutqMSy+Gbz2nYBwmVHth38K6z9IMhsOS5CPFGyhAMXFKFpkfR\nUtFVKvDrk+X3AZhPpbDa8p2SZWXfXf84hcDap2gGDxLwzUQ13S9Iytbpnz9/Hl1dXYXrzs5OnDlz\nBr29vejt7S1ZftQ/9gOxBH4iWlvNHhBbW40uGnTdqGkWCwD9BrPDKEx1yDQX4YAdw2qmGANFCMve\nAHnNlbYcVLbICUT4y0fxguZ1S95rFbQPiP6vg91s8jcQvp1uAdAKoaZhWnkgei1CzvQZghn7sFJe\n/iVBX/BmypqEqS5bGvePScMiEksL6X5BUlc3DKlUqqKGewBsRDxfJYmxzZLzeavN/R4EPv7nLeU9\noBBpqFzkJqZZ/5wD6IoqUAI1oDiHmMXL2SeTeZBMOMm8vIJ+SeLvty4N84/cP0a9bbRD+NFhhjwu\nigcAPV8WwNtKW9cgZv4uxHe2lAosKQzF6xPyWu8bQBY4y42GcLi2adMmTExMFK4nJibQ2Rk/mqlc\nmEvi02XBYtOvp6fm5435UvPzVhv+DgDDJcL9pRB2XwAEP753UinjBiQd/ZV/AGYTUw4xU93oX8uB\n4gKAsVQKG154Ae7GjUAuZ1zg64dwrrUCpV0HqFxD2HmYrEMu/KYhFmdNOCgWprapwfss6ba+qnUz\nS5646HVxQ505iNCGprb0NOlnBwisgjwAd0T0QdYh2+YoHYlLh1Ms3GVDQzhc2759O86dOwfOOW69\n9VYcO3YMR48ejV2eQQi1XAKVyXWLYL2uvXVMeJ5x4W/C89BlGRDUBUGbSmINwq4K+gHAV2Xlf/3r\nUEAYW2Ql06aftf5RTx9A8WyWAch6HtI3bgAXLxZZk0i2IrC+0euN4g4UC269fJL6bJgWYBnCbh3U\nDU9OjDoZhPCUTs7ieuIcNqQNIFlQF/3atEgOFFv/2OqJgy0WLtF8LPpMf//+/RgbG8OVK1fQ1dWF\nJ554AoODgzhy5Ah27tyJfD6PgwcPYsuWLbEbZv7x7xN09t1Vq4rcKMt0lQ/B/EPuB6yzcXUGzow5\nigXPBQC1rOaAAAAgAElEQVRrfHO6/zIzU6g76vVfXXhkCHTSHMVCQNdDy3xxF1FL4Srn1yBm306C\nsnLmPwsY9x1kAVzS8uYgnll/43H9a9WtQxbC8+YqhGfnHHYrqmHlnFnycJgHM5VKnBQwCHXhPhS/\n0eQAvIQgNgBBxGHRZ/q2GfyuXbuwa9eushpmEDP91QnWBlbl80ahNFmhnl1Hvq6rgiqN4vWHjYi2\nnWbww+zBvj2eRaR5EBYqqtIsbhwkjvCgca9WzyTEQugcxEJnCwweHS0wiJmwFNzyT6cbxQOUa8lr\najsNYf2zwnAvTj9tOAirVdS6ktQry8v6APEdYRBrGSnYJx9x30LkTlr5HZRt2PZi1CqkKFFfGkKn\nXynMP/6/CQR2x40bxh/lXs1S5xLMP95LEAvQTJvtc4QFss0zI9PSoxRTMl+/lmaq18YdED/6YUO9\npXAQFs66aaLelyT9AsxvIbVA92cj28wg7DNoxpBX5ucoHnAdVDaQMOVoKufCPkC3IL5HUpkvBwB+\n4JMCJn2+JZaEW0bcBqJxaAidfqVIe+y3E5SxOVDW020ml3sBTK1YUWTOadvgo8MhZsjyZ2dbjFRR\nzf4YKo/4ZIMjEHBAIOSkiaBb5fak2SZD8UJmC4A8zAuZSYP3vam0AwQCXN8J249iIcuUI7PcK4XM\n9zLEG9EcgN9CvMEB9oHfQfL/NYNYnM0pLiscpQ/36rEnOC+4aSglvJ3p6aLBgJmzEk1O3YT+BYgv\ndBzBKbH9wBJZzxts+OPOuhz/yLRjHGTeOJuP5KagtxFY75TCQSDcryFwt7AAodpZl7CfktMQ7n6l\nCmgNhHCVFkcZAMeVc7mWIYU704426xwOs9OyDbA7c1PTTGsKSZGqFNkHtX6G8ADzPUOeJEg1ECC8\ntabbxHvKVDoNrFuHzMwMOvz1K9fPt252FojYj0E0L02h3hn1j/tqULdhW1Yh/bofMUinD0Gg9nFj\njvhwhGem8pxBCEwXYT2wjmlTkERdbwCCNQf1jUmdYct64sBRLHTvRvCGc8pQRh3EotRiksswC3cH\n4WeVA8gC4q9jAOE3Glku7qKs3N+QbW1Fdz4fii2rvg3qb5YynSP8XLoprnymHMIL/bl8Hs/JtaHp\n6YIL60YJ+UnUn6ZQ7zD/aI53ZCZKmKustuRbDWCj5xl/OKrrXpudvglZRv1Bm7bfq9eTCOzgXSV9\nBqWtOuKuN+htloKhuA71/AOWe3FWZDgCQb8CdguhHAK//KpDNzdGG2qfTNey/XEUO1/Ty7htbSEv\nrXq9c0qf1EVV+X+V+fcg2IyXg7BCGja06RrSCKJW1E3oA2Jm/esE+dVXYj29Gsi6L0DMYKWduM0k\nUe3LAMQPmiFa2HYgIm6uVjeHeLYHEAxkSTZaca2ePj/9FgiVmBzc5L4Ek5mprZ+uf/6O0oYNB4F+\n+2alb8OW9hiC8IFRcCWPzQunvL/Xz+8hcITmQgjqPMRnkgfQn0rBe+897GlpwSXPw4eUQCk2X0WS\ndRC6/3shLKU6ALwK4PcRrUJ0Iu4lxuB+nGcyiZ22EY1FU6h3mH/8HzWo+wLMAuMC7Dpy1aOiemQI\nolABxYu+8tysNAqQevo87GsQqu5aHnXTvN+UaEfF1lemZ1TSXIR9AplQ32Rc/+hAuEU2odazGtED\no9zUZPosdBylLpPeXx0I2mAelIvakG6w4TuFi/CLo/ZBxUXwGea0ftYa04KuNOMsap928y4ZmkK9\nI0myCGvz2KirY+Ja+ZTiJYgfcJul3b1K+nsl6rqOQIBzmAWPo5VxDO3aVE+myEhJApVzv5zapt52\nFAz2Reokdco+6280spzc96BjmoXvhfi81kDM6MuBIXozmAkHdjVZpfD2duxJpwuLvACAbDbSiofM\nMgmVus70+2D3zWLCNpvW02+B+Ye2L6IObkjbhmg9t1puxr/3kqX+UuH3oJRlEMJtjaGMbeDSdxRz\nFM/0bXFsgbAKRpbhlrx6Ock8Att5abo5C+D/0sow2M0ZZyGeXZ0MMOXcRbQQVe+1QVgyOShPBSjV\nQMMl2tThCfNnW1vhtrWFoqLBUofjOMW2+H40syRtEkuLZaveuQDzotcF7TpqcJgybM4CkgkFZkiT\ns/57tfsc8XW26j/je4b7LoL1hgUEtvKm+oct5aNQd9iWQxvCqgyG4ti28nyPpY41EH13y+yD2oaK\n/P9wre5pCG+cJrpRnqdRB+EBWO8bR/h/1p3Pi4hoMd0kF22yIpqeplDvMP+YxHqnC2Zhtle7tsXd\nXQXE8oSpY5u9R8H9o5wt2tQLHMFnUSpWlgPz7lyG6sMQHgA44g1eXKvjOoSAV/8n70DM5PdCvE2p\nZSZRm+dpsdSbxGRYLc9hX+DnCH9uckLBlKOpLwSxGNRdp59Ez25yQSzTVWwqo1REe9OwC7g4H5Js\n07Yb2I1RxzyiVR+VMgWx2LoXxS4LHEN+pp2PQwjxDhQ/D9OupQmr1UcMwiEI9fIMYVWR3BT2tqHt\nOPp2m9XTOxAb0KTuX8VR6jb1L6Ocv4TgTcfR8g+0toIp7hN4JiNUMgRRB+om9Ll/TLJAFqWrr7SO\nAdhnYupgwxO0FQdHqVuqQtwqtyGRNv4ukvtuZxD9s+0T6AfwLoRJpgvxxUrahp5vDwIViy2IO0PY\nXfXLlvpsg32ff+QwD3xRi+HqZ9EPu8O6UJhOxxHtGDZeEcRiUDeh70D84P57gjIXYf5R6Sqidyz5\n3oEQ+iZUB2Jy16t0kaxuTDLVCwQzfduMUu7EnYTwnx/lNsAxpDEULxCq0Z3kwqnqgz4KtR5ZXpqt\nckuZ6xH3VJNXoDoDl+6G2sQ4xCI6969tG/hKMWxIYwhvoJOf8QKAd1euRNvatWDd3eCZDFoiZu7b\nbtwoLLyyMvsXwmCLL9OJ5qQpFnIB4YohyQLqWph/NPrGnCiHayaN/h4IV7jcv16JIEpV1IallxG4\nPHgTgQ7bhFwkdRFY+uhIaxvpQVKWk/0atvSFKflXIJhxcpgHEHlf7cN1iEFkAfYByaaqMZHkDa4S\n9MHmnlQK7vvfDwCYmZ1Fy40baIHwsJoUjvAzc/juItrbsbWnp2AiyVwXp7/7XXBTXGYTFQrtUiaY\n5F65+WiqhVybFYeJd2EWlnoAxCjdv2mmr275V3FRbEmi1vt/I5jhd0EIzVIxYXMI7P71PuobkriS\nj1vqG0XwGaYhNn+NQ3wmtv0Fav0OwgODqqfWHaipnjWj6jWd63UyiAFf7jvYAyFg5X3Hb1PWYdtx\nq/NfPK+gL1+BQCWk9jmLYjcMel/ht5+B+B//L/XG9DQwNgY3K1Zf2PAwGOexVDY8kxGfdy0FsMW9\nMqtNa8QSo+4LuUkCZ2+A+Ytr0+OaKOf1X77iy9d7qeftQrFpJUN4lq6TRrTdt0yXwk+aZ9pwYH8D\nGNeuJTJ9q+Geq5zbBkNmSCt1X6bZ/M8w7U/HtvHrJS2/+qajpqvnruWerV2bfyFTvGVZB4flLct3\ncWxqiyAWg7oL/Vpg2wi1B2b/PaVUEcf94z4IATzsX7uGvAzBeoDpnkznCASz1BNfgBhEGMQzvA7g\nVuW+2m/9GUxcQmCVErUQaoJF3OPKPdOsOaqMbiGj3ndh/1/MWfozb0h3LXlV1Pu2Xb5AiWdbUIZj\nXTWTyYSctulk5C5aHVLBEDWm7kI/udV8aWxvDyth1vfbZpG6ANKDgsxYys0g2qkaQ2C/vwbBjFAK\nHu4fb0X58XA5wqob6dOmVNnJGHWr9cvPaNyvu1xzUwfR/Wu33DOlOShtZcURhCC8gOQuOgBgQRH6\nuqBmrgvGudXZWUcuRyoYoi7UVaffh2RuGOK6Vq42zD9Kgeb6R1vfW2AfEGwzVlc5n0R00BMpwKNw\nDHlKlYHfrhOjbqljlwudcvC2bbaT/YnThyRwpU69bttbg+zLAIBcayuwciUWcmUsPUds9JODAOvr\nMwp3gkhCU1nvJLGqsHVWT48y7byOYrUMN6Sp6bMwu0Xos7RzCXbBafP1r28C6oZ91hzHA2W5xHHd\nzGGPZtWP8GfpWPIlxRbSUtbvaumTEO6NTelym1Q3ANx9N54bHYXb0RG5YUq2zRH8b99JJZmyLCJk\n0tl0NIX1jiTJz+ZNmFUxb2rXUV42VyHajBHKcRbFwlgeAeD9KD1r17HNJ9UF3g6IXaJrLXlnIdYX\nWgG8hcDCZDUCnzxORB9UOILg6Rx275ZAEBBEFZw6ugml7ntH7oGQi+Lw21RNXWV+jsARXKnANGoa\nh/jsHMO9AZgdok2l0yEdvLpWoeZVv3+tK2L8fGwCOJut6q5cMtMk4lJX9Q6QbDftB2HWce/Wrkvt\n3DXd09Pi5LHt1nRg1ynb/AKp/E+IRWfb4GXz/W86qpjUQg7CpqESPV8WQVSwOMjy8hk4gj0Kur2/\n2rbeBwb7wixH2G89Q3RwFomjXEuzy55PfAKZkRGk330XWFhALp8v/A9Vc9IUgv9tB1AyMLk1va8P\nuJjE81QJyEyTiEndhf47CcpEOlJTiNL9JzERlTD/yBH+EUWppmxRvko9r3QpIEPsuX66LZQiN/TR\nthg7B7N7ZVV1wlHapHLWcN8Es6QxQx6O4Fnl29AFiA1185b6HaV83IhiF/365P6AhYsX4XZ0YCqd\nFiaqhg1WTPsrcOOGMVB5LEgFQ9SJugv9JDN924/flm7C5KJBLkq6/lHdFKQfVUx9ZxACzGa9s9eQ\nzhF4l7Tpyl1LumNIN7UBiEAvtt3KAwgc17mIntXbzC4rwYH5bYsjnqpqAcGO4ihugWF9ZnoabHoa\nvL09RkvVgVQuRL2ou04/CVGmmCpRC74m1Q8zpMGSFofhiLIrEMzOZxF4rTxVYZsqqlqIWc51dP21\nG5Ffqnr09ZW3EEQ349FdLHp74ogeYE2oZeYh1GKlytbCRNgE6diJRqXqQv/EiRM4efIkrl69ioMH\nD+KP//iPjfnk6/UbCeq2/WD19CjXyqY6VC+NenoUF6CoCRAIW1NdEj3C1csoVq9UitT56/1QXSBU\nQifs4QltA44kC6Ejt+323QPxmcp65iA+k3shFo85hKprBcSbiVQ1rVDqNcH9/pQj9DkAV74FxF18\nJR070aBUXejff//9uP/++zE1NYW/+Iu/sAp9uRHptwnqjutPfwJmC5oJ2H3vmOrdi2g3CBthNuVk\n2hGw7/r8vqWsCYawukM9j0OUWwWunGcg3p5MfeGwf2laEHwesh59/aAbwjLJJpzTMA8oDyD8uZr6\nJttkKP5sHASLwklxAKDHf5Iq29xX7Y2A1giImMQS+gcOHMDJkyexfv16vPjii4X0kZERPProo8jn\n83jkkUfw2GOPFe49+eSTOHTokLVO5h+TRKWK8p6pYouwtQ/JvD+2QAip8VIZDTBD2wxhgcW1fHoZ\nFTmjdVBa+Ml74wirYGw+ZLhSt0QfjICwaacJfRYdtT7hGNJlGRO2/Q0qHOHBJAsxIVgH8eymvqh5\nL169iv5UqrDp6kYqhdYVKzC/di16IuLWViRYq/RGQCojIi6xhP7g4CAOHz6Mhx56qJCWz+dx6NAh\nnDp1Cps2bcKdd96J3bt3o7u7G1/84hexa9cu9PTo87wA5h9LhQhUiVLbqERZ77QmaA8Q/dyDsKB2\n/HObqoAb0t6n1WlqR+IibJ8OCEuh7/l9cf20SQTqJdNbhMkZnI1hrX0TcsHcZr2TKtFGpci6ueW+\nY2jfRfjNaHLlSuwFsEK6UGhpARYW0J3P4zlthy3zPOCuu2LFrY1C7RPPZAKfOzQLJ+pALKG/Y8cO\ncO0V9OzZs9i8eTMc/4s7MDCAEydO4NSpU/jZz36Gq1ev4rXXXsNnP/tZY53MP+qz9Cii1DYqUQu5\nJhWRzae/3DAkhy65Q1aWN/UFCEwtHcv9KJhSrxwQxxGooOQmLr0MM6RVmxWIVpHoVj1xzSjjUOku\nZOb/OQYhntRVQjkqGSZPfLfMoTSCWETK1umfP38eXV1dhevOzk6cOXMGX//613H48OGS5Zl/fAXC\nr0SfyeOgRpTaRiXqjcCkIuovzgrAvi9Alp+03E+jeNbJEKhc4sCUPxfBrN215JfeNGcRXlQu1QYA\nXNPqlQFhdJ289PT5upZf+tp/3b92/GOpPjAEjs8A8VaTJKhOHGRULdkW5Ey7EiuaOCoZRcdOMXGJ\ncqmmzx1J2UI/VQW/I30QgiWOwAfiO1yLsvIx9fpNhIWYtA65gEBYDBvK2ezhbW6VAbtHT67k0ctf\ns5RR6+73y61C8E9lWh6gWN0xbKmPIdpsU013EZip6ukmdGd0lxG8yeh1SBWX+j96HcULtfq1ZIWS\nPgwUZtrMkLeaqAMK6+ujmLhEWUifOw3hcG3Tpk2YmAgUKxMTE+js7IxdnvnHIwnajFLbqFyCWVhd\ngph96vc+iLDwYzALTFO7pnuuJX8UDsIDA0NgiaLqym0L0WtgN9XkSprajp5PJ+7O21LeLHWf+9v8\ndLkwLDeC6f15GcD/AfAbCNfKjtIv3QEbM5SHUq/pXkNAVjdEDBrC4dr27dtx7tw5cM5x66234tix\nYzh69Gjs8gxipp9kN60t/qyeHhVhyzPccxP0QUW+CehMIhB44wjrut9GMNvmCM9O9bqkJYrqa8eN\n2TdTv2S7cYkKeKLWb/sfyjyupT8Lfl02U9x+iJCPrQi7oNgK8x4KBiDb2oruvM1OyYLjCB88mnvl\nqXS6YLVTK8jqhojDos/09+/fj7GxMVy5cgVdXV144oknMDg4iCNHjmDnzp3I5/M4ePAgtmzZErth\nDjG7TmKnfxVm4XBVuza5WpDpNvvwcrCZkKpxf98C8Anl2kHYwZlaXr+WuAhm/isRrEFch1DnXEK8\nRdMZiFm3VDGVu9DqaP28F4rOPAHXYY9YBYgB2mb/lYaYNHS0tKBnx45Cei6bBUvoyMwkeAuLtZyH\nI1yVMwjQbJ6okEWf6dtm8Lt27cKuXbvKatiB+NH+jwRlosIgxs1ncoRW7hKbbVVD/VD1tw6mHLmh\nLyZk/6ICuDsoLXTnIWbMcoHWll/2K+6ehlJKPZuaSNXjm1A/X1OeAQBIpcBGRwtCuiOXq44lU9Ri\nbUIhTrN5olIaQqdfKcw/JhH6cXX6cX30lEpXYSgOKG5bWE4hcHmQQvEi8f80lHFhF45ydyyz3J/R\njtyST6IvmnKEHaw5yj0X4ZCLNidsan0SaeKaJFCOSqw3kRbfRqiERQ1vbwdT941UMMsmIU4sNg2h\n068HcTdnRVnvmByu7YGYNcrdoDP+9dsA/giBEL0E4VCsFDcg3ihU1cRLKN4nIAcGB0LoegiCqpt8\nztto8+9LNYtt8FjltylNUdU6GcJCW29PLaffn7H07zKEIzmb+mfST/85ggAtQBBgZdI/2lRACwBw\n4wbcjg7MzJqfWgp7hxydEQSABoiRWwuiwiWaolH1WPK7MM9gJbaNZStg1/erddnqNaWb8ul0+vnu\nseR/F2KAMNXHS7Q1AHuQdtUiSi27FuIzlPF+HUvdHHbna1FunAsL3NPTVlNYp6en4h21BFFvmkq9\n898TlJmLmb4RZuHycoK24hDlAE4uGOv31ehV5Xr3LIXNesmFeKMwCW+3zLbGEbZi4igduSouqyDe\nemyxglVsPnsIohloCvWO6x8vJChzCWbhpOuMo6x3TF42k6DWux7mdvYhrO835WGwLzi7Ee2r93II\nYtYyCLWIi2hduBSOertquMGocjpbYQ7fqMNhXxi2pbcAOI7wm4/NW2nVIYsbokmpm9B3kNx6pxP2\nnbEqtoDlu5HMn/oUzAusTDuaiBMLNwnqbNemZrkOoQ6xCf2cdo8Z8tjUJLVEuneQyPM5BIOcY7iv\nolplZVtb0X23H7q9TCFN+n+ikWgq9U4SoR/Xn77NDe9qmJ22yVmuXrf0ocMt7dq4A/FUEjZkuEa1\nzVI7S2V7w0o+ruXJIXi7+LmlniSb5aC0EdU3IHpTmO1NYTXMzuVMqOsyblubEPYWO3sS6MRSoynU\nOxKbnt5EXH/6UT56bE7bpOsDINgpq7sGMM2C5X1d7SAHKG7pyzjsO4zXIDzQyP5EkUZxyEP9qJK1\npN+DwExzHELFMu//uX4etS8OAtNOhtLmolL1tRKBu+nQLF1px6b2kf2TcBg+H4pcRRBG6ib0+yB+\nqEm8Kk7CrO+e1K5tm/DzEfd0173MkkdFzeNark31ANF6aQdCkEmTzmsIwgyugjCB7FTyAsVOzBjs\nApihWKUiWQfh3AwQwvmjEML/OKIHEZnOtPtc6aN6z1QHQ3E4yThtuFr+ybm5stxaE0Sj0hTqnT7/\nyBOUieta2bboWMrCg/lHbkmXC6fyWnU9PItoYSiR9uiOf+3CPIuXaaY6mZbOUPyPLLXr1hbm8Q5D\nuwzB/gHbZ3hNyavXqafZ4JY6SuFoZdzVceJsFUPBzIlGpSnUO8w/3p+gTJRVTtx8NusdB6VnmDr9\nEXllfXIXqzyaVDBx27ORhbCCUgVyKRcHUXBDuSmEZ/u1wFHqVzfLlYIjeLNxIByllQWphIhlQN2F\nfhK/N6bdtEDxTD9uvnJQ640KlzgFIXjlwqhUK0mVTSni5gOEcPwAwsI6SSxgHQfFn59bZl3jWtnX\nAdxqqF9HDo57IT67j8LuBsJBWO3Dum3OIgiCqLvQr4Yg1okKq2gK1q0Lk4xf/hrEYqM03ZyC0J1L\nfXpUvF2bc7QBhNU5si0gcD+w4LdjjzAsYP7xGsQ/cj2CASdtyAcUu3pW70vVkwz3qLqDkBZOpxH2\n0rkAscvX5tdStcyRbTFDPltaP8TAKcvJPprycz89k82i5xOfMA8sZGdPLHPqKvT7alR3Kd0/M/RF\nTbO5ZWCWo45j71rR5ihTzFugMr/5DHb9eD/M/bPVI/9sC7/y/g6ErW7k4KEHO+Ewwy39akEw8DAE\nayCmvsh0N5cjHTzRVDTFQi7zj0ns9G1BS/RdvVGRs66hOCaunp8n6JMJWZ5pR/jtq9eupQ6pnmFa\neta/1w/7JqzTEH5v9LJAYM7KlH7GUQWZ2pJ1MIjoY452jynnMJyrOCj+HzC/3dMABlpbkW5rw+tz\nc8KKa/VqXLt6FXd4BiVbuTp9gmhQmmIhV2KzqTdh8ynzknZtc4/wMsRirm4uqS8Eq7NTUz0q8v7L\nCCxfnIjyuo94m8C1qTCY/zcAu9BfB7t/HVPkMFM7OrYvyrAhLU59OhxhSyZHqWcAQLqtDU5Pj7jv\nW9Owvj7jwqtTrk6fXC8Qy4C6C/0k0Ztsu0X19LhWPhJp0aOXYdqRG/LJc5slTylsc9I4s+9y/nm6\nvx9mzhaCo7KF4Thcg1gjYYZ7aQDDfkBzWPJUA1IJEcuBuqt39FCHUURtuoqLzS+PKa3UfTXNNnip\nC7WA2Pg0ADFQtcBuATTv58tDDFbSj/9bCHbzxgn+ouMo52MQz6C/KUm4f5xC+b6ETO4oGAIdvrx/\nA8EOXR2bm4wkkA0+QQjqLvSTWO9E+dRRiTLZTOJwjSOwGDHxktKO7S1kJYp3pKp9k9fcv5Z5tyn5\nXARqFLW8es6Vsra+OFrb0vXEHpg/r0klXVrhcAgrowWU/vIwFK9ZMO0ICFccaQQbvHSsypoY6hgp\n7HkmI94WLP0hiOVC3dU7Nv8zJuKqd6KwOW2L8mFvU8GsUOqyeae8A9GzVGY5xoFp5/LaTVAHEG2t\npP6p6S/DLoxfRjBY6tY7QHiAUtceTH2IItYM3d9wlbRugmhW6jrT70MytYEpqLlMV0kaLhEIO1zT\n2zPlB4TqRQr7NxBWW1STpPr0SRTH5d1myJeH6PM4zM/IEbbQYcp1CsHiKxDY9Gch1FHcTz9uqFcG\nspf1SmRQGY7qf4YEsZRpKpPNIwnK2IKO7DGkJUX3qAmUnjHLWTwQLOQymPsovUeWsyCaQ7AOYPuH\nqd4p2yAE/zoUrzWo/XsfAksgwNxvBiGATyMQyMOG/Nw/piEGcwb7YLkqoi31T2WgtVW4TAaAbFa4\nSyZ9PLFMaAqTTeYfbTs5KyFKDRTHWoj5RzWalAuxCHsLApNPPdyhzKumSaSTtXJCIaq7e011y/qH\nlet7YVa/mNqP85ncDLMwNvVHTZM+dBwlLSrurY3ufB5M6uSnp4GLF0llQxBlUHehnyRurc1aRU+P\nUgOtM6SbMJXXZ8TMko9pRxUXYRXGSxAz8ynl/iyEYHT8NF19JetV9waoz8xh35ylmpbK+lV1jI7s\nZzmWQgzmHbSmfultqnlsPnfK6U+hjfZ2OD09ZINPLDvqvpBbC6LCJS4GWdhnzw6i1RiAEObDyjWz\nnLuW8gzRkbvkvTzEIPA72NVOjl+fvos5Cun7vwdiPUHdPsVgH2BeBrB35Upcam0FVq9GbmYG6Xw+\nVJYl6AcAq4WPQ6ohYplSN6Hv+seJBGXibrpqsZRvgdlFw0sInHoBgf5a5wLCu325pZ0oJwDqYqU8\nAmH1h9S1q2VMqIFT1Pwc0SobfbcuM2WqAN19tD5gOJZyra2t6P70p8GPH4djMK8sBxLsBBGmbkJ/\n2D8mWYSN6zLZ5v2yFWa1h4vwzNq1tKPns+0xcCzpgFCT2Op2YJ7Vm1wyMIhBzIWYpV8AMKrcLxXg\nXK2PG9oGwiaXpoDqvEQbEtMgrLfFAVxuaQE4j21Pz7OVRCImiOVJ1YX+G2+8gb/927/F9PQ0vv/9\n75fMn0RXHBX7VuUKzELiCsxBVJwEfVCx7SCNIlU6SwhmSBuAEPaq22J9Nm1bzFbVOEw5mtpxEcza\nTeskUwAegPB1H4V8i+L+9bAhDwOANdFLvFypx0GyUJsEQQiqLvRvu+02fPOb38S+fdF7bZl/tPnD\nMfFbmM0of6td/w7mWahM12fBSfz/qLyOsG95OaN9G8AfIbCB3wqhokkj2WY0G7Y4suq5FLQv++er\nIXa+XkWxldHpiLbUdtSFYw6hs7epnmT9DPaBTu3vOICFuTlkX3ihpO7ekXVTsBSCSEwsoX/gwAGc\nPLwiWlkAABJ0SURBVHkS69evx4svvlhIHxkZwaOPPop8Po9HHnkEjz32WOyGmX+8L0Fnb4J5Vs61\na5s//b0QO3J1UaGXN7VhSt8AszfLvYa0VX67xT4hk8O1a2bI4yrn7RCz4pUIq4Ak+2AeRK4p9bwO\nYCOEyse0X8DUB5mm7mdwLXkAALmctS5T3RlS7xBEYmIJ/cHBQRw+fBgPPfRQIS2fz+PQoUM4deoU\nNm3ahDvvvBO7d+/Gli1bEnUgiefzuDp928wyZanDVGccbO3ID9VUr80KJodoi5ukOLC/BehIoewi\n7N5YF9bDhrIuos099Xalmoi3t4uEmAu20mxTrc/N1dr3J0E0H7GE/o4dO8A1D4Vnz57F5s2b4fh2\nzgMDAzhx4gQ2bNiAL33pS8hkMnj66aets3/mH8chthj39fWV9QDVYh8CHT1HMCN1EAi1SYSFziVL\nXQsAfmO5Z/NFcxliE9aAlncxcfwjq0Hdss6ZlSuBu+4S3xvOC+6SVbihD6bPbbEgD51Evaim+wVJ\n2Tr98+fPo6urq3Dd2dmJM2fOYO3atXj22Wdj1dEHoSeOK/DjLuQmhUEI/WH/ug/FqhwHxbNn26y9\nBUKlYkJfquT+8ZQhr+79kvt9kIMPR/1xStxjWpq7Zg2Y/yVmlv+7qZx+DWDxImT5Ttt02OK0Tixj\npPuFhvC9k0oltUMJw/zjPyYoY3POpqdHqXdsqPf0yFsM1fuB6+aXUfXq3i+Z4W85U3aELIJYYjSE\n751NmzZhYiLYWjUxMYHOzs7Y5bdBCNf3ErRpU3fo6W/AbKf+BkQsV1sdzD9XQwoyBIFQrvn3pJXO\nGzAL3ssAoj4JtczLWpqpPh2OYKb/it8/deDzIHwa/RzCikiWcZRyjqHecSSLN6DSD2EdtApiP0Qe\n4g2M+fflEel0QV2SyWbhtivvRP7MfQphyxyeyVRtsxZBLEUaYqa/fft2nDt3Dpxz3HrrrTh27BiO\nHj0au7yM1vT/JGgzrmvljTBb1dwHoYd3tfQchCDX65bXNn/zaj3SNJMDuBvFztj0vprumdJUuJ/H\nUfL2A/iepS6OsND9DYTp5mUI/z6AEMwexCC5Vckr25N+gF5FsBFsHsJfECBUTR/w6/lvlmdQ05zu\nbru6pKenoPoJpfuDRFHd5DeHWCYs+kx///79GBsbw5UrV9DV1YUnnngCg4ODOHLkCHbu3Il8Po+D\nBw8mttwB7C4TKiFKvaNGspJwmB16ce1ogvlH3bUyQ3mCPQonYfkchKCWvu7fgxD6nQCmITZVcQjB\n3QExUHVo7XGY1zMkrt9OK+KtMfBstmCaGRdaLCWI6hFL6Ntm8Lt27cKuXbvKaphBLJgmYQPMgkf3\n1BkVVnGjpQ5T2qyfPmu4p6MuJjOEVULvg/CTI71pdsAMRzBDl64VXP9cd08siVLHSHfLDKUHoqj7\nA0pfoPVDDigMdtcPHIGZ5hSAnqWmqokRlpEgaklDqHcqhfnHJAu5cRdoy1nINSFVPqX82ADBQMO0\n9AEU27freSSOco8hLLRt5SpbTo9HN8xBZqCl2WxpHD8f6+kRCRdrEUWhdtCbBlFvGmIht1IYxEw/\niQuECzALngvaddRCru7BUhIV3MTkbAwIL/5yS70m8cYt7WQRvCV0QLxhuAhm+uo9SZyFV7m+IOtX\nYRH9UfsFhFVFej1xyWSz1ghoBEGYaYqZPoeYxU4mKBNXvePAvJC7D2HLHJV+Q7ocZNTBhsO8M/Ul\nS70mD0SFXakIq0rkbLrf0oaL4ghZDEHgFSBYI7mMYLHVJmRNaSZsvn6A8DPYnkulI5czPttAayu6\nSV1CEEaaYqY/7B+TBOeohnrHdm8NzKoZhmIrGBO2D9Lk5lnOlOWfjm1x24Hd747rH4ctZWuFo/Sh\n1HNFkW5rIzUKQSwCdVfvJLELt3mo1NNtdXqw2/r/HwQDkAchrN/2r/VQhCZsaqpZv17phdOL6F8p\noiJy5SBm90xL5xCWOia4n38S4cXnWYjPtMdS7iUAD6RSWJFK4bLnwfU8zEB8Tq5f30BrK9JtbUA6\nHWyichwgmzX721ms3bUEsQRpKvXOWwnK2AReua6RVVIIXCRMQ+j+ZcAVNULVDMzYBiT5ZiEF/QSE\nF9AobOIvD3v0sAsQ9vKme7aFaLlnQG4kk0J+DMB/RfjtRp5ziMFl44oVmAewYWEByOcLX6QpiH0K\nvK2tEINWncGzvj7jQi7triUIO8tWvWMThnp6lHrnJsu9DpitbBjCkbRcS3mbi4gPoDgql2PJW4o7\nIDaBMcO9ASTzWArYdf2uIY1p5+xGsccjph6np4GxMVqwJYgGo25CX5Jklm6KbyvTVcqx3okbhcm2\nK9hmWWSrlyFQpwDhna5ZiP6b5r62N4r5iL5xiAXlOYiBQb51RDmqexMiLkDVN8+RzTtB1JW66/ST\nCBXdEZpEf1u4DWbrHZnPVIcpTeIo5zZzRZ6gXlsa0+6b8rmGNED4MrKV0dMZwu3pOFBs61Fdb560\nWEsQyammTj/leV6564rlN5pKFWab+wB8P2YX9qRSOG5KB3BcqePjqVTB0ZjKzyEGmT7DvQxQVDfz\n/z4O4Pf8tDcB9BrKj0IIyisQ7g7W++nXIVQ/b0IMRqrfGkcpP45gAJRvP2oMXImLQF2k3uMQM/00\nih2+vQ7gF4bnkn86zO/PuwBu9p9J/Ty50ne1PEOxuSZvb4ezZw8Je4KoAqlUCpWK7Lqrd5LsKI3r\nWtnmamHcP5ruuRHt3opA0LoR5WUeZsij3leR+WxO00phyqO3H2dHsY7qfE2arsZp29HTp6fBTAFI\nCIKoC3UT+lIQvZGgTNwduXOW8nOwm0tKc0NAWO9Ib5QMwmqmVrwE+8D3MoRK6gaECek8xLPujShj\nYh7Bs81CWAC5EM/8AIQTOjW4i9x1q5Y39i+Vgvv+9wMAcrOzuDg/j97Ff3EkCCIBdRP63RBqliS+\nd0rN4CVRDtdsqDtdmdaOq5w7lvK29FK0we7QrRXB2oTeJ4b4rED4TaVXKa8fVWRam+EeANzx8Y+H\nXCGzvj6jy2SCICqjmjr9Wng2jgVDci+bBEEQy5G+vj4wxqpSV92EPkEQBLH41M16Rzoiex3A/47Z\nhe2pVMGKRkWvIyofgKJ7HoTVyW3+9SUI6xvp3fItAB+CsKq5jMAy53cINnvJMr+DsN65BULvrlrv\n/Fe/jhTE+sJqv+680rYK98tIG3zpXXMlhF5+AcIJXco/v+7X+Y7fvlrPbf5zvgVgTSqF9lQKKQBT\nnoe852GtUnfKz/uO3+YFAB9MpbC6pQXzLS1oW7tW7KDVd9v6Ea6K0PIRBFEeTWG90wjMQWzakkj/\nPNIqSO50vY6wW+MpBEL/d/7xPZiDpMjFYLmRbIN/7IIQ6h7EgHILxMLtCr9P8/55B8KCf51/fEe5\nt84vKwX+RYh1kJv8+icA3NbaCrS0IL1mDbJXr+JDnlcwEW3xn/Fd//wWv9yHAOQ8D2hpwbwv8Hk2\nC2Qy2LNxo+iTHw1rKp1Gj+Jrhw0Pg7mucL+gU+FgQIMMQZSBVwcAeJ7/ty9BF/Yp5byIOqLy2e4N\nKecPR9wrVcaWd59yP059ca7VtDj39OfSr9Vy1ejjUG+v53meOJrq8u+XS63qJYhGpRoim3T6BEEQ\ny4i6u2EgCIIgoiGTTYIgiGUEmWwSBEEQZVE39Y70ePl6ZK4wr8Psf1+vo1Q+/V4OwFUl/S2Eo11d\n8q9z/rW06vktgt3Ab/p5LkE4gJOmj9JkU/ZJuoyQcX09BBY6lyEc0N1Qrvf459JsU9YtjbZknikI\n9wzz/n1AWO+M+22+DGG9M6BZ7/RHWO/ImMIt/rOvWbkSubVrwaT1Ti6HqXRa7FhWrHeYGinLPzIY\nqNSdMrlpJojE1M/LJvloIQiCSEQ1ZCepdwiCIJYRJPQJgiCWEVUX+teuXcPDDz+Mz3zmM/jud79b\n7eqXBNUyrWpUmvn5mvnZAHo+ogYLuT/4wQ/Q39+PT33qUxgYGMCnP/1pY77+lPAI3wi+d6S/HOlT\nx4NYhH0DwmeN6otnNQJf9vK+er4A4BUAt0Ms4P4WwAcRLAr/FsJlcgeCBdkUxOLtFQi3CTcQuFVI\nQSzEtiJwvTCvlJWuGy4CWOunzwFY09oKALi4sIANLS14e+VK/NEf/iEy2Sw6cjnMzM2hbbVwNj05\nO4v5+XlsbGkpLPLmZmcxD2B+7drArQIgFkkdB6PDw8iMjKAjlyvcmpydxToAaembRykT1y2CybUC\nz2aFj/9uLWpwDdwtjI6Oos/kMqJJoOcjYu3pHRwc9NavX+9t27YtlP7jH//Y+4M/+ANv8+bN3lNP\nPeV5nuf93d/9nfeb3/zG8zzP+/SnP22sD1Vww/B8TDcMej7bvSFD2oP+9ZClzD2W8yHLeRIXBw/H\nKKP2Tb//vJb+sKFetQ7T8z9vSB/q7fWGhobEsYRLBptbhOeff976P1ZdK5j+P6G+Gdwt6HVHXZvO\nh4aGrH2LS9Tzxc1nulcqLc6zNuPzqelL6fmSfjc9bxHdMAwODmJkZCSUls/ncejQIYyMjGB8fBxH\njx7FK6+8gs7OTkxMTAAAFhYWTNVVhdEq5DPd09Oiri9ZzqsBj5FntIx7er16vtEY55US9xU8Kpft\nnl531LXtvFJiP19EPtO9Umlxn7VSGu35qq3SWaznq8d3E0hgssk5x3333YcXX3wRAPCrX/0Kjz/+\neGEweOqppwAAn//853Ho0CGk02ns2LED+/fvL240lSTYH0EQBCGJKbKtlK3TP3/+PLq6ugrXnZ2d\nOHPmDNasWYNvfetbkWUr7TRBEARRHmVb79BsnSAIYulRttDftGlTQXcPABMTE+js7KxKpwiCIIja\nULbQ3759O86dOwfOOa5fv45jx45h9+7d1ewbQRAEUWViCf39+/fjrrvuwquvvoquri788z//M1as\nWIEjR45g586d2Lp1Kx588EFs2bKl1v0lCIIgKqAuDtcIgiCI+tAQvnea3XXDG2+8gUceeQT79u2r\nd1dqwokTJ/CZz3wGAwMD+OlPf1rv7lSdbDaLz33uc+jv78c//dM/1bs7VefatWu48847cfLkyXp3\npeqMjo5ix44d+NznPoexsbF6d6fqeJ6HL3/5y/j85z+Pf/mXf4lVpiGEvnTd8I1vfAP/9m//Vu/u\nVJ3bbrsN3/zmN+vdjZpx//334xvf+AaeffZZHDt2rN7dqTrd3d145pln8Nxzz+EnP/lJvbtTdb76\n1a/iwQcfrHc3akJLSwtuvvlmzM3NNaWhyfHjx3H+/HmsWrUq9vPVTOgfOHAAGzZswEc+8pFQ+sjI\nCLq7u3H77bfj6aefBhC2+W/1fcY0OkmebylSzvM9+eSTOHTo0GJ2s2ySPt8Pf/jDgj+pRifJs/30\npz/F1q1bccstt9Sjq2WR5Pl27NiBH/3oR3jqqacwNDRUj+4mJsnzvfrqq/jYxz6Gr33ta3jmmWfi\nNVCxIwcLv/jFL7xf//rXIX898/Pz3oc//GHvjTfe8K5fv+599KMf9cbHx71//dd/9f793//d8zzP\nGxgYqFWXqkqS55Ps3bu3Hl0tiyTPt7Cw4P3lX/6ld+rUqTr2OBnl/P88z/N279692F1NTJJn+/KX\nv+w9+uij3p/8yZ94999/v7ewsFDHnsejnP/d3Nzckvn9JXm+73znO973vvc9z/M8r7+/P1b9NQuX\nuGPHDnDNW+LZs2exefNmOH44u4GBAZw4caLguuHkyZNLxuwzyfNt2LABX/rSl5DJZPD000/jscce\nW/wOJyTJ8506dQo/+9nPcPXqVbz22mv47Gc/u/gdTkiS57t06RJ+8IMfIJfL4Z577ln8ziYkybM9\n+eSTAIBvf/vbuOWWW5bEpsskz5fNZvGTn/wEU1NTOHz48OJ3tgySPN8XvvAFHD58GL/85S9jexdd\n1Bi5lbhuWArYnm/t2rV49tln69iz6mB7vq9//etL5gcVhe35ent70dvbW8eeVY7t2SQPP/xwPbpV\nNWzP98UvfhEPPPBAHXtWHWzPd9NNNyVeL1zUhdylMIuoBHq+pU0zP18zPxtAz5eERRX6ze66gZ5v\nadPMz9fMzwbQ8yVhUYV+s7tuoOdb2jTz8zXzswH0fImoyfKzJ6xwPvjBD3qrVq3yOjs7vW9961ue\n53nej370I+/3f//3vQ9/+MPeV77ylVo1X3Po+ej5GpVmfjbPo+er9PnIDQNBEMQyoiF25BIEQRCL\nAwl9giCIZQQJfYIgiGUECX2CIIhlBAl9giCIZQQJfYIgiGUECX2CIIhlBAl9giCIZQQJfYIgiGXE\n/w+Nu2GEOlQ+TwAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 78 }, { "cell_type": "code", "collapsed": false, "input": [ "fig = plt.semilogx(x, y, 'g^')\n" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAEECAYAAADQ7bj8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xt4VfWd7/H3DgkocglECDabGkl2hAAFrLk87VCCmQSF\nKehECTwK4Ta26WihelprWyTMaQX7lGmVkdZOYRpw2nBRCbYlgJym2vYYjij1HGjHVCPkJtOGQJFL\nbvt7/tgXEhJCVrJhJ/B5PU+e7L0uv/X7rez9+6zf2muvuMzMEBERcSAi3BUQEZG+R+EhIiKOKTxE\nRMQxhYeIiDim8BAREccUHiIi4lin4XH+/HnS0tKYPHkyycnJPPnkkwAUFBTgdruZMmUKU6ZMYffu\n3cF11qxZg8fjYezYsezduzc4/eDBg0ycOBGPx8Py5cuD0xsaGsjNzcXj8ZCens7Ro0eD8woLC0lK\nSiIpKYnNmzeHrNEiItJDdhlnzpwxM7OmpiZLS0uzN954wwoKCmzdunXtlj18+LBNmjTJGhsbraKi\nwhISEszr9ZqZWUpKipWVlZmZ2T333GO7d+82M7Pnn3/e8vPzzcysqKjIcnNzzcysrq7OxowZY/X1\n9VZfXx98LCIi4XfZ01YDBw4EoLGxkZaWFoYNGxYInXbLFhcXM3/+fKKiooiPjycxMZGysjJqa2s5\nffo0qampACxcuJCdO3cCsGvXLvLy8gDIyclh//79AOzZs4fs7Gyio6OJjo4mKyuLkpKSEMSliIj0\nVOTlFvB6vdxxxx28//775OfnM378eHbs2MH69evZvHkzd955J+vWrSM6OpqamhrS09OD67rdbqqr\nq4mKisLtdgenx8XFUV1dDUB1dTWjR4/2VSYykqFDh1JXV0dNTU2bdQJlteZyuXrWehGR61RHAwAn\nLjvyiIiI4NChQ1RVVfH6669TWlpKfn4+FRUVHDp0iFtuuYXHH3+8R5XoCTPr0c+qVat6vFxH8y6e\n1tnzrjy+ku1z2rbe0r4r9bcLRfuu1t+uO+1z+nq91tp3qbaG470XivZ1570XCl2+2mro0KHMmjWL\nt956i5EjR+JyuXC5XCxbtowDBw4AvhFFZWVlcJ2qqircbjdxcXFUVVW1mx5Y59ixYwA0Nzdz6tQp\nYmJi2pVVWVnZZiQSKhkZGT1erqN5F0/r7PmlHodCV8pz2raOpoejfVfqb9fR9GupfU5fr9da+y7V\n1nC89y63XG/uW7BO/OUvfwl+SH327FmbOnWqvfbaa1ZbWxtc5l//9V9t/vz5ZnbhA/OGhgb74IMP\nbMyYMcEPzFNTU+3NN980r9fb7gPzL37xi2Zm9vOf/7zNB+a33Xab1dfX24kTJ4KPW7tM9fu8VatW\nhbsKV9S13L5ruW1mal9fF4q+s9PPPGpra8nLy8Pr9eL1elmwYAGZmZksXLiQQ4cO4XK5uO2223jh\nhRcASE5OZu7cuSQnJxMZGcmGDRuCn0ts2LCBRYsWce7cOWbOnMndd98NwNKlS1mwYAEej4eYmBiK\niooAGD58OCtXriQlJQWAVatWER0dHdrk7OVCfqTQy1zL7buW2wZqn4DLn0J9ksvlCtn5OxGR60Uo\n+k59w1xERBxTeIiIiGMKDxERcUzhISIijik8RETEMYWHiIg4pvAQERHHFB4iIuKYwkNERBxTeIiI\niGMKDxERcUzhISIijik8RETEMYWHiIg4pvCQ64Zu3y8SOgoPuS6YGcv+eZkCRCREFB5yXXjp1ZfY\nfmQ7L//i5XBXReSaoP8kKNc8M+Mzcz/Dm+PfJP1wOr/f9vvgv0cWuR7pPwmKdMFLr77Eu4PfBRe8\nO/hdjT5EQkAjD7mmtR514AIMjT7kuqeRh8hltB51ABp9iISIRh5yTVvy6BI++NsHF8IDwGDMkDFs\nWr8pbPUSCadQ9J0KDxGR68wVP211/vx50tLSmDx5MsnJyTz55JMAnDhxgqysLJKSksjOzubkyZPB\nddasWYPH42Hs2LHs3bs3OP3gwYNMnDgRj8fD8uXLg9MbGhrIzc3F4/GQnp7O0aNHg/MKCwtJSkoi\nKSmJzZs396ihIiISQnYZZ86cMTOzpqYmS0tLszfeeMO++tWv2jPPPGNmZmvXrrUnnnjCzMwOHz5s\nkyZNssbGRquoqLCEhATzer1mZpaSkmJlZWVmZnbPPffY7t27zczs+eeft/z8fDMzKyoqstzcXDMz\nq6urszFjxlh9fb3V19cHH7fWheqLiMhFQtF3XvYD84EDBwLQ2NhIS0sLw4YNY9euXeTl5QGQl5fH\nzp07ASguLmb+/PlERUURHx9PYmIiZWVl1NbWcvr0aVJTUwFYuHBhcJ3WZeXk5LB//34A9uzZQ3Z2\nNtHR0URHR5OVlUVJSUlIg1NERLon8nILeL1e7rjjDt5//33y8/MZP348x48fJzY2FoDY2FiOHz8O\nQE1NDenp6cF13W431dXVREVF4Xa7g9Pj4uKorq4GoLq6mtGjR/sqExnJ0KFDqauro6amps06gbIu\nVlBQEHyckZFBRkaGg+aLiFz7SktLKS0tDWmZlw2PiIgIDh06xKlTp5gxYwa//vWv28x3uVxhvV6+\ndXiIiEh7Fx9Yr169usdldvl7HkOHDmXWrFkcPHiQ2NhYPvroIwBqa2sZOXIk4BtRVFZWBtepqqrC\n7XYTFxdHVVVVu+mBdY4dOwZAc3Mzp06dIiYmpl1ZlZWVbUYiIiISPp2Gx1//+tfglVTnzp1j3759\nTJkyhdmzZ1NYWAj4roi69957AZg9ezZFRUU0NjZSUVFBeXk5qampjBo1iiFDhlBWVoaZsWXLFubM\nmRNcJ1DWjh07yMzMBCA7O5u9e/dy8uRJ6uvr2bdvHzNmzLgye0FERBzp9LRVbW0teXl5eL1evF4v\nCxYsIDMzkylTpjB37lw2btxIfHw827ZtAyA5OZm5c+eSnJxMZGQkGzZsCJ7S2rBhA4sWLeLcuXPM\nnDmTu+++G4ClS5eyYMECPB4PMTExFBUVATB8+HBWrlxJSkoKAKtWrSI6OvqK7QgREek6fUlQROQ6\no3tbiYhIWCg8RETEMYWHiIg4pvAQERHHFB4iIuKYwkNERBxTeIiIiGMKDxERcUzhISIijik8RETE\nMYWHiIg4pvAQERHHFB4iIuKYwkNERBxTeIiIiGMKDxERcUzhISIijik8RETEMYWHiIg4pvAQERHH\nFB4iIuKYwkNERBzrNDwqKyuZPn0648ePZ8KECTz33HMAFBQU4Ha7mTJlClOmTGH37t3BddasWYPH\n42Hs2LHs3bs3OP3gwYNMnDgRj8fD8uXLg9MbGhrIzc3F4/GQnp7O0aNHg/MKCwtJSkoiKSmJzZs3\nh6zRIiLSQ9aJ2tpae+edd8zM7PTp05aUlGRHjhyxgoICW7duXbvlDx8+bJMmTbLGxkarqKiwhIQE\n83q9ZmaWkpJiZWVlZmZ2zz332O7du83M7Pnnn7f8/HwzMysqKrLc3FwzM6urq7MxY8ZYfX291dfX\nBx+3dpnqi4hIB0LRd3Y68hg1ahSTJ08GYNCgQYwbN47q6upA6LRbvri4mPnz5xMVFUV8fDyJiYmU\nlZVRW1vL6dOnSU1NBWDhwoXs3LkTgF27dpGXlwdATk4O+/fvB2DPnj1kZ2cTHR1NdHQ0WVlZlJSU\nhCQwRUSkZ7r8mceHH37IO++8Q3p6OgDr169n0qRJLF26lJMnTwJQU1OD2+0OruN2u6murm43PS4u\nLhhC1dXVjB49GoDIyEiGDh1KXV3dJcsSEZHwi+zKQh9//DH3338/zz77LIMGDSI/P5+nnnoKgJUr\nV/L444+zcePGK1rRSykoKAg+zsjIICMjIyz1EBHprUpLSyktLQ1pmZcNj6amJnJycnjooYe49957\nARg5cmRw/rJly/j85z8P+EYUlZWVwXlVVVW43W7i4uKoqqpqNz2wzrFjx/jEJz5Bc3Mzp06dIiYm\nhri4uDaNrays5K677mpXv9bhISIi7V18YL169eoel9npaSszY+nSpSQnJ7NixYrg9Nra2uDjV155\nhYkTJwIwe/ZsioqKaGxspKKigvLyclJTUxk1ahRDhgyhrKwMM2PLli3MmTMnuE5hYSEAO3bsIDMz\nE4Ds7Gz27t3LyZMnqa+vZ9++fcyYMaPHDRYRkZ7rdOTxu9/9jhdffJFPfepTTJkyBYCnn36an//8\n5xw6dAiXy8Vtt93GCy+8AEBycjJz584lOTmZyMhINmzYgMvlAmDDhg0sWrSIc+fOMXPmTO6++24A\nli5dyoIFC/B4PMTExFBUVATA8OHDWblyJSkpKQCsWrWK6OjoK7MXRETEEZd1dNlUH+FyuTq86ktE\nRC4tFH2nvmEuIiKOKTxERMQxhYeIiDim8BAREccUHiIi4pjCQ0REHFN4iIiIYwoPERFxTOEhIiKO\nKTxERMQxhYeIiDim8BAREccUHiIi4pjCQ0REHFN4iIiIYwoPERFxTOEhIiKOKTxERMQxhYeIiDim\n8BAREccUHiIi4pjCQ0REHFN4iIiIY52GR2VlJdOnT2f8+PFMmDCB5557DoATJ06QlZVFUlIS2dnZ\nnDx5MrjOmjVr8Hg8jB07lr179wanHzx4kIkTJ+LxeFi+fHlwekNDA7m5uXg8HtLT0zl69GhwXmFh\nIUlJSSQlJbF58+aQNVpERHrIOlFbW2vvvPOOmZmdPn3akpKS7MiRI/bVr37VnnnmGTMzW7t2rT3x\nxBNmZnb48GGbNGmSNTY2WkVFhSUkJJjX6zUzs5SUFCsrKzMzs3vuucd2795tZmbPP/+85efnm5lZ\nUVGR5ebmmplZXV2djRkzxurr662+vj74uLXLVF9ERDoQir6z05HHqFGjmDx5MgCDBg1i3LhxVFdX\ns2vXLvLy8gDIy8tj586dABQXFzN//nyioqKIj48nMTGRsrIyamtrOX36NKmpqQAsXLgwuE7rsnJy\ncti/fz8Ae/bsITs7m+joaKKjo8nKyqKkpCTk4SkiIs5FdnXBDz/8kHfeeYe0tDSOHz9ObGwsALGx\nsRw/fhyAmpoa0tPTg+u43W6qq6uJiorC7XYHp8fFxVFdXQ1AdXU1o0eP9lUmMpKhQ4dSV1dHTU1N\nm3UCZV2soKAg+DgjI4OMjIyuNklE5LpQWlpKaWlpSMvsUnh8/PHH5OTk8OyzzzJ48OA281wuFy6X\nK6SVcqJ1eIiISHsXH1ivXr26x2Ve9mqrpqYmcnJyWLBgAffeey/gG2189NFHANTW1jJy5EjAN6Ko\nrKwMrltVVYXb7SYuLo6qqqp20wPrHDt2DIDm5mZOnTpFTExMu7IqKyvbjERERCR8Og0PM2Pp0qUk\nJyezYsWK4PTZs2dTWFgI+K6ICoTK7NmzKSoqorGxkYqKCsrLy0lNTWXUqFEMGTKEsrIyzIwtW7Yw\nZ86cdmXt2LGDzMxMALKzs9m7dy8nT56kvr6effv2MWPGjNDvARERca6zT9PfeOMNc7lcNmnSJJs8\nebJNnjzZdu/ebXV1dZaZmWkej8eysrLaXAX1ne98xxISEuz222+3kpKS4PS33nrLJkyYYAkJCfbo\no48Gp58/f94eeOABS0xMtLS0NKuoqAjO27RpkyUmJlpiYqL99Kc/bVe/y1RfREQ6EIq+0+UvqE9y\nuVz04eqLiIRFKPpOfcNcREQcU3iIiIhjCg8REXFM4SEiIo4pPERExDGFh4iIOKbwEBERxxQeIiLi\nmMJDREQcU3iIiIhjCg8REXFM4SEiIo4pPERExDGFh4iIOKbwEBERxxQeIiLimMJDREQcU3iIiIhj\nCg8REXFM4SEiIo4pPERExDGFh4iIOKbwEBERxzoNjyVLlhAbG8vEiROD0woKCnC73UyZMoUpU6aw\ne/fu4Lw1a9bg8XgYO3Yse/fuDU4/ePAgEydOxOPxsHz58uD0hoYGcnNz8Xg8pKenc/To0eC8wsJC\nkpKSSEpKYvPmzSFprIiIhIh14vXXX7e3337bJkyYEJxWUFBg69ata7fs4cOHbdKkSdbY2GgVFRWW\nkJBgXq/XzMxSUlKsrKzMzMzuuece2717t5mZPf/885afn29mZkVFRZabm2tmZnV1dTZmzBirr6+3\n+vr64OOLXab6IiLSgVD0nZ2OPKZOncqwYcM6Cpx204qLi5k/fz5RUVHEx8eTmJhIWVkZtbW1nD59\nmtTUVAAWLlzIzp07Adi1axd5eXkA5OTksH//fgD27NlDdnY20dHRREdHk5WVRUlJSU8yUkREQiiy\nOyutX7+ezZs3c+edd7Ju3Tqio6OpqakhPT09uIzb7aa6upqoqCjcbndwelxcHNXV1QBUV1czevRo\nX0UiIxk6dCh1dXXU1NS0WSdQVkcKCgqCjzMyMsjIyOhOk0RErlmlpaWUlpaGtEzH4ZGfn89TTz0F\nwMqVK3n88cfZuHFjSCvlROvwEBGR9i4+sF69enWPy3R8tdXIkSNxuVy4XC6WLVvGgQMHAN+IorKy\nMrhcVVUVbrebuLg4qqqq2k0PrHPs2DEAmpubOXXqFDExMe3KqqysbDMSERGR8HIcHrW1tcHHr7zy\nSvBKrNmzZ1NUVERjYyMVFRWUl5eTmprKqFGjGDJkCGVlZZgZW7ZsYc6cOcF1CgsLAdixYweZmZkA\nZGdns3fvXk6ePEl9fT379u1jxowZPW6siIiERqenrebPn89vfvMb/vrXvzJ69GhWr15NaWkphw4d\nwuVycdttt/HCCy8AkJyczNy5c0lOTiYyMpINGzbgcrkA2LBhA4sWLeLcuXPMnDmTu+++G4ClS5ey\nYMECPB4PMTExFBUVATB8+HBWrlxJSkoKAKtWrSI6OvqK7QQREXHGZR1dOtVHuFyuDq/8EhGRSwtF\n36lvmIuIiGMKDxERcUzhISIijik8RETEMYWHiIg4pvAQERHHFB4iIuKYwkNERBxTeIiIiGMKDxER\ncUzhISIijik8RETEMYWHiPRJuilqeCk8RKTPMTOW/fMyBUgYKTxEpM956dWX2H5kOy//4uVwV+W6\npf/nISJ9ipnxmbmf4c3xb5J+OJ3fb/t98B/PSdfo/3mIyHXnpVdf4t3B74IL3h38rkYfYaKRh4j0\nGa1HHbgAQ6OPbtDIQ0SuK61HHYBGH2GkkYeI9BlLHl3CB3/74EJ4ABiMGTKGTes3ha1efU0o+k6F\nh4jIdUanrUREJCwUHiIi4lin4bFkyRJiY2OZOHFicNqJEyfIysoiKSmJ7OxsTp48GZy3Zs0aPB4P\nY8eOZe/evcHpBw8eZOLEiXg8HpYvXx6c3tDQQG5uLh6Ph/T0dI4ePRqcV1hYSFJSEklJSWzevDkk\njRURkdDoNDwWL15MSUlJm2lr164lKyuL9957j8zMTNauXQvAkSNH2Lp1K0eOHKGkpIQvfelLwXNq\n+fn5bNy4kfLycsrLy4Nlbty4kZiYGMrLy/nKV77CE088AfgC6l/+5V84cOAABw4cYPXq1W1CSkRE\nwqvT8Jg6dSrDhg1rM23Xrl3k5eUBkJeXx86dOwEoLi5m/vz5REVFER8fT2JiImVlZdTW1nL69GlS\nU1MBWLhwYXCd1mXl5OSwf/9+APbs2UN2djbR0dFER0eTlZXVLsRERCR8Ip2ucPz4cWJjYwGIjY3l\n+PHjANTU1JCenh5czu12U11dTVRUFG63Ozg9Li6O6upqAKqrqxk9erSvIpGRDB06lLq6Ompqatqs\nEyirIwUFBcHHGRkZZGRkOG2SiMg1rbS0lNLS0pCW6Tg8WnO5XGH/Vmfr8BARkfYuPrBevXp1j8t0\nfLVVbGwsH330EQC1tbWMHDkS8I0oKisrg8tVVVXhdruJi4ujqqqq3fTAOseOHQOgubmZU6dOERMT\n066sysrKNiMREREJL8fhMXv2bAoLCwHfFVH33ntvcHpRURGNjY1UVFRQXl5Oamoqo0aNYsiQIZSV\nlWFmbNmyhTlz5rQra8eOHWRmZgKQnZ3N3r17OXnyJPX19ezbt48ZM2aEpMEiIhIC1ol58+bZLbfc\nYlFRUeZ2u23Tpk1WV1dnmZmZ5vF4LCsry+rr64PLf+c737GEhAS7/fbbraSkJDj9rbfesgkTJlhC\nQoI9+uijwennz5+3Bx54wBITEy0tLc0qKiqC8zZt2mSJiYmWmJhoP/3pTzus32WqLyIiHQhF36nb\nk4iIXGd0exIREQkLhYeIiDim8BAREccUHiLSq+lzzd5J4SEivZaZseyflylAeiGFh4j0Wi+9+hLb\nj2zXv5nthXSproj0SmbGZ+Z+hjfHv0n64XR+v+33Yb8d0rVCl+qKyDXrpVdf4t3B74IL3h38rkYf\nvYxGHiLS67QedeACDI0+QkgjDxG5JrUedQAaffRCGnmISK+z5NElfPC3Dy6EB4DBmCFj2LR+U9jq\nda0IRd+p8BARuc7otJWIiISFwkNERBxTeIiIiGMKDxGRHrheP3dVeIiIdNP1fO8thYeISDddz/fe\n0qW6IiLd0JfvvaVLdUVEwuR6v/eWRh4iIg719XtvaeQhIhIGuveWRh4iIo719XtvhfXeVvHx8QwZ\nMoR+/foRFRXFgQMHOHHiBLm5uRw9epT4+Hi2bdtGdHQ0AGvWrGHTpk3069eP5557juzsbAAOHjzI\nokWLOH/+PDNnzuTZZ58FoKGhgYULF/L2228TExPD1q1bufXWW9tWXuEhIuJYWE9buVwuSktLeeed\ndzhw4AAAa9euJSsri/fee4/MzEzWrl0LwJEjR9i6dStHjhyhpKSEL33pS8GK5+fns3HjRsrLyykv\nL6ekpASAjRs3EhMTQ3l5OV/5yld44oknetRQEREJnR595nFxcu3atYu8vDwA8vLy2LlzJwDFxcXM\nnz+fqKgo4uPjSUxMpKysjNraWk6fPk1qaioACxcuDK7TuqycnBz279/fk6qKiITU9X7WI7K7K7pc\nLv7+7/+efv368YUvfIF/+qd/4vjx48TGxgIQGxvL8ePHAaipqSE9PT24rtvtprq6mqioKNxud3B6\nXFwc1dXVAFRXVzN69GhfJSMjGTp0KCdOnGD48OFt6lFQUBB8nJGRQUZGRnebJCLSJYFvlv/k+Z/0\niaurSktLKS0tDWmZ3Q6P3/3ud9xyyy385S9/ISsri7Fjx7aZ73K5rspObR0eIiJXQ+Cb5TN/MZOc\nz+eEuzqXdfGB9erVq3tcZrdPW91yyy0AjBgxgvvuu48DBw4QGxvLRx99BEBtbS0jR44EfCOKysrK\n4LpVVVW43W7i4uKoqqpqNz2wzrFjxwBobm7m1KlT7UYdIiJXm5mxbss6Tmec5nubv3fdnr7qVnic\nPXuW06dPA3DmzBn27t3LxIkTmT17NoWFhQAUFhZy7733AjB79myKiopobGykoqKC8vJyUlNTGTVq\nFEOGDKGsrAwzY8uWLcyZMye4TqCsHTt2kJmZ2ePGioj01PX+zfKAbl2qW1FRwX333Qf4RgUPPvgg\nTz75JCdOnGDu3LkcO3as3aW6Tz/9NJs2bSIyMpJnn32WGTNmABcu1T137hwzZ87kueeeA3yX6i5Y\nsIB33nmHmJgYioqKiI+Pb1t5XaorIldRX/9meYD+h7nCQ0Suoh27dpC3M4+zt54NTht4dCCb79vc\nJz77CFB4KDz6LDPrU0dq0veF4jXX179ZHqDwUHj0SX3tMkfp+/Saa0s3RpQ+6Xr+BzoSHnrNhZ5G\nHnJV9eV/oCN9k15z7WnkIX3O9XKZow5qeo/r5TV3tSk85KoJfLnq7Cd9V6qc/eTZa/JLVoHz69da\nu/qi6+U1Fw4KD7lqrpd/oHMlzq93pbNTh9je9fKaCwd95iFXzbVymWNnrsT59a5cKaSriTp2Pbzm\nukOX6rpceL1evVmk12j9JbJQfXlsx64dLPnXJfzH4/9xybK6soxIgD4whx4PP/twdkovcyXOr3fl\nJny6UZ+EQ58Pj568WfTBpoTSlTi/3pUrhXQ1kYRDnw+PnrxZ9MUhCaVf7fsVKS0pTKuYFvxJaUnh\nl3t/2a3yujKS0dVEEi59/jMPVnXvrpb64pD0dl25CV9vvVGf7l3Wu+kzD+j2UD2UQ/0+nL/Si3Vl\nJBPq0U4o6HTw9aHPjzym5U1zfOldKO/Jr0skr306inZGV371fhp5AI/mPEppYamja7ZD+cGmPjdx\n5lIv2HAdw1xuu909ir6S7elNx3sX1+VKXPnVm9orF/T58OjOCzRUQ/2rdYlkVzq4vuBSHXFXOuie\nXu7qpD6tdefgoKuB0502Bcr2er2Oyu/KY6f16qidlzsd3J0Q1imwXsr6MMAGLh5oO3btCMv2txdv\nt4GLBxoFV64eXq/XluQvMa/X2635vcn24u02eNrgdvvpUtMDWrexdTsvbnNH8zrbP5fbbktLi6Xf\nn26swtLvT7eWlpZO2xfYxuXK7axNl7O9eLsN+twgmz5zui3OX9ymnZ2Vf/HjlpaWDveNk9fTxe30\ner3B/UUBwf3Wlb9FV7choRGKrr/Ph8fFL9Cr5XJvlJ6U29rl3jxX6s0V6na03l8XdygXT7/UPthe\nvL3DDjFQTkfztu3c1uH+uTgYOqrv9FnT7cZFNxoF2I2LbrS7Zt0V7HQvbmPrTjn9/nTjqfbt7KxN\nlwumNvvqQcw11mU3fvZG27Frh7W0tLTbF61fF60fb9u5zW6aepNNnzXdtu+8UIeL63W511NHf7fW\nB1OBn9YHVU5fq5d6zUjPKTz84dGTo/6uvGk7crk3Sndc3CF2pYPr6pvrUkfsnXWEl+r4utqO1vt2\n+86OR2kXj962F29vd0SdlpNmrMI8mR4b9LlB7TrEQDkXd5bbdm6zEUkjOuzIWwdDoD6B4GlpabHt\nO7dbv5R+xoP+v++DWMSdETZ91nRb/MXFwQ470MbAiOCxbz3mKzfdFziBci/en4G/myfTYzdNvcnu\nmnXXZffz9uLtdkPeDUYavp9VWFpOmk2fNd0G/t2FNizOX2yp/5ganB/YVlpOmt3sudkYi3GHb9us\nwkZ8ekQwFC8X5K3rcvHfc/Eji23awmn2uYWfs2l503w/C6fZon9e1KWyOwrYKz2yv14pPMCSs5Nt\n2sJptvg5nnJ3AAAPdUlEQVSRxY7Xb2lpscRPJXYrQAJvlGl502zUZ0YF3yiBelxcZle20boD9Hq9\nNv0fpnf65unozXW5MOjocXNzsy3+4oUOO9ARdtTxddaei9cPdIgtLS024tMj2o3SWodjYHri3yVa\nVFqUbS/e7hs5vLLN+k3oZ6zCIuZHGA+27RAvLictJ80XNk9hIyeMND7t6/hb759tO7e1DQZ/OXkP\n59niLy62hIkJ5rnL17GSjrHS//tBjE9jAyYNsMdXPm6Dpw22jFkZ1tTU5NtmKtZvVD9j/oXl0+9P\nD45+Akf424u3B4MrYn6EcRcWkRIRnN/S0hLszNuF6DiMDHy/V2GRWZG+No7EUu5Lse07t9uNn73R\nXAkuYxUWlRVlAxYOMAow1zyXLzjGYtyFMc/f/nnY4ysfb1OvQJB39LfvbNTd0ShoSf6SdgcP23Zu\n63TkeKVG9uKj8MDXWXT3BfWVb37F+AT22Lceu+yyDQ0NHU7fvnO79buzX5uhf0tLi40ZPybYwba0\ntFjChIQ2HW5jY2PwcXNzc7sjs62vbLWIcRFt3jyBtgaWT7kv5cJ8/9F13hfyrKWlxfIezgu+CVuf\nIml9qiLQySenJVvkrZHBDjvlvhQj1be9ra9sDS7f3NxsLS0t1tTUZImfSrTGxkZramqypqYm83q9\nlveFPGtubvZ1cvN9HeLWV7baim+sMOZeaEfgNNBj33rswuhtlb+TvsX3O/GuRHto6UM2KGmQryN+\niGBn7prusqiHooyV2A15N9jj33rc1+k95e8gM3ydPLf69guf8pWZlpNmeQ/n2ZAJQy4Eg78+/R/q\nb9yMue50+dbLJdi5D4wb6HseOOKfhA2cOND4BkYMNuQTQ8w13d8xfxrD4y/3dizqoSgbPGFw8Ai/\nsbHRNyoIbHsVRqKvngkZCdbQ0GDxyfGWODHRFixbYIu+sMhaWlps6ytbrV9mP2OyLyiCYZZ2IdS4\nxbffWOUPiLkY8f594H+NkIaR4p/f6rUVeXuk7+/eql6BkV7gtR14PW19ZWu7UfeAhQNse/F2m/4P\n09uNCAd9blBwlBPYXlRiVPAAJbBc69N4rcMmVCN7uSAU4dHnv+cRcVcE21Zsc3w9udfrZUDCAJo/\nbCYyPpKG9xuIiOj44rPGxkYGDBhAQ0MD/fv3D043Mzx3eXi/9H0SpidQvr8cl8vF8ieX89za51j+\njeX84Ds/YMU3VvDsmmdZ8c0VfP/b36epqYn+/fvT2NhIREQEkQMi2bp9K4tfXczZHWcZeP9ABh8c\nzPGPj8NCoAb4BLAPih4vYt4/zuM/t/8nDz7wIDwFfA84C/0f6k/ji43B+hW9XETeojyik6I5Pv04\nUS9H0VTbBA8Cm4FGwA1UAgnAx/B34/+O3/6v30IkEAdUA8uBF3zzI4dE0vy3ZrgJOHNhH33tW1/j\nuz/4LgP6D6DhVAP0B74E/Mi/3I3ADUAEJE5J5M9v/RlOAoNg+O3DOfH+CTgPjPdthyZ/vT4JPATs\nAbKAI8CvgCjA/GV6gdHAcSAG+BBoBub4H9f7l/0U8Hvgc8DtwHvAn3x14iN/Qwb5t70w8EIBfgak\n+Zf5GzDQP+8d//I3AnXAYGCqf50bgHNACXA3kOzb1sBfD+Ts2LPQ4N/nif56uIDXfH9Hovxlfuzb\n1rQp03jr7bc4c/aMb95IIBco9rfjKDAD+E//32w68F/Aq/7XzXAg278P/wacAqb51w34k78e/f11\nSgT+CPSDm9+6mdrDtTz86MOYGS/ueJFxd4xj+CeG89s3f0vLgBZucN1A2ifTeOPPb+Cd6yX9cDq/\n2/o7Ppv7Wd4c8CYRFoE3yXthe38EIiG9odVy49/E81sPtU21pESn4I32hu1W6qWlpWRkZFzx7YTL\nNX9L9pKSElasWEFLSwvLli3jiSeeaDM/8MWttJw0/vf2/+3oi1yPfesxvv+d7wc7wce+9Rjr/ue6\nDpd1uVzB5Vrvrh3FO3jg3geC87YXb+cf/+Ef6devX3BaU1MTUVFRwectLS1t5uPC1xmdxTfN/NM+\nxtdRNHFh2Uv9Huhf75z/ceA9et7/vAlfGDS3KrOff50I/+8G/7T+rdZvAgb469Z6GzcQ/HIl1qqs\nQf76BJb3+rfXgq9jPeEv/0K+XQgAD/CBfzkXcNq/3CfxdZR/BHb717kN+LN/uy3+7biAaP96Sfg6\n+unAIX8bavCFSCyQB/wHsBjY6p832L+t/8LXEY/F1/FX+fezP/iYg6/TbvBvwwOU+/fvJ/37535g\nE7AEX4f+YKv9/DN/uxYBhf79NQ5fSEcD7/vrONTfxkAQDfXXowbI9G/3j8Bf8AVEIr4AeBd4AHgD\nX1DUALf463qTv521/m31a/W3PMmF0NwKkf0iaXY1+9pzHuKGxlHXv47GyEa8p70MrRvKkoeX8P2X\nvw+5ELErghEtIzg+4TgkwYAPBvDIrY/ww6ofcrb8LBFNEYwbMY6bY27m9QOvY4MN+sNAz0Dy3fm+\n5W49S8R/ReCN8JLeEN5bBhUUFFBQUBCWbV8N1/SXBFtaWnjkkUcoKSnhyJEj/PznP+ePf/xj+wWj\n4f/U/x9H1+F7vV7W/+d63xvI7SvjuRefw+v1tlu2sbGxzXKNjb6ez8z4+g++3mbe13/wdVZ8Y0Wb\naXdk3NHmeZv5Q4Fh/sfD8L2Z3fg6omhghP/3uC78vvWix4EOaLT/9634OpFkf7mj/MuM968X2F5y\nq3k3+9cPlPtJfz2T/fNi/T/j/M8j/POG+ttxK76j3jh/GYP80wf5lxmBrzOLA2b7l43F15kl+ucF\nOtRx/ue3+/fPOHwBNdZf39H+n2HAP/jr82f/OhOBIf6yP4Nv9DHK//tGfB3xMP82E4A38Y2YTraq\nz1n/tlz+3/3xjQAG+dsRqEssvhDE/zu11XOXv/3p/sep+EZKo/GF7O3+MocB9/jrdpu/vYF2DfHX\nE//yNfheN4Hn4BtB1ACz/NuZ6S8rUEa0f3sx+AK1Fl+gpfmXnwzNdc2Qgy9QhkL16WrODziPd7wX\nIuFU1CnWb1kPn/at4x3n5fjZ4759CTTc1sCPin/E2dFn4e/Be7eXwTcM5g73HdgA87Uv0XcjxxeK\nX/AtB77Ryfvwh0F/CL6nS0tLuZSO5l1u2sXzA8872053dLW8nravs+edtbunem14HDhwgMTEROLj\n44mKimLevHkUFxe3X/Bm8P4vL6+WvNrlsv/HU/+D5g+bfR3MKF8ZzR8289VVX2237IABA9osN2DA\nAABe2vUS75e+32be+79+n/XPrL8wLQb+7+/+b5tl1j+z3vemHYWvsxzeweMY/2OPbx2m+3/P8v++\np9XvGHwdwtBWj4fhO6IMbGcovk5uoL+Mm/B1OIHnzfg6pZv8z2/Ad7R+E77OsHW5gQ7I8HWgA1s9\nH4Svo7oZX1AF6vEJfx1at3WIv303AXdyoTMNdMwz/evM4MIR+CB//U775zf5nw/wl3mP//f7/m3f\niK+jPQpM9rfzD/7ns4Bj/p9kfAHzvr8tUfhGFsOAz/t/1+Eb0YCvk67E16kHjuyHARXA/wMO4guh\ng/hObe0HtgG/xDeSGesvJxD+b/rr9wd84ZLIhZFUsr8dY/CFUGBf4X9+h789+Kd/Cij1T3f56/V7\n/3qD/NMi/NP/iC+sovz7JhA+t3NhVDkMXzgO9/+tPvS3Owmab2z2tTOw7dZ1c8GZ8WfaBOcfBv2B\n9T9b7ztgcfnLcsHHEz5uG7CJcK7pXPCLtwqPzqeFKzx67WmrHTt2sGfPHv793/8dgBdffJGysjLW\nr18fXEb3GxIR6Z6edv2RIapHyHUlGHpp7omIXPN67WmruLg4Kisrg88rKytxu91hrJGIiAT02vC4\n8847KS8v58MPP6SxsZGtW7cye/bscFdLREToxaetIiMj+bd/+zdmzJhBS0sLS5cuZdy4ceGuloiI\n0Is/MBcRkd6r1562EhGR3uuaCo8zZ86Ql5fHww8/zM9+9rNwVyfkKioqWLZsGQ888EC4q3JFFBcX\n8/DDDzNv3jz27dsX7uqE3J/+9Cfy8/OZO3cuGzduDHd1Qu7MmTOkpKTwy1+G7/+nXymlpaVMnTqV\n/Px8fvOb34S7OiFnZnzzm9/ky1/+Mps3b+7SOtdUeLz88svMnTuXH//4x+zatSvc1Qm52267jZ/8\n5CfhrsYVM2fOHH784x/zox/9iK1bt4a7OiE3duxYfvjDH1JUVMSePXvCXZ2Q++53v0tubm64q3FF\nREREMHjwYBoaGq7Jqz537txJdXU1/fv373L7en14LFmyhNjYWCZOnNhmeklJCWPHjsXj8fDMM88A\nUF1dzejRowF894/qA5y0ry/qTvu+/e1v88gjj1zNanab0/a9+uqrzJo1i3nz5l3tqjrmpG379u0j\nOTmZESNGhKOq3eKkfVOnTuVXv/oVa9euZdWqVeGormNO2vfee+/x2c9+lu9973v88Ic/7NoGenxf\n3ivs9ddft7ffftsmTJgQnNbc3GwJCQlWUVFhjY2NNmnSJDty5Iht2bLFfvGLX5iZ2bx588JVZUec\ntC/g/vvvD0dVu8VJ+7xer33ta1+z1157LYw1dqY7fz8zs9mzZ1/tqjrmpG3f/OY3bcWKFZadnW1z\n5szpE/93ozt/u4aGhj7z/nPSvhdffNG2bdtmZmZz587tUvm99lLdgKlTp/Lhhx+2mdb6vldA8L5X\nX/7yl3nkkUf45S9/2We+E+KkfbGxsXzjG9/g0KFDPPPMM+3uMtwbOWnfa6+9xv79+/nb3/7Gn//8\nZ77whS9c/Qo75KR9//3f/83LL7/M+fPnmT59+tWvrENO2vbtb38bgMLCQkaMGNEnbh3kpH1/+tOf\n2LNnDydPnuTRRx+9+pXtBiftW758OY8++ihvvPFGl29F3+vDoyOtT08BuN1uysrKGDhwIJs2Xfl7\n/V9pl2rf8OHD+dGPfhTGmoXGpdq3fv36PvPG7Myl2jdt2jSmTZsWxpr13KXaFpCXlxeOaoXMpdr3\n9a9/nfvuuy+MNQuNS7XvxhtvdPx5aq//zKMjfeGopifUvr7tWm7ftdw2UPuc6JPhca3f90rt69uu\n5fZdy20Dtc+JPhke1/p9r9S+vu1abt+13DZQ+xy5Ih/zh9C8efPslltusf79+5vb7bZNmzaZmdmv\nfvUrS0pKsoSEBHv66afDXMvuU/vUvt7qWm6bmdrX0/bp3lYiIuJYnzxtJSIi4aXwEBERxxQeIiLi\nmMJDREQcU3iIiIhjCg8REXFM4SEiIo4pPERExDGFh4iIOPb/Acf6YHV/dFHoAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 79 }, { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "Annotation of genetable_clc in SQLShare" ] }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }