021.enhancer.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#!/usr/bin/env python\n",
    "# -*- encoding: utf-8 -*-\n",
    "'''\n",
    "@文件        :021.enhancer.ipynb\n",
    "@说明        :ABC enhancer for human\n",
    "@时间        :2023/08/21 17:20:05\n",
    "@作者        :Wu Baosheng\n",
    "@版本        :1.0\n",
    "'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os,re,sys\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "os.chdir('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer')\n",
    "df = pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer/AllPredictions.AvgHiC.ABC0.015.minus150.ForABCPaperV3.txt',sep='\\t',header=0,index_col=None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['fibroblast_of_lung-Roadmap']"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "[x for x in df['CellType'].unique() if 'lung' in x.lower()]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['A549_treated_with_ethanol_0.02_percent_for_1_hour-Roadmap']"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "[x for x in df['CellType'].unique() if 'A549' in x]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver = df[df['CellType']=='liver-ENCODE']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714276</td>\n",
       "      <td>promoter|chr1:713818-714426</td>\n",
       "      <td>promoter</td>\n",
       "      <td>3.897799</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.573692</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.485868</td>\n",
       "      <td>0.265022</td>\n",
       "      <td>0.481129</td>\n",
       "      <td>0.190119</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>369</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762795</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762645-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.909910</td>\n",
       "      <td>LINC00115</td>\n",
       "      <td>762902</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.506701</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.362726</td>\n",
       "      <td>0.126307</td>\n",
       "      <td>0.359188</td>\n",
       "      <td>0.132739</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>370</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762795</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762645-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.909910</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.506701</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.362726</td>\n",
       "      <td>0.126306</td>\n",
       "      <td>0.359188</td>\n",
       "      <td>0.132729</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>878</th>\n",
       "      <td>chr1</td>\n",
       "      <td>859329</td>\n",
       "      <td>859529</td>\n",
       "      <td>intergenic|chr1:859179-859679</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.211608</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.032714</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.033930</td>\n",
       "      <td>0.075040</td>\n",
       "      <td>0.024202</td>\n",
       "      <td>0.272992</td>\n",
       "      <td>0.055114</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>884</th>\n",
       "      <td>chr1</td>\n",
       "      <td>860018</td>\n",
       "      <td>860218</td>\n",
       "      <td>intergenic|chr1:859868-860368</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>3.423063</td>\n",
       "      <td>NOC2L</td>\n",
       "      <td>894679</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.533078</td>\n",
       "      <td>...</td>\n",
       "      <td>0.023099</td>\n",
       "      <td>0.022964</td>\n",
       "      <td>0.012852</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.014068</td>\n",
       "      <td>0.048156</td>\n",
       "      <td>0.015482</td>\n",
       "      <td>0.078606</td>\n",
       "      <td>0.014794</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716611</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>CMC4</td>\n",
       "      <td>154299547</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.653108</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232279</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.390153</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716612</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>MTCP1</td>\n",
       "      <td>154299547</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.653108</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232279</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.390153</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716613</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>BRCC3</td>\n",
       "      <td>154299694</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.646122</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232282</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.391263</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716826</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154444600</td>\n",
       "      <td>154444811</td>\n",
       "      <td>promoter|chrX:154444450-154444961</td>\n",
       "      <td>promoter</td>\n",
       "      <td>4.027038</td>\n",
       "      <td>VBP1</td>\n",
       "      <td>154444700</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.658709</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.037211</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.038427</td>\n",
       "      <td>0.154747</td>\n",
       "      <td>0.172342</td>\n",
       "      <td>0.497082</td>\n",
       "      <td>0.370809</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717096</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154842172</td>\n",
       "      <td>154842722</td>\n",
       "      <td>promoter|chrX:154842022-154842872</td>\n",
       "      <td>promoter</td>\n",
       "      <td>5.566828</td>\n",
       "      <td>TMLHE</td>\n",
       "      <td>154842622</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.582003</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.037942</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.039158</td>\n",
       "      <td>0.217984</td>\n",
       "      <td>0.264255</td>\n",
       "      <td>0.687147</td>\n",
       "      <td>0.545381</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>62122 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                               name  \\\n",
       "100      chr1     713968     714276        promoter|chr1:713818-714426   \n",
       "369      chr1     762795     763070        promoter|chr1:762645-763220   \n",
       "370      chr1     762795     763070        promoter|chr1:762645-763220   \n",
       "878      chr1     859329     859529      intergenic|chr1:859179-859679   \n",
       "884      chr1     860018     860218      intergenic|chr1:859868-860368   \n",
       "...       ...        ...        ...                                ...   \n",
       "7716611  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716612  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716613  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716826  chrX  154444600  154444811  promoter|chrX:154444450-154444961   \n",
       "7717096  chrX  154842172  154842722  promoter|chrX:154842022-154842872   \n",
       "\n",
       "              class  activity_base    TargetGene  TargetGeneTSS  \\\n",
       "100        promoter       3.897799  LOC100288069         714068   \n",
       "369        promoter       2.909910     LINC00115         762902   \n",
       "370        promoter       2.909910     LINC01128         762970   \n",
       "878      intergenic       2.211608        SAMD11         861120   \n",
       "884      intergenic       3.423063         NOC2L         894679   \n",
       "...             ...            ...           ...            ...   \n",
       "7716611    promoter       6.988388          CMC4      154299547   \n",
       "7716612    promoter       6.988388         MTCP1      154299547   \n",
       "7716613    promoter       6.988388         BRCC3      154299694   \n",
       "7716826    promoter       4.027038          VBP1      154444700   \n",
       "7717096    promoter       5.566828         TMLHE      154842622   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "100                       NaN                            0.573692  ...   \n",
       "369                       NaN                            0.506701  ...   \n",
       "370                       NaN                            0.506701  ...   \n",
       "878                       NaN                            0.425330  ...   \n",
       "884                       NaN                            0.533078  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7716611                   NaN                            0.653108  ...   \n",
       "7716612                   NaN                            0.653108  ...   \n",
       "7716613                   NaN                            0.646122  ...   \n",
       "7716826                   NaN                            0.658709  ...   \n",
       "7717096                   NaN                            0.582003  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "100              0.125372                    0.123436               0.123436   \n",
       "369              0.125372                    0.123436               0.123436   \n",
       "370              0.125372                    0.123436               0.123436   \n",
       "878              0.125372                    0.123436               0.032714   \n",
       "884              0.023099                    0.022964               0.012852   \n",
       "...                   ...                         ...                    ...   \n",
       "7716611          0.125372                    0.123436               0.041811   \n",
       "7716612          0.125372                    0.123436               0.041811   \n",
       "7716613          0.125372                    0.123436               0.041811   \n",
       "7716826          0.125372                    0.123436               0.037211   \n",
       "7717096          0.125372                    0.123436               0.037942   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "100             0.001216                   0.124652             0.485868   \n",
       "369             0.001216                   0.124652             0.362726   \n",
       "370             0.001216                   0.124652             0.362726   \n",
       "878             0.001216                   0.033930             0.075040   \n",
       "884             0.001216                   0.014068             0.048156   \n",
       "...                  ...                        ...                  ...   \n",
       "7716611         0.001216                   0.043026             0.300685   \n",
       "7716612         0.001216                   0.043026             0.300685   \n",
       "7716613         0.001216                   0.043026             0.300685   \n",
       "7716826         0.001216                   0.038427             0.154747   \n",
       "7717096         0.001216                   0.039158             0.217984   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "100       0.265022                  0.481129        0.190119  liver-ENCODE  \n",
       "369       0.126307                  0.359188        0.132739  liver-ENCODE  \n",
       "370       0.126306                  0.359188        0.132729  liver-ENCODE  \n",
       "878       0.024202                  0.272992        0.055114  liver-ENCODE  \n",
       "884       0.015482                  0.078606        0.014794  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7716611   0.232279                  0.862619        0.390153  liver-ENCODE  \n",
       "7716612   0.232279                  0.862619        0.390153  liver-ENCODE  \n",
       "7716613   0.232282                  0.862619        0.391263  liver-ENCODE  \n",
       "7716826   0.172342                  0.497082        0.370809  liver-ENCODE  \n",
       "7717096   0.264255                  0.687147        0.545381  liver-ENCODE  \n",
       "\n",
       "[62122 rows x 24 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714385</td>\n",
       "      <td>promoter|chr1:713818-714535</td>\n",
       "      <td>promoter</td>\n",
       "      <td>10.238452</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.811879</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>1.276242</td>\n",
       "      <td>0.659343</td>\n",
       "      <td>1.263794</td>\n",
       "      <td>0.521677</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>134</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714389</td>\n",
       "      <td>promoter|chr1:713818-714539</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.510178</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.686533</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.811506</td>\n",
       "      <td>0.586716</td>\n",
       "      <td>0.803591</td>\n",
       "      <td>0.431711</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>137</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066330</td>\n",
       "      <td>0.065544</td>\n",
       "      <td>0.067555</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068771</td>\n",
       "      <td>0.564927</td>\n",
       "      <td>0.251174</td>\n",
       "      <td>0.538417</td>\n",
       "      <td>0.242575</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>138</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC00115</td>\n",
       "      <td>762902</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066714</td>\n",
       "      <td>0.065921</td>\n",
       "      <td>0.067553</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068769</td>\n",
       "      <td>0.564908</td>\n",
       "      <td>0.251169</td>\n",
       "      <td>0.541513</td>\n",
       "      <td>0.243627</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>319</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762723</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762573-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.669664</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.528495</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.332779</td>\n",
       "      <td>0.197654</td>\n",
       "      <td>0.329533</td>\n",
       "      <td>0.197966</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717169</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154975511</td>\n",
       "      <td>154975711</td>\n",
       "      <td>intergenic|chrX:154975361-154975861</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>1.458903</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034592</td>\n",
       "      <td>0.034309</td>\n",
       "      <td>0.009497</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.010713</td>\n",
       "      <td>0.015630</td>\n",
       "      <td>0.018990</td>\n",
       "      <td>0.050053</td>\n",
       "      <td>0.044118</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717202</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154996931</td>\n",
       "      <td>154997496</td>\n",
       "      <td>promoter|chrX:154996781-154997646</td>\n",
       "      <td>promoter</td>\n",
       "      <td>3.122868</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.145535</td>\n",
       "      <td>0.176828</td>\n",
       "      <td>0.385474</td>\n",
       "      <td>0.339767</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717234</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154998175</td>\n",
       "      <td>154998375</td>\n",
       "      <td>genic|chrX:154998025-154998525</td>\n",
       "      <td>genic</td>\n",
       "      <td>0.309739</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.014435</td>\n",
       "      <td>0.017539</td>\n",
       "      <td>0.038233</td>\n",
       "      <td>0.033699</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717283</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155110783</td>\n",
       "      <td>155111567</td>\n",
       "      <td>promoter|chrX:155110633-155111717</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.397005</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.051684</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052900</td>\n",
       "      <td>0.338402</td>\n",
       "      <td>0.323493</td>\n",
       "      <td>0.789621</td>\n",
       "      <td>0.519858</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717325</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155117947</td>\n",
       "      <td>155118147</td>\n",
       "      <td>genic|chrX:155117797-155118297</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.121687</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.092194</td>\n",
       "      <td>0.090930</td>\n",
       "      <td>0.051775</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052991</td>\n",
       "      <td>0.059439</td>\n",
       "      <td>0.056820</td>\n",
       "      <td>0.101995</td>\n",
       "      <td>0.067150</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>113986 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "123      chr1     713968     714385          promoter|chr1:713818-714535   \n",
       "134      chr1     713968     714389          promoter|chr1:713818-714539   \n",
       "137      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "138      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "319      chr1     762723     763070          promoter|chr1:762573-763220   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7717169  chrX  154975511  154975711  intergenic|chrX:154975361-154975861   \n",
       "7717202  chrX  154996931  154997496    promoter|chrX:154996781-154997646   \n",
       "7717234  chrX  154998175  154998375       genic|chrX:154998025-154998525   \n",
       "7717283  chrX  155110783  155111567    promoter|chrX:155110633-155111717   \n",
       "7717325  chrX  155117947  155118147       genic|chrX:155117797-155118297   \n",
       "\n",
       "              class  activity_base    TargetGene  TargetGeneTSS  \\\n",
       "123        promoter      10.238452  LOC100288069         714068   \n",
       "134        promoter       6.510178  LOC100288069         714068   \n",
       "137        promoter       8.214603     LINC01128         762970   \n",
       "138        promoter       8.214603     LINC00115         762902   \n",
       "319        promoter       2.669664     LINC01128         762970   \n",
       "...             ...            ...           ...            ...   \n",
       "7717169  intergenic       1.458903         SPRY3      154997396   \n",
       "7717202    promoter       3.122868         SPRY3      154997396   \n",
       "7717234       genic       0.309739         SPRY3      154997396   \n",
       "7717283    promoter       6.397005         VAMP7      155110942   \n",
       "7717325       genic       1.121687         VAMP7      155110942   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "123                       NaN                            0.811879  ...   \n",
       "134                       NaN                            0.686533  ...   \n",
       "137                       NaN                            0.604347  ...   \n",
       "138                       NaN                            0.604347  ...   \n",
       "319                       NaN                            0.528495  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7717169                   NaN                            0.554913  ...   \n",
       "7717202                   NaN                            0.554913  ...   \n",
       "7717234                   NaN                            0.554913  ...   \n",
       "7717283                   NaN                            0.641967  ...   \n",
       "7717325                   NaN                            0.641967  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "123              0.125372                    0.123436               0.123436   \n",
       "134              0.125372                    0.123436               0.123436   \n",
       "137              0.066330                    0.065544               0.067555   \n",
       "138              0.066714                    0.065921               0.067553   \n",
       "319              0.125372                    0.123436               0.123436   \n",
       "...                   ...                         ...                    ...   \n",
       "7717169          0.034592                    0.034309               0.009497   \n",
       "7717202          0.125372                    0.123436               0.045387   \n",
       "7717234          0.125372                    0.123436               0.045387   \n",
       "7717283          0.125372                    0.123436               0.051684   \n",
       "7717325          0.092194                    0.090930               0.051775   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "123             0.001216                   0.124652             1.276242   \n",
       "134             0.001216                   0.124652             0.811506   \n",
       "137             0.001216                   0.068771             0.564927   \n",
       "138             0.001216                   0.068769             0.564908   \n",
       "319             0.001216                   0.124652             0.332779   \n",
       "...                  ...                        ...                  ...   \n",
       "7717169         0.001216                   0.010713             0.015630   \n",
       "7717202         0.001216                   0.046603             0.145535   \n",
       "7717234         0.001216                   0.046603             0.014435   \n",
       "7717283         0.001216                   0.052900             0.338402   \n",
       "7717325         0.001216                   0.052991             0.059439   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score  \\\n",
       "123       0.659343                  1.263794        0.521677   \n",
       "134       0.586716                  0.803591        0.431711   \n",
       "137       0.251174                  0.538417        0.242575   \n",
       "138       0.251169                  0.541513        0.243627   \n",
       "319       0.197654                  0.329533        0.197966   \n",
       "...            ...                       ...             ...   \n",
       "7717169   0.018990                  0.050053        0.044118   \n",
       "7717202   0.176828                  0.385474        0.339767   \n",
       "7717234   0.017539                  0.038233        0.033699   \n",
       "7717283   0.323493                  0.789621        0.519858   \n",
       "7717325   0.056820                  0.101995        0.067150   \n",
       "\n",
       "                                                  CellType  \n",
       "123                             fibroblast_of_lung-Roadmap  \n",
       "134      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "137                             fibroblast_of_lung-Roadmap  \n",
       "138                             fibroblast_of_lung-Roadmap  \n",
       "319      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "...                                                    ...  \n",
       "7717169  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717202  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717234  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717283  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717325  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "\n",
       "[113986 rows x 24 columns]"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lung"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver.to_csv('enhancer_liver.csv', index=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from Bio import SeqIO\n",
    "from Bio.Seq import Seq\n",
    "human_dir={}\n",
    "records = SeqIO.parse(\"./hg19.fa\", \"fasta\")\n",
    "for record in records:\n",
    "    id = record.id\n",
    "    seq = record.seq\n",
    "    human_dir[id] = seq"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "file = open('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer/hg19.fa.fai','r')\n",
    "chr_info = {}\n",
    "for line in file:\n",
    "    line = line.strip().split(',')[0].split(' ')\n",
    "    nc = line[0][1: ]\n",
    "    chr =line[-1]\n",
    "    if chr != 'mitochondrion':\n",
    "        id = 'chr' + chr\n",
    "        chr_info[id] = nc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'chr1': 'NC_000001.10', 'chrcontig': 'NT_167251.1', 'chr2': 'NC_000002.11', 'chr3': 'NC_000003.11', 'chr4': 'NC_000004.11', 'chr5': 'NC_000005.9', 'chr6': 'NC_000006.11', 'chr7': 'NC_000007.13', 'chr8': 'NC_000008.10', 'chr9': 'NC_000009.11', 'chr10': 'NC_000010.10', 'chr11': 'NC_000011.9', 'chr12': 'NC_000012.11', 'chr13': 'NC_000013.10', 'chr14': 'NC_000014.8', 'chr15': 'NC_000015.9', 'chr16': 'NC_000016.9', 'chr17': 'NC_000017.10', 'chr18': 'NC_000018.9', 'chr19': 'NC_000019.9', 'chr20': 'NC_000020.10', 'chr21': 'NC_000021.8', 'chr22': 'NC_000022.10', 'chrX': 'NC_000023.10', 'chrY': 'NC_000024.9', 'chrFIX': 'NW_003871103.3', 'chrNOVEL': 'NW_004504305.1'}\n"
     ]
    }
   ],
   "source": [
    "print(chr_info)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [],
   "source": [
    "tmp =[]\n",
    "table = open('./liver_enhancer.fa','w')\n",
    "for index,row in df_liver.iterrows():\n",
    "    name = row['name']\n",
    "    start = row['start']\n",
    "    end = row['end']\n",
    "    chr = row['chr']\n",
    "    id = name + '|' + 'chr'+ \"_\" + str(start) + '_' + str(end)\n",
    "    seq = human_dir[chr_info[chr]][start-1:end]\n",
    "    if id not in tmp:\n",
    "        table.write('>'+id+'\\n')\n",
    "        table.write(str(seq).upper()+'\\n')\n",
    "        if len(str(seq)) == 0:\n",
    "            print(name)\n",
    "        tmp.append(id)\n",
    "    else:\n",
    "        pass\n",
    "table.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "24662"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(tmp)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_lung = df[(df['CellType']=='fibroblast_of_lung-Roadmap') | (df['CellType']=='A549_treated_with_ethanol_0.02_percent_for_1_hour-Roadmap')]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>7302378</th>\n",
       "      <td>chr22</td>\n",
       "      <td>24256156</td>\n",
       "      <td>24256356</td>\n",
       "      <td>intergenic|chr22:24256006-24256506</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>5.227891</td>\n",
       "      <td>MIF</td>\n",
       "      <td>24236564</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.810616</td>\n",
       "      <td>...</td>\n",
       "      <td>0.037788</td>\n",
       "      <td>0.03746</td>\n",
       "      <td>0.014706</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.015921</td>\n",
       "      <td>0.083236</td>\n",
       "      <td>0.036618</td>\n",
       "      <td>0.195838</td>\n",
       "      <td>0.048013</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           chr     start       end                                name  \\\n",
       "7302378  chr22  24256156  24256356  intergenic|chr22:24256006-24256506   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "7302378  intergenic       5.227891        MIF       24236564   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "7302378                   NaN                            0.810616  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "7302378          0.037788                     0.03746               0.014706   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "7302378         0.001216                   0.015921             0.083236   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score  \\\n",
       "7302378   0.036618                  0.195838        0.048013   \n",
       "\n",
       "                           CellType  \n",
       "7302378  fibroblast_of_lung-Roadmap  \n",
       "\n",
       "[1 rows x 24 columns]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lung[df_lung['name']=='intergenic|chr22:24256006-24256506']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714385</td>\n",
       "      <td>promoter|chr1:713818-714535</td>\n",
       "      <td>promoter</td>\n",
       "      <td>10.238452</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.811879</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>1.276242</td>\n",
       "      <td>0.659343</td>\n",
       "      <td>1.263794</td>\n",
       "      <td>0.521677</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>134</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714389</td>\n",
       "      <td>promoter|chr1:713818-714539</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.510178</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.686533</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.811506</td>\n",
       "      <td>0.586716</td>\n",
       "      <td>0.803591</td>\n",
       "      <td>0.431711</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>137</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066330</td>\n",
       "      <td>0.065544</td>\n",
       "      <td>0.067555</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068771</td>\n",
       "      <td>0.564927</td>\n",
       "      <td>0.251174</td>\n",
       "      <td>0.538417</td>\n",
       "      <td>0.242575</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>138</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC00115</td>\n",
       "      <td>762902</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066714</td>\n",
       "      <td>0.065921</td>\n",
       "      <td>0.067553</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068769</td>\n",
       "      <td>0.564908</td>\n",
       "      <td>0.251169</td>\n",
       "      <td>0.541513</td>\n",
       "      <td>0.243627</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>319</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762723</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762573-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.669664</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.528495</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.332779</td>\n",
       "      <td>0.197654</td>\n",
       "      <td>0.329533</td>\n",
       "      <td>0.197966</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717169</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154975511</td>\n",
       "      <td>154975711</td>\n",
       "      <td>intergenic|chrX:154975361-154975861</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>1.458903</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034592</td>\n",
       "      <td>0.034309</td>\n",
       "      <td>0.009497</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.010713</td>\n",
       "      <td>0.015630</td>\n",
       "      <td>0.018990</td>\n",
       "      <td>0.050053</td>\n",
       "      <td>0.044118</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717202</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154996931</td>\n",
       "      <td>154997496</td>\n",
       "      <td>promoter|chrX:154996781-154997646</td>\n",
       "      <td>promoter</td>\n",
       "      <td>3.122868</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.145535</td>\n",
       "      <td>0.176828</td>\n",
       "      <td>0.385474</td>\n",
       "      <td>0.339767</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717234</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154998175</td>\n",
       "      <td>154998375</td>\n",
       "      <td>genic|chrX:154998025-154998525</td>\n",
       "      <td>genic</td>\n",
       "      <td>0.309739</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.014435</td>\n",
       "      <td>0.017539</td>\n",
       "      <td>0.038233</td>\n",
       "      <td>0.033699</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717283</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155110783</td>\n",
       "      <td>155111567</td>\n",
       "      <td>promoter|chrX:155110633-155111717</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.397005</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.051684</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052900</td>\n",
       "      <td>0.338402</td>\n",
       "      <td>0.323493</td>\n",
       "      <td>0.789621</td>\n",
       "      <td>0.519858</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717325</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155117947</td>\n",
       "      <td>155118147</td>\n",
       "      <td>genic|chrX:155117797-155118297</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.121687</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.092194</td>\n",
       "      <td>0.090930</td>\n",
       "      <td>0.051775</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052991</td>\n",
       "      <td>0.059439</td>\n",
       "      <td>0.056820</td>\n",
       "      <td>0.101995</td>\n",
       "      <td>0.067150</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>113986 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "123      chr1     713968     714385          promoter|chr1:713818-714535   \n",
       "134      chr1     713968     714389          promoter|chr1:713818-714539   \n",
       "137      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "138      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "319      chr1     762723     763070          promoter|chr1:762573-763220   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7717169  chrX  154975511  154975711  intergenic|chrX:154975361-154975861   \n",
       "7717202  chrX  154996931  154997496    promoter|chrX:154996781-154997646   \n",
       "7717234  chrX  154998175  154998375       genic|chrX:154998025-154998525   \n",
       "7717283  chrX  155110783  155111567    promoter|chrX:155110633-155111717   \n",
       "7717325  chrX  155117947  155118147       genic|chrX:155117797-155118297   \n",
       "\n",
       "              class  activity_base    TargetGene  TargetGeneTSS  \\\n",
       "123        promoter      10.238452  LOC100288069         714068   \n",
       "134        promoter       6.510178  LOC100288069         714068   \n",
       "137        promoter       8.214603     LINC01128         762970   \n",
       "138        promoter       8.214603     LINC00115         762902   \n",
       "319        promoter       2.669664     LINC01128         762970   \n",
       "...             ...            ...           ...            ...   \n",
       "7717169  intergenic       1.458903         SPRY3      154997396   \n",
       "7717202    promoter       3.122868         SPRY3      154997396   \n",
       "7717234       genic       0.309739         SPRY3      154997396   \n",
       "7717283    promoter       6.397005         VAMP7      155110942   \n",
       "7717325       genic       1.121687         VAMP7      155110942   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "123                       NaN                            0.811879  ...   \n",
       "134                       NaN                            0.686533  ...   \n",
       "137                       NaN                            0.604347  ...   \n",
       "138                       NaN                            0.604347  ...   \n",
       "319                       NaN                            0.528495  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7717169                   NaN                            0.554913  ...   \n",
       "7717202                   NaN                            0.554913  ...   \n",
       "7717234                   NaN                            0.554913  ...   \n",
       "7717283                   NaN                            0.641967  ...   \n",
       "7717325                   NaN                            0.641967  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "123              0.125372                    0.123436               0.123436   \n",
       "134              0.125372                    0.123436               0.123436   \n",
       "137              0.066330                    0.065544               0.067555   \n",
       "138              0.066714                    0.065921               0.067553   \n",
       "319              0.125372                    0.123436               0.123436   \n",
       "...                   ...                         ...                    ...   \n",
       "7717169          0.034592                    0.034309               0.009497   \n",
       "7717202          0.125372                    0.123436               0.045387   \n",
       "7717234          0.125372                    0.123436               0.045387   \n",
       "7717283          0.125372                    0.123436               0.051684   \n",
       "7717325          0.092194                    0.090930               0.051775   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "123             0.001216                   0.124652             1.276242   \n",
       "134             0.001216                   0.124652             0.811506   \n",
       "137             0.001216                   0.068771             0.564927   \n",
       "138             0.001216                   0.068769             0.564908   \n",
       "319             0.001216                   0.124652             0.332779   \n",
       "...                  ...                        ...                  ...   \n",
       "7717169         0.001216                   0.010713             0.015630   \n",
       "7717202         0.001216                   0.046603             0.145535   \n",
       "7717234         0.001216                   0.046603             0.014435   \n",
       "7717283         0.001216                   0.052900             0.338402   \n",
       "7717325         0.001216                   0.052991             0.059439   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score  \\\n",
       "123       0.659343                  1.263794        0.521677   \n",
       "134       0.586716                  0.803591        0.431711   \n",
       "137       0.251174                  0.538417        0.242575   \n",
       "138       0.251169                  0.541513        0.243627   \n",
       "319       0.197654                  0.329533        0.197966   \n",
       "...            ...                       ...             ...   \n",
       "7717169   0.018990                  0.050053        0.044118   \n",
       "7717202   0.176828                  0.385474        0.339767   \n",
       "7717234   0.017539                  0.038233        0.033699   \n",
       "7717283   0.323493                  0.789621        0.519858   \n",
       "7717325   0.056820                  0.101995        0.067150   \n",
       "\n",
       "                                                  CellType  \n",
       "123                             fibroblast_of_lung-Roadmap  \n",
       "134      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "137                             fibroblast_of_lung-Roadmap  \n",
       "138                             fibroblast_of_lung-Roadmap  \n",
       "319      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "...                                                    ...  \n",
       "7717169  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717202  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717234  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717283  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717325  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "\n",
       "[113986 rows x 24 columns]"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lung['type']=df['chr'].str+':'+df['start'].astype(str)+'-'+df['end'].astype(str)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "table = open('./lung_enhancer.fa','w')\n",
    "for index,row in df_lung.iterrows():\n",
    "    name = row['name']\n",
    "    start = row['start']\n",
    "    end = row['end']\n",
    "    chr = row['chr']\n",
    "    seq = human_dir[chr_info[chr]][start-1:end]\n",
    "    id = name + '|' + 'chr'+ \"_\" + str(start) + '_' + str(end)\n",
    "    table.write('>'+id+'\\n')\n",
    "    table.write(str(seq).upper()+'\\n')\n",
    "    if len(str(seq)) == 0:\n",
    "        print(name)\n",
    "table.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>123</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714385</td>\n",
       "      <td>promoter|chr1:713818-714535</td>\n",
       "      <td>promoter</td>\n",
       "      <td>10.238452</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.811879</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>1.276242</td>\n",
       "      <td>0.659343</td>\n",
       "      <td>1.263794</td>\n",
       "      <td>0.521677</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>134</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714389</td>\n",
       "      <td>promoter|chr1:713818-714539</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.510178</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.686533</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.811506</td>\n",
       "      <td>0.586716</td>\n",
       "      <td>0.803591</td>\n",
       "      <td>0.431711</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>137</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066330</td>\n",
       "      <td>0.065544</td>\n",
       "      <td>0.067555</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068771</td>\n",
       "      <td>0.564927</td>\n",
       "      <td>0.251174</td>\n",
       "      <td>0.538417</td>\n",
       "      <td>0.242575</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>138</th>\n",
       "      <td>chr1</td>\n",
       "      <td>752387</td>\n",
       "      <td>752850</td>\n",
       "      <td>promoter|chr1:752237-753000</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.214603</td>\n",
       "      <td>LINC00115</td>\n",
       "      <td>762902</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.604347</td>\n",
       "      <td>...</td>\n",
       "      <td>0.066714</td>\n",
       "      <td>0.065921</td>\n",
       "      <td>0.067553</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068769</td>\n",
       "      <td>0.564908</td>\n",
       "      <td>0.251169</td>\n",
       "      <td>0.541513</td>\n",
       "      <td>0.243627</td>\n",
       "      <td>fibroblast_of_lung-Roadmap</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>319</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762723</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762573-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.669664</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.528495</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.332779</td>\n",
       "      <td>0.197654</td>\n",
       "      <td>0.329533</td>\n",
       "      <td>0.197966</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717169</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154975511</td>\n",
       "      <td>154975711</td>\n",
       "      <td>intergenic|chrX:154975361-154975861</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>1.458903</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034592</td>\n",
       "      <td>0.034309</td>\n",
       "      <td>0.009497</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.010713</td>\n",
       "      <td>0.015630</td>\n",
       "      <td>0.018990</td>\n",
       "      <td>0.050053</td>\n",
       "      <td>0.044118</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717202</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154996931</td>\n",
       "      <td>154997496</td>\n",
       "      <td>promoter|chrX:154996781-154997646</td>\n",
       "      <td>promoter</td>\n",
       "      <td>3.122868</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.145535</td>\n",
       "      <td>0.176828</td>\n",
       "      <td>0.385474</td>\n",
       "      <td>0.339767</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717234</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154998175</td>\n",
       "      <td>154998375</td>\n",
       "      <td>genic|chrX:154998025-154998525</td>\n",
       "      <td>genic</td>\n",
       "      <td>0.309739</td>\n",
       "      <td>SPRY3</td>\n",
       "      <td>154997396</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.554913</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045387</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046603</td>\n",
       "      <td>0.014435</td>\n",
       "      <td>0.017539</td>\n",
       "      <td>0.038233</td>\n",
       "      <td>0.033699</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717283</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155110783</td>\n",
       "      <td>155111567</td>\n",
       "      <td>promoter|chrX:155110633-155111717</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.397005</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.051684</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052900</td>\n",
       "      <td>0.338402</td>\n",
       "      <td>0.323493</td>\n",
       "      <td>0.789621</td>\n",
       "      <td>0.519858</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717325</th>\n",
       "      <td>chrX</td>\n",
       "      <td>155117947</td>\n",
       "      <td>155118147</td>\n",
       "      <td>genic|chrX:155117797-155118297</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.121687</td>\n",
       "      <td>VAMP7</td>\n",
       "      <td>155110942</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.641967</td>\n",
       "      <td>...</td>\n",
       "      <td>0.092194</td>\n",
       "      <td>0.090930</td>\n",
       "      <td>0.051775</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.052991</td>\n",
       "      <td>0.059439</td>\n",
       "      <td>0.056820</td>\n",
       "      <td>0.101995</td>\n",
       "      <td>0.067150</td>\n",
       "      <td>A549_treated_with_ethanol_0.02_percent_for_1_h...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>113986 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "123      chr1     713968     714385          promoter|chr1:713818-714535   \n",
       "134      chr1     713968     714389          promoter|chr1:713818-714539   \n",
       "137      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "138      chr1     752387     752850          promoter|chr1:752237-753000   \n",
       "319      chr1     762723     763070          promoter|chr1:762573-763220   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7717169  chrX  154975511  154975711  intergenic|chrX:154975361-154975861   \n",
       "7717202  chrX  154996931  154997496    promoter|chrX:154996781-154997646   \n",
       "7717234  chrX  154998175  154998375       genic|chrX:154998025-154998525   \n",
       "7717283  chrX  155110783  155111567    promoter|chrX:155110633-155111717   \n",
       "7717325  chrX  155117947  155118147       genic|chrX:155117797-155118297   \n",
       "\n",
       "              class  activity_base    TargetGene  TargetGeneTSS  \\\n",
       "123        promoter      10.238452  LOC100288069         714068   \n",
       "134        promoter       6.510178  LOC100288069         714068   \n",
       "137        promoter       8.214603     LINC01128         762970   \n",
       "138        promoter       8.214603     LINC00115         762902   \n",
       "319        promoter       2.669664     LINC01128         762970   \n",
       "...             ...            ...           ...            ...   \n",
       "7717169  intergenic       1.458903         SPRY3      154997396   \n",
       "7717202    promoter       3.122868         SPRY3      154997396   \n",
       "7717234       genic       0.309739         SPRY3      154997396   \n",
       "7717283    promoter       6.397005         VAMP7      155110942   \n",
       "7717325       genic       1.121687         VAMP7      155110942   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "123                       NaN                            0.811879  ...   \n",
       "134                       NaN                            0.686533  ...   \n",
       "137                       NaN                            0.604347  ...   \n",
       "138                       NaN                            0.604347  ...   \n",
       "319                       NaN                            0.528495  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7717169                   NaN                            0.554913  ...   \n",
       "7717202                   NaN                            0.554913  ...   \n",
       "7717234                   NaN                            0.554913  ...   \n",
       "7717283                   NaN                            0.641967  ...   \n",
       "7717325                   NaN                            0.641967  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "123              0.125372                    0.123436               0.123436   \n",
       "134              0.125372                    0.123436               0.123436   \n",
       "137              0.066330                    0.065544               0.067555   \n",
       "138              0.066714                    0.065921               0.067553   \n",
       "319              0.125372                    0.123436               0.123436   \n",
       "...                   ...                         ...                    ...   \n",
       "7717169          0.034592                    0.034309               0.009497   \n",
       "7717202          0.125372                    0.123436               0.045387   \n",
       "7717234          0.125372                    0.123436               0.045387   \n",
       "7717283          0.125372                    0.123436               0.051684   \n",
       "7717325          0.092194                    0.090930               0.051775   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "123             0.001216                   0.124652             1.276242   \n",
       "134             0.001216                   0.124652             0.811506   \n",
       "137             0.001216                   0.068771             0.564927   \n",
       "138             0.001216                   0.068769             0.564908   \n",
       "319             0.001216                   0.124652             0.332779   \n",
       "...                  ...                        ...                  ...   \n",
       "7717169         0.001216                   0.010713             0.015630   \n",
       "7717202         0.001216                   0.046603             0.145535   \n",
       "7717234         0.001216                   0.046603             0.014435   \n",
       "7717283         0.001216                   0.052900             0.338402   \n",
       "7717325         0.001216                   0.052991             0.059439   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score  \\\n",
       "123       0.659343                  1.263794        0.521677   \n",
       "134       0.586716                  0.803591        0.431711   \n",
       "137       0.251174                  0.538417        0.242575   \n",
       "138       0.251169                  0.541513        0.243627   \n",
       "319       0.197654                  0.329533        0.197966   \n",
       "...            ...                       ...             ...   \n",
       "7717169   0.018990                  0.050053        0.044118   \n",
       "7717202   0.176828                  0.385474        0.339767   \n",
       "7717234   0.017539                  0.038233        0.033699   \n",
       "7717283   0.323493                  0.789621        0.519858   \n",
       "7717325   0.056820                  0.101995        0.067150   \n",
       "\n",
       "                                                  CellType  \n",
       "123                             fibroblast_of_lung-Roadmap  \n",
       "134      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "137                             fibroblast_of_lung-Roadmap  \n",
       "138                             fibroblast_of_lung-Roadmap  \n",
       "319      A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "...                                                    ...  \n",
       "7717169  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717202  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717234  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717283  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "7717325  A549_treated_with_ethanol_0.02_percent_for_1_h...  \n",
       "\n",
       "[113986 rows x 24 columns]"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_lung"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "human_high_varibale = pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/13.homoGenes/human_high_variable_genes.txt', index_col=None, header=None,sep='\\t',names=['gene'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_high_variable=df_liver[df_liver['TargetGene'].isin(human_high_varibale['gene'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>878</th>\n",
       "      <td>chr1</td>\n",
       "      <td>859329</td>\n",
       "      <td>859529</td>\n",
       "      <td>intergenic|chr1:859179-859679</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.211608</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.032714</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.033930</td>\n",
       "      <td>0.075040</td>\n",
       "      <td>0.024202</td>\n",
       "      <td>0.272992</td>\n",
       "      <td>0.055114</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>885</th>\n",
       "      <td>chr1</td>\n",
       "      <td>860018</td>\n",
       "      <td>860218</td>\n",
       "      <td>intergenic|chr1:859868-860368</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>3.423063</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.067077</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068293</td>\n",
       "      <td>0.233771</td>\n",
       "      <td>0.075396</td>\n",
       "      <td>0.422529</td>\n",
       "      <td>0.085304</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>952</th>\n",
       "      <td>chr1</td>\n",
       "      <td>872461</td>\n",
       "      <td>872661</td>\n",
       "      <td>genic|chr1:872311-872811</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.648334</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.060770</td>\n",
       "      <td>0.060079</td>\n",
       "      <td>0.028332</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.029548</td>\n",
       "      <td>0.048705</td>\n",
       "      <td>0.015708</td>\n",
       "      <td>0.099031</td>\n",
       "      <td>0.019993</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>960</th>\n",
       "      <td>chr1</td>\n",
       "      <td>873510</td>\n",
       "      <td>874703</td>\n",
       "      <td>genic|chr1:873360-874853</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.904584</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.054393</td>\n",
       "      <td>0.053809</td>\n",
       "      <td>0.028350</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.029566</td>\n",
       "      <td>0.145009</td>\n",
       "      <td>0.046768</td>\n",
       "      <td>0.263910</td>\n",
       "      <td>0.053281</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1129</th>\n",
       "      <td>chr1</td>\n",
       "      <td>876090</td>\n",
       "      <td>876290</td>\n",
       "      <td>genic|chr1:875940-876440</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.549120</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.047753</td>\n",
       "      <td>0.047276</td>\n",
       "      <td>0.019941</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.021157</td>\n",
       "      <td>0.096244</td>\n",
       "      <td>0.031041</td>\n",
       "      <td>0.215063</td>\n",
       "      <td>0.043419</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7705306</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153601838</td>\n",
       "      <td>153602038</td>\n",
       "      <td>genic|chrX:153601688-153602188</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.298913</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.085079</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.086294</td>\n",
       "      <td>0.198383</td>\n",
       "      <td>0.049837</td>\n",
       "      <td>0.283769</td>\n",
       "      <td>0.045245</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7705432</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153607015</td>\n",
       "      <td>153607696</td>\n",
       "      <td>promoter|chrX:153606865-153607846</td>\n",
       "      <td>promoter</td>\n",
       "      <td>15.685993</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.085079</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.086294</td>\n",
       "      <td>1.353612</td>\n",
       "      <td>0.340049</td>\n",
       "      <td>1.936217</td>\n",
       "      <td>0.308714</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7706120</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153627057</td>\n",
       "      <td>153627257</td>\n",
       "      <td>genic|chrX:153626907-153627407</td>\n",
       "      <td>genic</td>\n",
       "      <td>6.861326</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.031608</td>\n",
       "      <td>0.031366</td>\n",
       "      <td>0.013048</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.014264</td>\n",
       "      <td>0.097868</td>\n",
       "      <td>0.024586</td>\n",
       "      <td>0.215209</td>\n",
       "      <td>0.034313</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7711195</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153770112</td>\n",
       "      <td>153770312</td>\n",
       "      <td>genic|chrX:153769962-153770462</td>\n",
       "      <td>genic</td>\n",
       "      <td>7.777634</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005815</td>\n",
       "      <td>0.005826</td>\n",
       "      <td>0.006523</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.007739</td>\n",
       "      <td>0.060193</td>\n",
       "      <td>0.015121</td>\n",
       "      <td>0.045315</td>\n",
       "      <td>0.007225</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7712395</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153770892</td>\n",
       "      <td>153771374</td>\n",
       "      <td>genic|chrX:153770742-153771524</td>\n",
       "      <td>genic</td>\n",
       "      <td>12.060128</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005787</td>\n",
       "      <td>0.005798</td>\n",
       "      <td>0.006524</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.007739</td>\n",
       "      <td>0.093339</td>\n",
       "      <td>0.023448</td>\n",
       "      <td>0.069930</td>\n",
       "      <td>0.011150</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>4803 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                               name  \\\n",
       "878      chr1     859329     859529      intergenic|chr1:859179-859679   \n",
       "885      chr1     860018     860218      intergenic|chr1:859868-860368   \n",
       "952      chr1     872461     872661           genic|chr1:872311-872811   \n",
       "960      chr1     873510     874703           genic|chr1:873360-874853   \n",
       "1129     chr1     876090     876290           genic|chr1:875940-876440   \n",
       "...       ...        ...        ...                                ...   \n",
       "7705306  chrX  153601838  153602038     genic|chrX:153601688-153602188   \n",
       "7705432  chrX  153607015  153607696  promoter|chrX:153606865-153607846   \n",
       "7706120  chrX  153627057  153627257     genic|chrX:153626907-153627407   \n",
       "7711195  chrX  153770112  153770312     genic|chrX:153769962-153770462   \n",
       "7712395  chrX  153770892  153771374     genic|chrX:153770742-153771524   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "878      intergenic       2.211608     SAMD11         861120   \n",
       "885      intergenic       3.423063     SAMD11         861120   \n",
       "952           genic       1.648334     SAMD11         861120   \n",
       "960           genic       4.904584     SAMD11         861120   \n",
       "1129          genic       4.549120     SAMD11         861120   \n",
       "...             ...            ...        ...            ...   \n",
       "7705306       genic       2.298913       FLNA      153603006   \n",
       "7705432    promoter      15.685993       FLNA      153603006   \n",
       "7706120       genic       6.861326       FLNA      153603006   \n",
       "7711195       genic       7.777634       FLNA      153603006   \n",
       "7712395       genic      12.060128       FLNA      153603006   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "878                       NaN                            0.425330  ...   \n",
       "885                       NaN                            0.425330  ...   \n",
       "952                       NaN                            0.425330  ...   \n",
       "960                       NaN                            0.425330  ...   \n",
       "1129                      NaN                            0.425330  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7705306                   NaN                            0.437307  ...   \n",
       "7705432                   NaN                            0.437307  ...   \n",
       "7706120                   NaN                            0.437307  ...   \n",
       "7711195                   NaN                            0.437307  ...   \n",
       "7712395                   NaN                            0.437307  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "878              0.125372                    0.123436               0.032714   \n",
       "885              0.125372                    0.123436               0.067077   \n",
       "952              0.060770                    0.060079               0.028332   \n",
       "960              0.054393                    0.053809               0.028350   \n",
       "1129             0.047753                    0.047276               0.019941   \n",
       "...                   ...                         ...                    ...   \n",
       "7705306          0.125372                    0.123436               0.085079   \n",
       "7705432          0.125372                    0.123436               0.085079   \n",
       "7706120          0.031608                    0.031366               0.013048   \n",
       "7711195          0.005815                    0.005826               0.006523   \n",
       "7712395          0.005787                    0.005798               0.006524   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "878             0.001216                   0.033930             0.075040   \n",
       "885             0.001216                   0.068293             0.233771   \n",
       "952             0.001216                   0.029548             0.048705   \n",
       "960             0.001216                   0.029566             0.145009   \n",
       "1129            0.001216                   0.021157             0.096244   \n",
       "...                  ...                        ...                  ...   \n",
       "7705306         0.001216                   0.086294             0.198383   \n",
       "7705432         0.001216                   0.086294             1.353612   \n",
       "7706120         0.001216                   0.014264             0.097868   \n",
       "7711195         0.001216                   0.007739             0.060193   \n",
       "7712395         0.001216                   0.007739             0.093339   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "878       0.024202                  0.272992        0.055114  liver-ENCODE  \n",
       "885       0.075396                  0.422529        0.085304  liver-ENCODE  \n",
       "952       0.015708                  0.099031        0.019993  liver-ENCODE  \n",
       "960       0.046768                  0.263910        0.053281  liver-ENCODE  \n",
       "1129      0.031041                  0.215063        0.043419  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7705306   0.049837                  0.283769        0.045245  liver-ENCODE  \n",
       "7705432   0.340049                  1.936217        0.308714  liver-ENCODE  \n",
       "7706120   0.024586                  0.215209        0.034313  liver-ENCODE  \n",
       "7711195   0.015121                  0.045315        0.007225  liver-ENCODE  \n",
       "7712395   0.023448                  0.069930        0.011150  liver-ENCODE  \n",
       "\n",
       "[4803 rows x 24 columns]"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_high_variable.drop_duplicates(subset=['name'],keep='first')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [],
   "source": [
    "table = open('./liver_high_enhancer.fa','w')\n",
    "for index,row in df_liver_high_variable.iterrows():\n",
    "    name = row['name']\n",
    "    start = row['start']\n",
    "    end = row['end']\n",
    "    chr = row['chr']\n",
    "    seq = human_dir[chr_info[chr]][start-1:end]\n",
    "    table.write('>'+name+'\\n')\n",
    "    table.write(str(seq).upper()+'\\n')\n",
    "    if len(str(seq)) == 0:\n",
    "        print(name)\n",
    "table.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os,re,sys\n",
    "from Bio import SeqIO\n",
    "from Bio.Seq import Seq\n",
    "import pandas as pd\n",
    "human_dir={}\n",
    "records = SeqIO.parse(\"./mouse.fa\", \"fasta\")\n",
    "for record in records:\n",
    "    id = record.id\n",
    "    seq = record.seq\n",
    "    human_dir[id] = seq\n",
    "\n",
    "df = pd.read_csv('scACAT.txt',sep='\\t',header=None,names=['id','deepth'],index_col=None)\n",
    "set([x.split(\":\")[0] for x in df['id']])\n",
    "result = open('scACAT.fa','w')\n",
    "for i in df['id']:\n",
    "    if i=='MT':\n",
    "        continue\n",
    "    else:\n",
    "        chr = re.search(r'chr(\\w+)',i.split(':')[0]).group(1)\n",
    "        start = int(i.split(':')[1].split('-')[0])\n",
    "        end = int(i.split(':')[1].split('-')[1])\n",
    "        seq = human_dir[chr][start:end]\n",
    "        result.write('>'+i+'\\n'+str(seq)+'\\n')\n",
    "result.close()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2454812</th>\n",
       "      <td>chr5</td>\n",
       "      <td>180072333</td>\n",
       "      <td>180072533</td>\n",
       "      <td>genic|chr5:180072183-180072683</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.526214</td>\n",
       "      <td>FLT4</td>\n",
       "      <td>180076624</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.466637</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.043969</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.045185</td>\n",
       "      <td>0.068962</td>\n",
       "      <td>0.038165</td>\n",
       "      <td>0.188390</td>\n",
       "      <td>0.065526</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2457749</th>\n",
       "      <td>chr5</td>\n",
       "      <td>180236941</td>\n",
       "      <td>180237783</td>\n",
       "      <td>intergenic|chr5:180236791-180237933</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>50.443049</td>\n",
       "      <td>FLT4</td>\n",
       "      <td>180076624</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.466637</td>\n",
       "      <td>...</td>\n",
       "      <td>0.006019</td>\n",
       "      <td>0.006030</td>\n",
       "      <td>0.000615</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.001831</td>\n",
       "      <td>0.092355</td>\n",
       "      <td>0.051111</td>\n",
       "      <td>0.304157</td>\n",
       "      <td>0.105793</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "2454812  chr5  180072333  180072533       genic|chr5:180072183-180072683   \n",
       "2457749  chr5  180236941  180237783  intergenic|chr5:180236791-180237933   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "2454812       genic       1.526214       FLT4      180076624   \n",
       "2457749  intergenic      50.443049       FLT4      180076624   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "2454812                   NaN                            0.466637  ...   \n",
       "2457749                   NaN                            0.466637  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "2454812          0.125372                    0.123436               0.043969   \n",
       "2457749          0.006019                    0.006030               0.000615   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "2454812         0.001216                   0.045185             0.068962   \n",
       "2457749         0.001216                   0.001831             0.092355   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "2454812   0.038165                  0.188390        0.065526  liver-ENCODE  \n",
       "2457749   0.051111                  0.304157        0.105793  liver-ENCODE  \n",
       "\n",
       "[2 rows x 24 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver[df_liver['TargetGene']=='FLT4']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [],
   "source": [
    "mouse = pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/13.homoGenes/mouse.homoGenes.txt',sep='\\t',header=0,index_col=0)\n",
    "mouse.fillna('no value', inplace=True)\n",
    "mouse_wgd = []\n",
    "for index, row in mouse.iterrows():\n",
    "    if row['lancelet']== 'no value':\n",
    "        line = row['mouse'].split(',')\n",
    "        line2 = [x.strip() for x in line]\n",
    "        mouse_wgd = mouse_wgd + line2\n",
    "    else:\n",
    "        tmp =len(row['mouse'].split(','))-len(row['lancelet'].split(','))\n",
    "        line = row['mouse'].split(',')\n",
    "        line2 = [x.strip() for x in line]\n",
    "        mouse_wgd = mouse_wgd + line2[tmp: ]\n",
    "\n",
    "\n",
    "    # line = i.split(',')\n",
    "    # line2 = [x.strip() for x in line]\n",
    "    # human_wgd = human_wgd + line2[0: ]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_R1=pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/13.homoGenes/mouse_R2.txt',sep='\\t',header=0,index_col=0)\n",
    "df_R1.fillna('no value', inplace=True)\n",
    "R2 = []\n",
    "for index, row in df_R1.iterrows():\n",
    "    if row['lamprey']== 'no value':\n",
    "        line = row['mouse'].split(',')\n",
    "        line2 = [x.strip() for x in line]\n",
    "        R2 = R2 + line2\n",
    "    else:\n",
    "        tmp =len(row['mouse'].split(','))-len(row['lamprey'].split(','))\n",
    "        line = row['mouse'].split(',')\n",
    "        line2 = [x.strip() for x in line]\n",
    "        R2 = R2 + line2[tmp: ]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [],
   "source": [
    "R1=[x for x in mouse_wgd if x not in R2]\n",
    "R1=[x.upper() for x in R1]\n",
    "R2=[x.upper() for x in R2]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3751\n"
     ]
    }
   ],
   "source": [
    "print(len(R2))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_R1= df_liver_high_variable[df_liver_high_variable['TargetGene'].isin(R1)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_R2= df_liver_high_variable[df_liver_high_variable['TargetGene'].isin(R2)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 85,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>74521</th>\n",
       "      <td>chr1</td>\n",
       "      <td>16082992</td>\n",
       "      <td>16083192</td>\n",
       "      <td>intergenic|chr1:16082842-16083342</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.143405</td>\n",
       "      <td>FBLIM1</td>\n",
       "      <td>16085254</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.498900</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.028253</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.029469</td>\n",
       "      <td>0.063164</td>\n",
       "      <td>0.022091</td>\n",
       "      <td>0.264573</td>\n",
       "      <td>0.047719</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>74690</th>\n",
       "      <td>chr1</td>\n",
       "      <td>16107274</td>\n",
       "      <td>16107474</td>\n",
       "      <td>genic|chr1:16107124-16107624</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.253362</td>\n",
       "      <td>FBLIM1</td>\n",
       "      <td>16085254</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.498900</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034133</td>\n",
       "      <td>0.033856</td>\n",
       "      <td>0.010644</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.011860</td>\n",
       "      <td>0.050444</td>\n",
       "      <td>0.017643</td>\n",
       "      <td>0.144004</td>\n",
       "      <td>0.025973</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>104082</th>\n",
       "      <td>chr1</td>\n",
       "      <td>21671397</td>\n",
       "      <td>21672364</td>\n",
       "      <td>intergenic|chr1:21671247-21672514</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>26.304502</td>\n",
       "      <td>C1QA</td>\n",
       "      <td>22962998</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.411398</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000972</td>\n",
       "      <td>0.000984</td>\n",
       "      <td>0.000205</td>\n",
       "      <td>0.000972</td>\n",
       "      <td>0.001177</td>\n",
       "      <td>0.030965</td>\n",
       "      <td>0.015554</td>\n",
       "      <td>0.025889</td>\n",
       "      <td>0.010547</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>107613</th>\n",
       "      <td>chr1</td>\n",
       "      <td>22968540</td>\n",
       "      <td>22968740</td>\n",
       "      <td>intergenic|chr1:22968390-22968890</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>1.424871</td>\n",
       "      <td>C1QA</td>\n",
       "      <td>22962998</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.411398</td>\n",
       "      <td>...</td>\n",
       "      <td>0.112799</td>\n",
       "      <td>0.111124</td>\n",
       "      <td>0.042081</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043297</td>\n",
       "      <td>0.061692</td>\n",
       "      <td>0.030988</td>\n",
       "      <td>0.158337</td>\n",
       "      <td>0.064507</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123607</th>\n",
       "      <td>chr1</td>\n",
       "      <td>26606112</td>\n",
       "      <td>26606636</td>\n",
       "      <td>promoter|chr1:26605962-26606786</td>\n",
       "      <td>promoter</td>\n",
       "      <td>13.140121</td>\n",
       "      <td>SH3BGRL3</td>\n",
       "      <td>26606212</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.865977</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.042142</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043358</td>\n",
       "      <td>0.569724</td>\n",
       "      <td>0.158123</td>\n",
       "      <td>1.621965</td>\n",
       "      <td>0.246040</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7573598</th>\n",
       "      <td>chrX</td>\n",
       "      <td>70842985</td>\n",
       "      <td>70843185</td>\n",
       "      <td>genic|chrX:70842835-70843335</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.077731</td>\n",
       "      <td>CXCR3</td>\n",
       "      <td>70838367</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.423212</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.030543</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.031759</td>\n",
       "      <td>0.065987</td>\n",
       "      <td>0.076294</td>\n",
       "      <td>0.256467</td>\n",
       "      <td>0.136522</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7653580</th>\n",
       "      <td>chrX</td>\n",
       "      <td>128901553</td>\n",
       "      <td>128901753</td>\n",
       "      <td>genic|chrX:128901403-128901903</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.167014</td>\n",
       "      <td>SASH3</td>\n",
       "      <td>128913891</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.472747</td>\n",
       "      <td>...</td>\n",
       "      <td>0.057293</td>\n",
       "      <td>0.056661</td>\n",
       "      <td>0.014470</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.015686</td>\n",
       "      <td>0.018306</td>\n",
       "      <td>0.017490</td>\n",
       "      <td>0.066124</td>\n",
       "      <td>0.023872</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7657073</th>\n",
       "      <td>chrX</td>\n",
       "      <td>129091321</td>\n",
       "      <td>129091735</td>\n",
       "      <td>intergenic|chrX:129091171-129091885</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>6.635338</td>\n",
       "      <td>SASH3</td>\n",
       "      <td>128913891</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.472747</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005515</td>\n",
       "      <td>0.005527</td>\n",
       "      <td>0.001168</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002384</td>\n",
       "      <td>0.015820</td>\n",
       "      <td>0.015115</td>\n",
       "      <td>0.036677</td>\n",
       "      <td>0.013241</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7676688</th>\n",
       "      <td>chrX</td>\n",
       "      <td>135229458</td>\n",
       "      <td>135230013</td>\n",
       "      <td>promoter|chrX:135229308-135230163</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.571128</td>\n",
       "      <td>FHL1</td>\n",
       "      <td>135229558</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.507210</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.034564</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.035780</td>\n",
       "      <td>0.091995</td>\n",
       "      <td>0.130602</td>\n",
       "      <td>0.317370</td>\n",
       "      <td>0.221601</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7676913</th>\n",
       "      <td>chrX</td>\n",
       "      <td>135251521</td>\n",
       "      <td>135251721</td>\n",
       "      <td>genic|chrX:135251371-135251871</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.961939</td>\n",
       "      <td>FHL1</td>\n",
       "      <td>135229558</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.507210</td>\n",
       "      <td>...</td>\n",
       "      <td>0.034210</td>\n",
       "      <td>0.033933</td>\n",
       "      <td>0.005587</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.006802</td>\n",
       "      <td>0.013346</td>\n",
       "      <td>0.018947</td>\n",
       "      <td>0.066574</td>\n",
       "      <td>0.046484</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>761 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "74521    chr1   16082992   16083192    intergenic|chr1:16082842-16083342   \n",
       "74690    chr1   16107274   16107474         genic|chr1:16107124-16107624   \n",
       "104082   chr1   21671397   21672364    intergenic|chr1:21671247-21672514   \n",
       "107613   chr1   22968540   22968740    intergenic|chr1:22968390-22968890   \n",
       "123607   chr1   26606112   26606636      promoter|chr1:26605962-26606786   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7573598  chrX   70842985   70843185         genic|chrX:70842835-70843335   \n",
       "7653580  chrX  128901553  128901753       genic|chrX:128901403-128901903   \n",
       "7657073  chrX  129091321  129091735  intergenic|chrX:129091171-129091885   \n",
       "7676688  chrX  135229458  135230013    promoter|chrX:135229308-135230163   \n",
       "7676913  chrX  135251521  135251721       genic|chrX:135251371-135251871   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "74521    intergenic       2.143405     FBLIM1       16085254   \n",
       "74690         genic       4.253362     FBLIM1       16085254   \n",
       "104082   intergenic      26.304502       C1QA       22962998   \n",
       "107613   intergenic       1.424871       C1QA       22962998   \n",
       "123607     promoter      13.140121   SH3BGRL3       26606212   \n",
       "...             ...            ...        ...            ...   \n",
       "7573598       genic       2.077731      CXCR3       70838367   \n",
       "7653580       genic       1.167014      SASH3      128913891   \n",
       "7657073  intergenic       6.635338      SASH3      128913891   \n",
       "7676688    promoter       2.571128       FHL1      135229558   \n",
       "7676913       genic       1.961939       FHL1      135229558   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "74521                     NaN                            0.498900  ...   \n",
       "74690                     NaN                            0.498900  ...   \n",
       "104082                    NaN                            0.411398  ...   \n",
       "107613                    NaN                            0.411398  ...   \n",
       "123607                    NaN                            0.865977  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7573598                   NaN                            0.423212  ...   \n",
       "7653580                   NaN                            0.472747  ...   \n",
       "7657073                   NaN                            0.472747  ...   \n",
       "7676688                   NaN                            0.507210  ...   \n",
       "7676913                   NaN                            0.507210  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "74521            0.125372                    0.123436               0.028253   \n",
       "74690            0.034133                    0.033856               0.010644   \n",
       "104082           0.000972                    0.000984               0.000205   \n",
       "107613           0.112799                    0.111124               0.042081   \n",
       "123607           0.125372                    0.123436               0.042142   \n",
       "...                   ...                         ...                    ...   \n",
       "7573598          0.125372                    0.123436               0.030543   \n",
       "7653580          0.057293                    0.056661               0.014470   \n",
       "7657073          0.005515                    0.005527               0.001168   \n",
       "7676688          0.125372                    0.123436               0.034564   \n",
       "7676913          0.034210                    0.033933               0.005587   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "74521           0.001216                   0.029469             0.063164   \n",
       "74690           0.001216                   0.011860             0.050444   \n",
       "104082          0.000972                   0.001177             0.030965   \n",
       "107613          0.001216                   0.043297             0.061692   \n",
       "123607          0.001216                   0.043358             0.569724   \n",
       "...                  ...                        ...                  ...   \n",
       "7573598         0.001216                   0.031759             0.065987   \n",
       "7653580         0.001216                   0.015686             0.018306   \n",
       "7657073         0.001216                   0.002384             0.015820   \n",
       "7676688         0.001216                   0.035780             0.091995   \n",
       "7676913         0.001216                   0.006802             0.013346   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "74521     0.022091                  0.264573        0.047719  liver-ENCODE  \n",
       "74690     0.017643                  0.144004        0.025973  liver-ENCODE  \n",
       "104082    0.015554                  0.025889        0.010547  liver-ENCODE  \n",
       "107613    0.030988                  0.158337        0.064507  liver-ENCODE  \n",
       "123607    0.158123                  1.621965        0.246040  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7573598   0.076294                  0.256467        0.136522  liver-ENCODE  \n",
       "7653580   0.017490                  0.066124        0.023872  liver-ENCODE  \n",
       "7657073   0.015115                  0.036677        0.013241  liver-ENCODE  \n",
       "7676688   0.130602                  0.317370        0.221601  liver-ENCODE  \n",
       "7676913   0.018947                  0.066574        0.046484  liver-ENCODE  \n",
       "\n",
       "[761 rows x 24 columns]"
      ]
     },
     "execution_count": 85,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_R2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [],
   "source": [
    "blast=pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/014.mouse_liver_development/02.blastnParse.py.out',index_col=None,sep='\\t',header=None,names=['human','mouse'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 89,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_R1_mouse = df_liver_R1[df_liver_R1['name'].isin(blast['human'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_R2_mouse = df_liver_R2[df_liver_R2['name'].isin(blast['human'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 91,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>104082</th>\n",
       "      <td>chr1</td>\n",
       "      <td>21671397</td>\n",
       "      <td>21672364</td>\n",
       "      <td>intergenic|chr1:21671247-21672514</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>26.304502</td>\n",
       "      <td>C1QA</td>\n",
       "      <td>22962998</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.411398</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000972</td>\n",
       "      <td>0.000984</td>\n",
       "      <td>0.000205</td>\n",
       "      <td>0.000972</td>\n",
       "      <td>0.001177</td>\n",
       "      <td>0.030965</td>\n",
       "      <td>0.015554</td>\n",
       "      <td>0.025889</td>\n",
       "      <td>0.010547</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>123607</th>\n",
       "      <td>chr1</td>\n",
       "      <td>26606112</td>\n",
       "      <td>26606636</td>\n",
       "      <td>promoter|chr1:26605962-26606786</td>\n",
       "      <td>promoter</td>\n",
       "      <td>13.140121</td>\n",
       "      <td>SH3BGRL3</td>\n",
       "      <td>26606212</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.865977</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.042142</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043358</td>\n",
       "      <td>0.569724</td>\n",
       "      <td>0.158123</td>\n",
       "      <td>1.621965</td>\n",
       "      <td>0.246040</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>129947</th>\n",
       "      <td>chr1</td>\n",
       "      <td>26946849</td>\n",
       "      <td>26947665</td>\n",
       "      <td>intergenic|chr1:26946699-26947815</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>40.799530</td>\n",
       "      <td>PTAFR</td>\n",
       "      <td>28503455</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.436247</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000826</td>\n",
       "      <td>0.000837</td>\n",
       "      <td>0.000325</td>\n",
       "      <td>0.000826</td>\n",
       "      <td>0.001150</td>\n",
       "      <td>0.046930</td>\n",
       "      <td>0.019518</td>\n",
       "      <td>0.034133</td>\n",
       "      <td>0.009794</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>143628</th>\n",
       "      <td>chr1</td>\n",
       "      <td>28498699</td>\n",
       "      <td>28498899</td>\n",
       "      <td>genic|chr1:28498549-28499049</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.838767</td>\n",
       "      <td>PTAFR</td>\n",
       "      <td>28503455</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.436247</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.036465</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.037680</td>\n",
       "      <td>0.069286</td>\n",
       "      <td>0.028815</td>\n",
       "      <td>0.226970</td>\n",
       "      <td>0.065126</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>151851</th>\n",
       "      <td>chr1</td>\n",
       "      <td>31296977</td>\n",
       "      <td>31297510</td>\n",
       "      <td>intergenic|chr1:31296827-31297660</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>14.425468</td>\n",
       "      <td>SDC3</td>\n",
       "      <td>31381480</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.509492</td>\n",
       "      <td>...</td>\n",
       "      <td>0.010594</td>\n",
       "      <td>0.010579</td>\n",
       "      <td>0.007571</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.008787</td>\n",
       "      <td>0.126761</td>\n",
       "      <td>0.047650</td>\n",
       "      <td>0.152602</td>\n",
       "      <td>0.039463</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7494962</th>\n",
       "      <td>chrX</td>\n",
       "      <td>39968051</td>\n",
       "      <td>39968533</td>\n",
       "      <td>intergenic|chrX:39967901-39968683</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>12.617401</td>\n",
       "      <td>SRPX</td>\n",
       "      <td>38080177</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.447654</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000697</td>\n",
       "      <td>0.000707</td>\n",
       "      <td>0.000113</td>\n",
       "      <td>0.000697</td>\n",
       "      <td>0.000810</td>\n",
       "      <td>0.010226</td>\n",
       "      <td>0.017051</td>\n",
       "      <td>0.008922</td>\n",
       "      <td>0.009132</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7527818</th>\n",
       "      <td>chrX</td>\n",
       "      <td>49023441</td>\n",
       "      <td>49023641</td>\n",
       "      <td>genic|chrX:49023291-49023791</td>\n",
       "      <td>genic</td>\n",
       "      <td>6.256536</td>\n",
       "      <td>PLP2</td>\n",
       "      <td>49028183</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.512058</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.030186</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.031401</td>\n",
       "      <td>0.196464</td>\n",
       "      <td>0.081982</td>\n",
       "      <td>0.772282</td>\n",
       "      <td>0.138158</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7569958</th>\n",
       "      <td>chrX</td>\n",
       "      <td>70401856</td>\n",
       "      <td>70402056</td>\n",
       "      <td>intergenic|chrX:70401706-70402206</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>6.368728</td>\n",
       "      <td>CXCR3</td>\n",
       "      <td>70838367</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.423212</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002512</td>\n",
       "      <td>0.002529</td>\n",
       "      <td>0.000850</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002066</td>\n",
       "      <td>0.013157</td>\n",
       "      <td>0.015212</td>\n",
       "      <td>0.016105</td>\n",
       "      <td>0.008573</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7657073</th>\n",
       "      <td>chrX</td>\n",
       "      <td>129091321</td>\n",
       "      <td>129091735</td>\n",
       "      <td>intergenic|chrX:129091171-129091885</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>6.635338</td>\n",
       "      <td>SASH3</td>\n",
       "      <td>128913891</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.472747</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005515</td>\n",
       "      <td>0.005527</td>\n",
       "      <td>0.001168</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002384</td>\n",
       "      <td>0.015820</td>\n",
       "      <td>0.015115</td>\n",
       "      <td>0.036677</td>\n",
       "      <td>0.013241</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7676688</th>\n",
       "      <td>chrX</td>\n",
       "      <td>135229458</td>\n",
       "      <td>135230013</td>\n",
       "      <td>promoter|chrX:135229308-135230163</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.571128</td>\n",
       "      <td>FHL1</td>\n",
       "      <td>135229558</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.507210</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.034564</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.035780</td>\n",
       "      <td>0.091995</td>\n",
       "      <td>0.130602</td>\n",
       "      <td>0.317370</td>\n",
       "      <td>0.221601</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>240 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "104082   chr1   21671397   21672364    intergenic|chr1:21671247-21672514   \n",
       "123607   chr1   26606112   26606636      promoter|chr1:26605962-26606786   \n",
       "129947   chr1   26946849   26947665    intergenic|chr1:26946699-26947815   \n",
       "143628   chr1   28498699   28498899         genic|chr1:28498549-28499049   \n",
       "151851   chr1   31296977   31297510    intergenic|chr1:31296827-31297660   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7494962  chrX   39968051   39968533    intergenic|chrX:39967901-39968683   \n",
       "7527818  chrX   49023441   49023641         genic|chrX:49023291-49023791   \n",
       "7569958  chrX   70401856   70402056    intergenic|chrX:70401706-70402206   \n",
       "7657073  chrX  129091321  129091735  intergenic|chrX:129091171-129091885   \n",
       "7676688  chrX  135229458  135230013    promoter|chrX:135229308-135230163   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "104082   intergenic      26.304502       C1QA       22962998   \n",
       "123607     promoter      13.140121   SH3BGRL3       26606212   \n",
       "129947   intergenic      40.799530      PTAFR       28503455   \n",
       "143628        genic       1.838767      PTAFR       28503455   \n",
       "151851   intergenic      14.425468       SDC3       31381480   \n",
       "...             ...            ...        ...            ...   \n",
       "7494962  intergenic      12.617401       SRPX       38080177   \n",
       "7527818       genic       6.256536       PLP2       49028183   \n",
       "7569958  intergenic       6.368728      CXCR3       70838367   \n",
       "7657073  intergenic       6.635338      SASH3      128913891   \n",
       "7676688    promoter       2.571128       FHL1      135229558   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "104082                    NaN                            0.411398  ...   \n",
       "123607                    NaN                            0.865977  ...   \n",
       "129947                    NaN                            0.436247  ...   \n",
       "143628                    NaN                            0.436247  ...   \n",
       "151851                    NaN                            0.509492  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7494962                   NaN                            0.447654  ...   \n",
       "7527818                   NaN                            0.512058  ...   \n",
       "7569958                   NaN                            0.423212  ...   \n",
       "7657073                   NaN                            0.472747  ...   \n",
       "7676688                   NaN                            0.507210  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "104082           0.000972                    0.000984               0.000205   \n",
       "123607           0.125372                    0.123436               0.042142   \n",
       "129947           0.000826                    0.000837               0.000325   \n",
       "143628           0.125372                    0.123436               0.036465   \n",
       "151851           0.010594                    0.010579               0.007571   \n",
       "...                   ...                         ...                    ...   \n",
       "7494962          0.000697                    0.000707               0.000113   \n",
       "7527818          0.125372                    0.123436               0.030186   \n",
       "7569958          0.002512                    0.002529               0.000850   \n",
       "7657073          0.005515                    0.005527               0.001168   \n",
       "7676688          0.125372                    0.123436               0.034564   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "104082          0.000972                   0.001177             0.030965   \n",
       "123607          0.001216                   0.043358             0.569724   \n",
       "129947          0.000826                   0.001150             0.046930   \n",
       "143628          0.001216                   0.037680             0.069286   \n",
       "151851          0.001216                   0.008787             0.126761   \n",
       "...                  ...                        ...                  ...   \n",
       "7494962         0.000697                   0.000810             0.010226   \n",
       "7527818         0.001216                   0.031401             0.196464   \n",
       "7569958         0.001216                   0.002066             0.013157   \n",
       "7657073         0.001216                   0.002384             0.015820   \n",
       "7676688         0.001216                   0.035780             0.091995   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "104082    0.015554                  0.025889        0.010547  liver-ENCODE  \n",
       "123607    0.158123                  1.621965        0.246040  liver-ENCODE  \n",
       "129947    0.019518                  0.034133        0.009794  liver-ENCODE  \n",
       "143628    0.028815                  0.226970        0.065126  liver-ENCODE  \n",
       "151851    0.047650                  0.152602        0.039463  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7494962   0.017051                  0.008922        0.009132  liver-ENCODE  \n",
       "7527818   0.081982                  0.772282        0.138158  liver-ENCODE  \n",
       "7569958   0.015212                  0.016105        0.008573  liver-ENCODE  \n",
       "7657073   0.015115                  0.036677        0.013241  liver-ENCODE  \n",
       "7676688   0.130602                  0.317370        0.221601  liver-ENCODE  \n",
       "\n",
       "[240 rows x 24 columns]"
      ]
     },
     "execution_count": 91,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_R2_mouse"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 92,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(145, 24)"
      ]
     },
     "execution_count": 92,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_R1_mouse.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 93,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(240, 24)"
      ]
     },
     "execution_count": 93,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_R2_mouse.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 95,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>136939</th>\n",
       "      <td>chr1</td>\n",
       "      <td>27669515</td>\n",
       "      <td>27669715</td>\n",
       "      <td>genic|chr1:27669365-27669865</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.908763</td>\n",
       "      <td>SYTL1</td>\n",
       "      <td>27668482</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.498696</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.035278</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.036494</td>\n",
       "      <td>0.106152</td>\n",
       "      <td>0.034515</td>\n",
       "      <td>0.359046</td>\n",
       "      <td>0.075990</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>137083</th>\n",
       "      <td>chr1</td>\n",
       "      <td>27675815</td>\n",
       "      <td>27676462</td>\n",
       "      <td>genic|chr1:27675665-27676612</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.483270</td>\n",
       "      <td>SYTL1</td>\n",
       "      <td>27668482</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.498696</td>\n",
       "      <td>...</td>\n",
       "      <td>0.086357</td>\n",
       "      <td>0.085205</td>\n",
       "      <td>0.023264</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.024480</td>\n",
       "      <td>0.060791</td>\n",
       "      <td>0.019766</td>\n",
       "      <td>0.211587</td>\n",
       "      <td>0.044781</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>141823</th>\n",
       "      <td>chr1</td>\n",
       "      <td>28190667</td>\n",
       "      <td>28190867</td>\n",
       "      <td>intergenic|chr1:28190517-28191017</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.315472</td>\n",
       "      <td>THEMIS2</td>\n",
       "      <td>28199053</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.512669</td>\n",
       "      <td>...</td>\n",
       "      <td>0.080589</td>\n",
       "      <td>0.079546</td>\n",
       "      <td>0.036897</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.038113</td>\n",
       "      <td>0.088249</td>\n",
       "      <td>0.029721</td>\n",
       "      <td>0.184185</td>\n",
       "      <td>0.039754</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>142126</th>\n",
       "      <td>chr1</td>\n",
       "      <td>28200232</td>\n",
       "      <td>28200432</td>\n",
       "      <td>genic|chr1:28200082-28200582</td>\n",
       "      <td>genic</td>\n",
       "      <td>3.332576</td>\n",
       "      <td>THEMIS2</td>\n",
       "      <td>28199053</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.512669</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.044303</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.045519</td>\n",
       "      <td>0.151695</td>\n",
       "      <td>0.051089</td>\n",
       "      <td>0.411360</td>\n",
       "      <td>0.088786</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>142327</th>\n",
       "      <td>chr1</td>\n",
       "      <td>28208995</td>\n",
       "      <td>28209195</td>\n",
       "      <td>genic|chr1:28208845-28209345</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.527351</td>\n",
       "      <td>THEMIS2</td>\n",
       "      <td>28199053</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.512669</td>\n",
       "      <td>...</td>\n",
       "      <td>0.068115</td>\n",
       "      <td>0.067298</td>\n",
       "      <td>0.018275</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.019491</td>\n",
       "      <td>0.049261</td>\n",
       "      <td>0.016591</td>\n",
       "      <td>0.170085</td>\n",
       "      <td>0.036711</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691341</th>\n",
       "      <td>chrX</td>\n",
       "      <td>152753232</td>\n",
       "      <td>152753817</td>\n",
       "      <td>intergenic|chrX:152753082-152753967</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>6.322386</td>\n",
       "      <td>BGN</td>\n",
       "      <td>152760346</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.477025</td>\n",
       "      <td>...</td>\n",
       "      <td>0.095538</td>\n",
       "      <td>0.094209</td>\n",
       "      <td>0.018809</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.020025</td>\n",
       "      <td>0.126604</td>\n",
       "      <td>0.067286</td>\n",
       "      <td>0.595623</td>\n",
       "      <td>0.155744</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691431</th>\n",
       "      <td>chrX</td>\n",
       "      <td>152761003</td>\n",
       "      <td>152761203</td>\n",
       "      <td>genic|chrX:152760853-152761353</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.675277</td>\n",
       "      <td>BGN</td>\n",
       "      <td>152760346</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.477025</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045554</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046770</td>\n",
       "      <td>0.078353</td>\n",
       "      <td>0.041642</td>\n",
       "      <td>0.206790</td>\n",
       "      <td>0.054071</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691441</th>\n",
       "      <td>chrX</td>\n",
       "      <td>152761935</td>\n",
       "      <td>152762135</td>\n",
       "      <td>genic|chrX:152761785-152762285</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.468697</td>\n",
       "      <td>BGN</td>\n",
       "      <td>152760346</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.477025</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.045554</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.046770</td>\n",
       "      <td>0.068691</td>\n",
       "      <td>0.036507</td>\n",
       "      <td>0.181290</td>\n",
       "      <td>0.047404</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691458</th>\n",
       "      <td>chrX</td>\n",
       "      <td>152765208</td>\n",
       "      <td>152765408</td>\n",
       "      <td>genic|chrX:152765058-152765558</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.438929</td>\n",
       "      <td>BGN</td>\n",
       "      <td>152760346</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.477025</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.034625</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.035841</td>\n",
       "      <td>0.051572</td>\n",
       "      <td>0.027409</td>\n",
       "      <td>0.177616</td>\n",
       "      <td>0.046443</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7691473</th>\n",
       "      <td>chrX</td>\n",
       "      <td>152772847</td>\n",
       "      <td>152773047</td>\n",
       "      <td>genic|chrX:152772697-152773197</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.916468</td>\n",
       "      <td>BGN</td>\n",
       "      <td>152760346</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.477025</td>\n",
       "      <td>...</td>\n",
       "      <td>0.055846</td>\n",
       "      <td>0.055238</td>\n",
       "      <td>0.019684</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.020900</td>\n",
       "      <td>0.040054</td>\n",
       "      <td>0.021287</td>\n",
       "      <td>0.105862</td>\n",
       "      <td>0.027681</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>431 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                                 name  \\\n",
       "136939   chr1   27669515   27669715         genic|chr1:27669365-27669865   \n",
       "137083   chr1   27675815   27676462         genic|chr1:27675665-27676612   \n",
       "141823   chr1   28190667   28190867    intergenic|chr1:28190517-28191017   \n",
       "142126   chr1   28200232   28200432         genic|chr1:28200082-28200582   \n",
       "142327   chr1   28208995   28209195         genic|chr1:28208845-28209345   \n",
       "...       ...        ...        ...                                  ...   \n",
       "7691341  chrX  152753232  152753817  intergenic|chrX:152753082-152753967   \n",
       "7691431  chrX  152761003  152761203       genic|chrX:152760853-152761353   \n",
       "7691441  chrX  152761935  152762135       genic|chrX:152761785-152762285   \n",
       "7691458  chrX  152765208  152765408       genic|chrX:152765058-152765558   \n",
       "7691473  chrX  152772847  152773047       genic|chrX:152772697-152773197   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "136939        genic       2.908763      SYTL1       27668482   \n",
       "137083        genic       2.483270      SYTL1       27668482   \n",
       "141823   intergenic       2.315472    THEMIS2       28199053   \n",
       "142126        genic       3.332576    THEMIS2       28199053   \n",
       "142327        genic       2.527351    THEMIS2       28199053   \n",
       "...             ...            ...        ...            ...   \n",
       "7691341  intergenic       6.322386        BGN      152760346   \n",
       "7691431       genic       1.675277        BGN      152760346   \n",
       "7691441       genic       1.468697        BGN      152760346   \n",
       "7691458       genic       1.438929        BGN      152760346   \n",
       "7691473       genic       1.916468        BGN      152760346   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "136939                    NaN                            0.498696  ...   \n",
       "137083                    NaN                            0.498696  ...   \n",
       "141823                    NaN                            0.512669  ...   \n",
       "142126                    NaN                            0.512669  ...   \n",
       "142327                    NaN                            0.512669  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7691341                   NaN                            0.477025  ...   \n",
       "7691431                   NaN                            0.477025  ...   \n",
       "7691441                   NaN                            0.477025  ...   \n",
       "7691458                   NaN                            0.477025  ...   \n",
       "7691473                   NaN                            0.477025  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "136939           0.125372                    0.123436               0.035278   \n",
       "137083           0.086357                    0.085205               0.023264   \n",
       "141823           0.080589                    0.079546               0.036897   \n",
       "142126           0.125372                    0.123436               0.044303   \n",
       "142327           0.068115                    0.067298               0.018275   \n",
       "...                   ...                         ...                    ...   \n",
       "7691341          0.095538                    0.094209               0.018809   \n",
       "7691431          0.125372                    0.123436               0.045554   \n",
       "7691441          0.125372                    0.123436               0.045554   \n",
       "7691458          0.125372                    0.123436               0.034625   \n",
       "7691473          0.055846                    0.055238               0.019684   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "136939          0.001216                   0.036494             0.106152   \n",
       "137083          0.001216                   0.024480             0.060791   \n",
       "141823          0.001216                   0.038113             0.088249   \n",
       "142126          0.001216                   0.045519             0.151695   \n",
       "142327          0.001216                   0.019491             0.049261   \n",
       "...                  ...                        ...                  ...   \n",
       "7691341         0.001216                   0.020025             0.126604   \n",
       "7691431         0.001216                   0.046770             0.078353   \n",
       "7691441         0.001216                   0.046770             0.068691   \n",
       "7691458         0.001216                   0.035841             0.051572   \n",
       "7691473         0.001216                   0.020900             0.040054   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "136939    0.034515                  0.359046        0.075990  liver-ENCODE  \n",
       "137083    0.019766                  0.211587        0.044781  liver-ENCODE  \n",
       "141823    0.029721                  0.184185        0.039754  liver-ENCODE  \n",
       "142126    0.051089                  0.411360        0.088786  liver-ENCODE  \n",
       "142327    0.016591                  0.170085        0.036711  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7691341   0.067286                  0.595623        0.155744  liver-ENCODE  \n",
       "7691431   0.041642                  0.206790        0.054071  liver-ENCODE  \n",
       "7691441   0.036507                  0.181290        0.047404  liver-ENCODE  \n",
       "7691458   0.027409                  0.177616        0.046443  liver-ENCODE  \n",
       "7691473   0.021287                  0.105862        0.027681  liver-ENCODE  \n",
       "\n",
       "[431 rows x 24 columns]"
      ]
     },
     "execution_count": 95,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>5068390</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>NBEA</td>\n",
       "      <td>35516423</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.511365</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000721</td>\n",
       "      <td>0.000731</td>\n",
       "      <td>0.000061</td>\n",
       "      <td>0.000721</td>\n",
       "      <td>0.000782</td>\n",
       "      <td>0.014937</td>\n",
       "      <td>0.020186</td>\n",
       "      <td>0.013963</td>\n",
       "      <td>0.014332</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068391</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>RFC3</td>\n",
       "      <td>34392205</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.521387</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001678</td>\n",
       "      <td>0.001693</td>\n",
       "      <td>0.000259</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.001475</td>\n",
       "      <td>0.028162</td>\n",
       "      <td>0.026526</td>\n",
       "      <td>0.032328</td>\n",
       "      <td>0.020164</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068392</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>KL</td>\n",
       "      <td>33590570</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.433681</td>\n",
       "      <td>...</td>\n",
       "      <td>0.008408</td>\n",
       "      <td>0.008407</td>\n",
       "      <td>0.001714</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002930</td>\n",
       "      <td>0.055929</td>\n",
       "      <td>0.043236</td>\n",
       "      <td>0.160499</td>\n",
       "      <td>0.080289</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068393</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>PDS5B</td>\n",
       "      <td>33160563</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.863696</td>\n",
       "      <td>...</td>\n",
       "      <td>0.002086</td>\n",
       "      <td>0.002102</td>\n",
       "      <td>0.001469</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002685</td>\n",
       "      <td>0.051252</td>\n",
       "      <td>0.029175</td>\n",
       "      <td>0.040129</td>\n",
       "      <td>0.011981</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068394</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>N4BP2L2</td>\n",
       "      <td>33112932</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.853186</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001937</td>\n",
       "      <td>0.001953</td>\n",
       "      <td>0.000561</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.001776</td>\n",
       "      <td>0.033912</td>\n",
       "      <td>0.017204</td>\n",
       "      <td>0.037283</td>\n",
       "      <td>0.009275</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068395</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>BRCA2</td>\n",
       "      <td>32889616</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.802795</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001461</td>\n",
       "      <td>0.001475</td>\n",
       "      <td>0.000577</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.001793</td>\n",
       "      <td>0.034229</td>\n",
       "      <td>0.017171</td>\n",
       "      <td>0.028168</td>\n",
       "      <td>0.007686</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068397</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>FRY-AS1</td>\n",
       "      <td>32605776</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.651418</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001123</td>\n",
       "      <td>0.001136</td>\n",
       "      <td>0.000378</td>\n",
       "      <td>0.001123</td>\n",
       "      <td>0.001501</td>\n",
       "      <td>0.028659</td>\n",
       "      <td>0.017591</td>\n",
       "      <td>0.021693</td>\n",
       "      <td>0.008985</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5068398</th>\n",
       "      <td>chr13</td>\n",
       "      <td>33700000</td>\n",
       "      <td>33700529</td>\n",
       "      <td>genic|chr13:33699850-33700679</td>\n",
       "      <td>genic</td>\n",
       "      <td>19.090405</td>\n",
       "      <td>FRY</td>\n",
       "      <td>32605436</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.649910</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001123</td>\n",
       "      <td>0.001136</td>\n",
       "      <td>0.000378</td>\n",
       "      <td>0.001123</td>\n",
       "      <td>0.001501</td>\n",
       "      <td>0.028654</td>\n",
       "      <td>0.017587</td>\n",
       "      <td>0.021687</td>\n",
       "      <td>0.008980</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>8 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           chr     start       end                           name  class  \\\n",
       "5068390  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068391  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068392  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068393  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068394  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068395  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068397  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "5068398  chr13  33700000  33700529  genic|chr13:33699850-33700679  genic   \n",
       "\n",
       "         activity_base TargetGene  TargetGeneTSS  TargetGeneExpression  \\\n",
       "5068390      19.090405       NBEA       35516423                   NaN   \n",
       "5068391      19.090405       RFC3       34392205                   NaN   \n",
       "5068392      19.090405         KL       33590570                   NaN   \n",
       "5068393      19.090405      PDS5B       33160563                   NaN   \n",
       "5068394      19.090405    N4BP2L2       33112932                   NaN   \n",
       "5068395      19.090405      BRCA2       32889616                   NaN   \n",
       "5068397      19.090405    FRY-AS1       32605776                   NaN   \n",
       "5068398      19.090405        FRY       32605436                   NaN   \n",
       "\n",
       "         TargetGenePromoterActivityQuantile  ...  powerlaw_contact  \\\n",
       "5068390                            0.511365  ...          0.000721   \n",
       "5068391                            0.521387  ...          0.001678   \n",
       "5068392                            0.433681  ...          0.008408   \n",
       "5068393                            0.863696  ...          0.002086   \n",
       "5068394                            0.853186  ...          0.001937   \n",
       "5068395                            0.802795  ...          0.001461   \n",
       "5068397                            0.651418  ...          0.001123   \n",
       "5068398                            0.649910  ...          0.001123   \n",
       "\n",
       "         powerlaw_contact_reference  hic_contact_pl_scaled  hic_pseudocount  \\\n",
       "5068390                    0.000731               0.000061         0.000721   \n",
       "5068391                    0.001693               0.000259         0.001216   \n",
       "5068392                    0.008407               0.001714         0.001216   \n",
       "5068393                    0.002102               0.001469         0.001216   \n",
       "5068394                    0.001953               0.000561         0.001216   \n",
       "5068395                    0.001475               0.000577         0.001216   \n",
       "5068397                    0.001136               0.000378         0.001123   \n",
       "5068398                    0.001136               0.000378         0.001123   \n",
       "\n",
       "         hic_contact_pl_scaled_adj  ABC.Score.Numerator  ABC.Score  \\\n",
       "5068390                   0.000782             0.014937   0.020186   \n",
       "5068391                   0.001475             0.028162   0.026526   \n",
       "5068392                   0.002930             0.055929   0.043236   \n",
       "5068393                   0.002685             0.051252   0.029175   \n",
       "5068394                   0.001776             0.033912   0.017204   \n",
       "5068395                   0.001793             0.034229   0.017171   \n",
       "5068397                   0.001501             0.028659   0.017591   \n",
       "5068398                   0.001501             0.028654   0.017587   \n",
       "\n",
       "         powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "5068390                  0.013963        0.014332  liver-ENCODE  \n",
       "5068391                  0.032328        0.020164  liver-ENCODE  \n",
       "5068392                  0.160499        0.080289  liver-ENCODE  \n",
       "5068393                  0.040129        0.011981  liver-ENCODE  \n",
       "5068394                  0.037283        0.009275  liver-ENCODE  \n",
       "5068395                  0.028168        0.007686  liver-ENCODE  \n",
       "5068397                  0.021693        0.008985  liver-ENCODE  \n",
       "5068398                  0.021687        0.008980  liver-ENCODE  \n",
       "\n",
       "[8 rows x 24 columns]"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df[df['name']=='genic|chr13:33699850-33700679']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "import os,re,sys\n",
    "import pandas as pd\n",
    "def gene2names(gene_file,name_flie):\n",
    "    gene = pd.read_csv(gene_file,sep='\\t',header=0,index_col=None)\n",
    "    name = pd.read_csv(name_flie,sep='\\t',header=None,index_col=None,names=['gene_id','name','expression'])\n",
    "    name_dic = name.set_index('gene_id').to_dict()['name']\n",
    "    gene['gene_name'] = gene['gene_id'].map(name_dic)\n",
    "    gene.to_csv(f'{sys.argv[0]}.gene_name',sep='\\t',header=True,index=False)\n",
    "    return None\n",
    "gene2names('08.peak_annotation.py.out','lamprey.tsv')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [],
   "source": [
    "vertebrate_enhancers = pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer/vertebrate_enhancer_origin_info.txt', sep='\\t', header=0,index_col=None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(['HES4', 'PROX1', 'SLC7A11', 'PTGER4', 'DPYSL2', 'MYC', 'PHLDA1',\n",
       "       'SLCO3A1', 'PTGIR'], dtype=object)"
      ]
     },
     "execution_count": 69,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_high_variable[df_liver_high_variable['TargetGene'].isin(vertebrate_enhancers['gene_name'])]['TargetGene'].unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {},
   "outputs": [],
   "source": [
    "jisuo_enhancers = pd.read_csv('H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer/jisuo.origin.fa', sep='\\t', header=0,index_col=None)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_liver_high_variable[df_liver_high_variable['TargetGene'].isin(jisuo_enhancers['gene_name'])]['TargetGene'].unique()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2140780</th>\n",
       "      <td>chr5</td>\n",
       "      <td>39624987</td>\n",
       "      <td>39625187</td>\n",
       "      <td>intergenic|chr5:39624837-39625337</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>11.560091</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.001160</td>\n",
       "      <td>0.001173</td>\n",
       "      <td>0.000995</td>\n",
       "      <td>0.001160</td>\n",
       "      <td>0.002155</td>\n",
       "      <td>0.024912</td>\n",
       "      <td>0.021164</td>\n",
       "      <td>0.013564</td>\n",
       "      <td>0.004774</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2141531</th>\n",
       "      <td>chr5</td>\n",
       "      <td>40679683</td>\n",
       "      <td>40680131</td>\n",
       "      <td>promoter|chr5:40679533-40680281</td>\n",
       "      <td>promoter</td>\n",
       "      <td>8.104650</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.013519</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.014734</td>\n",
       "      <td>0.119417</td>\n",
       "      <td>0.101450</td>\n",
       "      <td>1.000406</td>\n",
       "      <td>0.352120</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2141616</th>\n",
       "      <td>chr5</td>\n",
       "      <td>40681303</td>\n",
       "      <td>40681658</td>\n",
       "      <td>genic|chr5:40681153-40681808</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.064353</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.027497</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.028712</td>\n",
       "      <td>0.116698</td>\n",
       "      <td>0.099140</td>\n",
       "      <td>0.501688</td>\n",
       "      <td>0.176583</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2141818</th>\n",
       "      <td>chr5</td>\n",
       "      <td>40688468</td>\n",
       "      <td>40688668</td>\n",
       "      <td>genic|chr5:40688318-40688818</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.658705</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.078513</td>\n",
       "      <td>0.077507</td>\n",
       "      <td>0.027570</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.028786</td>\n",
       "      <td>0.047748</td>\n",
       "      <td>0.040564</td>\n",
       "      <td>0.128561</td>\n",
       "      <td>0.045251</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2142167</th>\n",
       "      <td>chr5</td>\n",
       "      <td>40784222</td>\n",
       "      <td>40784615</td>\n",
       "      <td>genic|chr5:40784072-40784765</td>\n",
       "      <td>genic</td>\n",
       "      <td>7.167615</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.008781</td>\n",
       "      <td>0.008778</td>\n",
       "      <td>0.001662</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.002878</td>\n",
       "      <td>0.020630</td>\n",
       "      <td>0.017526</td>\n",
       "      <td>0.062917</td>\n",
       "      <td>0.022145</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr     start       end                               name  \\\n",
       "2140780  chr5  39624987  39625187  intergenic|chr5:39624837-39625337   \n",
       "2141531  chr5  40679683  40680131    promoter|chr5:40679533-40680281   \n",
       "2141616  chr5  40681303  40681658       genic|chr5:40681153-40681808   \n",
       "2141818  chr5  40688468  40688668       genic|chr5:40688318-40688818   \n",
       "2142167  chr5  40784222  40784615       genic|chr5:40784072-40784765   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "2140780  intergenic      11.560091     PTGER4       40680031   \n",
       "2141531    promoter       8.104650     PTGER4       40680031   \n",
       "2141616       genic       4.064353     PTGER4       40680031   \n",
       "2141818       genic       1.658705     PTGER4       40680031   \n",
       "2142167       genic       7.167615     PTGER4       40680031   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "2140780                   NaN                            0.750081  ...   \n",
       "2141531                   NaN                            0.750081  ...   \n",
       "2141616                   NaN                            0.750081  ...   \n",
       "2141818                   NaN                            0.750081  ...   \n",
       "2142167                   NaN                            0.750081  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "2140780          0.001160                    0.001173               0.000995   \n",
       "2141531          0.125372                    0.123436               0.013519   \n",
       "2141616          0.125372                    0.123436               0.027497   \n",
       "2141818          0.078513                    0.077507               0.027570   \n",
       "2142167          0.008781                    0.008778               0.001662   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "2140780         0.001160                   0.002155             0.024912   \n",
       "2141531         0.001216                   0.014734             0.119417   \n",
       "2141616         0.001216                   0.028712             0.116698   \n",
       "2141818         0.001216                   0.028786             0.047748   \n",
       "2142167         0.001216                   0.002878             0.020630   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "2140780   0.021164                  0.013564        0.004774  liver-ENCODE  \n",
       "2141531   0.101450                  1.000406        0.352120  liver-ENCODE  \n",
       "2141616   0.099140                  0.501688        0.176583  liver-ENCODE  \n",
       "2141818   0.040564                  0.128561        0.045251  liver-ENCODE  \n",
       "2142167   0.017526                  0.062917        0.022145  liver-ENCODE  \n",
       "\n",
       "[5 rows x 24 columns]"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver[df_liver['TargetGene']=='PTGER4']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>100</th>\n",
       "      <td>chr1</td>\n",
       "      <td>713968</td>\n",
       "      <td>714276</td>\n",
       "      <td>promoter|chr1:713818-714426</td>\n",
       "      <td>promoter</td>\n",
       "      <td>3.897799</td>\n",
       "      <td>LOC100288069</td>\n",
       "      <td>714068</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.573692</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.485868</td>\n",
       "      <td>0.265022</td>\n",
       "      <td>0.481129</td>\n",
       "      <td>0.190119</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>369</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762795</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762645-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.909910</td>\n",
       "      <td>LINC00115</td>\n",
       "      <td>762902</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.506701</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.362726</td>\n",
       "      <td>0.126307</td>\n",
       "      <td>0.359188</td>\n",
       "      <td>0.132739</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>370</th>\n",
       "      <td>chr1</td>\n",
       "      <td>762795</td>\n",
       "      <td>763070</td>\n",
       "      <td>promoter|chr1:762645-763220</td>\n",
       "      <td>promoter</td>\n",
       "      <td>2.909910</td>\n",
       "      <td>LINC01128</td>\n",
       "      <td>762970</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.506701</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.124652</td>\n",
       "      <td>0.362726</td>\n",
       "      <td>0.126306</td>\n",
       "      <td>0.359188</td>\n",
       "      <td>0.132729</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>878</th>\n",
       "      <td>chr1</td>\n",
       "      <td>859329</td>\n",
       "      <td>859529</td>\n",
       "      <td>intergenic|chr1:859179-859679</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.211608</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.032714</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.033930</td>\n",
       "      <td>0.075040</td>\n",
       "      <td>0.024202</td>\n",
       "      <td>0.272992</td>\n",
       "      <td>0.055114</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>884</th>\n",
       "      <td>chr1</td>\n",
       "      <td>860018</td>\n",
       "      <td>860218</td>\n",
       "      <td>intergenic|chr1:859868-860368</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>3.423063</td>\n",
       "      <td>NOC2L</td>\n",
       "      <td>894679</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.533078</td>\n",
       "      <td>...</td>\n",
       "      <td>0.023099</td>\n",
       "      <td>0.022964</td>\n",
       "      <td>0.012852</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.014068</td>\n",
       "      <td>0.048156</td>\n",
       "      <td>0.015482</td>\n",
       "      <td>0.078606</td>\n",
       "      <td>0.014794</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716611</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>CMC4</td>\n",
       "      <td>154299547</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.653108</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232279</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.390153</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716612</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>MTCP1</td>\n",
       "      <td>154299547</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.653108</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232279</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.390153</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716613</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154299412</td>\n",
       "      <td>154300250</td>\n",
       "      <td>promoter|chrX:154299262-154300400</td>\n",
       "      <td>promoter</td>\n",
       "      <td>6.988388</td>\n",
       "      <td>BRCC3</td>\n",
       "      <td>154299694</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.646122</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.041811</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.043026</td>\n",
       "      <td>0.300685</td>\n",
       "      <td>0.232282</td>\n",
       "      <td>0.862619</td>\n",
       "      <td>0.391263</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7716826</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154444600</td>\n",
       "      <td>154444811</td>\n",
       "      <td>promoter|chrX:154444450-154444961</td>\n",
       "      <td>promoter</td>\n",
       "      <td>4.027038</td>\n",
       "      <td>VBP1</td>\n",
       "      <td>154444700</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.658709</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.037211</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.038427</td>\n",
       "      <td>0.154747</td>\n",
       "      <td>0.172342</td>\n",
       "      <td>0.497082</td>\n",
       "      <td>0.370809</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7717096</th>\n",
       "      <td>chrX</td>\n",
       "      <td>154842172</td>\n",
       "      <td>154842722</td>\n",
       "      <td>promoter|chrX:154842022-154842872</td>\n",
       "      <td>promoter</td>\n",
       "      <td>5.566828</td>\n",
       "      <td>TMLHE</td>\n",
       "      <td>154842622</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.582003</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.037942</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.039158</td>\n",
       "      <td>0.217984</td>\n",
       "      <td>0.264255</td>\n",
       "      <td>0.687147</td>\n",
       "      <td>0.545381</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>62122 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                               name  \\\n",
       "100      chr1     713968     714276        promoter|chr1:713818-714426   \n",
       "369      chr1     762795     763070        promoter|chr1:762645-763220   \n",
       "370      chr1     762795     763070        promoter|chr1:762645-763220   \n",
       "878      chr1     859329     859529      intergenic|chr1:859179-859679   \n",
       "884      chr1     860018     860218      intergenic|chr1:859868-860368   \n",
       "...       ...        ...        ...                                ...   \n",
       "7716611  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716612  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716613  chrX  154299412  154300250  promoter|chrX:154299262-154300400   \n",
       "7716826  chrX  154444600  154444811  promoter|chrX:154444450-154444961   \n",
       "7717096  chrX  154842172  154842722  promoter|chrX:154842022-154842872   \n",
       "\n",
       "              class  activity_base    TargetGene  TargetGeneTSS  \\\n",
       "100        promoter       3.897799  LOC100288069         714068   \n",
       "369        promoter       2.909910     LINC00115         762902   \n",
       "370        promoter       2.909910     LINC01128         762970   \n",
       "878      intergenic       2.211608        SAMD11         861120   \n",
       "884      intergenic       3.423063         NOC2L         894679   \n",
       "...             ...            ...           ...            ...   \n",
       "7716611    promoter       6.988388          CMC4      154299547   \n",
       "7716612    promoter       6.988388         MTCP1      154299547   \n",
       "7716613    promoter       6.988388         BRCC3      154299694   \n",
       "7716826    promoter       4.027038          VBP1      154444700   \n",
       "7717096    promoter       5.566828         TMLHE      154842622   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "100                       NaN                            0.573692  ...   \n",
       "369                       NaN                            0.506701  ...   \n",
       "370                       NaN                            0.506701  ...   \n",
       "878                       NaN                            0.425330  ...   \n",
       "884                       NaN                            0.533078  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7716611                   NaN                            0.653108  ...   \n",
       "7716612                   NaN                            0.653108  ...   \n",
       "7716613                   NaN                            0.646122  ...   \n",
       "7716826                   NaN                            0.658709  ...   \n",
       "7717096                   NaN                            0.582003  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "100              0.125372                    0.123436               0.123436   \n",
       "369              0.125372                    0.123436               0.123436   \n",
       "370              0.125372                    0.123436               0.123436   \n",
       "878              0.125372                    0.123436               0.032714   \n",
       "884              0.023099                    0.022964               0.012852   \n",
       "...                   ...                         ...                    ...   \n",
       "7716611          0.125372                    0.123436               0.041811   \n",
       "7716612          0.125372                    0.123436               0.041811   \n",
       "7716613          0.125372                    0.123436               0.041811   \n",
       "7716826          0.125372                    0.123436               0.037211   \n",
       "7717096          0.125372                    0.123436               0.037942   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "100             0.001216                   0.124652             0.485868   \n",
       "369             0.001216                   0.124652             0.362726   \n",
       "370             0.001216                   0.124652             0.362726   \n",
       "878             0.001216                   0.033930             0.075040   \n",
       "884             0.001216                   0.014068             0.048156   \n",
       "...                  ...                        ...                  ...   \n",
       "7716611         0.001216                   0.043026             0.300685   \n",
       "7716612         0.001216                   0.043026             0.300685   \n",
       "7716613         0.001216                   0.043026             0.300685   \n",
       "7716826         0.001216                   0.038427             0.154747   \n",
       "7717096         0.001216                   0.039158             0.217984   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "100       0.265022                  0.481129        0.190119  liver-ENCODE  \n",
       "369       0.126307                  0.359188        0.132739  liver-ENCODE  \n",
       "370       0.126306                  0.359188        0.132729  liver-ENCODE  \n",
       "878       0.024202                  0.272992        0.055114  liver-ENCODE  \n",
       "884       0.015482                  0.078606        0.014794  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7716611   0.232279                  0.862619        0.390153  liver-ENCODE  \n",
       "7716612   0.232279                  0.862619        0.390153  liver-ENCODE  \n",
       "7716613   0.232282                  0.862619        0.391263  liver-ENCODE  \n",
       "7716826   0.172342                  0.497082        0.370809  liver-ENCODE  \n",
       "7717096   0.264255                  0.687147        0.545381  liver-ENCODE  \n",
       "\n",
       "[62122 rows x 24 columns]"
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>878</th>\n",
       "      <td>chr1</td>\n",
       "      <td>859329</td>\n",
       "      <td>859529</td>\n",
       "      <td>intergenic|chr1:859179-859679</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>2.211608</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.032714</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.033930</td>\n",
       "      <td>0.075040</td>\n",
       "      <td>0.024202</td>\n",
       "      <td>0.272992</td>\n",
       "      <td>0.055114</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>885</th>\n",
       "      <td>chr1</td>\n",
       "      <td>860018</td>\n",
       "      <td>860218</td>\n",
       "      <td>intergenic|chr1:859868-860368</td>\n",
       "      <td>intergenic</td>\n",
       "      <td>3.423063</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.067077</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.068293</td>\n",
       "      <td>0.233771</td>\n",
       "      <td>0.075396</td>\n",
       "      <td>0.422529</td>\n",
       "      <td>0.085304</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>952</th>\n",
       "      <td>chr1</td>\n",
       "      <td>872461</td>\n",
       "      <td>872661</td>\n",
       "      <td>genic|chr1:872311-872811</td>\n",
       "      <td>genic</td>\n",
       "      <td>1.648334</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.060770</td>\n",
       "      <td>0.060079</td>\n",
       "      <td>0.028332</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.029548</td>\n",
       "      <td>0.048705</td>\n",
       "      <td>0.015708</td>\n",
       "      <td>0.099031</td>\n",
       "      <td>0.019993</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>960</th>\n",
       "      <td>chr1</td>\n",
       "      <td>873510</td>\n",
       "      <td>874703</td>\n",
       "      <td>genic|chr1:873360-874853</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.904584</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.054393</td>\n",
       "      <td>0.053809</td>\n",
       "      <td>0.028350</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.029566</td>\n",
       "      <td>0.145009</td>\n",
       "      <td>0.046768</td>\n",
       "      <td>0.263910</td>\n",
       "      <td>0.053281</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1129</th>\n",
       "      <td>chr1</td>\n",
       "      <td>876090</td>\n",
       "      <td>876290</td>\n",
       "      <td>genic|chr1:875940-876440</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.549120</td>\n",
       "      <td>SAMD11</td>\n",
       "      <td>861120</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.425330</td>\n",
       "      <td>...</td>\n",
       "      <td>0.047753</td>\n",
       "      <td>0.047276</td>\n",
       "      <td>0.019941</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.021157</td>\n",
       "      <td>0.096244</td>\n",
       "      <td>0.031041</td>\n",
       "      <td>0.215063</td>\n",
       "      <td>0.043419</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7705306</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153601838</td>\n",
       "      <td>153602038</td>\n",
       "      <td>genic|chrX:153601688-153602188</td>\n",
       "      <td>genic</td>\n",
       "      <td>2.298913</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.085079</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.086294</td>\n",
       "      <td>0.198383</td>\n",
       "      <td>0.049837</td>\n",
       "      <td>0.283769</td>\n",
       "      <td>0.045245</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7705432</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153607015</td>\n",
       "      <td>153607696</td>\n",
       "      <td>promoter|chrX:153606865-153607846</td>\n",
       "      <td>promoter</td>\n",
       "      <td>15.685993</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.085079</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.086294</td>\n",
       "      <td>1.353612</td>\n",
       "      <td>0.340049</td>\n",
       "      <td>1.936217</td>\n",
       "      <td>0.308714</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7706120</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153627057</td>\n",
       "      <td>153627257</td>\n",
       "      <td>genic|chrX:153626907-153627407</td>\n",
       "      <td>genic</td>\n",
       "      <td>6.861326</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.031608</td>\n",
       "      <td>0.031366</td>\n",
       "      <td>0.013048</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.014264</td>\n",
       "      <td>0.097868</td>\n",
       "      <td>0.024586</td>\n",
       "      <td>0.215209</td>\n",
       "      <td>0.034313</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7711195</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153770112</td>\n",
       "      <td>153770312</td>\n",
       "      <td>genic|chrX:153769962-153770462</td>\n",
       "      <td>genic</td>\n",
       "      <td>7.777634</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005815</td>\n",
       "      <td>0.005826</td>\n",
       "      <td>0.006523</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.007739</td>\n",
       "      <td>0.060193</td>\n",
       "      <td>0.015121</td>\n",
       "      <td>0.045315</td>\n",
       "      <td>0.007225</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7712395</th>\n",
       "      <td>chrX</td>\n",
       "      <td>153770892</td>\n",
       "      <td>153771374</td>\n",
       "      <td>genic|chrX:153770742-153771524</td>\n",
       "      <td>genic</td>\n",
       "      <td>12.060128</td>\n",
       "      <td>FLNA</td>\n",
       "      <td>153603006</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.437307</td>\n",
       "      <td>...</td>\n",
       "      <td>0.005787</td>\n",
       "      <td>0.005798</td>\n",
       "      <td>0.006524</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.007739</td>\n",
       "      <td>0.093339</td>\n",
       "      <td>0.023448</td>\n",
       "      <td>0.069930</td>\n",
       "      <td>0.011150</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>6106 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr      start        end                               name  \\\n",
       "878      chr1     859329     859529      intergenic|chr1:859179-859679   \n",
       "885      chr1     860018     860218      intergenic|chr1:859868-860368   \n",
       "952      chr1     872461     872661           genic|chr1:872311-872811   \n",
       "960      chr1     873510     874703           genic|chr1:873360-874853   \n",
       "1129     chr1     876090     876290           genic|chr1:875940-876440   \n",
       "...       ...        ...        ...                                ...   \n",
       "7705306  chrX  153601838  153602038     genic|chrX:153601688-153602188   \n",
       "7705432  chrX  153607015  153607696  promoter|chrX:153606865-153607846   \n",
       "7706120  chrX  153627057  153627257     genic|chrX:153626907-153627407   \n",
       "7711195  chrX  153770112  153770312     genic|chrX:153769962-153770462   \n",
       "7712395  chrX  153770892  153771374     genic|chrX:153770742-153771524   \n",
       "\n",
       "              class  activity_base TargetGene  TargetGeneTSS  \\\n",
       "878      intergenic       2.211608     SAMD11         861120   \n",
       "885      intergenic       3.423063     SAMD11         861120   \n",
       "952           genic       1.648334     SAMD11         861120   \n",
       "960           genic       4.904584     SAMD11         861120   \n",
       "1129          genic       4.549120     SAMD11         861120   \n",
       "...             ...            ...        ...            ...   \n",
       "7705306       genic       2.298913       FLNA      153603006   \n",
       "7705432    promoter      15.685993       FLNA      153603006   \n",
       "7706120       genic       6.861326       FLNA      153603006   \n",
       "7711195       genic       7.777634       FLNA      153603006   \n",
       "7712395       genic      12.060128       FLNA      153603006   \n",
       "\n",
       "         TargetGeneExpression  TargetGenePromoterActivityQuantile  ...  \\\n",
       "878                       NaN                            0.425330  ...   \n",
       "885                       NaN                            0.425330  ...   \n",
       "952                       NaN                            0.425330  ...   \n",
       "960                       NaN                            0.425330  ...   \n",
       "1129                      NaN                            0.425330  ...   \n",
       "...                       ...                                 ...  ...   \n",
       "7705306                   NaN                            0.437307  ...   \n",
       "7705432                   NaN                            0.437307  ...   \n",
       "7706120                   NaN                            0.437307  ...   \n",
       "7711195                   NaN                            0.437307  ...   \n",
       "7712395                   NaN                            0.437307  ...   \n",
       "\n",
       "         powerlaw_contact  powerlaw_contact_reference  hic_contact_pl_scaled  \\\n",
       "878              0.125372                    0.123436               0.032714   \n",
       "885              0.125372                    0.123436               0.067077   \n",
       "952              0.060770                    0.060079               0.028332   \n",
       "960              0.054393                    0.053809               0.028350   \n",
       "1129             0.047753                    0.047276               0.019941   \n",
       "...                   ...                         ...                    ...   \n",
       "7705306          0.125372                    0.123436               0.085079   \n",
       "7705432          0.125372                    0.123436               0.085079   \n",
       "7706120          0.031608                    0.031366               0.013048   \n",
       "7711195          0.005815                    0.005826               0.006523   \n",
       "7712395          0.005787                    0.005798               0.006524   \n",
       "\n",
       "         hic_pseudocount  hic_contact_pl_scaled_adj  ABC.Score.Numerator  \\\n",
       "878             0.001216                   0.033930             0.075040   \n",
       "885             0.001216                   0.068293             0.233771   \n",
       "952             0.001216                   0.029548             0.048705   \n",
       "960             0.001216                   0.029566             0.145009   \n",
       "1129            0.001216                   0.021157             0.096244   \n",
       "...                  ...                        ...                  ...   \n",
       "7705306         0.001216                   0.086294             0.198383   \n",
       "7705432         0.001216                   0.086294             1.353612   \n",
       "7706120         0.001216                   0.014264             0.097868   \n",
       "7711195         0.001216                   0.007739             0.060193   \n",
       "7712395         0.001216                   0.007739             0.093339   \n",
       "\n",
       "         ABC.Score  powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "878       0.024202                  0.272992        0.055114  liver-ENCODE  \n",
       "885       0.075396                  0.422529        0.085304  liver-ENCODE  \n",
       "952       0.015708                  0.099031        0.019993  liver-ENCODE  \n",
       "960       0.046768                  0.263910        0.053281  liver-ENCODE  \n",
       "1129      0.031041                  0.215063        0.043419  liver-ENCODE  \n",
       "...            ...                       ...             ...           ...  \n",
       "7705306   0.049837                  0.283769        0.045245  liver-ENCODE  \n",
       "7705432   0.340049                  1.936217        0.308714  liver-ENCODE  \n",
       "7706120   0.024586                  0.215209        0.034313  liver-ENCODE  \n",
       "7711195   0.015121                  0.045315        0.007225  liver-ENCODE  \n",
       "7712395   0.023448                  0.069930        0.011150  liver-ENCODE  \n",
       "\n",
       "[6106 rows x 24 columns]"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_high_variable"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "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>chr</th>\n",
       "      <th>start</th>\n",
       "      <th>end</th>\n",
       "      <th>name</th>\n",
       "      <th>class</th>\n",
       "      <th>activity_base</th>\n",
       "      <th>TargetGene</th>\n",
       "      <th>TargetGeneTSS</th>\n",
       "      <th>TargetGeneExpression</th>\n",
       "      <th>TargetGenePromoterActivityQuantile</th>\n",
       "      <th>...</th>\n",
       "      <th>powerlaw_contact</th>\n",
       "      <th>powerlaw_contact_reference</th>\n",
       "      <th>hic_contact_pl_scaled</th>\n",
       "      <th>hic_pseudocount</th>\n",
       "      <th>hic_contact_pl_scaled_adj</th>\n",
       "      <th>ABC.Score.Numerator</th>\n",
       "      <th>ABC.Score</th>\n",
       "      <th>powerlaw.Score.Numerator</th>\n",
       "      <th>powerlaw.Score</th>\n",
       "      <th>CellType</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2141616</th>\n",
       "      <td>chr5</td>\n",
       "      <td>40681303</td>\n",
       "      <td>40681658</td>\n",
       "      <td>genic|chr5:40681153-40681808</td>\n",
       "      <td>genic</td>\n",
       "      <td>4.064353</td>\n",
       "      <td>PTGER4</td>\n",
       "      <td>40680031</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.750081</td>\n",
       "      <td>...</td>\n",
       "      <td>0.125372</td>\n",
       "      <td>0.123436</td>\n",
       "      <td>0.027497</td>\n",
       "      <td>0.001216</td>\n",
       "      <td>0.028712</td>\n",
       "      <td>0.116698</td>\n",
       "      <td>0.09914</td>\n",
       "      <td>0.501688</td>\n",
       "      <td>0.176583</td>\n",
       "      <td>liver-ENCODE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1 rows × 24 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "          chr     start       end                          name  class  \\\n",
       "2141616  chr5  40681303  40681658  genic|chr5:40681153-40681808  genic   \n",
       "\n",
       "         activity_base TargetGene  TargetGeneTSS  TargetGeneExpression  \\\n",
       "2141616       4.064353     PTGER4       40680031                   NaN   \n",
       "\n",
       "         TargetGenePromoterActivityQuantile  ...  powerlaw_contact  \\\n",
       "2141616                            0.750081  ...          0.125372   \n",
       "\n",
       "         powerlaw_contact_reference  hic_contact_pl_scaled  hic_pseudocount  \\\n",
       "2141616                    0.123436               0.027497         0.001216   \n",
       "\n",
       "         hic_contact_pl_scaled_adj  ABC.Score.Numerator  ABC.Score  \\\n",
       "2141616                   0.028712             0.116698    0.09914   \n",
       "\n",
       "         powerlaw.Score.Numerator  powerlaw.Score      CellType  \n",
       "2141616                  0.501688        0.176583  liver-ENCODE  \n",
       "\n",
       "[1 rows x 24 columns]"
      ]
     },
     "execution_count": 53,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_liver_high_variable[df_liver_high_variable['name']=='genic|chr5:40681153-40681808']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "names =[]\n",
    "from Bio import SeqIO\n",
    "records = SeqIO.parse(\"H:/002.singcell/01.liver/10.mutileSpecies/05.SAMap/015.enhancer/liver_high_enhancer.fa\", \"fasta\")\n",
    "for record in records:\n",
    "    id = record.id\n",
    "    seq = record.seq\n",
    "    names.append(id)\n",
    "    print(id)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "liver_high = df_liver_high_variable[df_liver_high_variable['name'].isin(names)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Users\\wubaosheng\\AppData\\Local\\Temp\\ipykernel_25376\\4060257618.py:1: SettingWithCopyWarning: \n",
      "A value is trying to be set on a copy of a slice from a DataFrame\n",
      "\n",
      "See the caveats in the documentation: https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#returning-a-view-versus-a-copy\n",
      "  liver_high.drop_duplicates(subset=['name'],keep='first',inplace=True)\n"
     ]
    }
   ],
   "source": [
    "liver_high.drop_duplicates(subset=['name'],keep='first',inplace=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "liver_high[['chr','start','end','name','TargetGene','TargetGeneTSS','CellType']].to_csv('liver_high_enhancer.csv',sep='\\t',header=True,index=False)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.7"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}