long_read_rna_seq_analysis.ipynb

{
  "cells": [
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "P5k_WL74Y2Ei"
      },
      "source": [
        "#$\\text{L-RAPiT: Long Read Analysis Pipeline for Transcriptomics}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "IaygDNNEZnph"
      },
      "source": [
        "$\\text{A pipeline to analyze Oxford Nanopore and PacBio third-generation long transcriptomic sequencing reads}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "hIjUpSQCZlw1"
      },
      "source": [
        "$\\text{Theodore Nelson}$\n",
        "\n",
        "$\\text{Columbia University Irving Medical Center}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_bUDAXz-7KYw"
      },
      "source": [
        "![flowchart.svg](data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0iMS4wIiBlbmNvZGluZz0iVVRGLTgiPz4KPCEtLSBEbyBub3QgZWRpdCB0aGlzIGZpbGUgd2l0aCBlZGl0b3JzIG90aGVyIHRoYW4gZGlhZ3JhbXMubmV0IC0tPgo8IURPQ1RZUEUgc3ZnIFBVQkxJQyAiLS8vVzNDLy9EVEQgU1ZHIDEuMS8vRU4iICJodHRwOi8vd3d3LnczLm9yZy9HcmFwaGljcy9TVkcvMS4xL0RURC9zdmcxMS5kdGQiPgo8c3ZnIHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8yMDAwL3N2ZyIgc3R5bGU9ImJhY2tncm91bmQtY29sb3I6IHJnYigyNTUsIDI1NSwgMjU1KTsiIHhtbG5zOnhsaW5rPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rIiB2ZXJzaW9uPSIxLjEiIHdpZHRoPSI1MjFweCIgaGVpZ2h0PSI4ODFweCIgdmlld0JveD0iLTAuNSAtMC41IDUyMSA4ODEiIGNvbnRlbnQ9IiZsdDtteGZpbGUgaG9zdD0mcXVvdDthcHAuZGlhZ3JhbXMubmV0JnF1b3Q7IG1vZGlmaWVkPSZxdW90OzIwMjItMTItMDdUMDU6NDY6MDMuMjE4WiZxdW90OyBhZ2VudD0mcXVvdDs1LjAgKE1hY2ludG9zaCkmcXVvdDsgZXRhZz0mcXVvdDthTU00Q1ZDWTZTMmlhbHoteTBQUiZxdW90OyB2ZXJzaW9uPSZxdW90OzIwLjYuMCZxdW90OyB0eXBlPSZxdW90O2dvb2dsZSZxdW90OyZndDsmbHQ7ZGlhZ3JhbSBpZD0mcXVvdDtDNVJCczQzb0RhLUtkelplTnR1eSZxdW90OyBuYW1lPSZxdW90O1BhZ2UtMSZxdW90OyZndDs3VjFkZTZJNEZQNDFYdHFIVDhYTDFtcTN1MU9mZG0wN08xZnpwQkExVXlRT3hMYk9yOThFUVJHd2ZCUUN0dW1GbFJDQ25Mem5QU2ZKeWFHakRwZHZWeTVZTFc2d0JlMk9JbGx2SGZXeW95aXliUFRvUDFheUNVcU12cjR0bWJ2SUNzcjJCVlAwQndhRlVsQzZSaGIwRGlvU2pHMkNWb2VGSm5ZY2FKS0RNdUM2K1BXdzJnemJoM2RkZ1RsTUZFeE5ZQ2RMdnlPTExMYWxodExmbC84RjBYd1IzbG51RGJabmxpQ3NIRHlKdHdBV2ZvMFVxYU9PT25ReEp0dHZ5N2NodEpuMFFybDh2OTU4dDc4OTk2Nyt2dk4rZzRlTGYrNG5qOTF0WStNaWwrd2V3WVVPS2QzMDdiMzZNUG45ZlA5dzkwdFQvek5kZlhMdWRCVTFlRGF5Q1FVR0xTcS80TkRCRHYxMzRlSzFZMEhXamtTUHNFc1dlSTRkWUgvRGVFVUxaVnI0Q3hLeUNYb2ZyQW1tUlF1eXRJT3o4QTJSLzlqbFozcHc5Q055NXZJdGFOay8ySVFIRG5FM2tZdlk0WS9vdWYxbC9sRjQzUXc3WkF5V3lHWUY1eTRDZGxBNHhXdlhaTDl2UVFnRm42S3I1L1NEeW90OXNBcmUyUnpqdVEzQkNubG5KbDc2SjB6UHJ6cWViWnVrWDNlTjV1eVdvUHU4OFBiMGNQSzBNYldOL0RnQkw1YzNQMyt0aWRGSDNRQ1JCTGh6U042cEYyb1c2NmdJeG9OZXY0SjRDYWs4YUFVWDJvQ2dsME5GQUlFK3pYZjFkcGZlWWtRZlJKRUM1VmQ3d1kwQzFlLzJKZW13amUwakJaZEZvUmR2U1pmUGVqMXA5eWZIMnUzSDJ0MktJTkV1L1JKNXZuMlJEL1FDb0ZmN2ZFQWZBbGcrZ084ZXpTY0c0THJ4cGlreHZHbGw4WmJkVW5VSUE3TUw3Vy9EVmZHak9mVW15Ky9qSDZ1ZjNZQlZYNEM5RHFTUlFCeTFKaXYyMVZ5Nzl1YkNCZVl6VS9xTDF3VWljTG9DZmorL1VwdDhpS2s5S0xkSGhBb2NPL1J3SUJYdXp4Zm9FdmoyYm8vdTlEY21UMWtaQkNXdmU1c2FDbjBSTWFlNkloMkh3WUhVaTRwWXl4WXhkS3h6NWtFd0tkdkE4NUI1VEppSDlxYVl1bGFzUkJHQjZta0NEY28rcUd0S1hFTUd2YkxjTG1XMVZMT3U5V29EUWtIZTNqazVFUmRuNy9Da096a25DUjQ5czh0TEV6VnY4UFFyQjA4QUFybVQxODh0Qm9GTTc3RFhLcWlvZzNnSGExVkJ4VkQ1UWtXcERTdFNicXdjb2FaS0xWUW13dlIySVV5dENtR3FiSnpwZlZsVGpPMW5yeDhEM0VBN0M4K3hUNVVyL0RROUczNFYrNVJkdms2bHp0T25URVcybENMaW5zMUVhS0dYQTFIM2ZxL1p0STgvTE90Ni9naVFqc1VrV1Z1OStTSUx6OU52OCtDL0hkYW45U3hBUUhkN3JXc2VOTGdkMnJIRzZMRHV5S0RPSDlMdEJuVGhjQzV5ejFJL1ZKb0I4L0NhUkx1cHoxTHFYalVMNFBDSFBvVUZGd2liMkxGQWVJSzI5QlN2SE53MGY3RVBqWk9SVE81ZXZuWThBbnlWdVhYeDNBVkw3NmdvWWl4MHlDY1o1RE5EdGozRU5uYjlhOVhSZUt5UFIzV3lUdHgzbFdVbFNUdHlHdThZZGRHT0xHZ25CeUFmVmpZR1ZvZjFUS0RXN29lVVZNaDIveXovUW1CNXpQZkRRcjVmMm1SZStaZlJ1MTI2ekNWL3Y5OUxBNlJDZ3pFZWpYdVhSbzBHUSt2SERZWWtOVzB3RkY0R280MUFMYXgwT3hmdkgrVE1aOGhiUVBjOUI3QmR0SGZTUFhDSlg1MVdHKzJUbHU3dCtzbEdKdU5xK29QS1VQWEpjckljbi9OcEFTZW5yVWQ5SGdmbTJsbXQyVU1Fam1MOUE4R3gvMWNqaHBJRHdlWXhsRGJGSjNndjA2NVBGOEI5RmlhZGkvREh5Q2JNZjhwZ2dpYXNqakc2cUpVeEVsWm4wRGhocEsxOUM4eG1Fc1lZZU9SdW1NNFl3bE90dmdNWVY5RHFkMnRnSTdKcEUybU1SdVBSNktKRzBraE1IelJQR21ucjJBS3ptYVN4UkE1YWdwWFNFWTRHRC9tZjIyanU1UFF6NmlRSTdndFNodEUwUVJncEJORlVvRXZCNFA5aXZaUVo2aEpRWlZ0Q1hlSUIrWU40SkVUK1VKZGpvS3MrbWlWVnNJTzZJQmFOdThzS3Z1UVQ4WnU1NzBNL1NURHU2R2NIb2JKQm9NbVc0ckN1R1l4cDRJdnNCOWx2L1JqdFN5dmNFMVVBc2Z5Mmh2aU5ibis2clBrYXNOUEdZSDhNTFJrakptai9sM20wajBpc2hsOFdxVk94M3ZEYUI1WEFaOTdkSkZRY1lCT3B0bUlWdkFJMzZxdFNNYzNSWlNtbUlOdmZVS202aEhJWEE0aGlBd2dQTE5uTzE2Z3pLK1lkdUhUQkI1NWwyOUFRTDFjdVpHd3ArZU9TSmZRVnVmWFJFZnhITDd1OW5NMUYwNG5vaUZMMFJPOUIwQkZ5RXN4VXgrd0dZeEZCL0hXSVY4dzB0NE9LandkRkNMUytROFgzTG5BODAwVXJNcVIzZElURHlDMmtDdHZJV3dqWjFtcnkyaFUxd0QxV1RkTWFwK1cwQ1g0QjJPd0IvQXQxa1RVMlNTc0ltVk1mK0ZuSmhHaHJFUzEwMkVhUElhWmN4bktvdFltVWpaRXg2dk1rNWVaOVpTVkh6cUJpUzJKTkxoYVV5NkJTckx1ekU2VVpPUmNJR2xveDY4Y3pXT1JldmszQWwvUHlyU0lHZHFVOGlCdndEUDlhUHduL2dhZVJjd0ZoaG00L2poYnlyWFVlODRnd2c4Zk4yL0lqZ3E4SEcvVytvRmZTNnpYdWxZZ3RLYVdJZmtwYzVNenZVZFN0ZnNxRlhzRkRaVHFBcXFrTlZoNXNwMmpGRk5tN3ZKZFlSRzRCNzFXZTdPM0xqOGJDTVVOYlIyTzlzc0cwU2p3c1Y2OHYyM0M2Wk1WOGJpa2pmVFdic2NBZUNuQmhwZmxZYVJkN1h2ZGZPSU9VVTh5c0tjZDY1Qmw2Q2xZbkxaYnJreHRXWFcvYXNLcHB5YzVPMWJCeWZLdEZNWWhreGs0cmVUY2RLTzJhUTAza0UwNGtBODJkVDFpTHQxVGZXeWpTTlNGdHdsL1lqRXlyUFlPQXJGMDRwQ1JBUkN3MXY5QVlGMXRyWnJPRmNHdGNpMTFuK1NKNTVIdEtMa29pYXJINXNiOHFGcmRLRWZPMzY3dHp3Y2VDanorQmNQY0J1SUtWMnpKd0ZHbExTckh5K052NXpValFNajlhaGhZeTJXTk04SXYvTHRYcHlrWW1jdWIwNitnbGM2NkpONVBVblRRdHdTVE5oeityYWNrbkJKS3pnNWZXTmtGSDBxYjVUMU1sbmFUSzg2UlRxTEtsQm1SWDdiQjlQamxOMmZ1SFBJTE1Kck0rMWIxWkx6RXpyL0o4QzhuOW55Y2RqbWN2cjdOWDYrZlQ1dkh5VHUySEx6OXFKZ2xLa2JROVRhNmtmOWtrS0ZyaVpkQmxnNWNURSs5MXBWTkozQ2dyblVweVNhRGlkQ3FwZXFla2VDT251aURHSmFGV2R0eS9rbE1uK0t4WkhjWGh4eU5OK3ZWbHlrckY2bWZjby9JVlh2ZWNnR0JQWVU5YWpyOFRNL2JKdG1xR1lkb0V2WURoQ2NLUUd0aXpjT05iQ1J5ZVNaRS9PYnZwbW1HWjR3MnRBcGFuQUV2cWtkWUZ5N1NtcTRQbDlmbGpkL0FIZVdpMCtEMnhMbjg4VEl4YmZtKzEvVnpUWGVPMWg3QXpoZkJadkMrTVV4OWNXNVEyMEd3alJGdFhoTWtXMU8xYWd1Q2VnWVhyWW1ZcUpYK21JRmd1dTB0YTRpa29SdExMTE9jbUtIcFdTelY3Qlo5bzJxa1lsTEtuajFxK1VTa3hYMW9hYzlwZ3dCVnphaEppa1NXR1lNSzh1c1RxSE1MeStZRytHTWh6cnh0a2J3RFFjMnJEUnlQN0Ura3JqVUU1bUNkYjBxdUNPVDEwTWZPcjl0VlpncW9iYkVGVzQzOD0mbHQ7L2RpYWdyYW0mZ3Q7Jmx0Oy9teGZpbGUmZ3Q7Ij48ZGVmcz48c3R5bGUgdHlwZT0idGV4dC9jc3MiPkBpbXBvcnQgdXJsKGh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCk7JiN4YTs8L3N0eWxlPjwvZGVmcz48Zz48cGF0aCBkPSJNIDI2MCAzMjAgTCAyNjAgMzQzLjYzIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAyNjAgMzQ4Ljg4IEwgMjU2LjUgMzQxLjg4IEwgMjYwIDM0My42MyBMIDI2My41IDM0MS44OCBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDI3MCAtMTcwIEwgMjY1IC0xNzAgUSAyNjAgLTE3MCAyNjAgLTE2MCBMIDI2MCA4MCBRIDI2MCA5MCAyNTUgOTAgTCAyNTIuNSA5MCBRIDI1MCA5MCAyNTUgOTAgTCAyNTcuNSA5MCBRIDI2MCA5MCAyNjAgMTAwIEwgMjYwIDM0MCBRIDI2MCAzNTAgMjY1IDM1MCBMIDI3MCAzNTAiIGZpbGw9Im5vbmUiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHRyYW5zZm9ybT0icm90YXRlKDkwLDI2MCw5MCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDEyMCAxNjAgTCAxOTMuNjMgMTYwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAxOTguODggMTYwIEwgMTkxLjg4IDE2My41IEwgMTkzLjYzIDE2MCBMIDE5MS44OCAxNTYuNSBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDQwMCAxNjAgTCAzMjYuMzcgMTYwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAzMjEuMTIgMTYwIEwgMzI4LjEyIDE1Ni41IEwgMzI2LjM3IDE2MCBMIDMyOC4xMiAxNjMuNSBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDI5MCAxODAgTCAzNjQuNDcgMjIyLjU1IiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAzNjkuMDMgMjI1LjE2IEwgMzYxLjIyIDIyNC43MyBMIDM2NC40NyAyMjIuNTUgTCAzNjQuNjkgMjE4LjY1IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxwYXRoIGQ9Ik0gMjMwIDE4MCBMIDE1NS43IDIxNy4xNSIgZmlsbD0ibm9uZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9InN0cm9rZSIvPjxwYXRoIGQ9Ik0gMTUxIDIxOS41IEwgMTU1LjcgMjEzLjI0IEwgMTU1LjcgMjE3LjE1IEwgMTU4LjgzIDIxOS41IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxwYXRoIGQ9Ik0gMjcwIDUzMCBMIDI2NSA1MzAgUSAyNjAgNTMwIDI2MCA1NDAgTCAyNjAgNzgwIFEgMjYwIDc5MCAyNTUgNzkwIEwgMjUyLjUgNzkwIFEgMjUwIDc5MCAyNTUgNzkwIEwgMjU3LjUgNzkwIFEgMjYwIDc5MCAyNjAgODAwIEwgMjYwIDEwNDAgUSAyNjAgMTA1MCAyNjUgMTA1MCBMIDI3MCAxMDUwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiB0cmFuc2Zvcm09InJvdGF0ZSgtOTAsMjYwLDc5MCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cmVjdCB4PSIyMDAiIHk9IjAiIHdpZHRoPSIxMjAiIGhlaWdodD0iODAiIHJ4PSIxMiIgcnk9IjEyIiBmaWxsPSIjZWZmNWZlIiBzdHJva2U9InJnYigwLCAwLCAwKSIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxnIHRyYW5zZm9ybT0idHJhbnNsYXRlKC0wLjUgLTAuNSkiPjxzd2l0Y2g+PGZvcmVpZ25PYmplY3Qgc3R5bGU9Im92ZXJmbG93OiB2aXNpYmxlOyB0ZXh0LWFsaWduOiBsZWZ0OyIgcG9pbnRlci1ldmVudHM9Im5vbmUiIHdpZHRoPSIxMDAlIiBoZWlnaHQ9IjEwMCUiIHJlcXVpcmVkRmVhdHVyZXM9Imh0dHA6Ly93d3cudzMub3JnL1RSL1NWRzExL2ZlYXR1cmUjRXh0ZW5zaWJpbGl0eSI+PGRpdiB4bWxucz0iaHR0cDovL3d3dy53My5vcmcvMTk5OS94aHRtbCIgc3R5bGU9ImRpc3BsYXk6IGZsZXg7IGFsaWduLWl0ZW1zOiB1bnNhZmUgY2VudGVyOyBqdXN0aWZ5LWNvbnRlbnQ6IHVuc2FmZSBjZW50ZXI7IHdpZHRoOiAxMThweDsgaGVpZ2h0OiAxcHg7IHBhZGRpbmctdG9wOiA0MHB4OyBtYXJnaW4tbGVmdDogMjAxcHg7Ij48ZGl2IHN0eWxlPSJib3gtc2l6aW5nOiBib3JkZXItYm94OyBmb250LXNpemU6IDBweDsgdGV4dC1hbGlnbjogY2VudGVyOyIgZGF0YS1kcmF3aW8tY29sb3JzPSJjb2xvcjogcmdiKDAsIDAsIDApOyAiPjxkaXYgc3R5bGU9ImRpc3BsYXk6IGlubGluZS1ibG9jazsgZm9udC1zaXplOiAxMnB4OyBmb250LWZhbWlseTogSGVsdmV0aWNhOyBjb2xvcjogcmdiKDAsIDAsIDApOyBsaW5lLWhlaWdodDogMS4yOyBwb2ludGVyLWV2ZW50czogYWxsOyB3aGl0ZS1zcGFjZTogbm9ybWFsOyBvdmVyZmxvdy13cmFwOiBub3JtYWw7Ij48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBmYWNlPSJBcmlhbCI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGI+QmlvY29uZGE8L2I+PC9mb250PjwvZm9udD48L2Rpdj48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+SW5zdGFsbCBQcm9ncmFtczwvZm9udD48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iMjYwIiB5PSI0NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPkJpb2NvbmRhSW5zdGFsbCBQcm9nLi4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cmVjdCB4PSIzNzAiIHk9IjEyMCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmZWY2ZDgiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDE2MHB4OyBtYXJnaW4tbGVmdDogMzcxcHg7Ij48ZGl2IHN0eWxlPSJib3gtc2l6aW5nOiBib3JkZXItYm94OyBmb250LXNpemU6IDBweDsgdGV4dC1hbGlnbjogY2VudGVyOyIgZGF0YS1kcmF3aW8tY29sb3JzPSJjb2xvcjogcmdiKDAsIDAsIDApOyAiPjxkaXYgc3R5bGU9ImRpc3BsYXk6IGlubGluZS1ibG9jazsgZm9udC1zaXplOiAxMnB4OyBmb250LWZhbWlseTogSGVsdmV0aWNhOyBjb2xvcjogcmdiKDAsIDAsIDApOyBsaW5lLWhlaWdodDogMS4yOyBwb2ludGVyLWV2ZW50czogYWxsOyB3aGl0ZS1zcGFjZTogbm9ybWFsOyBvdmVyZmxvdy13cmFwOiBub3JtYWw7Ij48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBmYWNlPSJBcmlhbCI+VXBsb2FkIDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBmYWNlPSJBcmlhbCI+UmVhZHMgdG8gPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIGZhY2U9IkFyaWFsIj48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij5Hb29nbGUgRHJpdmU8L2ZvbnQ+PGJyIC8+PC9mb250PjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSI0MzAiIHk9IjE2NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPlVwbG9hZC4uLjwvdGV4dD48L3N3aXRjaD48L2c+PHJlY3QgeD0iMzAiIHk9IjEyMCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmZWY2ZDgiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDE2MHB4OyBtYXJnaW4tbGVmdDogMzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5LaW5nZmlzaGVyPC9iPjxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+RG93bmxvYWQgPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5QdWJsaWMgRGF0YTwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iOTAiIHk9IjE2NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPktpbmdmaXNoZXIuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjIwMCIgeT0iMTIwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2ZmZmZmZiIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogMTYwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxmb250IHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIGZhY2U9IkFyaWFsIj5JbnB1dCBSZWFkczwvZm9udD48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iMjYwIiB5PSIxNjQiIGZpbGw9InJnYigwLCAwLCAwKSIgZm9udC1mYW1pbHk9IkhlbHZldGljYSIgZm9udC1zaXplPSIxMnB4IiB0ZXh0LWFuY2hvcj0ibWlkZGxlIj5JbnB1dCBSZWFkczwvdGV4dD48L3N3aXRjaD48L2c+PHJlY3QgeD0iMzAiIHk9IjIyMCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmOGViZmYiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDI2MHB4OyBtYXJnaW4tbGVmdDogMzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5TaGFyazwvYj48YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPkZpbHRlciBSZWFkczwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iOTAiIHk9IjI2NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPlNoYXJrLi4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cmVjdCB4PSIzNzAiIHk9IjIyMCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNlZWZlZWIiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDI2MHB4OyBtYXJnaW4tbGVmdDogMzcxcHg7Ij48ZGl2IHN0eWxlPSJib3gtc2l6aW5nOiBib3JkZXItYm94OyBmb250LXNpemU6IDBweDsgdGV4dC1hbGlnbjogY2VudGVyOyIgZGF0YS1kcmF3aW8tY29sb3JzPSJjb2xvcjogcmdiKDAsIDAsIDApOyAiPjxkaXYgc3R5bGU9ImRpc3BsYXk6IGlubGluZS1ibG9jazsgZm9udC1zaXplOiAxMnB4OyBmb250LWZhbWlseTogSGVsdmV0aWNhOyBjb2xvcjogcmdiKDAsIDAsIDApOyBsaW5lLWhlaWdodDogMS4yOyBwb2ludGVyLWV2ZW50czogYWxsOyB3aGl0ZS1zcGFjZTogbm9ybWFsOyBvdmVyZmxvdy13cmFwOiBub3JtYWw7Ij48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+PGI+RmFzdFFDPC9iPiA8YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPlJlYWQgUXVhbGl0eTwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iNDMwIiB5PSIyNjQiIGZpbGw9InJnYigwLCAwLCAwKSIgZm9udC1mYW1pbHk9IkhlbHZldGljYSIgZm9udC1zaXplPSIxMnB4IiB0ZXh0LWFuY2hvcj0ibWlkZGxlIj5GYXN0UUMuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjIwMCIgeT0iMjQwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2VmZjVmZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogMjgwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5taW5pbWFwMjwvYj48YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPkFsaWduIFJlYWRzPC9mb250PjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSIyNjAiIHk9IjI4NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPm1pbmltYXAyLi4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cGF0aCBkPSJNIDI2MCAyMDAgTCAyNjAgMjMzLjYzIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAyNjAgMjM4Ljg4IEwgMjU2LjUgMjMxLjg4IEwgMjYwIDIzMy42MyBMIDI2My41IDIzMS44OCBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDE1MCAyNjAgTCAxOTQuMDkgMjc3LjYzIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAxOTguOTYgMjc5LjU4IEwgMTkxLjE2IDI4MC4yMyBMIDE5NC4wOSAyNzcuNjMgTCAxOTMuNzYgMjczLjc0IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxwYXRoIGQ9Ik0gMjAwIDM5MCBMIDE4MCAzOTAgTCAxODAgNjEwIEwgMTgwIDcyMCIgZmlsbD0ibm9uZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9InN0cm9rZSIvPjxyZWN0IHg9IjIwMCIgeT0iMzUwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2VmZjVmZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogMzkwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5zYW10b29sczwvYj4gPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Zm9udCBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+Q29tcHJlc3MgQWxpZ25tZW50czwvZm9udD48YnIgLz48L2ZvbnQ+PC9kaXY+PC9kaXY+PC9kaXY+PC9kaXY+PC9mb3JlaWduT2JqZWN0Pjx0ZXh0IHg9IjI2MCIgeT0iMzk0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+c2FtdG9vbHMuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjM3MCIgeT0iMzUwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2VlZmVlYiIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogMzkwcHg7IG1hcmdpbi1sZWZ0OiAzNzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5waXN0aXM8L2I+PGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5BbGlnbm1lbnQgPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5RdWFsaXR5PC9mb250PjxiciAvPjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSI0MzAiIHk9IjM5NCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPnBpc3Rpcy4uLjwvdGV4dD48L3N3aXRjaD48L2c+PHJlY3QgeD0iMzAiIHk9IjY4MCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmZWY2ZDgiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDcyMHB4OyBtYXJnaW4tbGVmdDogMzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5UcmFuc2NyaXB0Q2xlYW48L2I+IDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+UG9saXNoIDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+QWxpZ25tZW50czwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iOTAiIHk9IjcyNCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPlRyYW5zY3JpcHRDbGVhbi4uLjwvdGV4dD48L3N3aXRjaD48L2c+PHJlY3QgeD0iMzAiIHk9IjM1MCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmOGU4ZTciIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDM5MHB4OyBtYXJnaW4tbGVmdDogMzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5zdmlzdDRnZXQ8L2I+IDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+R3JhcGggPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5HZW5lIENvdmVyYWdlPC9mb250PjxiciAvPjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSI5MCIgeT0iMzk0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+c3Zpc3Q0Z2V0Li4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cGF0aCBkPSJNIDE4MCAzOTAgTCAxNTYuMzcgMzkwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAxNTEuMTIgMzkwIEwgMTU4LjEyIDM4Ni41IEwgMTU2LjM3IDM5MCBMIDE1OC4xMiAzOTMuNSBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cmVjdCB4PSIzMCIgeT0iNDYwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2Y4ZThlNyIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogNTAwcHg7IG1hcmdpbi1sZWZ0OiAzMXB4OyI+PGRpdiBzdHlsZT0iYm94LXNpemluZzogYm9yZGVyLWJveDsgZm9udC1zaXplOiAwcHg7IHRleHQtYWxpZ246IGNlbnRlcjsiIGRhdGEtZHJhd2lvLWNvbG9ycz0iY29sb3I6IHJnYigwLCAwLCAwKTsgIj48ZGl2IHN0eWxlPSJkaXNwbGF5OiBpbmxpbmUtYmxvY2s7IGZvbnQtc2l6ZTogMTJweDsgZm9udC1mYW1pbHk6IEhlbHZldGljYTsgY29sb3I6IHJnYigwLCAwLCAwKTsgbGluZS1oZWlnaHQ6IDEuMjsgcG9pbnRlci1ldmVudHM6IGFsbDsgd2hpdGUtc3BhY2U6IG5vcm1hbDsgb3ZlcmZsb3ctd3JhcDogbm9ybWFsOyI+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPjxiPk1ha2VIdWI8L2I+IDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+R2VuZXJhdGUgYW48YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPkFsaWdubWVudCA8L2ZvbnQ+PGZvbnQgZmFjZT0iQXJpYWwiPlZpZXdlcjwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iOTAiIHk9IjUwNCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPk1ha2VIdWIuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjIwMCIgeT0iNDYwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2ZlZjZkOCIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogNTAwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5TdHJpbmdUaWU8L2I+PGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5Db2xsYXBzZTxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBmYWNlPSJBcmlhbCI+QWxpZ25tZW50czwvZm9udD48YnIgLz48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iMjYwIiB5PSI1MDQiIGZpbGw9InJnYigwLCAwLCAwKSIgZm9udC1mYW1pbHk9IkhlbHZldGljYSIgZm9udC1zaXplPSIxMnB4IiB0ZXh0LWFuY2hvcj0ibWlkZGxlIj5TdHJpbmdUaWUuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxwYXRoIGQ9Ik0gMTgwIDUwMCBMIDE1Ni4zNyA1MDAiIGZpbGw9Im5vbmUiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJzdHJva2UiLz48cGF0aCBkPSJNIDE1MS4xMiA1MDAgTCAxNTguMTIgNDk2LjUgTCAxNTYuMzcgNTAwIEwgMTU4LjEyIDUwMy41IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxyZWN0IHg9IjIwMCIgeT0iNTcwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2ZlZjZkOCIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogNjEwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5HZmZDb21wYXJlPC9iPjxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+Q3Jvc3MtUmVmZXJlbmNlPC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPkNvbGxhcHNlZCBBbGlnbm1lbnRzPGJyIC8+PC9kaXY+PC9kaXY+PC9kaXY+PC9kaXY+PC9mb3JlaWduT2JqZWN0Pjx0ZXh0IHg9IjI2MCIgeT0iNjE0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+R2ZmQ29tcGFyZS4uLjwvdGV4dD48L3N3aXRjaD48L2c+PHBhdGggZD0iTSAyNjAgNTQwIEwgMjYwIDU2My42MyIgZmlsbD0ibm9uZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9InN0cm9rZSIvPjxwYXRoIGQ9Ik0gMjYwIDU2OC44OCBMIDI1Ni41IDU2MS44OCBMIDI2MCA1NjMuNjMgTCAyNjMuNSA1NjEuODggWiIgZmlsbD0icmdiKDAsIDAsIDApIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PHJlY3QgeD0iMzcwIiB5PSI0NjAiIHdpZHRoPSIxMjAiIGhlaWdodD0iODAiIHJ4PSIxMiIgcnk9IjEyIiBmaWxsPSIjZmVmNmQ4IiBzdHJva2U9InJnYigwLCAwLCAwKSIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxnIHRyYW5zZm9ybT0idHJhbnNsYXRlKC0wLjUgLTAuNSkiPjxzd2l0Y2g+PGZvcmVpZ25PYmplY3Qgc3R5bGU9Im92ZXJmbG93OiB2aXNpYmxlOyB0ZXh0LWFsaWduOiBsZWZ0OyIgcG9pbnRlci1ldmVudHM9Im5vbmUiIHdpZHRoPSIxMDAlIiBoZWlnaHQ9IjEwMCUiIHJlcXVpcmVkRmVhdHVyZXM9Imh0dHA6Ly93d3cudzMub3JnL1RSL1NWRzExL2ZlYXR1cmUjRXh0ZW5zaWJpbGl0eSI+PGRpdiB4bWxucz0iaHR0cDovL3d3dy53My5vcmcvMTk5OS94aHRtbCIgc3R5bGU9ImRpc3BsYXk6IGZsZXg7IGFsaWduLWl0ZW1zOiB1bnNhZmUgY2VudGVyOyBqdXN0aWZ5LWNvbnRlbnQ6IHVuc2FmZSBjZW50ZXI7IHdpZHRoOiAxMThweDsgaGVpZ2h0OiAxcHg7IHBhZGRpbmctdG9wOiA1MDBweDsgbWFyZ2luLWxlZnQ6IDM3MXB4OyI+PGRpdiBzdHlsZT0iYm94LXNpemluZzogYm9yZGVyLWJveDsgZm9udC1zaXplOiAwcHg7IHRleHQtYWxpZ246IGNlbnRlcjsiIGRhdGEtZHJhd2lvLWNvbG9ycz0iY29sb3I6IHJnYigwLCAwLCAwKTsgIj48ZGl2IHN0eWxlPSJkaXNwbGF5OiBpbmxpbmUtYmxvY2s7IGZvbnQtc2l6ZTogMTJweDsgZm9udC1mYW1pbHk6IEhlbHZldGljYTsgY29sb3I6IHJnYigwLCAwLCAwKTsgbGluZS1oZWlnaHQ6IDEuMjsgcG9pbnRlci1ldmVudHM6IGFsbDsgd2hpdGUtc3BhY2U6IG5vcm1hbDsgb3ZlcmZsb3ctd3JhcDogbm9ybWFsOyI+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPjxiPmZlYXR1cmVDb3VudHM8L2I+IDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+UHJvZHVjZSA8YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPkdlbmUgQ291bnRzPGJyIC8+PC9mb250PjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSI0MzAiIHk9IjUwNCIgZmlsbD0icmdiKDAsIDAsIDApIiBmb250LWZhbWlseT0iSGVsdmV0aWNhIiBmb250LXNpemU9IjEycHgiIHRleHQtYW5jaG9yPSJtaWRkbGUiPmZlYXR1cmVDb3VudHMuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjM3MCIgeT0iNTcwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2ZlZjZkOCIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogNjEwcHg7IG1hcmdpbi1sZWZ0OiAzNzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5MSVFBPC9iPiA8YnIgLz48L2ZvbnQ+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyI+PGZvbnQgZGF0YS1mb250LXNyYz0iaHR0cHM6Ly9mb250cy5nb29nbGVhcGlzLmNvbS9jc3M/ZmFtaWx5PUFyaWFsIiBzdHlsZT0iZm9udC1zaXplOiAxNHB4OyIgZmFjZT0iQXJpYWwiPlByb2R1Y2UgPGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5UcmFuc2NyaXB0IENvdW50czxiciAvPjwvZm9udD48L2Rpdj48L2Rpdj48L2Rpdj48L2Rpdj48L2ZvcmVpZ25PYmplY3Q+PHRleHQgeD0iNDMwIiB5PSI2MTQiIGZpbGw9InJnYigwLCAwLCAwKSIgZm9udC1mYW1pbHk9IkhlbHZldGljYSIgZm9udC1zaXplPSIxMnB4IiB0ZXh0LWFuY2hvcj0ibWlkZGxlIj5MSVFBLi4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cmVjdCB4PSIzNzAiIHk9IjY4MCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmOGViZmYiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDcyMHB4OyBtYXJnaW4tbGVmdDogMzcxcHg7Ij48ZGl2IHN0eWxlPSJib3gtc2l6aW5nOiBib3JkZXItYm94OyBmb250LXNpemU6IDBweDsgdGV4dC1hbGlnbjogY2VudGVyOyIgZGF0YS1kcmF3aW8tY29sb3JzPSJjb2xvcjogcmdiKDAsIDAsIDApOyAiPjxkaXYgc3R5bGU9ImRpc3BsYXk6IGlubGluZS1ibG9jazsgZm9udC1zaXplOiAxMnB4OyBmb250LWZhbWlseTogSGVsdmV0aWNhOyBjb2xvcjogcmdiKDAsIDAsIDApOyBsaW5lLWhlaWdodDogMS4yOyBwb2ludGVyLWV2ZW50czogYWxsOyB3aGl0ZS1zcGFjZTogbm9ybWFsOyBvdmVyZmxvdy13cmFwOiBub3JtYWw7Ij48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+PGI+RkxBTUU8L2I+IDxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+UHJlZGljdCBOb3ZlbCBTcGxpY2luZyBFdmVudHM8YnIgLz48L2ZvbnQ+PC9kaXY+PC9kaXY+PC9kaXY+PC9kaXY+PC9mb3JlaWduT2JqZWN0Pjx0ZXh0IHg9IjQzMCIgeT0iNzI0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+RkxBTUUuLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxyZWN0IHg9IjIwMCIgeT0iODAwIiB3aWR0aD0iMTIwIiBoZWlnaHQ9IjgwIiByeD0iMTIiIHJ5PSIxMiIgZmlsbD0iI2VlZmVlYiIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48ZyB0cmFuc2Zvcm09InRyYW5zbGF0ZSgtMC41IC0wLjUpIj48c3dpdGNoPjxmb3JlaWduT2JqZWN0IHN0eWxlPSJvdmVyZmxvdzogdmlzaWJsZTsgdGV4dC1hbGlnbjogbGVmdDsiIHBvaW50ZXItZXZlbnRzPSJub25lIiB3aWR0aD0iMTAwJSIgaGVpZ2h0PSIxMDAlIiByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiPjxkaXYgeG1sbnM9Imh0dHA6Ly93d3cudzMub3JnLzE5OTkveGh0bWwiIHN0eWxlPSJkaXNwbGF5OiBmbGV4OyBhbGlnbi1pdGVtczogdW5zYWZlIGNlbnRlcjsganVzdGlmeS1jb250ZW50OiB1bnNhZmUgY2VudGVyOyB3aWR0aDogMTE4cHg7IGhlaWdodDogMXB4OyBwYWRkaW5nLXRvcDogODQwcHg7IG1hcmdpbi1sZWZ0OiAyMDFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5NdWx0aVFDPC9iPjogPGJyIC8+PC9mb250PjwvZGl2PjxkaXY+PGZvbnQgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgZmFjZT0iQXJpYWwiPkNvbXBpbGUgPGJyIC8+PC9mb250PjwvZGl2PjxkaXY+PGZvbnQgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgZmFjZT0iQXJpYWwiPlN0YXRpc3RpY3M8L2ZvbnQ+PC9kaXY+PC9kaXY+PC9kaXY+PC9kaXY+PC9mb3JlaWduT2JqZWN0Pjx0ZXh0IHg9IjI2MCIgeT0iODQ0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+TXVsdGlRQzouLi48L3RleHQ+PC9zd2l0Y2g+PC9nPjxwYXRoIGQ9Ik0gMzIwIDM5MCBMIDM0MCAzOTAgTCAzNDAgNjEwIEwgMzQwIDcyMCIgZmlsbD0ibm9uZSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9InN0cm9rZSIvPjxwYXRoIGQ9Ik0gMzQwIDM5MCBMIDM2My42MyAzOTAiIGZpbGw9Im5vbmUiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJzdHJva2UiLz48cGF0aCBkPSJNIDM2OC44OCAzOTAgTCAzNjEuODggMzkzLjUgTCAzNjMuNjMgMzkwIEwgMzYxLjg4IDM4Ni41IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxwYXRoIGQ9Ik0gMzQwIDQ5OS41IEwgMzYzLjYzIDQ5OS41IiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAzNjguODggNDk5LjUgTCAzNjEuODggNTAzIEwgMzYzLjYzIDQ5OS41IEwgMzYxLjg4IDQ5NiBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDM0MCA2MDkuNzcgTCAzNjMuNjMgNjA5Ljc3IiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAzNjguODggNjA5Ljc3IEwgMzYxLjg4IDYxMy4yNyBMIDM2My42MyA2MDkuNzcgTCAzNjEuODggNjA2LjI3IFoiIGZpbGw9InJnYigwLCAwLCAwKSIgc3Ryb2tlPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZS1taXRlcmxpbWl0PSIxMCIgcG9pbnRlci1ldmVudHM9ImFsbCIvPjxwYXRoIGQ9Ik0gMzQwIDcxOS43NyBMIDM2My42MyA3MTkuNzciIGZpbGw9Im5vbmUiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJzdHJva2UiLz48cGF0aCBkPSJNIDM2OC44OCA3MTkuNzcgTCAzNjEuODggNzIzLjI3IEwgMzYzLjYzIDcxOS43NyBMIDM2MS44OCA3MTYuMjcgWiIgZmlsbD0icmdiKDAsIDAsIDApIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PHJlY3QgeD0iMzAiIHk9IjU3MCIgd2lkdGg9IjEyMCIgaGVpZ2h0PSI4MCIgcng9IjEyIiByeT0iMTIiIGZpbGw9IiNmZWY2ZDgiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBwb2ludGVyLWV2ZW50cz0iYWxsIi8+PGcgdHJhbnNmb3JtPSJ0cmFuc2xhdGUoLTAuNSAtMC41KSI+PHN3aXRjaD48Zm9yZWlnbk9iamVjdCBzdHlsZT0ib3ZlcmZsb3c6IHZpc2libGU7IHRleHQtYWxpZ246IGxlZnQ7IiBwb2ludGVyLWV2ZW50cz0ibm9uZSIgd2lkdGg9IjEwMCUiIGhlaWdodD0iMTAwJSIgcmVxdWlyZWRGZWF0dXJlcz0iaHR0cDovL3d3dy53My5vcmcvVFIvU1ZHMTEvZmVhdHVyZSNFeHRlbnNpYmlsaXR5Ij48ZGl2IHhtbG5zPSJodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hodG1sIiBzdHlsZT0iZGlzcGxheTogZmxleDsgYWxpZ24taXRlbXM6IHVuc2FmZSBjZW50ZXI7IGp1c3RpZnktY29udGVudDogdW5zYWZlIGNlbnRlcjsgd2lkdGg6IDExOHB4OyBoZWlnaHQ6IDFweDsgcGFkZGluZy10b3A6IDYxMHB4OyBtYXJnaW4tbGVmdDogMzFweDsiPjxkaXYgc3R5bGU9ImJveC1zaXppbmc6IGJvcmRlci1ib3g7IGZvbnQtc2l6ZTogMHB4OyB0ZXh0LWFsaWduOiBjZW50ZXI7IiBkYXRhLWRyYXdpby1jb2xvcnM9ImNvbG9yOiByZ2IoMCwgMCwgMCk7ICI+PGRpdiBzdHlsZT0iZGlzcGxheTogaW5saW5lLWJsb2NrOyBmb250LXNpemU6IDEycHg7IGZvbnQtZmFtaWx5OiBIZWx2ZXRpY2E7IGNvbG9yOiByZ2IoMCwgMCwgMCk7IGxpbmUtaGVpZ2h0OiAxLjI7IHBvaW50ZXItZXZlbnRzOiBhbGw7IHdoaXRlLXNwYWNlOiBub3JtYWw7IG92ZXJmbG93LXdyYXA6IG5vcm1hbDsiPjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj48Yj5GdXNpb25TZWVrZXI8L2I+PGJyIC8+PC9mb250PjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiPjxmb250IGRhdGEtZm9udC1zcmM9Imh0dHBzOi8vZm9udHMuZ29vZ2xlYXBpcy5jb20vY3NzP2ZhbWlseT1BcmlhbCIgc3R5bGU9ImZvbnQtc2l6ZTogMTRweDsiIGZhY2U9IkFyaWFsIj5JZGVudGlmeTxiciAvPjwvZm9udD48L2Rpdj48ZGl2IHN0eWxlPSJmb250LXNpemU6IDE0cHg7Ij48Zm9udCBkYXRhLWZvbnQtc3JjPSJodHRwczovL2ZvbnRzLmdvb2dsZWFwaXMuY29tL2Nzcz9mYW1pbHk9QXJpYWwiIHN0eWxlPSJmb250LXNpemU6IDE0cHg7IiBmYWNlPSJBcmlhbCI+R2VuZSBGdXNpb25zPGJyIC8+PC9mb250PjwvZGl2PjwvZGl2PjwvZGl2PjwvZGl2PjwvZm9yZWlnbk9iamVjdD48dGV4dCB4PSI5MCIgeT0iNjE0IiBmaWxsPSJyZ2IoMCwgMCwgMCkiIGZvbnQtZmFtaWx5PSJIZWx2ZXRpY2EiIGZvbnQtc2l6ZT0iMTJweCIgdGV4dC1hbmNob3I9Im1pZGRsZSI+RnVzaW9uU2Vla2VyLi4uPC90ZXh0Pjwvc3dpdGNoPjwvZz48cGF0aCBkPSJNIDE4MSA2MTAgTCAxNTcuMzcgNjEwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAxNTIuMTIgNjEwIEwgMTU5LjEyIDYwNi41IEwgMTU3LjM3IDYxMCBMIDE1OS4xMiA2MTMuNSBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDE4MCA3MjAgTCAxNTYuMzcgNzIwIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAxNTEuMTIgNzIwIEwgMTU4LjEyIDcxNi41IEwgMTU2LjM3IDcyMCBMIDE1OC4xMiA3MjMuNSBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48cGF0aCBkPSJNIDI2MCA0MzAgTCAyNjAgNDUzLjYzIiBmaWxsPSJub25lIiBzdHJva2U9InJnYigwLCAwLCAwKSIgc3Ryb2tlLW1pdGVybGltaXQ9IjEwIiBwb2ludGVyLWV2ZW50cz0ic3Ryb2tlIi8+PHBhdGggZD0iTSAyNjAgNDU4Ljg4IEwgMjU2LjUgNDUxLjg4IEwgMjYwIDQ1My42MyBMIDI2My41IDQ1MS44OCBaIiBmaWxsPSJyZ2IoMCwgMCwgMCkiIHN0cm9rZT0icmdiKDAsIDAsIDApIiBzdHJva2UtbWl0ZXJsaW1pdD0iMTAiIHBvaW50ZXItZXZlbnRzPSJhbGwiLz48L2c+PHN3aXRjaD48ZyByZXF1aXJlZEZlYXR1cmVzPSJodHRwOi8vd3d3LnczLm9yZy9UUi9TVkcxMS9mZWF0dXJlI0V4dGVuc2liaWxpdHkiLz48YSB0cmFuc2Zvcm09InRyYW5zbGF0ZSgwLC01KSIgeGxpbms6aHJlZj0iaHR0cHM6Ly93d3cuZGlhZ3JhbXMubmV0L2RvYy9mYXEvc3ZnLWV4cG9ydC10ZXh0LXByb2JsZW1zIiB0YXJnZXQ9Il9ibGFuayI+PHRleHQgdGV4dC1hbmNob3I9Im1pZGRsZSIgZm9udC1zaXplPSIxMHB4IiB4PSI1MCUiIHk9IjEwMCUiPlRleHQgaXMgbm90IFNWRyAtIGNhbm5vdCBkaXNwbGF5PC90ZXh0PjwvYT48L3N3aXRjaD48L3N2Zz4=)"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "AUKHTlobUcf_"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{Parameter Input and User Instructions}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oM-NaZDK63YC"
      },
      "source": [
        "Please define where the file structure is within your Google Drive:"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "uc1rzY2x6tvu"
      },
      "source": [
        "<ul type=disc>\n",
        "<li><b>PIPELINE_FILE_PATH</b>: file path to location of long-read RNA sequencing analysis pipeline within your Google Colab/Google Drive/local file system - required for most applications. This must be defined and initialized first.</li>\n",
        "</ul>"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "cviQcgv7iqvo"
      },
      "source": [
        "Please note: when utilizing the default ```/content``` location on a colab machine, click on the folder icon (fourth icon from the top on the left menu) to see files listed in a graphical user interface.  "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_P4OAFNrzifE"
      },
      "outputs": [],
      "source": [
        "%env PIPELINE_FILE_PATH=/content"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_GxZ1UCRTDkt"
      },
      "source": [
        "Errors associated with the following command are most likely associated with the existence of the directory/folder where the pipeline will be installed. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Hpn2CRCkJTfT"
      },
      "outputs": [],
      "source": [
        "! mkdir $PIPELINE_FILE_PATH"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "I24mo0O8JF4P"
      },
      "outputs": [],
      "source": [
        "! cd $PIPELINE_FILE_PATH ; git clone https://github.com/Theo-Nelson/long-read-sequencing-pipeline"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "2USxlBU47A0L"
      },
      "source": [
        "Please modify the following parameters within the code box below to fit your own study requirements:  "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "f99vaAE0Z0AG"
      },
      "source": [
        "<li><b>ACC</b>: Run accession number for reads within the [European Nucleotide Archive](https://www.ebi.ac.uk/ena/browser/) (SRR...) or file path to location of long-read RNA sequencing data within your Google Drive/general file system - required for most applications</li>\n",
        "<li><b>INDEX_FILE_PATH</b>: file path to location of reference genome (e.g. .FASTA) within your Google Drive/general file system - required for most applications</li>\n",
        "<li><b>ANNOTATION_FILE_PATH</b>: file path to location of reference annotation (e.g. .GTF) within your Google Drive/general file system - required for most applications</li>\n",
        "<li><b>CHROMOSOME</b>: Name of Chromosome of Interest matching the name of the Chromosome within your Reference Annotation - required for svist4get</li>\n",
        "<li><b>CHROMOSOME_START</b>: Starting Location of Interest on the Chromosome - required for svist4get</li>\n",
        "<li><b>CHROMOSOME_FINISH</b>: Ending Location of Interest on the Chromosome - required for svist4get</li>\n",
        "<li><b>REGION_NAME</b>: Gene Name (should match annotation file) - required for FLAME and svist4get</li>\n",
        "<li><b>HUB_KEYWORD</b>: Short Keyword for your UCSC Track Hub - required for MakeHub</li>\n",
        "<li><b>HUB_NAME</b>: Longer Title for your UCSC Track Hub - required for MakeHub</li>\n",
        "<li><b>HUB_EMAIL</b>: Email for your UCSC Track Hub (if you publish your track hub then this email will be public) - required for MakeHub</li>\n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "tqGM1OFhUcoo"
      },
      "outputs": [],
      "source": [
        "%env ACC=DRR059311\n",
        "%env INDEX_FILE_PATH=${PIPELINE_FILE_PATH}/long-read-sequencing-pipeline/prebuilt_indices/hg38.fa\n",
        "%env ANNOTATION_FILE_PATH=${PIPELINE_FILE_PATH}/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001405.29_GRCh38.p14_genomic.gtf\n",
        "%env CHROMOSOME=chr12\n",
        "%env CHROMOSOME_START=116533422\n",
        "%env CHROMOSOME_FINISH=116536513\n",
        "%env REGION_NAME=LINC00173\n",
        "%env HUB_KEYWORD=LINC00173\n",
        "%env HUB_NAME=\"Human LINC00173\"\n",
        "%env HUB_EMAIL=your@email.address"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "RrTcAlHtbWfg"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{Mounting your Google Drive / Exporting to Your Local Hard Drive}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "C5_iJ_KkbYt3"
      },
      "source": [
        "This step allows for permanent storage of your bioinformatics analysis in Google Drive"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "rbS0P9JcJnnJ"
      },
      "source": [
        "<ul type=disc>\n",
        "<li><b>STORAGE_FILE_PATH</b>: file path to location where you would wish to store output of long-read RNA sequencing analysis pipeline - required to export data from Google Colab. </li>\n",
        "</ul>"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "igrdHc7NJ3Zr"
      },
      "outputs": [],
      "source": [
        "%env STORAGE_FILE_PATH=/content/drive/MyDrive/long-read-sequencing-pipeline"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "aQOezcA5ZrFA"
      },
      "outputs": [],
      "source": [
        "from google.colab import drive\n",
        "drive.mount('/content/drive')"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "e1fBbVvQgvYF"
      },
      "source": [
        "Errors associated with the following command are most likely associated with the existence of the directory/folder where the pipeline will be installed. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "aqCDej0lgrX7"
      },
      "outputs": [],
      "source": [
        "! mkdir $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "fZaStpCheIm0"
      },
      "source": [
        "Additionally, this command will export necessary parameter data to allow downloads to a local machine. Please ignore ```/bin/bash: line 0: export: `/content': not a valid identifier``` or similar errors. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "reOYVKmaeI81"
      },
      "outputs": [],
      "source": [
        "! export $ACC\n",
        "! export $PIPELINE_FILE_PATH\n",
        "! export $REGION_NAME"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "24vu4E_8vbKa"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{BioConda: Package Installations}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "lz3I3YU06E6z"
      },
      "source": [
        "BioConda is a software environment and package manager, providing acess to over 8,000 different software packages related to bioinformatics (documentation: [BioConda](https://bioconda.github.io/user/install.html) and [Managing Environments via Conda](https://docs.conda.io/projects/conda/en/latest/user-guide/tasks/manage-environments.html)). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "oxEPAO-_wpbz"
      },
      "outputs": [],
      "source": [
        "! wget https://repo.anaconda.com/miniconda/Miniconda3-py37_4.8.2-Linux-x86_64.sh\n",
        "! chmod +x Miniconda3-py37_4.8.2-Linux-x86_64.sh\n",
        "! bash ./Miniconda3-py37_4.8.2-Linux-x86_64.sh -b -f -p /usr/local"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "EZS5nfFJCNHz"
      },
      "outputs": [],
      "source": [
        "! export PYTHONHOME='/usr/local/lib/python3.7/site-packages/'"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "z1lysGn4hO-i"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{Kingfisher: procurement of sequence files - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "X4iuuuqPkKcT"
      },
      "source": [
        "The Kingfisher program allows for sequence files to be downloaded from the European Nucleotide Archive (documentation: https://github.com/wwood/kingfisher-download)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "l9lQPVf27Itd"
      },
      "outputs": [],
      "source": [
        "! conda create -n kingfisher -c conda-forge -c bioconda python=3.9 kingfisher -y"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "fDjCnO4C3NQu"
      },
      "outputs": [],
      "source": [
        "! conda install -c rpetit3 aspera-connect -y"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ksJvQOCtx8BJ"
      },
      "outputs": [],
      "source": [
        "! wget -qO- https://download.asperasoft.com/download/sw/connect/3.9.8/ibm-aspera-connect-3.9.8.176272-linux-g2.12-64.tar.gz | tar xvz"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "RHjWyX6Ei_jT"
      },
      "source": [
        "this command will pop up with an error message; please disregard"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "anXSgAYa3xQU"
      },
      "outputs": [],
      "source": [
        "! ./ibm-aspera-connect-3.9.8.176272-linux-g2.12-64.sh"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "F0hFNwmcv5qr"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{FastQC: quality control tool for high-throughput sequence data - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "68fqL9EGNOKE"
      },
      "source": [
        "FastQC is a program designed to spot potential problems in high througput sequencing datasets ([documentation](https://github.com/s-andrews/FastQC))."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_bRfEf9e3Sp-"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda fastqc -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "k_oe5__lB6Lz"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{Shark: fishing relevant reads in an RNA-Seq sample - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Pa8jDrcgB9xl"
      },
      "source": [
        "Shark is a tool to extract gene-specific reads from a RNA-seq sample (documentation: https://github.com/AlgoLab/shark)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "PEjwM6_7B96Z"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda shark -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "gyh52FxVTSQt"
      },
      "source": [
        "###$\\color{#e74b4b}{\\text{minimap2: A versatile pairwise aligner for spliced nucleotide sequences - installation}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "0bXGX4mEeegN"
      },
      "source": [
        "minimap2 is a long-read sequencing aligner (documentation: https://github.com/lh3/minimap2). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "yZz8ktFYTQ1F"
      },
      "outputs": [],
      "source": [
        "! conda create -n minimap2 -c conda-forge -c bioconda minimap2=2.23 -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "dxspwFM3_1qN"
      },
      "source": [
        "###$\\color{#e74b4b}{\\text{samtools: Write/Index SAM to BAM - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "SXphSNhPA9di"
      },
      "source": [
        "samtools allows for manipulation of high-throughput sequencing data (documentation: http://www.htslib.org/) \n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "v0_UcWdBHDGQ"
      },
      "outputs": [],
      "source": [
        "! conda create -n samtools -c conda-forge -c bioconda python=3.7 samtools -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Y9qrP4Ip3xwz"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{TranscriptClean: correct mismatches, microindels, and noncanonical splice junctions - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "bB4Zr5Ef38qi"
      },
      "source": [
        "TranscriptClean is a command-line program which corrects long-read mismatches, microindels and noncanonical splice junctions (documentation: https://github.com/mortazavilab/TranscriptClean). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "oXIztTkj3zbu"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda pyfasta pyranges samtools -y"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "BfsI59QG4BH_"
      },
      "outputs": [],
      "source": [
        "! wget https://github.com/mortazavilab/TranscriptClean/archive/refs/tags/v2.0.3.tar.gz\n",
        "! tar xvf v2.0.3.tar.gz "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "3U_BpHSJA9Xe"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{FLAME: long-read splice variant annotation - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7_98iJtTy1wr"
      },
      "source": [
        "Full-Length Adjacency Matrix Enumeration (FLAME) is a program which can detect and quantify novel splice junctions on an annotated gene locus (documentation: https://github.com/marabouboy/FLAME). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "cyzQjKxz7H1v"
      },
      "outputs": [],
      "source": [
        "! sudo apt-get install bedtools"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "rvgXztd99NMC"
      },
      "outputs": [],
      "source": [
        "! pip install pysam"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "pKvCCPWvBJ2b"
      },
      "outputs": [],
      "source": [
        "! git clone https://github.com/marabouboy/FLAME"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "TZaGFSOoVUCo"
      },
      "outputs": [],
      "source": [
        "! chmod 755 /content/FLAME/FLAME/FLAME.py"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "7Erq_G3QVlz9"
      },
      "outputs": [],
      "source": [
        "! chmod 755 /content/FLAME/setup.py"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "grBakyZXVv3s"
      },
      "outputs": [],
      "source": [
        "! cd /content/FLAME/ ; python3 setup.py install"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "yhR96h7O5CHN"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{featureCounts: assign sequence reads to genomic features - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "VnlTGqwG5JYz"
      },
      "source": [
        "featureCounts is a read summarization program suitable for counting reads generated from either RNA or genomic DNA sequencing experiments (documentation: http://subread.sourceforge.net/). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Vd48uWN45VkA"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda subread -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "vZPSGeXlh7kh"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{LIQA: transcript quantification}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "vMq1SRmih_SS"
      },
      "source": [
        "LIQA is a program which quantifies isoform/transcript expression based on long read RNA sequencing data (documentation: https://github.com/WGLab/LIQA)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ggyRF2F953H8"
      },
      "outputs": [],
      "source": [
        "! pip install liqa"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oOS4NvQ8o5Q7"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{FusionSeeker: detect gene fusions - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "mJKun3uw-Ihy"
      },
      "source": [
        "FussionCaller is a gene fusion caller for long-read single-molecular sequencing data (documentation: https://github.com/Maggi-Chen/FusionSeeker). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "fvgPhbBwmkzP"
      },
      "outputs": [],
      "source": [
        "! git clone https://github.com/Theo-Nelson/FusionSeeker.git"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "bEPIzmwa2dFO"
      },
      "outputs": [],
      "source": [
        "! git clone https://github.com/ruanjue/bsalign.git"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "mlblnJtF2eyk"
      },
      "outputs": [],
      "source": [
        "! cd /content/bsalign && make"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EFHE0ybf9rRI"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{StringTie: transcript assembly - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "WCeEyrb4S0xW"
      },
      "source": [
        "StringTie is a program which can produce transcriptomes specific to the sample input (documentation: http://ccb.jhu.edu/software/stringtie/index.shtml?t=manual#run)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ljsRmfFksCki"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda stringtie -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "CId-MKlu9uMZ"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{GffCompare: transcript assembly statistics - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Hy2ietg1UPWj"
      },
      "source": [
        "GffCompare is a program which can evaluate the specificity and novelty of transcripts within a sample-specific transcriptome to a reference transcriptome (documentation: http://ccb.jhu.edu/software/stringtie/gffcompare.shtml)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "iKx4Bn4wSyLb"
      },
      "outputs": [],
      "source": [
        "! conda install -c bioconda gffcompare -y"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "npZwllsupfEJ"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{svist4get: visualize genomic tracks from sequencing experiments - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "niesqrQ1nHMM"
      },
      "source": [
        "svist4get allows you to view read coverage at a defined region on a chromosome (documentation: https://bitbucket.org/artegorov/svist4get/src/master/)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "iaRNkYPEoa0T"
      },
      "outputs": [],
      "source": [
        "! apt-get update"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "-gNi5jxXA5oq"
      },
      "outputs": [],
      "source": [
        "! apt-get install libmagickwand-dev"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "3HS96PJ45f26"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/svist4get/policy_revised.xml /etc/ImageMagick-6/policy.xml"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ZNVRJydu1wIV"
      },
      "outputs": [],
      "source": [
        "! python3 -m pip install svist4get"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "5CwywN4tc74s"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{Pistis: quality control plotting for long reads - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ckgZcyKcpFyC"
      },
      "source": [
        "Pistis generates long-read specific quality control graphs, including a plot demonstrating read alignment percentage to the reference genome (documentation: https://github.com/mbhall88/pistis)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "S5mHHUr1d50q"
      },
      "outputs": [],
      "source": [
        "! pip3 install pistis"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EQn9Ikb3YkFC"
      },
      "source": [
        "###$\\color{#0072ff}{\\text{MakeHub: generate UCSC assembly hubs - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "X7J1qXU1bmXd"
      },
      "source": [
        "MakeHub is a command line tool for the fully automatic generation of of track data hubs for visualizing genomes with the UCSC genome browser (documentation: https://github.com/Gaius-Augustus/MakeHub)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_eYWwdsmdsBa"
      },
      "outputs": [],
      "source": [
        "! python3.7 -m pip install biopython"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "yudzuOTW2IQ-"
      },
      "outputs": [],
      "source": [
        "! sudo apt install augustus augustus-data augustus-doc"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "kpfzCCTKuqMt"
      },
      "source": [
        "###$\\color{#ff00d5}{\\text{MultiQC: aggregate bioinformatics analysis - installation}}$\n",
        "\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Jluk65BgNoEJ"
      },
      "source": [
        "MultiQC is a program which allows you to combine reports for as many samples as you wish ([documentation](https://multiqc.info/docs/))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "B6OuaTChNs4J"
      },
      "outputs": [],
      "source": [
        "! pip install multiqc"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oBdBA_CDdE0z"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{Kingfisher: procurement of sequence files - usage}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "IE4T9swGdQY-"
      },
      "source": [
        "The Kingfisher program allows for sequence files to be downloaded from the European Nucleotide Archive (documentation: https://github.com/wwood/kingfisher-download)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ReZdoTfnd5uQ"
      },
      "outputs": [],
      "source": [
        "! cd $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq ; source activate kingfisher ; kingfisher get -r $ACC -m ena-ascp aws-http prefetch"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "V5WMZKjYx6YU"
      },
      "source": [
        "The next command will unzip a fastq file, if neccessary. Please do not be concerned if this command throws an error. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "wwpiZ9OU6WPn"
      },
      "outputs": [],
      "source": [
        "! cd $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq ; gunzip *.gz"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "izrcNnG0REyD"
      },
      "source": [
        "The next two commands simply standardize the file name away from uncommon variants provided by depositors in the European Nucleotide Archive. In rare instances pipeline users may need to directly manipulate these commands to standardize the filename. Please do not be concerned if either command throws an error. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "4pmkDR2x6SzT"
      },
      "outputs": [],
      "source": [
        "! cd $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq ; mv ${ACC}_1.fastq $ACC.fastq"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "5_r7Ed86Q-cI"
      },
      "outputs": [],
      "source": [
        "! cd $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq ; mv ${ACC}_subreads.fastq $ACC.fastq"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_MislmkXZcz-"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{Kingfisher: procurement of sequence files - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "-mPcrcf-Kmfl"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "U0pvTRDjKtqG"
      },
      "outputs": [],
      "source": [
        "! eval cp $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$ACC.fastq $STORAGE_FILE_PATH/$ACC.fastq"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "txLoObMrKrnJ"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "sWZhEEl3KuIu"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/fastq/\",os.path.expandvars(os.path.expandvars(\"$ACC\")),\".fastq\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "SQb5HTiwtfMn"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{FastQC: quality control tool for high-throughput sequence data - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZoutKrvVtjAY"
      },
      "source": [
        "FastQC is a program designed to spot potential problems in high througput sequencing datasets ([documentation](https://github.com/s-andrews/FastQC))."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "bTYA5m273eVN"
      },
      "outputs": [],
      "source": [
        "! eval fastqc -t 2 $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$ACC.fastq --outdir $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastqc"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "nHGHIAKCSUen"
      },
      "source": [
        "For the best viewing experience, please download the HTML output to your Hard Drive and open in your local browser, such as Google Chrome, Firefox or Edge. "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "NOrgEp-fbtBt"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{FastQC: quality control tool for high-throughput sequence data - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "bGhAXFSfbtBx"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "z2dsXIK4btBy"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastqc/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "6CZoOeptbtB0"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "H3dCOerJqUoS"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastqc.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastqc/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "LqiW6frsbtB1"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/fastqc.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7uHwZ4Zj1vgC"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{Reference Genome - installation}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "5FvumDB7c6bH"
      },
      "source": [
        "The Reference Genome provides a scaffold to align long-read data. These commands install the hg38 genome and ensembl annotation availiable from UCSC. You can download more current genomes by utilizing the appropriate links from NCBI RefSeq, Ensembl, or other reference genome providers. If you are unsure of how to find other species, we recommend checking out the list of species available in the ```current_fasta``` and ```current_gtf``` folders: http://ftp.ensembl.org/pub/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_D5cG18VURVJ"
      },
      "outputs": [],
      "source": [
        "! wget -P $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Jy4bAiZ6z3Qb"
      },
      "outputs": [],
      "source": [
        "! wget -P $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices https://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/405/GCA_000001405.29_GRCh38.p14/GCA_000001405.29_GRCh38.p14_genomic.gtf.gz"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "bz8hNYEvUXOE"
      },
      "outputs": [],
      "source": [
        "! gunzip -c $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/hg38.fa.gz > $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/hg38.fa"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "z2Hqu_wb0NOv"
      },
      "outputs": [],
      "source": [
        "! gunzip -c $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001405.29_GRCh38.p14_genomic.gtf.gz > $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001405.29_GRCh38.p14_genomic.gtf"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "iulCArtGA6sw"
      },
      "source": [
        "Another example is provided for which installs the UCSC murine mm39 genome. Note that the reference annotation is from RefSeq. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "rf6DUsMHBUmp"
      },
      "outputs": [],
      "source": [
        "! wget -P $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices https://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/635/GCA_000001635.9_GRCm39/GCA_000001635.9_GRCm39_genomic.fna.gz"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "IL949aPzBUsA"
      },
      "outputs": [],
      "source": [
        "! wget -P $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices https://ftp.ncbi.nlm.nih.gov/genomes/all/GCA/000/001/635/GCA_000001635.9_GRCm39/GCA_000001635.9_GRCm39_genomic.gtf.gz"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "vAqaExDBBUxi"
      },
      "outputs": [],
      "source": [
        "! gunzip -c $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001635.9_GRCm39_genomic.fna.gz > $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001635.9_GRCm39_genomic.fa"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "suKULpveBU20"
      },
      "outputs": [],
      "source": [
        "! gunzip -c $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001635.9_GRCm39_genomic.gtf.gz > $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/GCA_000001635.9_GRCm39_genomic.gtf"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "D6iKV4YDiJQk"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{Shark: fishing relevant reads in an RNA-Seq sample - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "9BtvLo6dHRxt"
      },
      "source": [
        "Shark is a tool to extract gene-specific reads from a RNA-seq sample (documentation: https://github.com/AlgoLab/shark)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Jjx8LrclHVCH"
      },
      "outputs": [],
      "source": [
        "! echo ${CHROMOSOME}$'\\t'${CHROMOSOME_START}$'\\t'${CHROMOSOME_FINISH} > $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$REGION_NAME.bed"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "BxccuEsbFpyJ"
      },
      "outputs": [],
      "source": [
        "! eval bedtools getfasta -fi $INDEX_FILE_PATH -bed $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$REGION_NAME.bed -fo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$REGION_NAME.fasta"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "FB1MpOyLnRxD"
      },
      "outputs": [],
      "source": [
        "! eval shark -c .40 -k 10 -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$REGION_NAME.fasta -1 $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$ACC.fastq -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$REGION_NAME$ACC.fastq"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "cDYzmwgwdZtX"
      },
      "source": [
        "In order to continue analysis with just the filtered reads, please utilize the following commands:"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "1uQp4_LvYtPi"
      },
      "outputs": [],
      "source": [
        "import os\n",
        "os.environ[\"ACC_cp\"] = os.environ[\"ACC\"]"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "mZGMFpxiRROC"
      },
      "outputs": [],
      "source": [
        "%env ACC = $REGION_NAME$ACC_cp"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "EJ0KI6ZxJInS"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{Shark: fishing relevant reads in an RNA-Seq sample - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Q3quwJI4JInW"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "11khL1RgJInX"
      },
      "outputs": [],
      "source": [
        "! eval cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$REGION_NAME.fastq $STORAGE_FILE_PATH/$REGION_NAME$ACC.fastq"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ePctHvJkJInY"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "a1IIimzeJInb"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/fastq/\",os.environ[\"REGION_NAME\"],os.environ[\"ACC\"],\".fastq\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "oBPQyTcKOPDF"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{minimap2: A versatile pairwise aligner for spliced nucleotide sequences - index minimization}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "of339UK9OT_0"
      },
      "source": [
        "minimap2 is a long-read sequencing aligner (documentation: https://github.com/lh3/minimap2). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ZFQz8s-vFrNG"
      },
      "outputs": [],
      "source": [
        "! source activate minimap2 ; eval minimap2 -I2G -k15 -w5 -d $INDEX_FILE_PATH.mmi $INDEX_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "f8EfwQZze29B"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{minimap2: A versatile pairwise aligner for spliced nucleotide sequences - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "xlj72R1ofOAz"
      },
      "source": [
        "minimap2 is a long-read sequencing aligner (documentation: https://github.com/lh3/minimap2). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "YjPnVj5PTRtK"
      },
      "outputs": [],
      "source": [
        "! source activate minimap2 ; eval minimap2 --split-prefix -k15 -ax splice $INDEX_FILE_PATH.mmi $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$ACC.fastq > $(eval echo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/sam/$ACC.sam)"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "J_OH_OItuBjO"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{minimap2: A versatile pairwise aligner for spliced nucleotide sequences - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tCKOQBzbuBjR"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "xEM-Bg_nuBjS"
      },
      "outputs": [],
      "source": [
        "! eval cp $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/sam/$ACC.sam $STORAGE_FILE_PATH/$ACC.sam"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sfmGI2VruBjS"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "zpy-NJKXuBjU"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/sam/\",os.path.expandvars(os.path.expandvars(\"$ACC\")),\".sam\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sjns1XbAkILI"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{samtools: Write/Index SAM to BAM - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "KfgkJ0bJkYDv"
      },
      "source": [
        "samtools allows for manipulation of high-throughput sequencing data (documentation: http://www.htslib.org/) \n"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "fAFTMfy2HDOs"
      },
      "outputs": [],
      "source": [
        "! source activate samtools ; eval samtools view -S -b $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/sam/$ACC.sam > $(eval echo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.bam)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "qYOZbl8PJo-h"
      },
      "outputs": [],
      "source": [
        "! source activate samtools ; eval samtools sort $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.bam -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam  "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ypjDsfP7JpM6"
      },
      "outputs": [],
      "source": [
        "! source activate samtools ; eval samtools index $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "swcyphLAvE52"
      },
      "source": [
        "##$\\color{#e74b4b}{\\text{samtools: Write/Index SAM to BAM - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "GbZSSgbkvE54"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "qi5xpmLnvE55"
      },
      "outputs": [],
      "source": [
        "! eval cp $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam $STORAGE_FILE_PATH/$ACC.sorted.bam"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_QPGPrNnXNi_"
      },
      "outputs": [],
      "source": [
        "! eval cp $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam.bai $STORAGE_FILE_PATH/$ACC.sorted.bam.bai"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "psCo62ZRvE57"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "TSjtffSHvE59"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/bam/\",os.path.expandvars(os.path.expandvars(\"$ACC\")),\".sorted.bam\"]))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "sD5KruF8XV86"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/bam/\",os.path.expandvars(os.path.expandvars(\"$ACC\")),\".sorted.bam.bai\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "zqQcBdxa4X7q"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{TranscriptClean: correct mismatches, microindels, and noncanonical splice junctions - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "r5CNnQ5w4k_K"
      },
      "source": [
        "TranscriptClean is a command-line program which corrects long-read mismatches, microindels and noncanonical splice junctions (documentation: https://github.com/mortazavilab/TranscriptClean). "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZU_kVYSvWXdN"
      },
      "source": [
        "Please note that the corrected reads should not be utilized for downstream-level base-calling analysis such as variant calling. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "sJPb4PSU5LmU"
      },
      "outputs": [],
      "source": [
        "! eval python /content/TranscriptClean-2.0.3/TranscriptClean.py --sam $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/sam/$ACC.sam --genome $INDEX_FILE_PATH --outprefix $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/transcriptclean/$ACC"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "11Qa37-tea-K"
      },
      "source": [
        "Convert the resulting sam file to a bam file:"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "bn_4YOoJd68N"
      },
      "outputs": [],
      "source": [
        "! eval samtools view -S -b $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/transcriptclean/$(eval echo $ACC)_clean.sam > $(eval echo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$(eval echo $ACC)_clean.bam)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "kVdkAQw0eStU"
      },
      "outputs": [],
      "source": [
        "! eval samtools sort $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$(eval echo $ACC)_clean.bam -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$(eval echo $ACC)_clean.sorted.bam  "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "G8A4OL32eXcI"
      },
      "outputs": [],
      "source": [
        "! eval samtools index $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$(eval echo $ACC)_clean.sorted.bam "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "iUwbIIRkwWrL"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{TranscriptClean: correct mismatches, microindels, and noncanonical splice junctions - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "1K-zlCYDwWrL"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "j_1xnKJJwWrM"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/transcriptclean/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "tcK1cY7pwWrN"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_t5YZRckdd_T"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/transcriptclean.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/transcriptclean/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "RW7dqOgpwWrP"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/transcriptclean.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "BKjT28HjsBkW"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{StringTie: transcript assembly - usage}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "FuoJtl20TNAx"
      },
      "source": [
        "StringTie is a program which can collapse sample reads into transcriptomes specific to the sample input (documentation: http://ccb.jhu.edu/software/stringtie/index.shtml?t=manual#run)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_gRdTG0lsDQb"
      },
      "outputs": [],
      "source": [
        "! eval stringtie -L -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/stringtie/${ACC}.gtf -G $ANNOTATION_FILE_PATH $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "xiPsVkfrUw6S"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{StringTie: transcript assembly - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "or7JeoLlUw6U"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "EWrA0eBPUw6V"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/stringtie $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Y7ib8gNMUw6Y"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "MglPWDSTUw6b"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/stringtie.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/stringtie/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "EvsEJOAdUw6c"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/stringtie.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "TAN1_OR6TfkL"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{GffCompare: transcript assembly statistics - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "wjUMZNQ0TwR9"
      },
      "source": [
        "GffCompare is a program which can evaluate the specificity and novelty of transcripts within a sample-specific transcriptome to a reference transcriptome (documentation: http://ccb.jhu.edu/software/stringtie/gffcompare.shtml)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "G3BM5t4ZzVHD"
      },
      "outputs": [],
      "source": [
        "! eval gffcompare -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/gffcompare/$ACC -r $ANNOTATION_FILE_PATH -R -Q $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/stringtie/${ACC}.gtf"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "W8SCLFLZUyVs"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{GffCompare: transcript assembly statistics - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "TEKaFxwJUyVs"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "zsj9Lt__UyVt"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/gffcompare $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "kC2Gvg3nUyVu"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "U1eeIJFhUyVv"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/gffcompare.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/gffcompare/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ic6DjrjFUyVv"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/gffcompare.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "y3bt2k_tBGnT"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{FLAME: gene-specific long-read splice variant annotation - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "kXs3BJZqygs7"
      },
      "source": [
        "Full-Length Adjacency Matrix Enumeration (FLAME) is a program which can detect and quantify novel splice junctions on an annotated gene locus (documentation: https://github.com/marabouboy/FLAME). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "JB86-GcXZHIT"
      },
      "outputs": [],
      "source": [
        "! eval bedtools bamtobed -bed12 -i $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam > $(eval echo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bed/$ACC.sorted.bed12)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_Y4DKdo5LhBB"
      },
      "outputs": [],
      "source": [
        "! eval python3 /content/FLAME/FLAME/FLAME.py -I $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bed/$ACC.sorted.bed12 -GTF $ANNOTATION_FILE_PATH -G $REGION_NAME -O $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/flame/$ACC"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "yNrDbLpg4H6_"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{FLAME: gene-specific long-read splice variant annotation - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "gtxRHyEO4H7A"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "cLxi8NXw4H7C"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/flame $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "lM7C8DCy4H7D"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "KdiyIHxkwdFT"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/flame.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/flame/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "T3GgoMYP4H7F"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/flame.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "7yVbaSBc5cu7"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{featureCounts: assign sequence reads to genomic features - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "-2kbK23S5oi5"
      },
      "source": [
        "featureCounts is a read summarization program suitable for counting reads generated from either RNA or genomic DNA sequencing experiments (documentation: http://subread.sourceforge.net/). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "1ycdSYt_5plb"
      },
      "outputs": [],
      "source": [
        "! eval featureCounts -O -L -a $ANNOTATION_FILE_PATH -t exon -g gene_id -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/featureCounts/$ACC.txt $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "a2s7KTWRw3ry"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{featureCounts: assign sequence reads to genomic features - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZsTqqT1aw3r0"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "StNofMQFw3r1"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/featureCounts/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "jqBqz7Ztw3r2"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "CpJsMgUNVogc"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/featureCounts.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/featureCounts/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "udK69J1FpWIF"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/featureCounts.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "v8jXqioM50fY"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{LIQA: transcript quantification - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "MrRLmJ1Ojbcc"
      },
      "source": [
        "LIQA is a program which quantifies isoform/transcript expression based on long read RNA sequencing data (documentation: https://github.com/WGLab/LIQA)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "2_zvCVli55EX"
      },
      "outputs": [],
      "source": [
        "! eval liqa -task refgene -ref $ANNOTATION_FILE_PATH -format gtf -out $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$ACC.refgene"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "Hy-E4YBT7iCr"
      },
      "outputs": [],
      "source": [
        "! eval liqa -task quantify -refgene $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/prebuilt_indices/$ACC.refgene -bam $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam -out $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/liqa/${ACC}_quantification.txt -max_distance 10 -f_weight 1"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "fQy3NibXkZww"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{LIQA: transcript quantification - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "zfLiaJmBkZwx"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "lNQuGyYNkZwy"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/liqa/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZJ50Bs_RkZwz"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "4c9Rnp_HkZwz"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/liqa.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/liqa/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "TPLA1-_6kZw0"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/liqa.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "VXojzd-9mjPc"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{FusionSeeker: detect gene fusions - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "kAb1EVhSpyhS"
      },
      "source": [
        "FussionCaller is a gene fusion caller for long-read single-molecular sequencing data (documentation: https://github.com/Maggi-Chen/FusionSeeker). "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "_x9eJ_oDmmU8"
      },
      "outputs": [],
      "source": [
        "! eval /content/FusionSeeker/fusionseeker --bam $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam --outpath $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fusionseeker/ --ref $INDEX_FILE_PATH --gtf $ANNOTATION_FILE_PATH "
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "veJUugPx8-Q_"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{FusionSeeker: detect gene fusions - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "A9IgEX-88-RE"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "nx5uVzJX8-RF"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fusionseeker $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "6JwKYuga8-RG"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "aO9VI51c8-RI"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fusionseeker.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fusionseeker/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ge2aCdde8-RJ"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/fusionseeker.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "S9r4MCB0lvAf"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{svist4get: visualize genomic tracks from sequencing experiments - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "3T2Y3dcAnaWu"
      },
      "source": [
        "svist4get allows you to view read coverage at a defined region on a chromosome (documentation: https://bitbucket.org/artegorov/svist4get/src/master/)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ehHeatJQCZ42"
      },
      "outputs": [],
      "source": [
        "! eval bedtools genomecov -split -ibam $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam -bg > $(eval echo $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bed/$ACC.sorted.bedgraph)"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "OpQMW3gI15ZI"
      },
      "outputs": [],
      "source": [
        "! eval svist4get -bg $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bed/$ACC.sorted.bedgraph -gtf $ANNOTATION_FILE_PATH -fa $INDEX_FILE_PATH -bl Long-Read Coverage -w $CHROMOSOME $CHROMOSOME_START $CHROMOSOME_FINISH -it \"$REGION_NAME\" -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/svist4get/$ACC"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "emZMFcIZ4Z7q"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{svist4get: visualize genomic tracks from sequencing experiments - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "rB3h07vK4Z7s"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "HY784C9N4Z7t"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/svist4get/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "w6hQ-NhH4Z7u"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "c9bDP-P1vvA4"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/svist4get.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/svist4get/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "XyJTQFQZ4Z7w"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/svist4get.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "MjZ_tYz4pEHX"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{Pistis: quality control plotting for long reads - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "RFhhJITCphA_"
      },
      "source": [
        "Pistis generates long-read specific quality control graphs, including a plot demonstrating read alignment percentage to the reference genome (documentation: https://github.com/mbhall88/pistis)"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "B5Zv-aSzVPwK"
      },
      "source": [
        "Please note that the report generates assuming alignment of more than 50,000 reads. If this is not the case, please add the flag `--downsample INTEGER` replacing `INTEGER` with a number less than the number of aligned reads. Additionally, please ignore deprecation warnings. "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ik-boaBxeNjP"
      },
      "outputs": [],
      "source": [
        "! eval pistis --downsample 4000 -f $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/fastq/$ACC.fastq -b $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam  -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/pistis/$ACC.pdf"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZzePzpXU48T0"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{Pistis: quality control plotting for long reads - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "Ev9sTiLI48T2"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "rUoqVlq548T3"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/pistis/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "_XM7DHSH48T4"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "uPjj7MrHa23h"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/pistis.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/pistis/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "ToqqgNwD48T7"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/pistis.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "POLtVVCEKndf"
      },
      "source": [
        "##$\\color{#2977d6}{\\text{MakeHub: generate UCSC assembly hubs - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "ZcW7zXOoKrCy"
      },
      "source": [
        "MakeHub is a command line tool for the fully automatic generation of of track data hubs for visualizing genomes with the UCSC genome browser (documentation: https://github.com/Gaius-Augustus/MakeHub)."
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "fMN2TYdniufR"
      },
      "outputs": [],
      "source": [
        "! chmod 755 $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub/make_hub.py"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "pVA2mrrjZ2yM"
      },
      "outputs": [],
      "source": [
        "! eval $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub/make_hub.py -l $HUB_KEYWORD -L $HUB_NAME -g $INDEX_FILE_PATH -e \\\n",
        "  $HUB_EMAIL -a $ANNOTATION_FILE_PATH -b $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/bam/$ACC.sorted.bam -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub/"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "E2nQAifZ6Fn9"
      },
      "source": [
        "##$\\color{#4f7bb0}{\\text{MakeHub: generate UCSC assembly hubs - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "NTlyCHAj6Fn_"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "OyvrqB4n6FoA"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "AB1LYUCh6FoA"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "u_2BDvY0bz8Y"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/makehub/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "-vJjpEWh6FoD"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/makehub.zip\"]))"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "m5pU9i4DvGxQ"
      },
      "source": [
        "##$\\color{#d42bb4}{\\text{MultiQC: aggregate bioinformatics analysis - usage}}$"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "sPB5udNzvVdJ"
      },
      "source": [
        "MultiQC is a program which allows you to combine reports for as many samples as you wish ([documentation](https://multiqc.info/docs/))"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "SEBsMVQ5OG4U"
      },
      "outputs": [],
      "source": [
        "! eval multiqc $PIPELINE_FILE_PATH/long-read-sequencing-pipeline -o $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/multiqc/$ACC.pdf"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "n-WhLFAz5e0s"
      },
      "source": [
        "##$\\color{#a7588f}{\\text{MultiQC: aggregate bioinformatics analysis - export}}$\n"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "n98P06aH5e0t"
      },
      "source": [
        "To Store Resulting Files in your Google Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "28HahdKB5e0u"
      },
      "outputs": [],
      "source": [
        "! cp -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/multiqc/ $STORAGE_FILE_PATH"
      ]
    },
    {
      "cell_type": "markdown",
      "metadata": {
        "id": "FmRDMa3Q5e0v"
      },
      "source": [
        "To Store Resulting Files in your Local Hard Drive: "
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "mj5NMZTucrsV"
      },
      "outputs": [],
      "source": [
        "! zip -r $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/multiqc.zip $PIPELINE_FILE_PATH/long-read-sequencing-pipeline/multiqc/"
      ]
    },
    {
      "cell_type": "code",
      "execution_count": null,
      "metadata": {
        "id": "5o0K3bsz5e0x"
      },
      "outputs": [],
      "source": [
        "from google.colab import files\n",
        "import os\n",
        "\n",
        "files.download(\"\".join([os.environ[\"PIPELINE_FILE_PATH\"],\"/long-read-sequencing-pipeline/multiqc.zip\"]))"
      ]
    }
  ],
  "metadata": {
    "colab": {
      "collapsed_sections": [
        "AUKHTlobUcf_",
        "RrTcAlHtbWfg",
        "z1lysGn4hO-i",
        "F0hFNwmcv5qr",
        "k_oe5__lB6Lz",
        "gyh52FxVTSQt",
        "dxspwFM3_1qN",
        "Y9qrP4Ip3xwz",
        "3U_BpHSJA9Xe",
        "yhR96h7O5CHN",
        "vZPSGeXlh7kh",
        "oOS4NvQ8o5Q7",
        "EFHE0ybf9rRI",
        "CId-MKlu9uMZ",
        "npZwllsupfEJ",
        "5CwywN4tc74s",
        "EQn9Ikb3YkFC",
        "kpfzCCTKuqMt",
        "oBdBA_CDdE0z",
        "_MislmkXZcz-",
        "SQb5HTiwtfMn",
        "NOrgEp-fbtBt",
        "7uHwZ4Zj1vgC",
        "D6iKV4YDiJQk",
        "EJ0KI6ZxJInS",
        "oBPQyTcKOPDF",
        "f8EfwQZze29B",
        "J_OH_OItuBjO",
        "sjns1XbAkILI",
        "swcyphLAvE52",
        "zqQcBdxa4X7q",
        "iUwbIIRkwWrL",
        "BKjT28HjsBkW",
        "xiPsVkfrUw6S",
        "TAN1_OR6TfkL",
        "W8SCLFLZUyVs",
        "y3bt2k_tBGnT",
        "yNrDbLpg4H6_",
        "7yVbaSBc5cu7",
        "a2s7KTWRw3ry",
        "v8jXqioM50fY",
        "fQy3NibXkZww",
        "VXojzd-9mjPc",
        "veJUugPx8-Q_",
        "S9r4MCB0lvAf",
        "emZMFcIZ4Z7q",
        "MjZ_tYz4pEHX",
        "ZzePzpXU48T0",
        "POLtVVCEKndf",
        "E2nQAifZ6Fn9",
        "m5pU9i4DvGxQ",
        "n-WhLFAz5e0s"
      ],
      "machine_shape": "hm",
      "provenance": []
    },
    "gpuClass": "standard",
    "kernelspec": {
      "display_name": "Python 3",
      "name": "python3"
    },
    "language_info": {
      "name": "python"
    }
  },
  "nbformat": 4,
  "nbformat_minor": 0
}