diff --git a/pysequila/aoh.ipynb b/pysequila/aoh.ipynb
index fa3ee7a..47647aa 100644
--- a/pysequila/aoh.ipynb
+++ b/pysequila/aoh.ipynb
@@ -4,7 +4,7 @@
    "cell_type": "code",
    "execution_count": null,
    "metadata": {
-    "collapsed": false,
+    "collapsed": true,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -18,27 +18,38 @@
     "mkdir data\n",
     "cd data\n",
     "    \n",
-    "#wget http://biodatageeks.ii.pw.edu.pl/sequila/data/WES/NA12878.proper.wes.chr1.md.bam\n",
     "wget http://biodatageeks.ii.pw.edu.pl/sequila/data/WES/NA12878.proper.wes.chr1.md.bam\n",
     "wget http://biodatageeks.ii.pw.edu.pl/sequila/data/WES/Homo_sapiens_assembly18.chr1.fasta\n",
-    "wget http://biodatageeks.ii.pw.edu.pl/sequila/data/WES/Homo_sapiens_assembly18.chr1.fasta.fai\n"
+    "wget http://biodatageeks.ii.pw.edu.pl/sequila/data/WES/Homo_sapiens_assembly18.chr1.fasta.fai"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": [
+    "pip install pandas"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": true
+   },
    "outputs": [],
    "source": [
-    "pip install pandas\n",
     "pip install matplotlib"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {
-    "collapsed": false,
+    "collapsed": true,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -59,9 +70,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {
-    "collapsed": false,
+    "collapsed": true,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -75,14 +86,14 @@
     "spark = SparkSession \\\n",
     ".builder \\\n",
     ".appName(f'{app_name}') \\\n",
+    ".config('spark.driver.memory','16g') \\\n",
     ".getOrCreate()"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
    "metadata": {
-    "collapsed": false,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -90,18 +101,7 @@
      "name": "#%%\n"
     }
    },
-   "outputs": [
-    {
-     "data": {
-      "text/plain": [
-       "DataFrame[]"
-      ]
-     },
-     "execution_count": 5,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "from pysequila import SequilaSession\n",
     "ss = SequilaSession(spark)\n",
@@ -113,72 +113,8 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
-   "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Collecting pandas\n",
-      "  Downloading pandas-1.1.4-cp37-cp37m-manylinux1_x86_64.whl (9.5 MB)\n",
-      "\u001b[K     |████████████████████████████████| 9.5 MB 1.6 MB/s eta 0:00:01     |▉                               | 235 kB 451 kB/s eta 0:00:21\n",
-      "\u001b[?25hCollecting numpy>=1.15.4\n",
-      "  Downloading numpy-1.19.3-cp37-cp37m-manylinux2010_x86_64.whl (14.9 MB)\n",
-      "\u001b[K     |████████████████████████████████| 14.9 MB 1.2 MB/s eta 0:00:01\n",
-      "\u001b[?25hRequirement already satisfied: python-dateutil>=2.7.3 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from pandas) (2.8.1)\n",
-      "Collecting pytz>=2017.2\n",
-      "  Downloading pytz-2020.1-py2.py3-none-any.whl (510 kB)\n",
-      "\u001b[K     |████████████████████████████████| 510 kB 3.3 MB/s eta 0:00:01\n",
-      "\u001b[?25hRequirement already satisfied: six>=1.5 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from python-dateutil>=2.7.3->pandas) (1.15.0)\n",
-      "Installing collected packages: numpy, pytz, pandas\n",
-      "Successfully installed numpy-1.19.3 pandas-1.1.4 pytz-2020.1\n",
-      "Note: you may need to restart the kernel to use updated packages.\n"
-     ]
-    }
-   ],
-   "source": []
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 14,
-   "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Collecting matplotlib\n",
-      "  Downloading matplotlib-3.3.2-cp37-cp37m-manylinux1_x86_64.whl (11.6 MB)\n",
-      "\u001b[K     |████████████████████████████████| 11.6 MB 6.2 MB/s eta 0:00:01    |██████████████████████          | 8.0 MB 6.2 MB/s eta 0:00:01\n",
-      "\u001b[?25hCollecting pyparsing!=2.0.4,!=2.1.2,!=2.1.6,>=2.0.3\n",
-      "  Downloading pyparsing-2.4.7-py2.py3-none-any.whl (67 kB)\n",
-      "\u001b[K     |████████████████████████████████| 67 kB 2.1 MB/s eta 0:00:01\n",
-      "\u001b[?25hCollecting kiwisolver>=1.0.1\n",
-      "  Downloading kiwisolver-1.2.0-cp37-cp37m-manylinux1_x86_64.whl (88 kB)\n",
-      "\u001b[K     |████████████████████████████████| 88 kB 2.6 MB/s eta 0:00:01\n",
-      "\u001b[?25hRequirement already satisfied: python-dateutil>=2.1 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from matplotlib) (2.8.1)\n",
-      "Collecting cycler>=0.10\n",
-      "  Downloading cycler-0.10.0-py2.py3-none-any.whl (6.5 kB)\n",
-      "Collecting pillow>=6.2.0\n",
-      "  Downloading Pillow-8.0.1-cp37-cp37m-manylinux1_x86_64.whl (2.2 MB)\n",
-      "\u001b[K     |████████████████████████████████| 2.2 MB 3.9 MB/s eta 0:00:01\n",
-      "\u001b[?25hRequirement already satisfied: numpy>=1.15 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from matplotlib) (1.19.3)\n",
-      "Requirement already satisfied: certifi>=2020.06.20 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from matplotlib) (2020.6.20)\n",
-      "Requirement already satisfied: six>=1.5 in /home/jovyan/venv/pysequila/lib/python3.7/site-packages (from python-dateutil>=2.1->matplotlib) (1.15.0)\n",
-      "Installing collected packages: pyparsing, kiwisolver, cycler, pillow, matplotlib\n",
-      "Successfully installed cycler-0.10.0 kiwisolver-1.2.0 matplotlib-3.3.2 pillow-8.0.1 pyparsing-2.4.7\n",
-      "Note: you may need to restart the kernel to use updated packages.\n"
-     ]
-    }
-   ],
-   "source": []
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 23,
+   "execution_count": null,
    "metadata": {
-    "collapsed": false,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -186,530 +122,7 @@
      "name": "#%%\n"
     }
    },
-   "outputs": [
-    {
-     "data": {
-      "text/html": [
-       "<div>\n",
-       "<style scoped>\n",
-       "    .dataframe tbody tr th:only-of-type {\n",
-       "        vertical-align: middle;\n",
-       "    }\n",
-       "\n",
-       "    .dataframe tbody tr th {\n",
-       "        vertical-align: top;\n",
-       "    }\n",
-       "\n",
-       "    .dataframe thead th {\n",
-       "        text-align: right;\n",
-       "    }\n",
-       "</style>\n",
-       "<table border=\"1\" class=\"dataframe\">\n",
-       "  <thead>\n",
-       "    <tr style=\"text-align: right;\">\n",
-       "      <th></th>\n",
-       "      <th>sample_id</th>\n",
-       "      <th>qname</th>\n",
-       "      <th>flag</th>\n",
-       "      <th>contig</th>\n",
-       "      <th>pos</th>\n",
-       "      <th>pos_start</th>\n",
-       "      <th>pos_end</th>\n",
-       "      <th>mapq</th>\n",
-       "      <th>cigar</th>\n",
-       "      <th>rnext</th>\n",
-       "      <th>pnext</th>\n",
-       "      <th>tlen</th>\n",
-       "      <th>seq</th>\n",
-       "      <th>qual</th>\n",
-       "      <th>tag_AM</th>\n",
-       "      <th>tag_AS</th>\n",
-       "      <th>tag_BC</th>\n",
-       "      <th>tag_BQ</th>\n",
-       "      <th>tag_BZ</th>\n",
-       "      <th>tag_CB</th>\n",
-       "      <th>tag_CC</th>\n",
-       "      <th>tag_CG</th>\n",
-       "      <th>tag_CM</th>\n",
-       "      <th>tag_CO</th>\n",
-       "      <th>tag_CP</th>\n",
-       "      <th>tag_CQ</th>\n",
-       "      <th>tag_CR</th>\n",
-       "      <th>tag_CS</th>\n",
-       "      <th>tag_CT</th>\n",
-       "      <th>tag_CY</th>\n",
-       "      <th>tag_E2</th>\n",
-       "      <th>tag_FI</th>\n",
-       "      <th>tag_FS</th>\n",
-       "      <th>tag_FZ</th>\n",
-       "      <th>tag_H0</th>\n",
-       "      <th>tag_H1</th>\n",
-       "      <th>tag_H2</th>\n",
-       "      <th>tag_HI</th>\n",
-       "      <th>tag_IH</th>\n",
-       "      <th>tag_LB</th>\n",
-       "      <th>tag_MC</th>\n",
-       "      <th>tag_MD</th>\n",
-       "      <th>tag_MI</th>\n",
-       "      <th>tag_MQ</th>\n",
-       "      <th>tag_NH</th>\n",
-       "      <th>tag_NM</th>\n",
-       "      <th>tag_OA</th>\n",
-       "      <th>tag_OC</th>\n",
-       "      <th>tag_OP</th>\n",
-       "      <th>tag_OQ</th>\n",
-       "      <th>tag_OX</th>\n",
-       "      <th>tag_PG</th>\n",
-       "      <th>tag_PQ</th>\n",
-       "      <th>tag_PT</th>\n",
-       "      <th>tag_PU</th>\n",
-       "      <th>tag_Q2</th>\n",
-       "      <th>tag_QT</th>\n",
-       "      <th>tag_QX</th>\n",
-       "      <th>tag_R2</th>\n",
-       "      <th>tag_RG</th>\n",
-       "      <th>tag_RX</th>\n",
-       "      <th>tag_SA</th>\n",
-       "      <th>tag_SM</th>\n",
-       "      <th>tag_TC</th>\n",
-       "      <th>tag_U2</th>\n",
-       "      <th>tag_UQ</th>\n",
-       "    </tr>\n",
-       "  </thead>\n",
-       "  <tbody>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>NA12878</td>\n",
-       "      <td>61DC0AAXX100127:8:58:8295:16397</td>\n",
-       "      <td>99</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>7</td>\n",
-       "      <td>7</td>\n",
-       "      <td>107</td>\n",
-       "      <td>99</td>\n",
-       "      <td>101M</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>27</td>\n",
-       "      <td>101</td>\n",
-       "      <td>AGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTT...</td>\n",
-       "      <td>ADEAEDBCFCEFFFCDEDFGHEGGHEGBIIFIHGHHIHFFIIFFGG...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>CAAGAGCGACAAAGCGGCCAACAGACAATACTAGACAACGAACAGG...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>66G34</td>\n",
-       "      <td>None</td>\n",
-       "      <td>99.0</td>\n",
-       "      <td>None</td>\n",
-       "      <td>1</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>36</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>1</th>\n",
-       "      <td>NA12878</td>\n",
-       "      <td>61DC0AAXX100127:8:35:1455:1743</td>\n",
-       "      <td>1113</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>9</td>\n",
-       "      <td>9</td>\n",
-       "      <td>109</td>\n",
-       "      <td>90</td>\n",
-       "      <td>101M</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>9</td>\n",
-       "      <td>101</td>\n",
-       "      <td>GTCTGTCACCCTTGTAGCCGCTCACGGGAGCTCTCCATGCATTTGG...</td>\n",
-       "      <td>##########################G1/7=@AEHB5=F=DFC@FE...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>CAAACAATGTAGAAATAAACACACAAACCTAAAAATGGTTGAACTA...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>4A7A0T2A2A44G36</td>\n",
-       "      <td>None</td>\n",
-       "      <td>NaN</td>\n",
-       "      <td>None</td>\n",
-       "      <td>6</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>##########################B1+8;&gt;&gt;BB&gt;1&lt;B9BB@&gt;BB...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>34</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>2</th>\n",
-       "      <td>NA12878</td>\n",
-       "      <td>61DC0AAXX100127:8:1:5111:5929</td>\n",
-       "      <td>1177</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>10</td>\n",
-       "      <td>10</td>\n",
-       "      <td>110</td>\n",
-       "      <td>99</td>\n",
-       "      <td>101M</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>10</td>\n",
-       "      <td>101</td>\n",
-       "      <td>TCTATCCCCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGT...</td>\n",
-       "      <td>B?=&lt;=:\"$.,/;;@?&gt;CI@CAIEDIIGGGHDEFBIEGCIFHHGIHE...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>GTCTAAAAAACTAATTAACAAATATTTGTTATAAATGTTTAGGTTG...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>6A56G37</td>\n",
-       "      <td>None</td>\n",
-       "      <td>NaN</td>\n",
-       "      <td>None</td>\n",
-       "      <td>2</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>B&lt;?79:&amp;*0+.8:@=:@C=@=CBCDCBBBB@AD@CDCADCCCCCCC...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>38</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>3</th>\n",
-       "      <td>NA12878</td>\n",
-       "      <td>61DC0AAXX100127:8:47:17866:4686</td>\n",
-       "      <td>153</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>20</td>\n",
-       "      <td>20</td>\n",
-       "      <td>120</td>\n",
-       "      <td>99</td>\n",
-       "      <td>101M</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>20</td>\n",
-       "      <td>101</td>\n",
-       "      <td>TATTATCCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCT...</td>\n",
-       "      <td>###########ICA@GGG=DFDFHC=@GH&gt;HG&lt;IAFBHHFDCHCFC...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>ATGATAAACTAACACTTGTTATATATGTATGGGTTGTGGGCTTCTG...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>5A47G47</td>\n",
-       "      <td>None</td>\n",
-       "      <td>NaN</td>\n",
-       "      <td>None</td>\n",
-       "      <td>2</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>###########B@B@BBB=AB@BB&gt;;?BB:CB&lt;CAC@CCBACCABA...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61DC0.8</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>37</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>4</th>\n",
-       "      <td>NA12878</td>\n",
-       "      <td>61CC3AAXX100125:5:79:19223:6261</td>\n",
-       "      <td>99</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>25</td>\n",
-       "      <td>25</td>\n",
-       "      <td>100</td>\n",
-       "      <td>99</td>\n",
-       "      <td>76M</td>\n",
-       "      <td>MT</td>\n",
-       "      <td>30</td>\n",
-       "      <td>76</td>\n",
-       "      <td>ACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGG...</td>\n",
-       "      <td>ABCADDFCE@GGEHGFHFHHEEIHEEFFII&gt;DEFFFIBDHGIIIII...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>CAACAGACATTTCTAGAGAATGAACGGGTTGCAGGGAAGAGTTTTT...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>48G27</td>\n",
-       "      <td>None</td>\n",
-       "      <td>99.0</td>\n",
-       "      <td>None</td>\n",
-       "      <td>1</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC@CCCCCCCCCCCCCCC...</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61CC3.5</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>61CC3.5</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>None</td>\n",
-       "      <td>29</td>\n",
-       "    </tr>\n",
-       "  </tbody>\n",
-       "</table>\n",
-       "</div>"
-      ],
-      "text/plain": [
-       "  sample_id                            qname  flag contig  pos  pos_start  \\\n",
-       "0   NA12878  61DC0AAXX100127:8:58:8295:16397    99     MT    7          7   \n",
-       "1   NA12878   61DC0AAXX100127:8:35:1455:1743  1113     MT    9          9   \n",
-       "2   NA12878    61DC0AAXX100127:8:1:5111:5929  1177     MT   10         10   \n",
-       "3   NA12878  61DC0AAXX100127:8:47:17866:4686   153     MT   20         20   \n",
-       "4   NA12878  61CC3AAXX100125:5:79:19223:6261    99     MT   25         25   \n",
-       "\n",
-       "   pos_end  mapq cigar rnext  pnext  tlen  \\\n",
-       "0      107    99  101M    MT     27   101   \n",
-       "1      109    90  101M    MT      9   101   \n",
-       "2      110    99  101M    MT     10   101   \n",
-       "3      120    99  101M    MT     20   101   \n",
-       "4      100    99   76M    MT     30    76   \n",
-       "\n",
-       "                                                 seq  \\\n",
-       "0  AGGTCTATCACCCTATTAACCACTCACGGGAGCTCTCCATGCATTT...   \n",
-       "1  GTCTGTCACCCTTGTAGCCGCTCACGGGAGCTCTCCATGCATTTGG...   \n",
-       "2  TCTATCCCCCTATTAACCACTCACGGGAGCTCTCCATGCATTTGGT...   \n",
-       "3  TATTATCCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCT...   \n",
-       "4  ACCACTCACGGGAGCTCTCCATGCATTTGGTATTTTCGTCTGGGGG...   \n",
-       "\n",
-       "                                                qual tag_AM tag_AS tag_BC  \\\n",
-       "0  ADEAEDBCFCEFFFCDEDFGHEGGHEGBIIFIHGHHIHFFIIFFGG...   None   None   None   \n",
-       "1  ##########################G1/7=@AEHB5=F=DFC@FE...   None   None   None   \n",
-       "2  B?=<=:\"$.,/;;@?>CI@CAIEDIIGGGHDEFBIEGCIFHHGIHE...   None   None   None   \n",
-       "3  ###########ICA@GGG=DFDFHC=@GH>HG<IAFBHHFDCHCFC...   None   None   None   \n",
-       "4  ABCADDFCE@GGEHGFHFHHEEIHEEFFII>DEFFFIBDHGIIIII...   None   None   None   \n",
-       "\n",
-       "  tag_BQ tag_BZ tag_CB tag_CC tag_CG tag_CM tag_CO tag_CP tag_CQ tag_CR  \\\n",
-       "0   None   None   None   None   None   None   None   None   None   None   \n",
-       "1   None   None   None   None   None   None   None   None   None   None   \n",
-       "2   None   None   None   None   None   None   None   None   None   None   \n",
-       "3   None   None   None   None   None   None   None   None   None   None   \n",
-       "4   None   None   None   None   None   None   None   None   None   None   \n",
-       "\n",
-       "  tag_CS tag_CT tag_CY                                             tag_E2  \\\n",
-       "0   None   None   None  CAAGAGCGACAAAGCGGCCAACAGACAATACTAGACAACGAACAGG...   \n",
-       "1   None   None   None  CAAACAATGTAGAAATAAACACACAAACCTAAAAATGGTTGAACTA...   \n",
-       "2   None   None   None  GTCTAAAAAACTAATTAACAAATATTTGTTATAAATGTTTAGGTTG...   \n",
-       "3   None   None   None  ATGATAAACTAACACTTGTTATATATGTATGGGTTGTGGGCTTCTG...   \n",
-       "4   None   None   None  CAACAGACATTTCTAGAGAATGAACGGGTTGCAGGGAAGAGTTTTT...   \n",
-       "\n",
-       "  tag_FI tag_FS tag_FZ tag_H0 tag_H1 tag_H2 tag_HI tag_IH tag_LB tag_MC  \\\n",
-       "0   None   None   None   None   None   None   None   None   None   None   \n",
-       "1   None   None   None   None   None   None   None   None   None   None   \n",
-       "2   None   None   None   None   None   None   None   None   None   None   \n",
-       "3   None   None   None   None   None   None   None   None   None   None   \n",
-       "4   None   None   None   None   None   None   None   None   None   None   \n",
-       "\n",
-       "            tag_MD tag_MI  tag_MQ tag_NH  tag_NM tag_OA tag_OC tag_OP  \\\n",
-       "0            66G34   None    99.0   None       1   None   None   None   \n",
-       "1  4A7A0T2A2A44G36   None     NaN   None       6   None   None   None   \n",
-       "2          6A56G37   None     NaN   None       2   None   None   None   \n",
-       "3          5A47G47   None     NaN   None       2   None   None   None   \n",
-       "4            48G27   None    99.0   None       1   None   None   None   \n",
-       "\n",
-       "                                              tag_OQ tag_OX   tag_PG tag_PQ  \\\n",
-       "0  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC...   None  61DC0.8   None   \n",
-       "1  ##########################B1+8;>>BB>1<B9BB@>BB...   None  61DC0.8   None   \n",
-       "2  B<?79:&*0+.8:@=:@C=@=CBCDCBBBB@AD@CDCADCCCCCCC...   None  61DC0.8   None   \n",
-       "3  ###########B@B@BBB=AB@BB>;?BB:CB<CAC@CCBACCABA...   None  61DC0.8   None   \n",
-       "4  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC@CCCCCCCCCCCCCCC...   None  61CC3.5   None   \n",
-       "\n",
-       "  tag_PT tag_PU tag_Q2 tag_QT tag_QX tag_R2   tag_RG tag_RX tag_SA tag_SM  \\\n",
-       "0   None   None   None   None   None   None  61DC0.8   None   None   None   \n",
-       "1   None   None   None   None   None   None  61DC0.8   None   None   None   \n",
-       "2   None   None   None   None   None   None  61DC0.8   None   None   None   \n",
-       "3   None   None   None   None   None   None  61DC0.8   None   None   None   \n",
-       "4   None   None   None   None   None   None  61CC3.5   None   None   None   \n",
-       "\n",
-       "  tag_TC tag_U2  tag_UQ  \n",
-       "0   None   None      36  \n",
-       "1   None   None      34  \n",
-       "2   None   None      38  \n",
-       "3   None   None      37  \n",
-       "4   None   None      29  "
-      ]
-     },
-     "execution_count": 23,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "import pandas as pd\n",
     "pd.options.display.max_columns = None\n",
@@ -720,7 +133,7 @@
    "cell_type": "code",
    "execution_count": null,
    "metadata": {
-    "collapsed": false,
+    "collapsed": true,
     "jupyter": {
      "outputs_hidden": false
     },
@@ -730,27 +143,44 @@
    },
    "outputs": [],
    "source": [
-    "import matplotlib.pyplot as plt\n",
-    "df = ss.sql(f\"SELECT * FROM  pileup('{table_name}', '{sample_id}', '{ref_path}', true) WHERE contig='1' \").toPandas()\n",
+    "df = ss.sql(f\"SELECT * FROM  pileup('{table_name}', '{sample_id}', '{ref_path}', true) WHERE contig='1'\").toPandas()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": [
     "df['baf'] = df.countNonRef/df.coverage\n",
-    "df2 = df[(((df.baf > 0.35) & (df.baf< 0.65)) | (df.baf > 0.9)) & (df.coverage>20) & (df.coverage<300)]\n",
-    "plt.rcParams['figure.figsize'] = [15,3]; plt.scatter(df2.pos_start, df2.baf, s=30,alpha=0.1);\n",
-    "plt.xlabel('Chromsome 1');plt.ylabel('B-allele frequency'); plt.ylim(0,1); plt.show()\n"
+    "df2 = df[(((df.baf > 0.35) & (df.baf< 0.65)) | (df.baf > 0.9)) & (df.coverage>20) & (df.coverage<300)]"
    ]
   },
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {},
+   "metadata": {
+    "collapsed": true
+   },
    "outputs": [],
-   "source": []
+   "source": [
+    "import matplotlib.pyplot as plt\n",
+    "plt.rcParams['figure.figsize'] = [15,3]\n",
+    "plt.scatter(df2.pos_start, df2.baf, s=30,alpha=0.1)\n",
+    "plt.xlabel('Chromosome 1')\n",
+    "plt.ylabel('B-allele frequency')\n",
+    "plt.ylim(0,1)\n",
+    "plt.show()"
+   ]
   }
  ],
  "metadata": {
   "kernelspec": {
-   "display_name": "pysequila",
+   "display_name": "Python 3",
    "language": "python",
-   "name": "pysequila"
+   "name": "python3"
   },
   "language_info": {
    "codemirror_mode": {
@@ -762,7 +192,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.7.8"
+   "version": "3.6.7"
   }
  },
  "nbformat": 4,