scripts/gnomad-variants-af.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```\n",
    "dx download Annotation/gnomad_Freq/gnomad_v2.1_sv.sites.sorted.tsv\n",
    "dx download chr-subset-genes-variants/all_variants_chr21_grch37.tsv\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Loading required package: stats4\n",
      "\n",
      "Loading required package: BiocGenerics\n",
      "\n",
      "Loading required package: parallel\n",
      "\n",
      "\n",
      "Attaching package: ‘BiocGenerics’\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:parallel’:\n",
      "\n",
      "    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,\n",
      "    clusterExport, clusterMap, parApply, parCapply, parLapply,\n",
      "    parLapplyLB, parRapply, parSapply, parSapplyLB\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:dplyr’:\n",
      "\n",
      "    combine, intersect, setdiff, union\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:stats’:\n",
      "\n",
      "    IQR, mad, sd, var, xtabs\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:base’:\n",
      "\n",
      "    anyDuplicated, append, as.data.frame, basename, cbind, colnames,\n",
      "    dirname, do.call, duplicated, eval, evalq, Filter, Find, get, grep,\n",
      "    grepl, intersect, is.unsorted, lapply, Map, mapply, match, mget,\n",
      "    order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank,\n",
      "    rbind, Reduce, rownames, sapply, setdiff, sort, table, tapply,\n",
      "    union, unique, unsplit, which, which.max, which.min\n",
      "\n",
      "\n",
      "Loading required package: S4Vectors\n",
      "\n",
      "\n",
      "Attaching package: ‘S4Vectors’\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:dplyr’:\n",
      "\n",
      "    first, rename\n",
      "\n",
      "\n",
      "The following object is masked from ‘package:base’:\n",
      "\n",
      "    expand.grid\n",
      "\n",
      "\n",
      "Loading required package: IRanges\n",
      "\n",
      "\n",
      "Attaching package: ‘IRanges’\n",
      "\n",
      "\n",
      "The following objects are masked from ‘package:dplyr’:\n",
      "\n",
      "    collapse, desc, slice\n",
      "\n",
      "\n",
      "Loading required package: GenomeInfoDb\n",
      "\n"
     ]
    }
   ],
   "source": [
    "library(dplyr)\n",
    "library(GenomicRanges)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 6</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>V1</th><th scope=col>V2</th><th scope=col>V3</th><th scope=col>V4</th><th scope=col>V5</th><th scope=col>V6</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>chr1</td><td>10642</td><td> 10643</td><td><span style=white-space:pre-wrap>&lt;BND&gt;                                                         </span></td><td>gnomAD-SV_v2.1_BND_1_1 </td><td><span style=white-space:pre-wrap>0.006786                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>2</th><td>chr1</td><td>21000</td><td> 26000</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_1 </td><td><span style=white-space:pre-wrap>0.016073                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>3</th><td>chr1</td><td>40000</td><td> 47200</td><td><span style=white-space:pre-wrap>&lt;DUP&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DUP_1_1 </td><td><span style=white-space:pre-wrap>0.068963                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>4</th><td>chr1</td><td>46000</td><td>111250</td><td>&lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt;</td><td>gnomAD-SV_v2.1_MCNV_1_1</td><td>0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>chr1</td><td>50941</td><td> 51053</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_2 </td><td><span style=white-space:pre-wrap>0.004052                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>6</th><td>chr1</td><td>54665</td><td> 54666</td><td><span style=white-space:pre-wrap>&lt;INS&gt;                                                         </span></td><td>gnomAD-SV_v2.1_INS_1_1 </td><td><span style=white-space:pre-wrap>9.4e-05                                                                       </span></td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 6\n",
       "\\begin{tabular}{r|llllll}\n",
       "  & V1 & V2 & V3 & V4 & V5 & V6\\\\\n",
       "  & <chr> & <int> & <int> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & chr1 & 10642 &  10643 & <BND>                                                          & gnomAD-SV\\_v2.1\\_BND\\_1\\_1  & 0.006786                                                                      \\\\\n",
       "\t2 & chr1 & 21000 &  26000 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_1  & 0.016073                                                                      \\\\\n",
       "\t3 & chr1 & 40000 &  47200 & <DUP>                                                          & gnomAD-SV\\_v2.1\\_DUP\\_1\\_1  & 0.068963                                                                      \\\\\n",
       "\t4 & chr1 & 46000 & 111250 & <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8> & gnomAD-SV\\_v2.1\\_MCNV\\_1\\_1 & 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\\\\\n",
       "\t5 & chr1 & 50941 &  51053 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_2  & 0.004052                                                                      \\\\\n",
       "\t6 & chr1 & 54665 &  54666 & <INS>                                                          & gnomAD-SV\\_v2.1\\_INS\\_1\\_1  & 9.4e-05                                                                       \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 6\n",
       "\n",
       "| <!--/--> | V1 &lt;chr&gt; | V2 &lt;int&gt; | V3 &lt;int&gt; | V4 &lt;chr&gt; | V5 &lt;chr&gt; | V6 &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|\n",
       "| 1 | chr1 | 10642 |  10643 | &lt;BND&gt;                                                          | gnomAD-SV_v2.1_BND_1_1  | 0.006786                                                                       |\n",
       "| 2 | chr1 | 21000 |  26000 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_1  | 0.016073                                                                       |\n",
       "| 3 | chr1 | 40000 |  47200 | &lt;DUP&gt;                                                          | gnomAD-SV_v2.1_DUP_1_1  | 0.068963                                                                       |\n",
       "| 4 | chr1 | 46000 | 111250 | &lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt; | gnomAD-SV_v2.1_MCNV_1_1 | 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 |\n",
       "| 5 | chr1 | 50941 |  51053 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_2  | 0.004052                                                                       |\n",
       "| 6 | chr1 | 54665 |  54666 | &lt;INS&gt;                                                          | gnomAD-SV_v2.1_INS_1_1  | 9.4e-05                                                                        |\n",
       "\n"
      ],
      "text/plain": [
       "  V1   V2    V3    \n",
       "1 chr1 10642  10643\n",
       "2 chr1 21000  26000\n",
       "3 chr1 40000  47200\n",
       "4 chr1 46000 111250\n",
       "5 chr1 50941  51053\n",
       "6 chr1 54665  54666\n",
       "  V4                                                            \n",
       "1 <BND>                                                         \n",
       "2 <DEL>                                                         \n",
       "3 <DUP>                                                         \n",
       "4 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>\n",
       "5 <DEL>                                                         \n",
       "6 <INS>                                                         \n",
       "  V5                     \n",
       "1 gnomAD-SV_v2.1_BND_1_1 \n",
       "2 gnomAD-SV_v2.1_DEL_1_1 \n",
       "3 gnomAD-SV_v2.1_DUP_1_1 \n",
       "4 gnomAD-SV_v2.1_MCNV_1_1\n",
       "5 gnomAD-SV_v2.1_DEL_1_2 \n",
       "6 gnomAD-SV_v2.1_INS_1_1 \n",
       "  V6                                                                            \n",
       "1 0.006786                                                                      \n",
       "2 0.016073                                                                      \n",
       "3 0.068963                                                                      \n",
       "4 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\n",
       "5 0.004052                                                                      \n",
       "6 9.4e-05                                                                       "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gnomad = read.table('gnomad_v2.1_sv.sites.sorted.tsv', as.is=TRUE)\n",
    "head(gnomad)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 6</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>chr</th><th scope=col>start</th><th scope=col>end</th><th scope=col>type</th><th scope=col>id</th><th scope=col>af</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>chr1</td><td>10642</td><td> 10643</td><td><span style=white-space:pre-wrap>&lt;BND&gt;                                                         </span></td><td>gnomAD-SV_v2.1_BND_1_1 </td><td><span style=white-space:pre-wrap>0.006786                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>2</th><td>chr1</td><td>21000</td><td> 26000</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_1 </td><td><span style=white-space:pre-wrap>0.016073                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>3</th><td>chr1</td><td>40000</td><td> 47200</td><td><span style=white-space:pre-wrap>&lt;DUP&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DUP_1_1 </td><td><span style=white-space:pre-wrap>0.068963                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>4</th><td>chr1</td><td>46000</td><td>111250</td><td>&lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt;</td><td>gnomAD-SV_v2.1_MCNV_1_1</td><td>0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>chr1</td><td>50941</td><td> 51053</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_2 </td><td><span style=white-space:pre-wrap>0.004052                                                                      </span></td></tr>\n",
       "\t<tr><th scope=row>6</th><td>chr1</td><td>54665</td><td> 54666</td><td><span style=white-space:pre-wrap>&lt;INS&gt;                                                         </span></td><td>gnomAD-SV_v2.1_INS_1_1 </td><td><span style=white-space:pre-wrap>9.4e-05                                                                       </span></td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 6\n",
       "\\begin{tabular}{r|llllll}\n",
       "  & chr & start & end & type & id & af\\\\\n",
       "  & <chr> & <int> & <int> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & chr1 & 10642 &  10643 & <BND>                                                          & gnomAD-SV\\_v2.1\\_BND\\_1\\_1  & 0.006786                                                                      \\\\\n",
       "\t2 & chr1 & 21000 &  26000 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_1  & 0.016073                                                                      \\\\\n",
       "\t3 & chr1 & 40000 &  47200 & <DUP>                                                          & gnomAD-SV\\_v2.1\\_DUP\\_1\\_1  & 0.068963                                                                      \\\\\n",
       "\t4 & chr1 & 46000 & 111250 & <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8> & gnomAD-SV\\_v2.1\\_MCNV\\_1\\_1 & 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\\\\\n",
       "\t5 & chr1 & 50941 &  51053 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_2  & 0.004052                                                                      \\\\\n",
       "\t6 & chr1 & 54665 &  54666 & <INS>                                                          & gnomAD-SV\\_v2.1\\_INS\\_1\\_1  & 9.4e-05                                                                       \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 6\n",
       "\n",
       "| <!--/--> | chr &lt;chr&gt; | start &lt;int&gt; | end &lt;int&gt; | type &lt;chr&gt; | id &lt;chr&gt; | af &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|\n",
       "| 1 | chr1 | 10642 |  10643 | &lt;BND&gt;                                                          | gnomAD-SV_v2.1_BND_1_1  | 0.006786                                                                       |\n",
       "| 2 | chr1 | 21000 |  26000 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_1  | 0.016073                                                                       |\n",
       "| 3 | chr1 | 40000 |  47200 | &lt;DUP&gt;                                                          | gnomAD-SV_v2.1_DUP_1_1  | 0.068963                                                                       |\n",
       "| 4 | chr1 | 46000 | 111250 | &lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt; | gnomAD-SV_v2.1_MCNV_1_1 | 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 |\n",
       "| 5 | chr1 | 50941 |  51053 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_2  | 0.004052                                                                       |\n",
       "| 6 | chr1 | 54665 |  54666 | &lt;INS&gt;                                                          | gnomAD-SV_v2.1_INS_1_1  | 9.4e-05                                                                        |\n",
       "\n"
      ],
      "text/plain": [
       "  chr  start end   \n",
       "1 chr1 10642  10643\n",
       "2 chr1 21000  26000\n",
       "3 chr1 40000  47200\n",
       "4 chr1 46000 111250\n",
       "5 chr1 50941  51053\n",
       "6 chr1 54665  54666\n",
       "  type                                                          \n",
       "1 <BND>                                                         \n",
       "2 <DEL>                                                         \n",
       "3 <DUP>                                                         \n",
       "4 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>\n",
       "5 <DEL>                                                         \n",
       "6 <INS>                                                         \n",
       "  id                     \n",
       "1 gnomAD-SV_v2.1_BND_1_1 \n",
       "2 gnomAD-SV_v2.1_DEL_1_1 \n",
       "3 gnomAD-SV_v2.1_DUP_1_1 \n",
       "4 gnomAD-SV_v2.1_MCNV_1_1\n",
       "5 gnomAD-SV_v2.1_DEL_1_2 \n",
       "6 gnomAD-SV_v2.1_INS_1_1 \n",
       "  af                                                                            \n",
       "1 0.006786                                                                      \n",
       "2 0.016073                                                                      \n",
       "3 0.068963                                                                      \n",
       "4 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\n",
       "5 0.004052                                                                      \n",
       "6 9.4e-05                                                                       "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "colnames(gnomad) = c('chr', 'start', 'end', 'type', 'id', 'af')\n",
    "head(gnomad)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "\n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <BND> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         52604 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          <CN=0>,<CN=1>,<CN=2> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            87 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   <CN=0>,<CN=1>,<CN=2>,<CN=3> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            71 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            95 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           206 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                              <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           183 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           132 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            83 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            76 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            44 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            22 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            13 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            15 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            11 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            10 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             6 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             9 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             4 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             5 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             4 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             2 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             3 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             2 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             2 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             3 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             2 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54>,<CN=55>,<CN=56>,<CN=57> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54>,<CN=55>,<CN=56>,<CN=57>,<CN=58>,<CN=59>,<CN=60>,<CN=61>,<CN=62>,<CN=63>,<CN=64>,<CN=65>,<CN=66>,<CN=67>,<CN=68>,<CN=69>,<CN=70>,<CN=71>,<CN=72>,<CN=73>,<CN=74>,<CN=75> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54>,<CN=55>,<CN=56>,<CN=57>,<CN=58>,<CN=59>,<CN=60>,<CN=61>,<CN=62>,<CN=63>,<CN=64>,<CN=65>,<CN=66>,<CN=67>,<CN=68>,<CN=69>,<CN=70>,<CN=71>,<CN=72>,<CN=73>,<CN=74>,<CN=75>,<CN=76>,<CN=77>,<CN=78>,<CN=79>,<CN=80>,<CN=81>,<CN=82>,<CN=83>,<CN=84>,<CN=85>,<CN=86>,<CN=87>,<CN=88>,<CN=89>,<CN=90>,<CN=91>,<CN=92>,<CN=93>,<CN=94>,<CN=95>,<CN=96>,<CN=97>,<CN=98>,<CN=99> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                   <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54>,<CN=55>,<CN=56>,<CN=57>,<CN=58>,<CN=59>,<CN=60>,<CN=61>,<CN=62>,<CN=63>,<CN=64>,<CN=65>,<CN=66>,<CN=67>,<CN=68>,<CN=69>,<CN=70>,<CN=71>,<CN=72>,<CN=73>,<CN=74>,<CN=75>,<CN=76>,<CN=77>,<CN=78>,<CN=79>,<CN=80>,<CN=81>,<CN=82>,<CN=83>,<CN=84>,<CN=85>,<CN=86>,<CN=87>,<CN=88>,<CN=89>,<CN=90>,<CN=91>,<CN=92>,<CN=93>,<CN=94>,<CN=95>,<CN=96>,<CN=97>,<CN=98>,<CN=99>,<CN=100>,<CN=101>,<CN=102>,<CN=103>,<CN=104>,<CN=105>,<CN=106>,<CN=107>,<CN=108>,<CN=109>,<CN=110>,<CN=111>,<CN=112>,<CN=113>,<CN=114>,<CN=115>,<CN=116>,<CN=117>,<CN=118>,<CN=119>,<CN=120>,<CN=121>,<CN=122>,<CN=123>,<CN=124>,<CN=125>,<CN=126>,<CN=127>,<CN=128>,<CN=129>,<CN=130>,<CN=131>,<CN=132>,<CN=133>,<CN=134>,<CN=135>,<CN=136>,<CN=137>,<CN=138>,<CN=139>,<CN=140>,<CN=141>,<CN=142>,<CN=143>,<CN=144>,<CN=145>,<CN=146>,<CN=147>,<CN=148>,<CN=149>,<CN=150>,<CN=151>,<CN=152>,<CN=153>,<CN=154>,<CN=155>,<CN=156>,<CN=157>,<CN=158>,<CN=159>,<CN=160>,<CN=161>,<CN=162>,<CN=163>,<CN=164>,<CN=165>,<CN=166>,<CN=167>,<CN=168>,<CN=169>,<CN=170>,<CN=171>,<CN=172>,<CN=173>,<CN=174>,<CN=175>,<CN=176>,<CN=177>,<CN=178>,<CN=179>,<CN=180>,<CN=181>,<CN=182>,<CN=183>,<CN=184>,<CN=185> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "<CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>,<CN=9>,<CN=10>,<CN=11>,<CN=12>,<CN=13>,<CN=14>,<CN=15>,<CN=16>,<CN=17>,<CN=18>,<CN=19>,<CN=20>,<CN=21>,<CN=22>,<CN=23>,<CN=24>,<CN=25>,<CN=26>,<CN=27>,<CN=28>,<CN=29>,<CN=30>,<CN=31>,<CN=32>,<CN=33>,<CN=34>,<CN=35>,<CN=36>,<CN=37>,<CN=38>,<CN=39>,<CN=40>,<CN=41>,<CN=42>,<CN=43>,<CN=44>,<CN=45>,<CN=46>,<CN=47>,<CN=48>,<CN=49>,<CN=50>,<CN=51>,<CN=52>,<CN=53>,<CN=54>,<CN=55>,<CN=56>,<CN=57>,<CN=58>,<CN=59>,<CN=60>,<CN=61>,<CN=62>,<CN=63>,<CN=64>,<CN=65>,<CN=66>,<CN=67>,<CN=68>,<CN=69>,<CN=70>,<CN=71>,<CN=72>,<CN=73>,<CN=74>,<CN=75>,<CN=76>,<CN=77>,<CN=78>,<CN=79>,<CN=80>,<CN=81>,<CN=82>,<CN=83>,<CN=84>,<CN=85>,<CN=86>,<CN=87>,<CN=88>,<CN=89>,<CN=90>,<CN=91>,<CN=92>,<CN=93>,<CN=94>,<CN=95>,<CN=96>,<CN=97>,<CN=98>,<CN=99>,<CN=100>,<CN=101>,<CN=102>,<CN=103>,<CN=104>,<CN=105>,<CN=106>,<CN=107>,<CN=108>,<CN=109>,<CN=110>,<CN=111>,<CN=112>,<CN=113>,<CN=114>,<CN=115>,<CN=116>,<CN=117>,<CN=118>,<CN=119>,<CN=120>,<CN=121>,<CN=122>,<CN=123>,<CN=124>,<CN=125>,<CN=126>,<CN=127>,<CN=128>,<CN=129>,<CN=130>,<CN=131>,<CN=132>,<CN=133>,<CN=134>,<CN=135>,<CN=136>,<CN=137>,<CN=138>,<CN=139>,<CN=140>,<CN=141>,<CN=142>,<CN=143>,<CN=144>,<CN=145>,<CN=146>,<CN=147>,<CN=148>,<CN=149>,<CN=150>,<CN=151>,<CN=152>,<CN=153>,<CN=154>,<CN=155>,<CN=156>,<CN=157>,<CN=158>,<CN=159>,<CN=160>,<CN=161>,<CN=162>,<CN=163>,<CN=164>,<CN=165>,<CN=166>,<CN=167>,<CN=168>,<CN=169>,<CN=170>,<CN=171>,<CN=172>,<CN=173>,<CN=174>,<CN=175>,<CN=176>,<CN=177>,<CN=178>,<CN=179>,<CN=180>,<CN=181>,<CN=182>,<CN=183>,<CN=184>,<CN=185>,<CN=186>,<CN=187>,<CN=188>,<CN=189>,<CN=190>,<CN=191>,<CN=192>,<CN=193>,<CN=194>,<CN=195>,<CN=196>,<CN=197>,<CN=198>,<CN=199>,<CN=200>,<CN=201>,<CN=202>,<CN=203>,<CN=204>,<CN=205>,<CN=206>,<CN=207>,<CN=208>,<CN=209>,<CN=210>,<CN=211>,<CN=212>,<CN=213>,<CN=214>,<CN=215>,<CN=216>,<CN=217>,<CN=218>,<CN=219>,<CN=220>,<CN=221>,<CN=222>,<CN=223>,<CN=224>,<CN=225>,<CN=226>,<CN=227>,<CN=228>,<CN=229>,<CN=230>,<CN=231>,<CN=232>,<CN=233>,<CN=234>,<CN=235>,<CN=236>,<CN=237>,<CN=238>,<CN=239>,<CN=240>,<CN=241>,<CN=242>,<CN=243>,<CN=244>,<CN=245>,<CN=246>,<CN=247>,<CN=248>,<CN=249>,<CN=250>,<CN=251>,<CN=252>,<CN=253>,<CN=254>,<CN=255>,<CN=256>,<CN=257>,<CN=258>,<CN=259>,<CN=260>,<CN=261>,<CN=262>,<CN=263>,<CN=264>,<CN=265>,<CN=266>,<CN=267>,<CN=268>,<CN=269>,<CN=270>,<CN=271>,<CN=272>,<CN=273>,<CN=274>,<CN=275>,<CN=276>,<CN=277>,<CN=278>,<CN=279>,<CN=280>,<CN=281>,<CN=282>,<CN=283>,<CN=284>,<CN=285>,<CN=286>,<CN=287>,<CN=288>,<CN=289>,<CN=290>,<CN=291>,<CN=292> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             1 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CPX> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          4778 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <CTX> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                             8 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <DEL> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                        169635 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <DUP> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         49571 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  <INS:ME:ALU> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         60475 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                <INS:ME:LINE1> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         10018 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  <INS:ME:SVA> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                          6417 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      <INS:ME> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           672 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <INS> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         31443 \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                         <INV> \n",
       "                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                           748 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "table(gnomad$type)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 7</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>chr</th><th scope=col>start</th><th scope=col>end</th><th scope=col>type</th><th scope=col>id</th><th scope=col>af</th><th scope=col>svtype</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>chr1</td><td>10642</td><td> 10643</td><td><span style=white-space:pre-wrap>&lt;BND&gt;                                                         </span></td><td>gnomAD-SV_v2.1_BND_1_1 </td><td><span style=white-space:pre-wrap>0.006786                                                                      </span></td><td>&lt;BND&gt;</td></tr>\n",
       "\t<tr><th scope=row>2</th><td>chr1</td><td>21000</td><td> 26000</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_1 </td><td><span style=white-space:pre-wrap>0.016073                                                                      </span></td><td><span style=white-space:pre-wrap>DEL  </span></td></tr>\n",
       "\t<tr><th scope=row>3</th><td>chr1</td><td>40000</td><td> 47200</td><td><span style=white-space:pre-wrap>&lt;DUP&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DUP_1_1 </td><td><span style=white-space:pre-wrap>0.068963                                                                      </span></td><td><span style=white-space:pre-wrap>DUP  </span></td></tr>\n",
       "\t<tr><th scope=row>4</th><td>chr1</td><td>46000</td><td>111250</td><td>&lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt;</td><td>gnomAD-SV_v2.1_MCNV_1_1</td><td>0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05</td><td><span style=white-space:pre-wrap>CNV  </span></td></tr>\n",
       "\t<tr><th scope=row>5</th><td>chr1</td><td>50941</td><td> 51053</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_2 </td><td><span style=white-space:pre-wrap>0.004052                                                                      </span></td><td><span style=white-space:pre-wrap>DEL  </span></td></tr>\n",
       "\t<tr><th scope=row>6</th><td>chr1</td><td>54665</td><td> 54666</td><td><span style=white-space:pre-wrap>&lt;INS&gt;                                                         </span></td><td>gnomAD-SV_v2.1_INS_1_1 </td><td><span style=white-space:pre-wrap>9.4e-05                                                                       </span></td><td><span style=white-space:pre-wrap>INS  </span></td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 7\n",
       "\\begin{tabular}{r|lllllll}\n",
       "  & chr & start & end & type & id & af & svtype\\\\\n",
       "  & <chr> & <int> & <int> & <chr> & <chr> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & chr1 & 10642 &  10643 & <BND>                                                          & gnomAD-SV\\_v2.1\\_BND\\_1\\_1  & 0.006786                                                                       & <BND>\\\\\n",
       "\t2 & chr1 & 21000 &  26000 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_1  & 0.016073                                                                       & DEL  \\\\\n",
       "\t3 & chr1 & 40000 &  47200 & <DUP>                                                          & gnomAD-SV\\_v2.1\\_DUP\\_1\\_1  & 0.068963                                                                       & DUP  \\\\\n",
       "\t4 & chr1 & 46000 & 111250 & <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8> & gnomAD-SV\\_v2.1\\_MCNV\\_1\\_1 & 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 & CNV  \\\\\n",
       "\t5 & chr1 & 50941 &  51053 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_2  & 0.004052                                                                       & DEL  \\\\\n",
       "\t6 & chr1 & 54665 &  54666 & <INS>                                                          & gnomAD-SV\\_v2.1\\_INS\\_1\\_1  & 9.4e-05                                                                        & INS  \\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 7\n",
       "\n",
       "| <!--/--> | chr &lt;chr&gt; | start &lt;int&gt; | end &lt;int&gt; | type &lt;chr&gt; | id &lt;chr&gt; | af &lt;chr&gt; | svtype &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|---|---|\n",
       "| 1 | chr1 | 10642 |  10643 | &lt;BND&gt;                                                          | gnomAD-SV_v2.1_BND_1_1  | 0.006786                                                                       | &lt;BND&gt; |\n",
       "| 2 | chr1 | 21000 |  26000 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_1  | 0.016073                                                                       | DEL   |\n",
       "| 3 | chr1 | 40000 |  47200 | &lt;DUP&gt;                                                          | gnomAD-SV_v2.1_DUP_1_1  | 0.068963                                                                       | DUP   |\n",
       "| 4 | chr1 | 46000 | 111250 | &lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt; | gnomAD-SV_v2.1_MCNV_1_1 | 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 | CNV   |\n",
       "| 5 | chr1 | 50941 |  51053 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_2  | 0.004052                                                                       | DEL   |\n",
       "| 6 | chr1 | 54665 |  54666 | &lt;INS&gt;                                                          | gnomAD-SV_v2.1_INS_1_1  | 9.4e-05                                                                        | INS   |\n",
       "\n"
      ],
      "text/plain": [
       "  chr  start end   \n",
       "1 chr1 10642  10643\n",
       "2 chr1 21000  26000\n",
       "3 chr1 40000  47200\n",
       "4 chr1 46000 111250\n",
       "5 chr1 50941  51053\n",
       "6 chr1 54665  54666\n",
       "  type                                                          \n",
       "1 <BND>                                                         \n",
       "2 <DEL>                                                         \n",
       "3 <DUP>                                                         \n",
       "4 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>\n",
       "5 <DEL>                                                         \n",
       "6 <INS>                                                         \n",
       "  id                     \n",
       "1 gnomAD-SV_v2.1_BND_1_1 \n",
       "2 gnomAD-SV_v2.1_DEL_1_1 \n",
       "3 gnomAD-SV_v2.1_DUP_1_1 \n",
       "4 gnomAD-SV_v2.1_MCNV_1_1\n",
       "5 gnomAD-SV_v2.1_DEL_1_2 \n",
       "6 gnomAD-SV_v2.1_INS_1_1 \n",
       "  af                                                                            \n",
       "1 0.006786                                                                      \n",
       "2 0.016073                                                                      \n",
       "3 0.068963                                                                      \n",
       "4 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\n",
       "5 0.004052                                                                      \n",
       "6 9.4e-05                                                                       \n",
       "  svtype\n",
       "1 <BND> \n",
       "2 DEL   \n",
       "3 DUP   \n",
       "4 CNV   \n",
       "5 DEL   \n",
       "6 INS   "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gnomad = gnomad %>% mutate(svtype=ifelse(grepl('CN', type), 'CNV', type),\n",
    "                           svtype=ifelse(grepl('INS', type), 'INS', svtype),\n",
    "                           svtype=ifelse(grepl('DEL', type), 'DEL', svtype),\n",
    "                           svtype=ifelse(grepl('DUP', type), 'DUP', svtype),\n",
    "                           svtype=ifelse(grepl('INV', type), 'INV', svtype))\n",
    "head(gnomad)                "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 8</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>chr</th><th scope=col>start</th><th scope=col>end</th><th scope=col>type</th><th scope=col>id</th><th scope=col>af</th><th scope=col>svtype</th><th scope=col>af.max</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>chr1</td><td>10642</td><td> 10643</td><td><span style=white-space:pre-wrap>&lt;BND&gt;                                                         </span></td><td>gnomAD-SV_v2.1_BND_1_1 </td><td><span style=white-space:pre-wrap>0.006786                                                                      </span></td><td>&lt;BND&gt;</td><td>0.006786</td></tr>\n",
       "\t<tr><th scope=row>2</th><td>chr1</td><td>21000</td><td> 26000</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_1 </td><td><span style=white-space:pre-wrap>0.016073                                                                      </span></td><td><span style=white-space:pre-wrap>DEL  </span></td><td>0.016073</td></tr>\n",
       "\t<tr><th scope=row>3</th><td>chr1</td><td>40000</td><td> 47200</td><td><span style=white-space:pre-wrap>&lt;DUP&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DUP_1_1 </td><td><span style=white-space:pre-wrap>0.068963                                                                      </span></td><td><span style=white-space:pre-wrap>DUP  </span></td><td>0.068963</td></tr>\n",
       "\t<tr><th scope=row>4</th><td>chr1</td><td>46000</td><td>111250</td><td>&lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt;</td><td>gnomAD-SV_v2.1_MCNV_1_1</td><td>0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05</td><td><span style=white-space:pre-wrap>CNV  </span></td><td>0.349762</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>chr1</td><td>50941</td><td> 51053</td><td><span style=white-space:pre-wrap>&lt;DEL&gt;                                                         </span></td><td>gnomAD-SV_v2.1_DEL_1_2 </td><td><span style=white-space:pre-wrap>0.004052                                                                      </span></td><td><span style=white-space:pre-wrap>DEL  </span></td><td>0.004052</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>chr1</td><td>54665</td><td> 54666</td><td><span style=white-space:pre-wrap>&lt;INS&gt;                                                         </span></td><td>gnomAD-SV_v2.1_INS_1_1 </td><td><span style=white-space:pre-wrap>9.4e-05                                                                       </span></td><td><span style=white-space:pre-wrap>INS  </span></td><td>0.000094</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 8\n",
       "\\begin{tabular}{r|llllllll}\n",
       "  & chr & start & end & type & id & af & svtype & af.max\\\\\n",
       "  & <chr> & <int> & <int> & <chr> & <chr> & <chr> & <chr> & <dbl>\\\\\n",
       "\\hline\n",
       "\t1 & chr1 & 10642 &  10643 & <BND>                                                          & gnomAD-SV\\_v2.1\\_BND\\_1\\_1  & 0.006786                                                                       & <BND> & 0.006786\\\\\n",
       "\t2 & chr1 & 21000 &  26000 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_1  & 0.016073                                                                       & DEL   & 0.016073\\\\\n",
       "\t3 & chr1 & 40000 &  47200 & <DUP>                                                          & gnomAD-SV\\_v2.1\\_DUP\\_1\\_1  & 0.068963                                                                       & DUP   & 0.068963\\\\\n",
       "\t4 & chr1 & 46000 & 111250 & <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8> & gnomAD-SV\\_v2.1\\_MCNV\\_1\\_1 & 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 & CNV   & 0.349762\\\\\n",
       "\t5 & chr1 & 50941 &  51053 & <DEL>                                                          & gnomAD-SV\\_v2.1\\_DEL\\_1\\_2  & 0.004052                                                                       & DEL   & 0.004052\\\\\n",
       "\t6 & chr1 & 54665 &  54666 & <INS>                                                          & gnomAD-SV\\_v2.1\\_INS\\_1\\_1  & 9.4e-05                                                                        & INS   & 0.000094\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 8\n",
       "\n",
       "| <!--/--> | chr &lt;chr&gt; | start &lt;int&gt; | end &lt;int&gt; | type &lt;chr&gt; | id &lt;chr&gt; | af &lt;chr&gt; | svtype &lt;chr&gt; | af.max &lt;dbl&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|\n",
       "| 1 | chr1 | 10642 |  10643 | &lt;BND&gt;                                                          | gnomAD-SV_v2.1_BND_1_1  | 0.006786                                                                       | &lt;BND&gt; | 0.006786 |\n",
       "| 2 | chr1 | 21000 |  26000 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_1  | 0.016073                                                                       | DEL   | 0.016073 |\n",
       "| 3 | chr1 | 40000 |  47200 | &lt;DUP&gt;                                                          | gnomAD-SV_v2.1_DUP_1_1  | 0.068963                                                                       | DUP   | 0.068963 |\n",
       "| 4 | chr1 | 46000 | 111250 | &lt;CN=0&gt;,&lt;CN=1&gt;,&lt;CN=2&gt;,&lt;CN=3&gt;,&lt;CN=4&gt;,&lt;CN=5&gt;,&lt;CN=6&gt;,&lt;CN=7&gt;,&lt;CN=8&gt; | gnomAD-SV_v2.1_MCNV_1_1 | 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05 | CNV   | 0.349762 |\n",
       "| 5 | chr1 | 50941 |  51053 | &lt;DEL&gt;                                                          | gnomAD-SV_v2.1_DEL_1_2  | 0.004052                                                                       | DEL   | 0.004052 |\n",
       "| 6 | chr1 | 54665 |  54666 | &lt;INS&gt;                                                          | gnomAD-SV_v2.1_INS_1_1  | 9.4e-05                                                                        | INS   | 0.000094 |\n",
       "\n"
      ],
      "text/plain": [
       "  chr  start end   \n",
       "1 chr1 10642  10643\n",
       "2 chr1 21000  26000\n",
       "3 chr1 40000  47200\n",
       "4 chr1 46000 111250\n",
       "5 chr1 50941  51053\n",
       "6 chr1 54665  54666\n",
       "  type                                                          \n",
       "1 <BND>                                                         \n",
       "2 <DEL>                                                         \n",
       "3 <DUP>                                                         \n",
       "4 <CN=0>,<CN=1>,<CN=2>,<CN=3>,<CN=4>,<CN=5>,<CN=6>,<CN=7>,<CN=8>\n",
       "5 <DEL>                                                         \n",
       "6 <INS>                                                         \n",
       "  id                     \n",
       "1 gnomAD-SV_v2.1_BND_1_1 \n",
       "2 gnomAD-SV_v2.1_DEL_1_1 \n",
       "3 gnomAD-SV_v2.1_DUP_1_1 \n",
       "4 gnomAD-SV_v2.1_MCNV_1_1\n",
       "5 gnomAD-SV_v2.1_DEL_1_2 \n",
       "6 gnomAD-SV_v2.1_INS_1_1 \n",
       "  af                                                                            \n",
       "1 0.006786                                                                      \n",
       "2 0.016073                                                                      \n",
       "3 0.068963                                                                      \n",
       "4 0.09479,0.270108,0.349762,0.177452,0.083632,0.019986,0.003784,0.000388,9.7e-05\n",
       "5 0.004052                                                                      \n",
       "6 9.4e-05                                                                       \n",
       "  svtype af.max  \n",
       "1 <BND>  0.006786\n",
       "2 DEL    0.016073\n",
       "3 DUP    0.068963\n",
       "4 CNV    0.349762\n",
       "5 DEL    0.004052\n",
       "6 INS    0.000094"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "gnomad$af.max = unlist(lapply(strsplit(gnomad$af, split=','), function(x) max(as.numeric(x))))\n",
    "head(gnomad)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 5</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>chr</th><th scope=col>start</th><th scope=col>end</th><th scope=col>type</th><th scope=col>variant_id</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>21</td><td>      1</td><td>23875167</td><td>DEL</td><td>nssv1613109 </td></tr>\n",
       "\t<tr><th scope=row>2</th><td>21</td><td>      1</td><td> 9411300</td><td>DEL</td><td>nssv3000025 </td></tr>\n",
       "\t<tr><th scope=row>3</th><td>21</td><td>9410851</td><td> 9483551</td><td>DUP</td><td>essv25939335</td></tr>\n",
       "\t<tr><th scope=row>4</th><td>21</td><td>9411194</td><td> 9841977</td><td>DEL</td><td>nssv14267892</td></tr>\n",
       "\t<tr><th scope=row>5</th><td>21</td><td>9411194</td><td> 9841977</td><td>DEL</td><td>nssv14267893</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>21</td><td>9411194</td><td> 9841977</td><td>DEL</td><td>nssv14267894</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 5\n",
       "\\begin{tabular}{r|lllll}\n",
       "  & chr & start & end & type & variant\\_id\\\\\n",
       "  & <int> & <int> & <int> & <chr> & <chr>\\\\\n",
       "\\hline\n",
       "\t1 & 21 &       1 & 23875167 & DEL & nssv1613109 \\\\\n",
       "\t2 & 21 &       1 &  9411300 & DEL & nssv3000025 \\\\\n",
       "\t3 & 21 & 9410851 &  9483551 & DUP & essv25939335\\\\\n",
       "\t4 & 21 & 9411194 &  9841977 & DEL & nssv14267892\\\\\n",
       "\t5 & 21 & 9411194 &  9841977 & DEL & nssv14267893\\\\\n",
       "\t6 & 21 & 9411194 &  9841977 & DEL & nssv14267894\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 5\n",
       "\n",
       "| <!--/--> | chr &lt;int&gt; | start &lt;int&gt; | end &lt;int&gt; | type &lt;chr&gt; | variant_id &lt;chr&gt; |\n",
       "|---|---|---|---|---|---|\n",
       "| 1 | 21 |       1 | 23875167 | DEL | nssv1613109  |\n",
       "| 2 | 21 |       1 |  9411300 | DEL | nssv3000025  |\n",
       "| 3 | 21 | 9410851 |  9483551 | DUP | essv25939335 |\n",
       "| 4 | 21 | 9411194 |  9841977 | DEL | nssv14267892 |\n",
       "| 5 | 21 | 9411194 |  9841977 | DEL | nssv14267893 |\n",
       "| 6 | 21 | 9411194 |  9841977 | DEL | nssv14267894 |\n",
       "\n"
      ],
      "text/plain": [
       "  chr start   end      type variant_id  \n",
       "1 21        1 23875167 DEL  nssv1613109 \n",
       "2 21        1  9411300 DEL  nssv3000025 \n",
       "3 21  9410851  9483551 DUP  essv25939335\n",
       "4 21  9411194  9841977 DEL  nssv14267892\n",
       "5 21  9411194  9841977 DEL  nssv14267893\n",
       "6 21  9411194  9841977 DEL  nssv14267894"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "vars.df = read.table('all_variants_chr21_grch37.tsv', as.is=TRUE, header=TRUE)\n",
    "head(vars.df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "\n",
       "   CNV    DEL    DUP    INS    INV \n",
       " 24671 219242  48684 118236   9147 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "table(vars.df$type)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n",
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n",
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n",
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n",
      "`summarise()` ungrouping output (override with `.groups` argument)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "vars.af = lapply(unique(vars.df$type), function(tt){\n",
    "    vars.gr = vars.df %>% mutate(chr=paste0('chr', chr)) %>% filter(type==tt) %>% makeGRangesFromDataFrame(keep.extra.columns = TRUE)\n",
    "    gnomad.gr = gnomad %>% filter(svtype==tt) %>% makeGRangesFromDataFrame(keep.extra.columns = TRUE)\n",
    "    findOverlaps(vars.gr, gnomad.gr) %>% as.data.frame %>%\n",
    "    mutate(af=gnomad.gr$af.max[subjectHits], variant_id=vars.gr$variant_id[queryHits],\n",
    "          qw=width(vars.gr)[queryHits], sw=width(gnomad.gr)[subjectHits],\n",
    "          qsw=width(pintersect(vars.gr[queryHits], gnomad.gr[subjectHits]))) %>%\n",
    "    filter(qsw>.9 * qw, qsw > .9 * sw) %>% \n",
    "    group_by(variant_id) %>% summarize(af=max(af))\n",
    "}) %>% bind_rows()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A tibble: 6 × 2</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>variant_id</th><th scope=col>af</th></tr>\n",
       "\t<tr><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td>essv100052</td><td>0.298128</td></tr>\n",
       "\t<tr><td>essv100397</td><td>0.298128</td></tr>\n",
       "\t<tr><td>essv100590</td><td>0.298128</td></tr>\n",
       "\t<tr><td>essv101052</td><td>0.298128</td></tr>\n",
       "\t<tr><td>essv101179</td><td>0.298128</td></tr>\n",
       "\t<tr><td>essv101524</td><td>0.298128</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A tibble: 6 × 2\n",
       "\\begin{tabular}{ll}\n",
       " variant\\_id & af\\\\\n",
       " <chr> & <dbl>\\\\\n",
       "\\hline\n",
       "\t essv100052 & 0.298128\\\\\n",
       "\t essv100397 & 0.298128\\\\\n",
       "\t essv100590 & 0.298128\\\\\n",
       "\t essv101052 & 0.298128\\\\\n",
       "\t essv101179 & 0.298128\\\\\n",
       "\t essv101524 & 0.298128\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A tibble: 6 × 2\n",
       "\n",
       "| variant_id &lt;chr&gt; | af &lt;dbl&gt; |\n",
       "|---|---|\n",
       "| essv100052 | 0.298128 |\n",
       "| essv100397 | 0.298128 |\n",
       "| essv100590 | 0.298128 |\n",
       "| essv101052 | 0.298128 |\n",
       "| essv101179 | 0.298128 |\n",
       "| essv101524 | 0.298128 |\n",
       "\n"
      ],
      "text/plain": [
       "  variant_id af      \n",
       "1 essv100052 0.298128\n",
       "2 essv100397 0.298128\n",
       "3 essv100590 0.298128\n",
       "4 essv101052 0.298128\n",
       "5 essv101179 0.298128\n",
       "6 essv101524 0.298128"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "head(vars.af)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 43,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. \n",
       "0.000046 0.128653 0.384190 0.350648 0.493359 0.999515 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "summary(vars.af$af)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": [
    "write.table(vars.af, file='gnomad-af-variants-chr21.tsv', quote=FALSE, row.names=FALSE, sep='\\t')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```\n",
    "dx upload gnomad-af-variants-chr21.tsv --path chr-subset-genes-variants/\n",
    "dx upload gnomad-variants-af.ipynb --path chr-subset-genes-variants/\n",
    "```"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}