analysis_round_2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 200,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The watermark extension is already loaded. To reload it, use:\n",
      "  %reload_ext watermark\n",
      "numpy   1.18.4\n",
      "seaborn 0.10.1\n",
      "pandas  1.0.4\n",
      "\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "import warnings\n",
    "import seaborn as sns\n",
    "\n",
    "from scipy.stats import pearsonr\n",
    "\n",
    "sns.set_style(\"darkgrid\")\n",
    "np.random.seed(930525)\n",
    "pd.set_option('display.max_columns', 20)\n",
    "pd.set_option('display.max_rows', 200)\n",
    "\n",
    "warnings.simplefilter('once')\n",
    "\n",
    "%matplotlib inline\n",
    "%load_ext watermark\n",
    "%watermark --iversions"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 201,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_folder = \"/mnt/btrfs/data/italy\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 202,
   "metadata": {},
   "outputs": [],
   "source": [
    "metadata_file = \"data/metadata.csv\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 203,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_metadata_file = pd.read_csv(metadata_file, sep=\",\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 204,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['sample_id', 'Sample_Number', 'Sample_Code', 'Shotgun', '16s',\n",
       "       'viral_load', 'nCoVN1', 'nCoVN2', 'RNase P', 'ConcN1', 'ConcN2', 'N1',\n",
       "       'Sex', 'Age', 'Age_range', 'Age_average', 'Blood_sampling_date',\n",
       "       'Swab date', 'Severe_Case', 'REUMA', 'PNEUMO', 'ONCO', 'HEMATO',\n",
       "       'DIABETES', 'METABOLIC', 'CARDIO', 'HYPERTENSION', 'NEPHRO', 'TB',\n",
       "       'HIV', 'IMMUNOSUPPRESSION', 'HEPATOPATHY', 'HCQ', 'LPV', 'HCQ_LPV',\n",
       "       'HCQ_LPV_01', 'HCQ_LPV_02', 'RDV', 'TCZ', 'BCT', 'Steroide',\n",
       "       'SYMPTOMS_ONSET', 'SYMPTOMS_TO_SWAB_DAYS', 'Swab_day_quartile',\n",
       "       'Set_severe', 'swab_days_avg', 'SD_q1_severe', 'SD_q1_HB02', 'swab_out',\n",
       "       'OUTCOME', 'Hospitalization_Days', 'CRP', 'CRP_range', 'LDH',\n",
       "       'LDH_range', 'Ferritin', 'CPK', 'CPK_range', 'WBC', 'WBC_range',\n",
       "       'Neutrophil_Count', 'AST', 'ALT', 'Albumin', 'D.Dimer', 'Lowest_SpO2',\n",
       "       'Hospital_bin', 'Hospital_bin_01', 'Hospital_bin_02', 'q1', 'bin_02',\n",
       "       'HCQ.1', 'LPV.1', 'HB_02_HCQ', 'HB_02_LVD', 'Out_age_range'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 204,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_metadata_file.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 205,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "No     68\n",
       "Yes    32\n",
       "Name: RDV, dtype: int64"
      ]
     },
     "execution_count": 205,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_metadata_file[\"RDV\"].value_counts()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 206,
   "metadata": {},
   "outputs": [],
   "source": [
    "taxatable = \"/mnt/btrfs/data/italy/taxatable/filter_qc_taxatable.txt\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 207,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_raw = pd.read_csv(taxatable, sep=\"\\t\", index_col=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 208,
   "metadata": {},
   "outputs": [],
   "source": [
    "columns = df_taxa_raw.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 209,
   "metadata": {},
   "outputs": [],
   "source": [
    "sample_ids = [\"_\".join(_.split(\".\")[:2]) if not _.startswith(\"pos\") else \"_\".join(_.split(\".\")[:3]) for _ in columns]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 210,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa = df_taxa_raw.copy().T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 211,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pos.dna.02.S35.001.fa    1.0\n",
      "CAPMA.83.S12.001.fa      1.0\n",
      "FOGI.63.S13.001.fa       1.0\n",
      "COCL.11.S7.001.fa        1.0\n",
      "WACH.38.S1.001.fa        1.0\n",
      "PRPI.42.S31.001.fa       1.0\n",
      "RODE.80.S8.001.fa        1.0\n",
      "pos.dna.01.S17.001.fa    1.0\n",
      "MEEN.93.S11.001.fa       1.0\n",
      "DETH.41.S27.001.fa       1.0\n",
      "TECA.45.S15.001.fa       1.0\n",
      "BODA.52.S9.001.fa        1.0\n",
      "CHMA.61.S3.001.fa        1.0\n",
      "ASTER.02.S34.001.fa      1.0\n",
      "BIRO.74.S19.001.fa       1.0\n",
      "CAAU.07.S2.001.fa        1.0\n",
      "DIGMA.71.S21.001.fa      1.0\n",
      "FAFE.16.S10.001.fa       1.0\n",
      "OPSE.53.S5.001.fa        1.0\n",
      "LAES.91.S29.001.fa       1.0\n",
      "DOEC.81.S30.001.fa       1.0\n",
      "EKJO.64.S32.001.fa       1.0\n",
      "DIAL.15.S23.001.fa       1.0\n",
      "GIFI.21.S22.001.fa       1.0\n",
      "TAOT.65.S24.001.fa       1.0\n",
      "DEFI.14.S6.001.fa        1.0\n",
      "MUNA.98.S28.001.fa       1.0\n",
      "ROAL.75.S33.001.fa       1.0\n",
      "MERI.28.S20.001.fa       1.0\n",
      "BUMA.05.S4.001.fa        1.0\n",
      "MOEM.48.S16.001.fa       1.0\n",
      "PAMA.46.S18.001.fa       1.0\n",
      "VUGE.37.S14.001.fa       1.0\n",
      "TOSA.76.S25.001.fa       1.0\n",
      "DEIS.70.S26.001.fa       1.0\n",
      "dtype: float64\n"
     ]
    }
   ],
   "source": [
    "df_taxa_ra = df_taxa.apply(lambda x: x / x.sum(), axis=1)\n",
    "\n",
    "print(df_taxa_ra.sum(axis=1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 212,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "#OTU ID\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_atypica             2.851332e-03\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_rogosae             3.755483e-03\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_parvula             1.260326e-03\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_dispar              2.185924e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_tobetsuensis        2.696719e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_parvula_A           2.443783e-03\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900552445         3.482613e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900549845         3.694330e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_dispar_A            6.185859e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp003463825         1.172892e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_montpellierensis    4.341006e-06\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_seminalis       2.523844e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900550175         8.987449e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp002959755         2.335864e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900556785         4.074784e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900552715         1.032947e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_sp000431435     3.440359e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_rodentium           1.286452e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_infantium           1.252720e-04\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_criceti         3.058097e-07\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_ratti           4.754968e-07\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900538355         6.638122e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_caviae              3.746169e-06\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_magna           8.018181e-08\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_denticariosi        9.687082e-06\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900545205         3.817587e-05\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_sp900545795     1.864693e-09\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900549805         2.610571e-08\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900550455         1.139328e-06\n",
       "dtype: float64"
      ]
     },
     "execution_count": 212,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(df_taxa.mean(axis=0) / df_taxa.mean().sum())[genera]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 213,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "#OTU ID\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_atypica             0.971429\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_rogosae             0.971429\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_parvula             0.914286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_dispar              0.914286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_tobetsuensis        0.771429\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_parvula_A           0.942857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900552445         0.857143\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900549845         0.828571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_dispar_A            0.914286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp003463825         0.828571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_montpellierensis    0.742857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_seminalis       0.828571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900550175         0.771429\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp002959755         0.942857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900556785         0.914286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900552715         0.742857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_sp000431435     0.828571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_rodentium           0.628571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_infantium           0.857143\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_criceti         0.142857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_ratti           0.171429\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900538355         0.742857\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_caviae              0.428571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_magna           0.314286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_denticariosi        0.514286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900545205         0.628571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella_A;s__Veillonella_A_sp900545795     0.028571\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900549805         0.114286\n",
       "k__Bacteria;p__Firmicutes_C;c__Negativicutes;o__Veillonellales;f__Veillonellaceae;g__Veillonella;s__Veillonella_sp900550455         0.228571\n",
       "dtype: float64"
      ]
     },
     "execution_count": 213,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(((df_taxa > 0).sum(axis=0)) / df_taxa.shape[0])[genera]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 214,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.0036772251058606437\n",
      "0.0056389236160887075\n",
      "0.002186069515610472\n",
      "0.00041927854251735307\n",
      "0.0003574163793665351\n",
      "0.0040247657932257945\n",
      "0.0005146253983254719\n",
      "6.105916470448266e-05\n",
      "0.0011684249338513872\n",
      "0.00020648962167528102\n",
      "4.72115014996214e-06\n",
      "0.00029510066075474765\n",
      "0.00015088308544970556\n",
      "0.0004143575493339903\n",
      "0.0007157524851392719\n",
      "0.00015809464732849724\n",
      "3.595941007504857e-05\n",
      "2.2662311391491095e-05\n",
      "0.00021912281026804502\n",
      "4.14592475827635e-07\n",
      "7.859542155777695e-07\n",
      "0.00011046211409347521\n",
      "7.075057341688369e-06\n",
      "9.283664667603348e-08\n",
      "1.7363778223646576e-05\n",
      "6.744169158907973e-05\n",
      "1.5316721072936311e-09\n",
      "3.123461795789539e-08\n",
      "1.5660990319765224e-06\n"
     ]
    }
   ],
   "source": [
    "genus = \"g__Veillonella\"\n",
    "\n",
    "genera = [col for col in df_taxa.columns if genus in col]\n",
    "\n",
    "for species in genera:\n",
    "    print(df_taxa_ra.loc[:, species].mean())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 215,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(35, 7513)"
      ]
     },
     "execution_count": 215,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_taxa.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 216,
   "metadata": {},
   "outputs": [],
   "source": [
    "# we only want species present in greater than 95% of the samples\n",
    "mask_low_prevalence = ((df_taxa > 0).sum(axis=0)) / df_taxa.shape[0] >= .8"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 217,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "455"
      ]
     },
     "execution_count": 217,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# mask low abundance\n",
    "mask_low_abundance = (df_taxa.median(axis=0) / df_taxa.median().sum()) >= .0001\n",
    "\n",
    "mask_low_abundance.sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 218,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa = df_taxa.loc[:, mask_low_prevalence & mask_low_abundance]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 219,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_export = df_taxa.reset_index().copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 221,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Clostridium_Q;s__Clostridium_Q_saccharolyticum_A'"
      ]
     },
     "execution_count": 221,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(df_taxa_export.columns)[-1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 222,
   "metadata": {},
   "outputs": [],
   "source": [
    "sample_ids = [\"_\".join(_.split(\".\")[:2]) if not _.startswith(\"pos\") else \"_\".join(_.split(\".\")[:3]) for _ in df_taxa_export[\"index\"]]\n",
    "\n",
    "df_taxa_export[\"index\"] = sample_ids\n",
    "df_taxa_export.loc[df_taxa_export[\"index\"] == \"MERI_28\", \"index\"] = \"MEIR_28\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 223,
   "metadata": {},
   "outputs": [],
   "source": [
    "cols = list(df_taxa_export.columns)\n",
    "\n",
    "cols[0] = \"SampleID\"\n",
    "\n",
    "df_taxa_export.columns = cols"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 230,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_export.T.to_csv(\"./data/shogun.98.p8.la0001.csv\", index=True, header=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/bhillmann/.conda/envs/type_1/lib/python3.8/site-packages/skbio/sequence/_sequence.py:28: DeprecationWarning: Using or importing the ABCs from 'collections' instead of from 'collections.abc' is deprecated since Python 3.3, and in 3.9 it will stop working\n",
      "  collections.Sequence, SkbioObject):\n",
      "/home/bhillmann/.conda/envs/type_1/lib/python3.8/site-packages/skbio/io/format/gff3.py:217: DeprecationWarning: Using or importing the ABCs from 'collections' instead of from 'collections.abc' is deprecated since Python 3.3, and in 3.9 it will stop working\n",
      "  from collections import Iterable\n",
      "<frozen importlib._bootstrap>:219: RuntimeWarning: numpy.ufunc size changed, may indicate binary incompatibility. Expected 192 from C header, got 216 from PyObject\n",
      "<frozen importlib._bootstrap>:219: RuntimeWarning: numpy.ufunc size changed, may indicate binary incompatibility. Expected 192 from C header, got 216 from PyObject\n"
     ]
    }
   ],
   "source": [
    "from skbio.stats.composition import clr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_clr = df_taxa.copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_clr = pd.DataFrame(clr(df_taxa + 1))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_clr.columns = df_taxa.columns\n",
    "df_taxa_clr.index = df_taxa.index"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "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>#OTU ID</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Phocaeicola;s__Phocaeicola_vulgatus</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_fragilis</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli_D</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_flexneri</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_bouchesdurhonensis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Tannerellaceae;g__Parabacteroides;s__Parabacteroides_distasonis_A</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_uniformis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_faecis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_thetaiotaomicron</th>\n",
       "      <th>...</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Porphyromonadaceae;g__Porphyromonas;s__Porphyromonas_uenonis_A</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_sp001553265</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_oris</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp900066825</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Faecalimonas;s__Faecalimonas_sp900556835</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Hungatella_A;s__Hungatella_A_sp003478355</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp000177015</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Paraprevotella;s__Paraprevotella_sp003477995</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_oleiciplenus</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__CAG-81;s__CAG-81_sp900066785</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>pos.dna.02.S35.001.fa</th>\n",
       "      <td>12.564805</td>\n",
       "      <td>12.596202</td>\n",
       "      <td>10.432703</td>\n",
       "      <td>10.382420</td>\n",
       "      <td>10.717464</td>\n",
       "      <td>-0.529343</td>\n",
       "      <td>-0.683494</td>\n",
       "      <td>5.178716</td>\n",
       "      <td>1.714401</td>\n",
       "      <td>2.002083</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CAPMA.83.S12.001.fa</th>\n",
       "      <td>1.502473</td>\n",
       "      <td>0.108900</td>\n",
       "      <td>8.026002</td>\n",
       "      <td>6.006419</td>\n",
       "      <td>4.688253</td>\n",
       "      <td>-3.136293</td>\n",
       "      <td>0.814951</td>\n",
       "      <td>1.198817</td>\n",
       "      <td>3.197282</td>\n",
       "      <td>2.198838</td>\n",
       "      <td>...</td>\n",
       "      <td>3.138940</td>\n",
       "      <td>4.841704</td>\n",
       "      <td>3.443421</td>\n",
       "      <td>-2.730828</td>\n",
       "      <td>-2.124692</td>\n",
       "      <td>-0.929852</td>\n",
       "      <td>-4.522588</td>\n",
       "      <td>-2.443146</td>\n",
       "      <td>-3.675290</td>\n",
       "      <td>-3.829440</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>FOGI.63.S13.001.fa</th>\n",
       "      <td>3.594296</td>\n",
       "      <td>-0.611240</td>\n",
       "      <td>0.322862</td>\n",
       "      <td>-0.505445</td>\n",
       "      <td>-2.583584</td>\n",
       "      <td>-0.913521</td>\n",
       "      <td>-0.274768</td>\n",
       "      <td>2.714890</td>\n",
       "      <td>0.887994</td>\n",
       "      <td>2.404863</td>\n",
       "      <td>...</td>\n",
       "      <td>2.407628</td>\n",
       "      <td>1.353132</td>\n",
       "      <td>-0.121505</td>\n",
       "      <td>-2.117494</td>\n",
       "      <td>-3.004797</td>\n",
       "      <td>-1.820056</td>\n",
       "      <td>-3.537014</td>\n",
       "      <td>-2.626143</td>\n",
       "      <td>-1.696893</td>\n",
       "      <td>-2.913825</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>COCL.11.S7.001.fa</th>\n",
       "      <td>1.536088</td>\n",
       "      <td>1.074962</td>\n",
       "      <td>4.663100</td>\n",
       "      <td>3.188711</td>\n",
       "      <td>1.774582</td>\n",
       "      <td>-2.081677</td>\n",
       "      <td>-2.133863</td>\n",
       "      <td>1.997295</td>\n",
       "      <td>2.537797</td>\n",
       "      <td>4.726931</td>\n",
       "      <td>...</td>\n",
       "      <td>2.444156</td>\n",
       "      <td>4.645267</td>\n",
       "      <td>1.707124</td>\n",
       "      <td>-2.901118</td>\n",
       "      <td>-3.214776</td>\n",
       "      <td>-2.330573</td>\n",
       "      <td>-4.213305</td>\n",
       "      <td>-3.594266</td>\n",
       "      <td>-0.489334</td>\n",
       "      <td>-3.856630</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>WACH.38.S1.001.fa</th>\n",
       "      <td>4.240597</td>\n",
       "      <td>-0.421227</td>\n",
       "      <td>2.999156</td>\n",
       "      <td>1.341303</td>\n",
       "      <td>0.023436</td>\n",
       "      <td>-1.574201</td>\n",
       "      <td>-1.208627</td>\n",
       "      <td>2.685130</td>\n",
       "      <td>0.050759</td>\n",
       "      <td>1.814724</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.117196</td>\n",
       "      <td>3.141126</td>\n",
       "      <td>1.670049</td>\n",
       "      <td>-2.016654</td>\n",
       "      <td>-1.179972</td>\n",
       "      <td>-1.957565</td>\n",
       "      <td>-4.815676</td>\n",
       "      <td>-1.953475</td>\n",
       "      <td>-2.093441</td>\n",
       "      <td>-1.995148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PRPI.42.S31.001.fa</th>\n",
       "      <td>3.169797</td>\n",
       "      <td>-1.438640</td>\n",
       "      <td>7.172793</td>\n",
       "      <td>4.708191</td>\n",
       "      <td>3.824696</td>\n",
       "      <td>-0.532085</td>\n",
       "      <td>0.522833</td>\n",
       "      <td>2.376027</td>\n",
       "      <td>1.803052</td>\n",
       "      <td>1.172110</td>\n",
       "      <td>...</td>\n",
       "      <td>2.369289</td>\n",
       "      <td>-5.137255</td>\n",
       "      <td>-6.075525</td>\n",
       "      <td>0.473331</td>\n",
       "      <td>-1.156355</td>\n",
       "      <td>-0.917108</td>\n",
       "      <td>-1.554945</td>\n",
       "      <td>-1.227844</td>\n",
       "      <td>-0.051808</td>\n",
       "      <td>-1.882509</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>RODE.80.S8.001.fa</th>\n",
       "      <td>5.918568</td>\n",
       "      <td>-1.343871</td>\n",
       "      <td>4.211689</td>\n",
       "      <td>4.242758</td>\n",
       "      <td>4.215763</td>\n",
       "      <td>-0.235021</td>\n",
       "      <td>1.030012</td>\n",
       "      <td>2.403992</td>\n",
       "      <td>1.124451</td>\n",
       "      <td>2.924393</td>\n",
       "      <td>...</td>\n",
       "      <td>3.704593</td>\n",
       "      <td>1.527456</td>\n",
       "      <td>-1.924759</td>\n",
       "      <td>-0.191832</td>\n",
       "      <td>-1.928109</td>\n",
       "      <td>-1.326253</td>\n",
       "      <td>-1.318927</td>\n",
       "      <td>-1.467165</td>\n",
       "      <td>-2.030491</td>\n",
       "      <td>-1.091414</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>pos.dna.01.S17.001.fa</th>\n",
       "      <td>10.758325</td>\n",
       "      <td>10.466699</td>\n",
       "      <td>8.427199</td>\n",
       "      <td>8.411471</td>\n",
       "      <td>8.758813</td>\n",
       "      <td>-1.616304</td>\n",
       "      <td>-0.923157</td>\n",
       "      <td>3.486100</td>\n",
       "      <td>2.323658</td>\n",
       "      <td>1.673341</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.923157</td>\n",
       "      <td>-1.010168</td>\n",
       "      <td>-1.798626</td>\n",
       "      <td>-0.667223</td>\n",
       "      <td>-0.635475</td>\n",
       "      <td>-1.010168</td>\n",
       "      <td>-2.309451</td>\n",
       "      <td>-2.714916</td>\n",
       "      <td>-2.155300</td>\n",
       "      <td>-1.703315</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MEEN.93.S11.001.fa</th>\n",
       "      <td>4.110142</td>\n",
       "      <td>6.091029</td>\n",
       "      <td>4.004449</td>\n",
       "      <td>3.807656</td>\n",
       "      <td>2.422215</td>\n",
       "      <td>-0.461099</td>\n",
       "      <td>0.197382</td>\n",
       "      <td>3.913378</td>\n",
       "      <td>1.410490</td>\n",
       "      <td>4.790671</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.942135</td>\n",
       "      <td>-2.301677</td>\n",
       "      <td>-2.368428</td>\n",
       "      <td>0.770305</td>\n",
       "      <td>1.697052</td>\n",
       "      <td>0.457184</td>\n",
       "      <td>-0.315932</td>\n",
       "      <td>-3.596317</td>\n",
       "      <td>-1.389233</td>\n",
       "      <td>-1.218305</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DETH.41.S27.001.fa</th>\n",
       "      <td>3.183976</td>\n",
       "      <td>-3.761992</td>\n",
       "      <td>-0.928224</td>\n",
       "      <td>-2.920100</td>\n",
       "      <td>-3.462587</td>\n",
       "      <td>0.307518</td>\n",
       "      <td>-3.401551</td>\n",
       "      <td>1.633435</td>\n",
       "      <td>0.197562</td>\n",
       "      <td>-0.076657</td>\n",
       "      <td>...</td>\n",
       "      <td>1.070765</td>\n",
       "      <td>0.845063</td>\n",
       "      <td>0.288315</td>\n",
       "      <td>1.436088</td>\n",
       "      <td>-1.135138</td>\n",
       "      <td>-0.484491</td>\n",
       "      <td>-1.364955</td>\n",
       "      <td>-0.362683</td>\n",
       "      <td>-2.735072</td>\n",
       "      <td>0.163711</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TECA.45.S15.001.fa</th>\n",
       "      <td>3.134289</td>\n",
       "      <td>4.631672</td>\n",
       "      <td>2.388486</td>\n",
       "      <td>3.961840</td>\n",
       "      <td>-0.605582</td>\n",
       "      <td>-0.851349</td>\n",
       "      <td>0.515884</td>\n",
       "      <td>3.626052</td>\n",
       "      <td>-0.435956</td>\n",
       "      <td>-0.125998</td>\n",
       "      <td>...</td>\n",
       "      <td>1.255953</td>\n",
       "      <td>-2.019154</td>\n",
       "      <td>-2.957723</td>\n",
       "      <td>0.560763</td>\n",
       "      <td>-1.981214</td>\n",
       "      <td>0.301699</td>\n",
       "      <td>-1.118817</td>\n",
       "      <td>-2.342227</td>\n",
       "      <td>-1.738692</td>\n",
       "      <td>-0.965607</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BODA.52.S9.001.fa</th>\n",
       "      <td>4.091882</td>\n",
       "      <td>1.095134</td>\n",
       "      <td>0.574600</td>\n",
       "      <td>-0.936630</td>\n",
       "      <td>-1.537898</td>\n",
       "      <td>0.710902</td>\n",
       "      <td>-0.204385</td>\n",
       "      <td>2.611502</td>\n",
       "      <td>1.548782</td>\n",
       "      <td>-0.011895</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.812368</td>\n",
       "      <td>-2.907754</td>\n",
       "      <td>-2.278180</td>\n",
       "      <td>0.806878</td>\n",
       "      <td>-1.435304</td>\n",
       "      <td>0.612192</td>\n",
       "      <td>-0.467510</td>\n",
       "      <td>-0.727855</td>\n",
       "      <td>-1.865448</td>\n",
       "      <td>0.195402</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CHMA.61.S3.001.fa</th>\n",
       "      <td>6.751274</td>\n",
       "      <td>1.924646</td>\n",
       "      <td>3.074376</td>\n",
       "      <td>2.536282</td>\n",
       "      <td>1.136575</td>\n",
       "      <td>1.722341</td>\n",
       "      <td>1.558442</td>\n",
       "      <td>6.117136</td>\n",
       "      <td>3.203964</td>\n",
       "      <td>5.417908</td>\n",
       "      <td>...</td>\n",
       "      <td>-6.522518</td>\n",
       "      <td>-5.223235</td>\n",
       "      <td>-4.676691</td>\n",
       "      <td>0.755190</td>\n",
       "      <td>0.602765</td>\n",
       "      <td>0.508339</td>\n",
       "      <td>-0.216243</td>\n",
       "      <td>0.590895</td>\n",
       "      <td>0.945615</td>\n",
       "      <td>-0.637804</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ASTER.02.S34.001.fa</th>\n",
       "      <td>7.511852</td>\n",
       "      <td>-0.131282</td>\n",
       "      <td>-3.570534</td>\n",
       "      <td>-3.458830</td>\n",
       "      <td>-3.434732</td>\n",
       "      <td>0.951255</td>\n",
       "      <td>1.994439</td>\n",
       "      <td>3.549830</td>\n",
       "      <td>2.011697</td>\n",
       "      <td>4.680516</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.781043</td>\n",
       "      <td>-1.652275</td>\n",
       "      <td>-2.837292</td>\n",
       "      <td>0.171671</td>\n",
       "      <td>-0.867344</td>\n",
       "      <td>-0.210790</td>\n",
       "      <td>-1.297471</td>\n",
       "      <td>1.110924</td>\n",
       "      <td>-0.424424</td>\n",
       "      <td>0.490222</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BIRO.74.S19.001.fa</th>\n",
       "      <td>6.603156</td>\n",
       "      <td>4.111350</td>\n",
       "      <td>-1.495615</td>\n",
       "      <td>-1.615759</td>\n",
       "      <td>-2.541918</td>\n",
       "      <td>-0.057356</td>\n",
       "      <td>-0.108232</td>\n",
       "      <td>4.973692</td>\n",
       "      <td>1.036538</td>\n",
       "      <td>4.116220</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.852802</td>\n",
       "      <td>-4.475508</td>\n",
       "      <td>-3.103602</td>\n",
       "      <td>1.275340</td>\n",
       "      <td>-0.624752</td>\n",
       "      <td>-0.903965</td>\n",
       "      <td>1.473348</td>\n",
       "      <td>1.265571</td>\n",
       "      <td>1.526685</td>\n",
       "      <td>-1.350001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CAAU.07.S2.001.fa</th>\n",
       "      <td>5.567414</td>\n",
       "      <td>-1.544217</td>\n",
       "      <td>-4.652267</td>\n",
       "      <td>-5.025247</td>\n",
       "      <td>-6.750757</td>\n",
       "      <td>-0.619196</td>\n",
       "      <td>1.828631</td>\n",
       "      <td>5.706051</td>\n",
       "      <td>1.065846</td>\n",
       "      <td>3.684255</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.491175</td>\n",
       "      <td>-0.625905</td>\n",
       "      <td>-4.495425</td>\n",
       "      <td>0.780260</td>\n",
       "      <td>0.033700</td>\n",
       "      <td>-1.256430</td>\n",
       "      <td>-2.549514</td>\n",
       "      <td>0.328233</td>\n",
       "      <td>0.564564</td>\n",
       "      <td>-0.054863</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DIGMA.71.S21.001.fa</th>\n",
       "      <td>5.374936</td>\n",
       "      <td>0.454759</td>\n",
       "      <td>-1.006752</td>\n",
       "      <td>0.006649</td>\n",
       "      <td>-1.588852</td>\n",
       "      <td>-0.308543</td>\n",
       "      <td>-0.375151</td>\n",
       "      <td>6.540831</td>\n",
       "      <td>0.876363</td>\n",
       "      <td>1.213984</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.997182</td>\n",
       "      <td>-5.422697</td>\n",
       "      <td>-2.208633</td>\n",
       "      <td>4.121940</td>\n",
       "      <td>-0.361749</td>\n",
       "      <td>1.308846</td>\n",
       "      <td>-0.674639</td>\n",
       "      <td>0.604328</td>\n",
       "      <td>-0.262493</td>\n",
       "      <td>0.456438</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>FAFE.16.S10.001.fa</th>\n",
       "      <td>5.584302</td>\n",
       "      <td>-0.067135</td>\n",
       "      <td>-1.416951</td>\n",
       "      <td>-2.181225</td>\n",
       "      <td>-2.503694</td>\n",
       "      <td>1.577985</td>\n",
       "      <td>0.749284</td>\n",
       "      <td>5.106302</td>\n",
       "      <td>6.007668</td>\n",
       "      <td>3.278542</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.479546</td>\n",
       "      <td>-5.285171</td>\n",
       "      <td>-5.385255</td>\n",
       "      <td>1.517647</td>\n",
       "      <td>-0.069201</td>\n",
       "      <td>1.602478</td>\n",
       "      <td>0.442297</td>\n",
       "      <td>-0.797606</td>\n",
       "      <td>-1.121834</td>\n",
       "      <td>0.994868</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>OPSE.53.S5.001.fa</th>\n",
       "      <td>7.120223</td>\n",
       "      <td>1.635535</td>\n",
       "      <td>1.917310</td>\n",
       "      <td>-0.146372</td>\n",
       "      <td>-1.667841</td>\n",
       "      <td>1.399707</td>\n",
       "      <td>3.893911</td>\n",
       "      <td>4.110070</td>\n",
       "      <td>1.131665</td>\n",
       "      <td>4.004854</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.781932</td>\n",
       "      <td>-4.310000</td>\n",
       "      <td>-4.449762</td>\n",
       "      <td>-0.532751</td>\n",
       "      <td>-0.095620</td>\n",
       "      <td>-1.119344</td>\n",
       "      <td>-2.203747</td>\n",
       "      <td>0.762726</td>\n",
       "      <td>-0.606017</td>\n",
       "      <td>-1.502694</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>LAES.91.S29.001.fa</th>\n",
       "      <td>3.183280</td>\n",
       "      <td>-0.423308</td>\n",
       "      <td>5.162297</td>\n",
       "      <td>3.618158</td>\n",
       "      <td>2.340659</td>\n",
       "      <td>-0.533104</td>\n",
       "      <td>-2.047337</td>\n",
       "      <td>4.861458</td>\n",
       "      <td>-0.527922</td>\n",
       "      <td>1.678388</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.659404</td>\n",
       "      <td>-3.815749</td>\n",
       "      <td>-2.204774</td>\n",
       "      <td>1.430708</td>\n",
       "      <td>-2.543399</td>\n",
       "      <td>0.343855</td>\n",
       "      <td>0.111238</td>\n",
       "      <td>-2.009986</td>\n",
       "      <td>0.495566</td>\n",
       "      <td>-1.028367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DOEC.81.S30.001.fa</th>\n",
       "      <td>4.092462</td>\n",
       "      <td>-0.679445</td>\n",
       "      <td>3.202142</td>\n",
       "      <td>3.084641</td>\n",
       "      <td>3.229331</td>\n",
       "      <td>-0.325333</td>\n",
       "      <td>1.312362</td>\n",
       "      <td>3.325615</td>\n",
       "      <td>4.108773</td>\n",
       "      <td>1.343689</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.364446</td>\n",
       "      <td>-2.452479</td>\n",
       "      <td>-3.483327</td>\n",
       "      <td>0.449270</td>\n",
       "      <td>-1.440480</td>\n",
       "      <td>-0.737190</td>\n",
       "      <td>-1.007263</td>\n",
       "      <td>-2.681997</td>\n",
       "      <td>-0.129765</td>\n",
       "      <td>-0.393402</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EKJO.64.S32.001.fa</th>\n",
       "      <td>2.758383</td>\n",
       "      <td>1.899094</td>\n",
       "      <td>3.268204</td>\n",
       "      <td>1.459727</td>\n",
       "      <td>1.149572</td>\n",
       "      <td>2.914325</td>\n",
       "      <td>-1.524576</td>\n",
       "      <td>1.952459</td>\n",
       "      <td>1.117723</td>\n",
       "      <td>-0.472484</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.004701</td>\n",
       "      <td>2.080180</td>\n",
       "      <td>5.461192</td>\n",
       "      <td>-0.952057</td>\n",
       "      <td>-2.130712</td>\n",
       "      <td>-0.723798</td>\n",
       "      <td>-2.536177</td>\n",
       "      <td>-3.517006</td>\n",
       "      <td>-1.907569</td>\n",
       "      <td>-1.843030</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DIAL.15.S23.001.fa</th>\n",
       "      <td>4.740969</td>\n",
       "      <td>4.392127</td>\n",
       "      <td>1.668718</td>\n",
       "      <td>0.472244</td>\n",
       "      <td>-1.063350</td>\n",
       "      <td>-0.037714</td>\n",
       "      <td>1.751842</td>\n",
       "      <td>3.866595</td>\n",
       "      <td>-0.562924</td>\n",
       "      <td>2.009991</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.826564</td>\n",
       "      <td>-1.847991</td>\n",
       "      <td>-2.486585</td>\n",
       "      <td>0.842740</td>\n",
       "      <td>-1.756216</td>\n",
       "      <td>1.252917</td>\n",
       "      <td>-0.041654</td>\n",
       "      <td>-1.053558</td>\n",
       "      <td>0.031977</td>\n",
       "      <td>-0.745484</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>GIFI.21.S22.001.fa</th>\n",
       "      <td>3.184368</td>\n",
       "      <td>2.202335</td>\n",
       "      <td>-0.188357</td>\n",
       "      <td>-0.344703</td>\n",
       "      <td>-0.019281</td>\n",
       "      <td>2.197504</td>\n",
       "      <td>-1.151795</td>\n",
       "      <td>2.699819</td>\n",
       "      <td>1.567684</td>\n",
       "      <td>0.284930</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.117893</td>\n",
       "      <td>3.923379</td>\n",
       "      <td>5.297858</td>\n",
       "      <td>0.153738</td>\n",
       "      <td>-2.727331</td>\n",
       "      <td>-1.341037</td>\n",
       "      <td>-3.132796</td>\n",
       "      <td>-2.909653</td>\n",
       "      <td>-2.573181</td>\n",
       "      <td>-3.132796</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TAOT.65.S24.001.fa</th>\n",
       "      <td>2.689188</td>\n",
       "      <td>1.152850</td>\n",
       "      <td>0.797345</td>\n",
       "      <td>-1.012114</td>\n",
       "      <td>-2.398408</td>\n",
       "      <td>-0.606649</td>\n",
       "      <td>-1.571730</td>\n",
       "      <td>3.131588</td>\n",
       "      <td>-0.432295</td>\n",
       "      <td>-0.878582</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.552559</td>\n",
       "      <td>3.106923</td>\n",
       "      <td>5.611904</td>\n",
       "      <td>0.568337</td>\n",
       "      <td>-2.958024</td>\n",
       "      <td>-1.399879</td>\n",
       "      <td>-2.734880</td>\n",
       "      <td>-2.264877</td>\n",
       "      <td>-2.552559</td>\n",
       "      <td>-1.208824</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DEFI.14.S6.001.fa</th>\n",
       "      <td>3.475730</td>\n",
       "      <td>5.737301</td>\n",
       "      <td>3.139103</td>\n",
       "      <td>5.188512</td>\n",
       "      <td>0.212743</td>\n",
       "      <td>-1.604646</td>\n",
       "      <td>0.610018</td>\n",
       "      <td>3.058043</td>\n",
       "      <td>-0.740189</td>\n",
       "      <td>2.086746</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.038194</td>\n",
       "      <td>1.291514</td>\n",
       "      <td>-0.052186</td>\n",
       "      <td>-0.119842</td>\n",
       "      <td>1.011532</td>\n",
       "      <td>-1.295560</td>\n",
       "      <td>-2.430593</td>\n",
       "      <td>-3.090619</td>\n",
       "      <td>-1.534714</td>\n",
       "      <td>-1.908701</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MUNA.98.S28.001.fa</th>\n",
       "      <td>2.218994</td>\n",
       "      <td>-1.491165</td>\n",
       "      <td>-1.683421</td>\n",
       "      <td>-2.960675</td>\n",
       "      <td>-3.069716</td>\n",
       "      <td>-2.292634</td>\n",
       "      <td>-0.164856</td>\n",
       "      <td>4.237870</td>\n",
       "      <td>-0.495943</td>\n",
       "      <td>1.740706</td>\n",
       "      <td>...</td>\n",
       "      <td>0.302494</td>\n",
       "      <td>0.470776</td>\n",
       "      <td>-0.653966</td>\n",
       "      <td>-0.514026</td>\n",
       "      <td>-0.816727</td>\n",
       "      <td>-0.916989</td>\n",
       "      <td>-1.873143</td>\n",
       "      <td>-4.406562</td>\n",
       "      <td>-1.715021</td>\n",
       "      <td>-0.154564</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ROAL.75.S33.001.fa</th>\n",
       "      <td>2.617716</td>\n",
       "      <td>0.620977</td>\n",
       "      <td>1.251969</td>\n",
       "      <td>0.998512</td>\n",
       "      <td>1.236305</td>\n",
       "      <td>-0.857894</td>\n",
       "      <td>-0.566466</td>\n",
       "      <td>2.787953</td>\n",
       "      <td>-0.119041</td>\n",
       "      <td>0.706430</td>\n",
       "      <td>...</td>\n",
       "      <td>2.554822</td>\n",
       "      <td>1.559750</td>\n",
       "      <td>0.630927</td>\n",
       "      <td>0.524060</td>\n",
       "      <td>-2.359089</td>\n",
       "      <td>-2.246501</td>\n",
       "      <td>-3.651487</td>\n",
       "      <td>-2.515834</td>\n",
       "      <td>-1.401189</td>\n",
       "      <td>-2.638173</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MERI.28.S20.001.fa</th>\n",
       "      <td>1.645391</td>\n",
       "      <td>-2.516686</td>\n",
       "      <td>6.880540</td>\n",
       "      <td>5.298956</td>\n",
       "      <td>5.544301</td>\n",
       "      <td>-2.042753</td>\n",
       "      <td>0.192413</td>\n",
       "      <td>2.897733</td>\n",
       "      <td>0.579711</td>\n",
       "      <td>-1.315868</td>\n",
       "      <td>...</td>\n",
       "      <td>0.439290</td>\n",
       "      <td>2.206092</td>\n",
       "      <td>-0.342471</td>\n",
       "      <td>-0.432335</td>\n",
       "      <td>-3.599673</td>\n",
       "      <td>-0.984146</td>\n",
       "      <td>-0.108111</td>\n",
       "      <td>-3.217738</td>\n",
       "      <td>-3.448850</td>\n",
       "      <td>-2.688044</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BUMA.05.S4.001.fa</th>\n",
       "      <td>5.184688</td>\n",
       "      <td>0.726147</td>\n",
       "      <td>-1.573127</td>\n",
       "      <td>-1.091536</td>\n",
       "      <td>-2.000538</td>\n",
       "      <td>-1.530629</td>\n",
       "      <td>0.112467</td>\n",
       "      <td>4.452070</td>\n",
       "      <td>-0.874077</td>\n",
       "      <td>1.570038</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.666852</td>\n",
       "      <td>-1.504957</td>\n",
       "      <td>-2.146903</td>\n",
       "      <td>-0.352957</td>\n",
       "      <td>-1.530629</td>\n",
       "      <td>-0.746264</td>\n",
       "      <td>-1.368205</td>\n",
       "      <td>-2.212961</td>\n",
       "      <td>-1.443736</td>\n",
       "      <td>-1.600636</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MOEM.48.S16.001.fa</th>\n",
       "      <td>5.631256</td>\n",
       "      <td>1.947917</td>\n",
       "      <td>-1.263681</td>\n",
       "      <td>-1.753768</td>\n",
       "      <td>-2.266186</td>\n",
       "      <td>0.196565</td>\n",
       "      <td>2.102101</td>\n",
       "      <td>4.604422</td>\n",
       "      <td>2.247472</td>\n",
       "      <td>3.520333</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.459434</td>\n",
       "      <td>-5.226689</td>\n",
       "      <td>-7.098492</td>\n",
       "      <td>0.976780</td>\n",
       "      <td>2.587275</td>\n",
       "      <td>1.486921</td>\n",
       "      <td>-1.131064</td>\n",
       "      <td>-0.922624</td>\n",
       "      <td>0.683690</td>\n",
       "      <td>0.319989</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PAMA.46.S18.001.fa</th>\n",
       "      <td>2.826640</td>\n",
       "      <td>1.425925</td>\n",
       "      <td>0.920859</td>\n",
       "      <td>-0.381243</td>\n",
       "      <td>-0.867528</td>\n",
       "      <td>-1.068199</td>\n",
       "      <td>-2.690059</td>\n",
       "      <td>2.445372</td>\n",
       "      <td>-0.587451</td>\n",
       "      <td>-0.174381</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.335513</td>\n",
       "      <td>1.219461</td>\n",
       "      <td>-0.602835</td>\n",
       "      <td>-0.996740</td>\n",
       "      <td>-0.908771</td>\n",
       "      <td>-0.037698</td>\n",
       "      <td>-2.690059</td>\n",
       "      <td>-1.335513</td>\n",
       "      <td>-2.977741</td>\n",
       "      <td>-2.327154</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>VUGE.37.S14.001.fa</th>\n",
       "      <td>4.700142</td>\n",
       "      <td>-1.786927</td>\n",
       "      <td>4.158713</td>\n",
       "      <td>1.681301</td>\n",
       "      <td>0.816839</td>\n",
       "      <td>-1.146786</td>\n",
       "      <td>0.389677</td>\n",
       "      <td>2.790468</td>\n",
       "      <td>2.500457</td>\n",
       "      <td>0.155391</td>\n",
       "      <td>...</td>\n",
       "      <td>0.423770</td>\n",
       "      <td>-1.154333</td>\n",
       "      <td>-3.232563</td>\n",
       "      <td>0.427640</td>\n",
       "      <td>-0.313027</td>\n",
       "      <td>0.374880</td>\n",
       "      <td>-0.862795</td>\n",
       "      <td>-2.197253</td>\n",
       "      <td>-1.232786</td>\n",
       "      <td>-0.935806</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TOSA.76.S25.001.fa</th>\n",
       "      <td>2.213781</td>\n",
       "      <td>6.897451</td>\n",
       "      <td>-1.166429</td>\n",
       "      <td>-0.946432</td>\n",
       "      <td>-1.813533</td>\n",
       "      <td>-1.793730</td>\n",
       "      <td>1.110166</td>\n",
       "      <td>2.397261</td>\n",
       "      <td>2.802775</td>\n",
       "      <td>0.494040</td>\n",
       "      <td>...</td>\n",
       "      <td>0.618325</td>\n",
       "      <td>1.556518</td>\n",
       "      <td>1.321092</td>\n",
       "      <td>0.577063</td>\n",
       "      <td>-2.275568</td>\n",
       "      <td>-0.003606</td>\n",
       "      <td>-2.781117</td>\n",
       "      <td>-3.086498</td>\n",
       "      <td>-2.393351</td>\n",
       "      <td>-2.393351</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DEIS.70.S26.001.fa</th>\n",
       "      <td>2.781823</td>\n",
       "      <td>1.181892</td>\n",
       "      <td>1.358823</td>\n",
       "      <td>0.169754</td>\n",
       "      <td>-0.047490</td>\n",
       "      <td>4.855276</td>\n",
       "      <td>-1.869103</td>\n",
       "      <td>2.909639</td>\n",
       "      <td>-1.153483</td>\n",
       "      <td>-1.433785</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.406268</td>\n",
       "      <td>3.801718</td>\n",
       "      <td>10.340782</td>\n",
       "      <td>1.103640</td>\n",
       "      <td>-3.573851</td>\n",
       "      <td>0.557308</td>\n",
       "      <td>-3.350707</td>\n",
       "      <td>-1.246573</td>\n",
       "      <td>0.276297</td>\n",
       "      <td>-3.350707</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>35 rows × 181 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Phocaeicola;s__Phocaeicola_vulgatus  \\\n",
       "pos.dna.02.S35.001.fa                                          12.564805                                                                      \n",
       "CAPMA.83.S12.001.fa                                             1.502473                                                                      \n",
       "FOGI.63.S13.001.fa                                              3.594296                                                                      \n",
       "COCL.11.S7.001.fa                                               1.536088                                                                      \n",
       "WACH.38.S1.001.fa                                               4.240597                                                                      \n",
       "PRPI.42.S31.001.fa                                              3.169797                                                                      \n",
       "RODE.80.S8.001.fa                                               5.918568                                                                      \n",
       "pos.dna.01.S17.001.fa                                          10.758325                                                                      \n",
       "MEEN.93.S11.001.fa                                              4.110142                                                                      \n",
       "DETH.41.S27.001.fa                                              3.183976                                                                      \n",
       "TECA.45.S15.001.fa                                              3.134289                                                                      \n",
       "BODA.52.S9.001.fa                                               4.091882                                                                      \n",
       "CHMA.61.S3.001.fa                                               6.751274                                                                      \n",
       "ASTER.02.S34.001.fa                                             7.511852                                                                      \n",
       "BIRO.74.S19.001.fa                                              6.603156                                                                      \n",
       "CAAU.07.S2.001.fa                                               5.567414                                                                      \n",
       "DIGMA.71.S21.001.fa                                             5.374936                                                                      \n",
       "FAFE.16.S10.001.fa                                              5.584302                                                                      \n",
       "OPSE.53.S5.001.fa                                               7.120223                                                                      \n",
       "LAES.91.S29.001.fa                                              3.183280                                                                      \n",
       "DOEC.81.S30.001.fa                                              4.092462                                                                      \n",
       "EKJO.64.S32.001.fa                                              2.758383                                                                      \n",
       "DIAL.15.S23.001.fa                                              4.740969                                                                      \n",
       "GIFI.21.S22.001.fa                                              3.184368                                                                      \n",
       "TAOT.65.S24.001.fa                                              2.689188                                                                      \n",
       "DEFI.14.S6.001.fa                                               3.475730                                                                      \n",
       "MUNA.98.S28.001.fa                                              2.218994                                                                      \n",
       "ROAL.75.S33.001.fa                                              2.617716                                                                      \n",
       "MERI.28.S20.001.fa                                              1.645391                                                                      \n",
       "BUMA.05.S4.001.fa                                               5.184688                                                                      \n",
       "MOEM.48.S16.001.fa                                              5.631256                                                                      \n",
       "PAMA.46.S18.001.fa                                              2.826640                                                                      \n",
       "VUGE.37.S14.001.fa                                              4.700142                                                                      \n",
       "TOSA.76.S25.001.fa                                              2.213781                                                                      \n",
       "DEIS.70.S26.001.fa                                              2.781823                                                                      \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_fragilis  \\\n",
       "pos.dna.02.S35.001.fa                                          12.596202                                                                      \n",
       "CAPMA.83.S12.001.fa                                             0.108900                                                                      \n",
       "FOGI.63.S13.001.fa                                             -0.611240                                                                      \n",
       "COCL.11.S7.001.fa                                               1.074962                                                                      \n",
       "WACH.38.S1.001.fa                                              -0.421227                                                                      \n",
       "PRPI.42.S31.001.fa                                             -1.438640                                                                      \n",
       "RODE.80.S8.001.fa                                              -1.343871                                                                      \n",
       "pos.dna.01.S17.001.fa                                          10.466699                                                                      \n",
       "MEEN.93.S11.001.fa                                              6.091029                                                                      \n",
       "DETH.41.S27.001.fa                                             -3.761992                                                                      \n",
       "TECA.45.S15.001.fa                                              4.631672                                                                      \n",
       "BODA.52.S9.001.fa                                               1.095134                                                                      \n",
       "CHMA.61.S3.001.fa                                               1.924646                                                                      \n",
       "ASTER.02.S34.001.fa                                            -0.131282                                                                      \n",
       "BIRO.74.S19.001.fa                                              4.111350                                                                      \n",
       "CAAU.07.S2.001.fa                                              -1.544217                                                                      \n",
       "DIGMA.71.S21.001.fa                                             0.454759                                                                      \n",
       "FAFE.16.S10.001.fa                                             -0.067135                                                                      \n",
       "OPSE.53.S5.001.fa                                               1.635535                                                                      \n",
       "LAES.91.S29.001.fa                                             -0.423308                                                                      \n",
       "DOEC.81.S30.001.fa                                             -0.679445                                                                      \n",
       "EKJO.64.S32.001.fa                                              1.899094                                                                      \n",
       "DIAL.15.S23.001.fa                                              4.392127                                                                      \n",
       "GIFI.21.S22.001.fa                                              2.202335                                                                      \n",
       "TAOT.65.S24.001.fa                                              1.152850                                                                      \n",
       "DEFI.14.S6.001.fa                                               5.737301                                                                      \n",
       "MUNA.98.S28.001.fa                                             -1.491165                                                                      \n",
       "ROAL.75.S33.001.fa                                              0.620977                                                                      \n",
       "MERI.28.S20.001.fa                                             -2.516686                                                                      \n",
       "BUMA.05.S4.001.fa                                               0.726147                                                                      \n",
       "MOEM.48.S16.001.fa                                              1.947917                                                                      \n",
       "PAMA.46.S18.001.fa                                              1.425925                                                                      \n",
       "VUGE.37.S14.001.fa                                             -1.786927                                                                      \n",
       "TOSA.76.S25.001.fa                                              6.897451                                                                      \n",
       "DEIS.70.S26.001.fa                                              1.181892                                                                      \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli  \\\n",
       "pos.dna.02.S35.001.fa                                          10.432703                                                                                   \n",
       "CAPMA.83.S12.001.fa                                             8.026002                                                                                   \n",
       "FOGI.63.S13.001.fa                                              0.322862                                                                                   \n",
       "COCL.11.S7.001.fa                                               4.663100                                                                                   \n",
       "WACH.38.S1.001.fa                                               2.999156                                                                                   \n",
       "PRPI.42.S31.001.fa                                              7.172793                                                                                   \n",
       "RODE.80.S8.001.fa                                               4.211689                                                                                   \n",
       "pos.dna.01.S17.001.fa                                           8.427199                                                                                   \n",
       "MEEN.93.S11.001.fa                                              4.004449                                                                                   \n",
       "DETH.41.S27.001.fa                                             -0.928224                                                                                   \n",
       "TECA.45.S15.001.fa                                              2.388486                                                                                   \n",
       "BODA.52.S9.001.fa                                               0.574600                                                                                   \n",
       "CHMA.61.S3.001.fa                                               3.074376                                                                                   \n",
       "ASTER.02.S34.001.fa                                            -3.570534                                                                                   \n",
       "BIRO.74.S19.001.fa                                             -1.495615                                                                                   \n",
       "CAAU.07.S2.001.fa                                              -4.652267                                                                                   \n",
       "DIGMA.71.S21.001.fa                                            -1.006752                                                                                   \n",
       "FAFE.16.S10.001.fa                                             -1.416951                                                                                   \n",
       "OPSE.53.S5.001.fa                                               1.917310                                                                                   \n",
       "LAES.91.S29.001.fa                                              5.162297                                                                                   \n",
       "DOEC.81.S30.001.fa                                              3.202142                                                                                   \n",
       "EKJO.64.S32.001.fa                                              3.268204                                                                                   \n",
       "DIAL.15.S23.001.fa                                              1.668718                                                                                   \n",
       "GIFI.21.S22.001.fa                                             -0.188357                                                                                   \n",
       "TAOT.65.S24.001.fa                                              0.797345                                                                                   \n",
       "DEFI.14.S6.001.fa                                               3.139103                                                                                   \n",
       "MUNA.98.S28.001.fa                                             -1.683421                                                                                   \n",
       "ROAL.75.S33.001.fa                                              1.251969                                                                                   \n",
       "MERI.28.S20.001.fa                                              6.880540                                                                                   \n",
       "BUMA.05.S4.001.fa                                              -1.573127                                                                                   \n",
       "MOEM.48.S16.001.fa                                             -1.263681                                                                                   \n",
       "PAMA.46.S18.001.fa                                              0.920859                                                                                   \n",
       "VUGE.37.S14.001.fa                                              4.158713                                                                                   \n",
       "TOSA.76.S25.001.fa                                             -1.166429                                                                                   \n",
       "DEIS.70.S26.001.fa                                              1.358823                                                                                   \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli_D  \\\n",
       "pos.dna.02.S35.001.fa                                          10.382420                                                                                     \n",
       "CAPMA.83.S12.001.fa                                             6.006419                                                                                     \n",
       "FOGI.63.S13.001.fa                                             -0.505445                                                                                     \n",
       "COCL.11.S7.001.fa                                               3.188711                                                                                     \n",
       "WACH.38.S1.001.fa                                               1.341303                                                                                     \n",
       "PRPI.42.S31.001.fa                                              4.708191                                                                                     \n",
       "RODE.80.S8.001.fa                                               4.242758                                                                                     \n",
       "pos.dna.01.S17.001.fa                                           8.411471                                                                                     \n",
       "MEEN.93.S11.001.fa                                              3.807656                                                                                     \n",
       "DETH.41.S27.001.fa                                             -2.920100                                                                                     \n",
       "TECA.45.S15.001.fa                                              3.961840                                                                                     \n",
       "BODA.52.S9.001.fa                                              -0.936630                                                                                     \n",
       "CHMA.61.S3.001.fa                                               2.536282                                                                                     \n",
       "ASTER.02.S34.001.fa                                            -3.458830                                                                                     \n",
       "BIRO.74.S19.001.fa                                             -1.615759                                                                                     \n",
       "CAAU.07.S2.001.fa                                              -5.025247                                                                                     \n",
       "DIGMA.71.S21.001.fa                                             0.006649                                                                                     \n",
       "FAFE.16.S10.001.fa                                             -2.181225                                                                                     \n",
       "OPSE.53.S5.001.fa                                              -0.146372                                                                                     \n",
       "LAES.91.S29.001.fa                                              3.618158                                                                                     \n",
       "DOEC.81.S30.001.fa                                              3.084641                                                                                     \n",
       "EKJO.64.S32.001.fa                                              1.459727                                                                                     \n",
       "DIAL.15.S23.001.fa                                              0.472244                                                                                     \n",
       "GIFI.21.S22.001.fa                                             -0.344703                                                                                     \n",
       "TAOT.65.S24.001.fa                                             -1.012114                                                                                     \n",
       "DEFI.14.S6.001.fa                                               5.188512                                                                                     \n",
       "MUNA.98.S28.001.fa                                             -2.960675                                                                                     \n",
       "ROAL.75.S33.001.fa                                              0.998512                                                                                     \n",
       "MERI.28.S20.001.fa                                              5.298956                                                                                     \n",
       "BUMA.05.S4.001.fa                                              -1.091536                                                                                     \n",
       "MOEM.48.S16.001.fa                                             -1.753768                                                                                     \n",
       "PAMA.46.S18.001.fa                                             -0.381243                                                                                     \n",
       "VUGE.37.S14.001.fa                                              1.681301                                                                                     \n",
       "TOSA.76.S25.001.fa                                             -0.946432                                                                                     \n",
       "DEIS.70.S26.001.fa                                              0.169754                                                                                     \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_flexneri  \\\n",
       "pos.dna.02.S35.001.fa                                          10.717464                                                                                       \n",
       "CAPMA.83.S12.001.fa                                             4.688253                                                                                       \n",
       "FOGI.63.S13.001.fa                                             -2.583584                                                                                       \n",
       "COCL.11.S7.001.fa                                               1.774582                                                                                       \n",
       "WACH.38.S1.001.fa                                               0.023436                                                                                       \n",
       "PRPI.42.S31.001.fa                                              3.824696                                                                                       \n",
       "RODE.80.S8.001.fa                                               4.215763                                                                                       \n",
       "pos.dna.01.S17.001.fa                                           8.758813                                                                                       \n",
       "MEEN.93.S11.001.fa                                              2.422215                                                                                       \n",
       "DETH.41.S27.001.fa                                             -3.462587                                                                                       \n",
       "TECA.45.S15.001.fa                                             -0.605582                                                                                       \n",
       "BODA.52.S9.001.fa                                              -1.537898                                                                                       \n",
       "CHMA.61.S3.001.fa                                               1.136575                                                                                       \n",
       "ASTER.02.S34.001.fa                                            -3.434732                                                                                       \n",
       "BIRO.74.S19.001.fa                                             -2.541918                                                                                       \n",
       "CAAU.07.S2.001.fa                                              -6.750757                                                                                       \n",
       "DIGMA.71.S21.001.fa                                            -1.588852                                                                                       \n",
       "FAFE.16.S10.001.fa                                             -2.503694                                                                                       \n",
       "OPSE.53.S5.001.fa                                              -1.667841                                                                                       \n",
       "LAES.91.S29.001.fa                                              2.340659                                                                                       \n",
       "DOEC.81.S30.001.fa                                              3.229331                                                                                       \n",
       "EKJO.64.S32.001.fa                                              1.149572                                                                                       \n",
       "DIAL.15.S23.001.fa                                             -1.063350                                                                                       \n",
       "GIFI.21.S22.001.fa                                             -0.019281                                                                                       \n",
       "TAOT.65.S24.001.fa                                             -2.398408                                                                                       \n",
       "DEFI.14.S6.001.fa                                               0.212743                                                                                       \n",
       "MUNA.98.S28.001.fa                                             -3.069716                                                                                       \n",
       "ROAL.75.S33.001.fa                                              1.236305                                                                                       \n",
       "MERI.28.S20.001.fa                                              5.544301                                                                                       \n",
       "BUMA.05.S4.001.fa                                              -2.000538                                                                                       \n",
       "MOEM.48.S16.001.fa                                             -2.266186                                                                                       \n",
       "PAMA.46.S18.001.fa                                             -0.867528                                                                                       \n",
       "VUGE.37.S14.001.fa                                              0.816839                                                                                       \n",
       "TOSA.76.S25.001.fa                                             -1.813533                                                                                       \n",
       "DEIS.70.S26.001.fa                                             -0.047490                                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_bouchesdurhonensis  \\\n",
       "pos.dna.02.S35.001.fa                                          -0.529343                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -3.136293                                                                                \n",
       "FOGI.63.S13.001.fa                                             -0.913521                                                                                \n",
       "COCL.11.S7.001.fa                                              -2.081677                                                                                \n",
       "WACH.38.S1.001.fa                                              -1.574201                                                                                \n",
       "PRPI.42.S31.001.fa                                             -0.532085                                                                                \n",
       "RODE.80.S8.001.fa                                              -0.235021                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -1.616304                                                                                \n",
       "MEEN.93.S11.001.fa                                             -0.461099                                                                                \n",
       "DETH.41.S27.001.fa                                              0.307518                                                                                \n",
       "TECA.45.S15.001.fa                                             -0.851349                                                                                \n",
       "BODA.52.S9.001.fa                                               0.710902                                                                                \n",
       "CHMA.61.S3.001.fa                                               1.722341                                                                                \n",
       "ASTER.02.S34.001.fa                                             0.951255                                                                                \n",
       "BIRO.74.S19.001.fa                                             -0.057356                                                                                \n",
       "CAAU.07.S2.001.fa                                              -0.619196                                                                                \n",
       "DIGMA.71.S21.001.fa                                            -0.308543                                                                                \n",
       "FAFE.16.S10.001.fa                                              1.577985                                                                                \n",
       "OPSE.53.S5.001.fa                                               1.399707                                                                                \n",
       "LAES.91.S29.001.fa                                             -0.533104                                                                                \n",
       "DOEC.81.S30.001.fa                                             -0.325333                                                                                \n",
       "EKJO.64.S32.001.fa                                              2.914325                                                                                \n",
       "DIAL.15.S23.001.fa                                             -0.037714                                                                                \n",
       "GIFI.21.S22.001.fa                                              2.197504                                                                                \n",
       "TAOT.65.S24.001.fa                                             -0.606649                                                                                \n",
       "DEFI.14.S6.001.fa                                              -1.604646                                                                                \n",
       "MUNA.98.S28.001.fa                                             -2.292634                                                                                \n",
       "ROAL.75.S33.001.fa                                             -0.857894                                                                                \n",
       "MERI.28.S20.001.fa                                             -2.042753                                                                                \n",
       "BUMA.05.S4.001.fa                                              -1.530629                                                                                \n",
       "MOEM.48.S16.001.fa                                              0.196565                                                                                \n",
       "PAMA.46.S18.001.fa                                             -1.068199                                                                                \n",
       "VUGE.37.S14.001.fa                                             -1.146786                                                                                \n",
       "TOSA.76.S25.001.fa                                             -1.793730                                                                                \n",
       "DEIS.70.S26.001.fa                                              4.855276                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Tannerellaceae;g__Parabacteroides;s__Parabacteroides_distasonis_A  \\\n",
       "pos.dna.02.S35.001.fa                                          -0.683494                                                                                  \n",
       "CAPMA.83.S12.001.fa                                             0.814951                                                                                  \n",
       "FOGI.63.S13.001.fa                                             -0.274768                                                                                  \n",
       "COCL.11.S7.001.fa                                              -2.133863                                                                                  \n",
       "WACH.38.S1.001.fa                                              -1.208627                                                                                  \n",
       "PRPI.42.S31.001.fa                                              0.522833                                                                                  \n",
       "RODE.80.S8.001.fa                                               1.030012                                                                                  \n",
       "pos.dna.01.S17.001.fa                                          -0.923157                                                                                  \n",
       "MEEN.93.S11.001.fa                                              0.197382                                                                                  \n",
       "DETH.41.S27.001.fa                                             -3.401551                                                                                  \n",
       "TECA.45.S15.001.fa                                              0.515884                                                                                  \n",
       "BODA.52.S9.001.fa                                              -0.204385                                                                                  \n",
       "CHMA.61.S3.001.fa                                               1.558442                                                                                  \n",
       "ASTER.02.S34.001.fa                                             1.994439                                                                                  \n",
       "BIRO.74.S19.001.fa                                             -0.108232                                                                                  \n",
       "CAAU.07.S2.001.fa                                               1.828631                                                                                  \n",
       "DIGMA.71.S21.001.fa                                            -0.375151                                                                                  \n",
       "FAFE.16.S10.001.fa                                              0.749284                                                                                  \n",
       "OPSE.53.S5.001.fa                                               3.893911                                                                                  \n",
       "LAES.91.S29.001.fa                                             -2.047337                                                                                  \n",
       "DOEC.81.S30.001.fa                                              1.312362                                                                                  \n",
       "EKJO.64.S32.001.fa                                             -1.524576                                                                                  \n",
       "DIAL.15.S23.001.fa                                              1.751842                                                                                  \n",
       "GIFI.21.S22.001.fa                                             -1.151795                                                                                  \n",
       "TAOT.65.S24.001.fa                                             -1.571730                                                                                  \n",
       "DEFI.14.S6.001.fa                                               0.610018                                                                                  \n",
       "MUNA.98.S28.001.fa                                             -0.164856                                                                                  \n",
       "ROAL.75.S33.001.fa                                             -0.566466                                                                                  \n",
       "MERI.28.S20.001.fa                                              0.192413                                                                                  \n",
       "BUMA.05.S4.001.fa                                               0.112467                                                                                  \n",
       "MOEM.48.S16.001.fa                                              2.102101                                                                                  \n",
       "PAMA.46.S18.001.fa                                             -2.690059                                                                                  \n",
       "VUGE.37.S14.001.fa                                              0.389677                                                                                  \n",
       "TOSA.76.S25.001.fa                                              1.110166                                                                                  \n",
       "DEIS.70.S26.001.fa                                             -1.869103                                                                                  \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_uniformis  \\\n",
       "pos.dna.02.S35.001.fa                                           5.178716                                                                       \n",
       "CAPMA.83.S12.001.fa                                             1.198817                                                                       \n",
       "FOGI.63.S13.001.fa                                              2.714890                                                                       \n",
       "COCL.11.S7.001.fa                                               1.997295                                                                       \n",
       "WACH.38.S1.001.fa                                               2.685130                                                                       \n",
       "PRPI.42.S31.001.fa                                              2.376027                                                                       \n",
       "RODE.80.S8.001.fa                                               2.403992                                                                       \n",
       "pos.dna.01.S17.001.fa                                           3.486100                                                                       \n",
       "MEEN.93.S11.001.fa                                              3.913378                                                                       \n",
       "DETH.41.S27.001.fa                                              1.633435                                                                       \n",
       "TECA.45.S15.001.fa                                              3.626052                                                                       \n",
       "BODA.52.S9.001.fa                                               2.611502                                                                       \n",
       "CHMA.61.S3.001.fa                                               6.117136                                                                       \n",
       "ASTER.02.S34.001.fa                                             3.549830                                                                       \n",
       "BIRO.74.S19.001.fa                                              4.973692                                                                       \n",
       "CAAU.07.S2.001.fa                                               5.706051                                                                       \n",
       "DIGMA.71.S21.001.fa                                             6.540831                                                                       \n",
       "FAFE.16.S10.001.fa                                              5.106302                                                                       \n",
       "OPSE.53.S5.001.fa                                               4.110070                                                                       \n",
       "LAES.91.S29.001.fa                                              4.861458                                                                       \n",
       "DOEC.81.S30.001.fa                                              3.325615                                                                       \n",
       "EKJO.64.S32.001.fa                                              1.952459                                                                       \n",
       "DIAL.15.S23.001.fa                                              3.866595                                                                       \n",
       "GIFI.21.S22.001.fa                                              2.699819                                                                       \n",
       "TAOT.65.S24.001.fa                                              3.131588                                                                       \n",
       "DEFI.14.S6.001.fa                                               3.058043                                                                       \n",
       "MUNA.98.S28.001.fa                                              4.237870                                                                       \n",
       "ROAL.75.S33.001.fa                                              2.787953                                                                       \n",
       "MERI.28.S20.001.fa                                              2.897733                                                                       \n",
       "BUMA.05.S4.001.fa                                               4.452070                                                                       \n",
       "MOEM.48.S16.001.fa                                              4.604422                                                                       \n",
       "PAMA.46.S18.001.fa                                              2.445372                                                                       \n",
       "VUGE.37.S14.001.fa                                              2.790468                                                                       \n",
       "TOSA.76.S25.001.fa                                              2.397261                                                                       \n",
       "DEIS.70.S26.001.fa                                              2.909639                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_faecis  \\\n",
       "pos.dna.02.S35.001.fa                                           1.714401                                                                    \n",
       "CAPMA.83.S12.001.fa                                             3.197282                                                                    \n",
       "FOGI.63.S13.001.fa                                              0.887994                                                                    \n",
       "COCL.11.S7.001.fa                                               2.537797                                                                    \n",
       "WACH.38.S1.001.fa                                               0.050759                                                                    \n",
       "PRPI.42.S31.001.fa                                              1.803052                                                                    \n",
       "RODE.80.S8.001.fa                                               1.124451                                                                    \n",
       "pos.dna.01.S17.001.fa                                           2.323658                                                                    \n",
       "MEEN.93.S11.001.fa                                              1.410490                                                                    \n",
       "DETH.41.S27.001.fa                                              0.197562                                                                    \n",
       "TECA.45.S15.001.fa                                             -0.435956                                                                    \n",
       "BODA.52.S9.001.fa                                               1.548782                                                                    \n",
       "CHMA.61.S3.001.fa                                               3.203964                                                                    \n",
       "ASTER.02.S34.001.fa                                             2.011697                                                                    \n",
       "BIRO.74.S19.001.fa                                              1.036538                                                                    \n",
       "CAAU.07.S2.001.fa                                               1.065846                                                                    \n",
       "DIGMA.71.S21.001.fa                                             0.876363                                                                    \n",
       "FAFE.16.S10.001.fa                                              6.007668                                                                    \n",
       "OPSE.53.S5.001.fa                                               1.131665                                                                    \n",
       "LAES.91.S29.001.fa                                             -0.527922                                                                    \n",
       "DOEC.81.S30.001.fa                                              4.108773                                                                    \n",
       "EKJO.64.S32.001.fa                                              1.117723                                                                    \n",
       "DIAL.15.S23.001.fa                                             -0.562924                                                                    \n",
       "GIFI.21.S22.001.fa                                              1.567684                                                                    \n",
       "TAOT.65.S24.001.fa                                             -0.432295                                                                    \n",
       "DEFI.14.S6.001.fa                                              -0.740189                                                                    \n",
       "MUNA.98.S28.001.fa                                             -0.495943                                                                    \n",
       "ROAL.75.S33.001.fa                                             -0.119041                                                                    \n",
       "MERI.28.S20.001.fa                                              0.579711                                                                    \n",
       "BUMA.05.S4.001.fa                                              -0.874077                                                                    \n",
       "MOEM.48.S16.001.fa                                              2.247472                                                                    \n",
       "PAMA.46.S18.001.fa                                             -0.587451                                                                    \n",
       "VUGE.37.S14.001.fa                                              2.500457                                                                    \n",
       "TOSA.76.S25.001.fa                                              2.802775                                                                    \n",
       "DEIS.70.S26.001.fa                                             -1.153483                                                                    \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_thetaiotaomicron  \\\n",
       "pos.dna.02.S35.001.fa                                           2.002083                                                                              \n",
       "CAPMA.83.S12.001.fa                                             2.198838                                                                              \n",
       "FOGI.63.S13.001.fa                                              2.404863                                                                              \n",
       "COCL.11.S7.001.fa                                               4.726931                                                                              \n",
       "WACH.38.S1.001.fa                                               1.814724                                                                              \n",
       "PRPI.42.S31.001.fa                                              1.172110                                                                              \n",
       "RODE.80.S8.001.fa                                               2.924393                                                                              \n",
       "pos.dna.01.S17.001.fa                                           1.673341                                                                              \n",
       "MEEN.93.S11.001.fa                                              4.790671                                                                              \n",
       "DETH.41.S27.001.fa                                             -0.076657                                                                              \n",
       "TECA.45.S15.001.fa                                             -0.125998                                                                              \n",
       "BODA.52.S9.001.fa                                              -0.011895                                                                              \n",
       "CHMA.61.S3.001.fa                                               5.417908                                                                              \n",
       "ASTER.02.S34.001.fa                                             4.680516                                                                              \n",
       "BIRO.74.S19.001.fa                                              4.116220                                                                              \n",
       "CAAU.07.S2.001.fa                                               3.684255                                                                              \n",
       "DIGMA.71.S21.001.fa                                             1.213984                                                                              \n",
       "FAFE.16.S10.001.fa                                              3.278542                                                                              \n",
       "OPSE.53.S5.001.fa                                               4.004854                                                                              \n",
       "LAES.91.S29.001.fa                                              1.678388                                                                              \n",
       "DOEC.81.S30.001.fa                                              1.343689                                                                              \n",
       "EKJO.64.S32.001.fa                                             -0.472484                                                                              \n",
       "DIAL.15.S23.001.fa                                              2.009991                                                                              \n",
       "GIFI.21.S22.001.fa                                              0.284930                                                                              \n",
       "TAOT.65.S24.001.fa                                             -0.878582                                                                              \n",
       "DEFI.14.S6.001.fa                                               2.086746                                                                              \n",
       "MUNA.98.S28.001.fa                                              1.740706                                                                              \n",
       "ROAL.75.S33.001.fa                                              0.706430                                                                              \n",
       "MERI.28.S20.001.fa                                             -1.315868                                                                              \n",
       "BUMA.05.S4.001.fa                                               1.570038                                                                              \n",
       "MOEM.48.S16.001.fa                                              3.520333                                                                              \n",
       "PAMA.46.S18.001.fa                                             -0.174381                                                                              \n",
       "VUGE.37.S14.001.fa                                              0.155391                                                                              \n",
       "TOSA.76.S25.001.fa                                              0.494040                                                                              \n",
       "DEIS.70.S26.001.fa                                             -1.433785                                                                              \n",
       "\n",
       "#OTU ID                ...  \\\n",
       "pos.dna.02.S35.001.fa  ...   \n",
       "CAPMA.83.S12.001.fa    ...   \n",
       "FOGI.63.S13.001.fa     ...   \n",
       "COCL.11.S7.001.fa      ...   \n",
       "WACH.38.S1.001.fa      ...   \n",
       "PRPI.42.S31.001.fa     ...   \n",
       "RODE.80.S8.001.fa      ...   \n",
       "pos.dna.01.S17.001.fa  ...   \n",
       "MEEN.93.S11.001.fa     ...   \n",
       "DETH.41.S27.001.fa     ...   \n",
       "TECA.45.S15.001.fa     ...   \n",
       "BODA.52.S9.001.fa      ...   \n",
       "CHMA.61.S3.001.fa      ...   \n",
       "ASTER.02.S34.001.fa    ...   \n",
       "BIRO.74.S19.001.fa     ...   \n",
       "CAAU.07.S2.001.fa      ...   \n",
       "DIGMA.71.S21.001.fa    ...   \n",
       "FAFE.16.S10.001.fa     ...   \n",
       "OPSE.53.S5.001.fa      ...   \n",
       "LAES.91.S29.001.fa     ...   \n",
       "DOEC.81.S30.001.fa     ...   \n",
       "EKJO.64.S32.001.fa     ...   \n",
       "DIAL.15.S23.001.fa     ...   \n",
       "GIFI.21.S22.001.fa     ...   \n",
       "TAOT.65.S24.001.fa     ...   \n",
       "DEFI.14.S6.001.fa      ...   \n",
       "MUNA.98.S28.001.fa     ...   \n",
       "ROAL.75.S33.001.fa     ...   \n",
       "MERI.28.S20.001.fa     ...   \n",
       "BUMA.05.S4.001.fa      ...   \n",
       "MOEM.48.S16.001.fa     ...   \n",
       "PAMA.46.S18.001.fa     ...   \n",
       "VUGE.37.S14.001.fa     ...   \n",
       "TOSA.76.S25.001.fa     ...   \n",
       "DEIS.70.S26.001.fa     ...   \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Porphyromonadaceae;g__Porphyromonas;s__Porphyromonas_uenonis_A  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                               \n",
       "CAPMA.83.S12.001.fa                                             3.138940                                                                               \n",
       "FOGI.63.S13.001.fa                                              2.407628                                                                               \n",
       "COCL.11.S7.001.fa                                               2.444156                                                                               \n",
       "WACH.38.S1.001.fa                                              -2.117196                                                                               \n",
       "PRPI.42.S31.001.fa                                              2.369289                                                                               \n",
       "RODE.80.S8.001.fa                                               3.704593                                                                               \n",
       "pos.dna.01.S17.001.fa                                          -0.923157                                                                               \n",
       "MEEN.93.S11.001.fa                                             -2.942135                                                                               \n",
       "DETH.41.S27.001.fa                                              1.070765                                                                               \n",
       "TECA.45.S15.001.fa                                              1.255953                                                                               \n",
       "BODA.52.S9.001.fa                                              -1.812368                                                                               \n",
       "CHMA.61.S3.001.fa                                              -6.522518                                                                               \n",
       "ASTER.02.S34.001.fa                                            -0.781043                                                                               \n",
       "BIRO.74.S19.001.fa                                             -4.852802                                                                               \n",
       "CAAU.07.S2.001.fa                                              -0.491175                                                                               \n",
       "DIGMA.71.S21.001.fa                                            -3.997182                                                                               \n",
       "FAFE.16.S10.001.fa                                             -4.479546                                                                               \n",
       "OPSE.53.S5.001.fa                                              -3.781932                                                                               \n",
       "LAES.91.S29.001.fa                                             -0.659404                                                                               \n",
       "DOEC.81.S30.001.fa                                             -2.364446                                                                               \n",
       "EKJO.64.S32.001.fa                                             -1.004701                                                                               \n",
       "DIAL.15.S23.001.fa                                             -2.826564                                                                               \n",
       "GIFI.21.S22.001.fa                                             -1.117893                                                                               \n",
       "TAOT.65.S24.001.fa                                             -2.552559                                                                               \n",
       "DEFI.14.S6.001.fa                                              -3.038194                                                                               \n",
       "MUNA.98.S28.001.fa                                              0.302494                                                                               \n",
       "ROAL.75.S33.001.fa                                              2.554822                                                                               \n",
       "MERI.28.S20.001.fa                                              0.439290                                                                               \n",
       "BUMA.05.S4.001.fa                                              -0.666852                                                                               \n",
       "MOEM.48.S16.001.fa                                             -4.459434                                                                               \n",
       "PAMA.46.S18.001.fa                                             -1.335513                                                                               \n",
       "VUGE.37.S14.001.fa                                              0.423770                                                                               \n",
       "TOSA.76.S25.001.fa                                              0.618325                                                                               \n",
       "DEIS.70.S26.001.fa                                             -0.406268                                                                               \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_sp001553265  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                       \n",
       "CAPMA.83.S12.001.fa                                             4.841704                                                                       \n",
       "FOGI.63.S13.001.fa                                              1.353132                                                                       \n",
       "COCL.11.S7.001.fa                                               4.645267                                                                       \n",
       "WACH.38.S1.001.fa                                               3.141126                                                                       \n",
       "PRPI.42.S31.001.fa                                             -5.137255                                                                       \n",
       "RODE.80.S8.001.fa                                               1.527456                                                                       \n",
       "pos.dna.01.S17.001.fa                                          -1.010168                                                                       \n",
       "MEEN.93.S11.001.fa                                             -2.301677                                                                       \n",
       "DETH.41.S27.001.fa                                              0.845063                                                                       \n",
       "TECA.45.S15.001.fa                                             -2.019154                                                                       \n",
       "BODA.52.S9.001.fa                                              -2.907754                                                                       \n",
       "CHMA.61.S3.001.fa                                              -5.223235                                                                       \n",
       "ASTER.02.S34.001.fa                                            -1.652275                                                                       \n",
       "BIRO.74.S19.001.fa                                             -4.475508                                                                       \n",
       "CAAU.07.S2.001.fa                                              -0.625905                                                                       \n",
       "DIGMA.71.S21.001.fa                                            -5.422697                                                                       \n",
       "FAFE.16.S10.001.fa                                             -5.285171                                                                       \n",
       "OPSE.53.S5.001.fa                                              -4.310000                                                                       \n",
       "LAES.91.S29.001.fa                                             -3.815749                                                                       \n",
       "DOEC.81.S30.001.fa                                             -2.452479                                                                       \n",
       "EKJO.64.S32.001.fa                                              2.080180                                                                       \n",
       "DIAL.15.S23.001.fa                                             -1.847991                                                                       \n",
       "GIFI.21.S22.001.fa                                              3.923379                                                                       \n",
       "TAOT.65.S24.001.fa                                              3.106923                                                                       \n",
       "DEFI.14.S6.001.fa                                               1.291514                                                                       \n",
       "MUNA.98.S28.001.fa                                              0.470776                                                                       \n",
       "ROAL.75.S33.001.fa                                              1.559750                                                                       \n",
       "MERI.28.S20.001.fa                                              2.206092                                                                       \n",
       "BUMA.05.S4.001.fa                                              -1.504957                                                                       \n",
       "MOEM.48.S16.001.fa                                             -5.226689                                                                       \n",
       "PAMA.46.S18.001.fa                                              1.219461                                                                       \n",
       "VUGE.37.S14.001.fa                                             -1.154333                                                                       \n",
       "TOSA.76.S25.001.fa                                              1.556518                                                                       \n",
       "DEIS.70.S26.001.fa                                              3.801718                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_oris  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                \n",
       "CAPMA.83.S12.001.fa                                             3.443421                                                                \n",
       "FOGI.63.S13.001.fa                                             -0.121505                                                                \n",
       "COCL.11.S7.001.fa                                               1.707124                                                                \n",
       "WACH.38.S1.001.fa                                               1.670049                                                                \n",
       "PRPI.42.S31.001.fa                                             -6.075525                                                                \n",
       "RODE.80.S8.001.fa                                              -1.924759                                                                \n",
       "pos.dna.01.S17.001.fa                                          -1.798626                                                                \n",
       "MEEN.93.S11.001.fa                                             -2.368428                                                                \n",
       "DETH.41.S27.001.fa                                              0.288315                                                                \n",
       "TECA.45.S15.001.fa                                             -2.957723                                                                \n",
       "BODA.52.S9.001.fa                                              -2.278180                                                                \n",
       "CHMA.61.S3.001.fa                                              -4.676691                                                                \n",
       "ASTER.02.S34.001.fa                                            -2.837292                                                                \n",
       "BIRO.74.S19.001.fa                                             -3.103602                                                                \n",
       "CAAU.07.S2.001.fa                                              -4.495425                                                                \n",
       "DIGMA.71.S21.001.fa                                            -2.208633                                                                \n",
       "FAFE.16.S10.001.fa                                             -5.385255                                                                \n",
       "OPSE.53.S5.001.fa                                              -4.449762                                                                \n",
       "LAES.91.S29.001.fa                                             -2.204774                                                                \n",
       "DOEC.81.S30.001.fa                                             -3.483327                                                                \n",
       "EKJO.64.S32.001.fa                                              5.461192                                                                \n",
       "DIAL.15.S23.001.fa                                             -2.486585                                                                \n",
       "GIFI.21.S22.001.fa                                              5.297858                                                                \n",
       "TAOT.65.S24.001.fa                                              5.611904                                                                \n",
       "DEFI.14.S6.001.fa                                              -0.052186                                                                \n",
       "MUNA.98.S28.001.fa                                             -0.653966                                                                \n",
       "ROAL.75.S33.001.fa                                              0.630927                                                                \n",
       "MERI.28.S20.001.fa                                             -0.342471                                                                \n",
       "BUMA.05.S4.001.fa                                              -2.146903                                                                \n",
       "MOEM.48.S16.001.fa                                             -7.098492                                                                \n",
       "PAMA.46.S18.001.fa                                             -0.602835                                                                \n",
       "VUGE.37.S14.001.fa                                             -3.232563                                                                \n",
       "TOSA.76.S25.001.fa                                              1.321092                                                                \n",
       "DEIS.70.S26.001.fa                                             10.340782                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp900066825  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -2.730828                                                                                \n",
       "FOGI.63.S13.001.fa                                             -2.117494                                                                                \n",
       "COCL.11.S7.001.fa                                              -2.901118                                                                                \n",
       "WACH.38.S1.001.fa                                              -2.016654                                                                                \n",
       "PRPI.42.S31.001.fa                                              0.473331                                                                                \n",
       "RODE.80.S8.001.fa                                              -0.191832                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -0.667223                                                                                \n",
       "MEEN.93.S11.001.fa                                              0.770305                                                                                \n",
       "DETH.41.S27.001.fa                                              1.436088                                                                                \n",
       "TECA.45.S15.001.fa                                              0.560763                                                                                \n",
       "BODA.52.S9.001.fa                                               0.806878                                                                                \n",
       "CHMA.61.S3.001.fa                                               0.755190                                                                                \n",
       "ASTER.02.S34.001.fa                                             0.171671                                                                                \n",
       "BIRO.74.S19.001.fa                                              1.275340                                                                                \n",
       "CAAU.07.S2.001.fa                                               0.780260                                                                                \n",
       "DIGMA.71.S21.001.fa                                             4.121940                                                                                \n",
       "FAFE.16.S10.001.fa                                              1.517647                                                                                \n",
       "OPSE.53.S5.001.fa                                              -0.532751                                                                                \n",
       "LAES.91.S29.001.fa                                              1.430708                                                                                \n",
       "DOEC.81.S30.001.fa                                              0.449270                                                                                \n",
       "EKJO.64.S32.001.fa                                             -0.952057                                                                                \n",
       "DIAL.15.S23.001.fa                                              0.842740                                                                                \n",
       "GIFI.21.S22.001.fa                                              0.153738                                                                                \n",
       "TAOT.65.S24.001.fa                                              0.568337                                                                                \n",
       "DEFI.14.S6.001.fa                                              -0.119842                                                                                \n",
       "MUNA.98.S28.001.fa                                             -0.514026                                                                                \n",
       "ROAL.75.S33.001.fa                                              0.524060                                                                                \n",
       "MERI.28.S20.001.fa                                             -0.432335                                                                                \n",
       "BUMA.05.S4.001.fa                                              -0.352957                                                                                \n",
       "MOEM.48.S16.001.fa                                              0.976780                                                                                \n",
       "PAMA.46.S18.001.fa                                             -0.996740                                                                                \n",
       "VUGE.37.S14.001.fa                                              0.427640                                                                                \n",
       "TOSA.76.S25.001.fa                                              0.577063                                                                                \n",
       "DEIS.70.S26.001.fa                                              1.103640                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Faecalimonas;s__Faecalimonas_sp900556835  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                            \n",
       "CAPMA.83.S12.001.fa                                            -2.124692                                                                            \n",
       "FOGI.63.S13.001.fa                                             -3.004797                                                                            \n",
       "COCL.11.S7.001.fa                                              -3.214776                                                                            \n",
       "WACH.38.S1.001.fa                                              -1.179972                                                                            \n",
       "PRPI.42.S31.001.fa                                             -1.156355                                                                            \n",
       "RODE.80.S8.001.fa                                              -1.928109                                                                            \n",
       "pos.dna.01.S17.001.fa                                          -0.635475                                                                            \n",
       "MEEN.93.S11.001.fa                                              1.697052                                                                            \n",
       "DETH.41.S27.001.fa                                             -1.135138                                                                            \n",
       "TECA.45.S15.001.fa                                             -1.981214                                                                            \n",
       "BODA.52.S9.001.fa                                              -1.435304                                                                            \n",
       "CHMA.61.S3.001.fa                                               0.602765                                                                            \n",
       "ASTER.02.S34.001.fa                                            -0.867344                                                                            \n",
       "BIRO.74.S19.001.fa                                             -0.624752                                                                            \n",
       "CAAU.07.S2.001.fa                                               0.033700                                                                            \n",
       "DIGMA.71.S21.001.fa                                            -0.361749                                                                            \n",
       "FAFE.16.S10.001.fa                                             -0.069201                                                                            \n",
       "OPSE.53.S5.001.fa                                              -0.095620                                                                            \n",
       "LAES.91.S29.001.fa                                             -2.543399                                                                            \n",
       "DOEC.81.S30.001.fa                                             -1.440480                                                                            \n",
       "EKJO.64.S32.001.fa                                             -2.130712                                                                            \n",
       "DIAL.15.S23.001.fa                                             -1.756216                                                                            \n",
       "GIFI.21.S22.001.fa                                             -2.727331                                                                            \n",
       "TAOT.65.S24.001.fa                                             -2.958024                                                                            \n",
       "DEFI.14.S6.001.fa                                               1.011532                                                                            \n",
       "MUNA.98.S28.001.fa                                             -0.816727                                                                            \n",
       "ROAL.75.S33.001.fa                                             -2.359089                                                                            \n",
       "MERI.28.S20.001.fa                                             -3.599673                                                                            \n",
       "BUMA.05.S4.001.fa                                              -1.530629                                                                            \n",
       "MOEM.48.S16.001.fa                                              2.587275                                                                            \n",
       "PAMA.46.S18.001.fa                                             -0.908771                                                                            \n",
       "VUGE.37.S14.001.fa                                             -0.313027                                                                            \n",
       "TOSA.76.S25.001.fa                                             -2.275568                                                                            \n",
       "DEIS.70.S26.001.fa                                             -3.573851                                                                            \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Hungatella_A;s__Hungatella_A_sp003478355  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                            \n",
       "CAPMA.83.S12.001.fa                                            -0.929852                                                                            \n",
       "FOGI.63.S13.001.fa                                             -1.820056                                                                            \n",
       "COCL.11.S7.001.fa                                              -2.330573                                                                            \n",
       "WACH.38.S1.001.fa                                              -1.957565                                                                            \n",
       "PRPI.42.S31.001.fa                                             -0.917108                                                                            \n",
       "RODE.80.S8.001.fa                                              -1.326253                                                                            \n",
       "pos.dna.01.S17.001.fa                                          -1.010168                                                                            \n",
       "MEEN.93.S11.001.fa                                              0.457184                                                                            \n",
       "DETH.41.S27.001.fa                                             -0.484491                                                                            \n",
       "TECA.45.S15.001.fa                                              0.301699                                                                            \n",
       "BODA.52.S9.001.fa                                               0.612192                                                                            \n",
       "CHMA.61.S3.001.fa                                               0.508339                                                                            \n",
       "ASTER.02.S34.001.fa                                            -0.210790                                                                            \n",
       "BIRO.74.S19.001.fa                                             -0.903965                                                                            \n",
       "CAAU.07.S2.001.fa                                              -1.256430                                                                            \n",
       "DIGMA.71.S21.001.fa                                             1.308846                                                                            \n",
       "FAFE.16.S10.001.fa                                              1.602478                                                                            \n",
       "OPSE.53.S5.001.fa                                              -1.119344                                                                            \n",
       "LAES.91.S29.001.fa                                              0.343855                                                                            \n",
       "DOEC.81.S30.001.fa                                             -0.737190                                                                            \n",
       "EKJO.64.S32.001.fa                                             -0.723798                                                                            \n",
       "DIAL.15.S23.001.fa                                              1.252917                                                                            \n",
       "GIFI.21.S22.001.fa                                             -1.341037                                                                            \n",
       "TAOT.65.S24.001.fa                                             -1.399879                                                                            \n",
       "DEFI.14.S6.001.fa                                              -1.295560                                                                            \n",
       "MUNA.98.S28.001.fa                                             -0.916989                                                                            \n",
       "ROAL.75.S33.001.fa                                             -2.246501                                                                            \n",
       "MERI.28.S20.001.fa                                             -0.984146                                                                            \n",
       "BUMA.05.S4.001.fa                                              -0.746264                                                                            \n",
       "MOEM.48.S16.001.fa                                              1.486921                                                                            \n",
       "PAMA.46.S18.001.fa                                             -0.037698                                                                            \n",
       "VUGE.37.S14.001.fa                                              0.374880                                                                            \n",
       "TOSA.76.S25.001.fa                                             -0.003606                                                                            \n",
       "DEIS.70.S26.001.fa                                              0.557308                                                                            \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp000177015  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -4.522588                                                                                \n",
       "FOGI.63.S13.001.fa                                             -3.537014                                                                                \n",
       "COCL.11.S7.001.fa                                              -4.213305                                                                                \n",
       "WACH.38.S1.001.fa                                              -4.815676                                                                                \n",
       "PRPI.42.S31.001.fa                                             -1.554945                                                                                \n",
       "RODE.80.S8.001.fa                                              -1.318927                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -2.309451                                                                                \n",
       "MEEN.93.S11.001.fa                                             -0.315932                                                                                \n",
       "DETH.41.S27.001.fa                                             -1.364955                                                                                \n",
       "TECA.45.S15.001.fa                                             -1.118817                                                                                \n",
       "BODA.52.S9.001.fa                                              -0.467510                                                                                \n",
       "CHMA.61.S3.001.fa                                              -0.216243                                                                                \n",
       "ASTER.02.S34.001.fa                                            -1.297471                                                                                \n",
       "BIRO.74.S19.001.fa                                              1.473348                                                                                \n",
       "CAAU.07.S2.001.fa                                              -2.549514                                                                                \n",
       "DIGMA.71.S21.001.fa                                            -0.674639                                                                                \n",
       "FAFE.16.S10.001.fa                                              0.442297                                                                                \n",
       "OPSE.53.S5.001.fa                                              -2.203747                                                                                \n",
       "LAES.91.S29.001.fa                                              0.111238                                                                                \n",
       "DOEC.81.S30.001.fa                                             -1.007263                                                                                \n",
       "EKJO.64.S32.001.fa                                             -2.536177                                                                                \n",
       "DIAL.15.S23.001.fa                                             -0.041654                                                                                \n",
       "GIFI.21.S22.001.fa                                             -3.132796                                                                                \n",
       "TAOT.65.S24.001.fa                                             -2.734880                                                                                \n",
       "DEFI.14.S6.001.fa                                              -2.430593                                                                                \n",
       "MUNA.98.S28.001.fa                                             -1.873143                                                                                \n",
       "ROAL.75.S33.001.fa                                             -3.651487                                                                                \n",
       "MERI.28.S20.001.fa                                             -0.108111                                                                                \n",
       "BUMA.05.S4.001.fa                                              -1.368205                                                                                \n",
       "MOEM.48.S16.001.fa                                             -1.131064                                                                                \n",
       "PAMA.46.S18.001.fa                                             -2.690059                                                                                \n",
       "VUGE.37.S14.001.fa                                             -0.862795                                                                                \n",
       "TOSA.76.S25.001.fa                                             -2.781117                                                                                \n",
       "DEIS.70.S26.001.fa                                             -3.350707                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Paraprevotella;s__Paraprevotella_sp003477995  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                               \n",
       "CAPMA.83.S12.001.fa                                            -2.443146                                                                               \n",
       "FOGI.63.S13.001.fa                                             -2.626143                                                                               \n",
       "COCL.11.S7.001.fa                                              -3.594266                                                                               \n",
       "WACH.38.S1.001.fa                                              -1.953475                                                                               \n",
       "PRPI.42.S31.001.fa                                             -1.227844                                                                               \n",
       "RODE.80.S8.001.fa                                              -1.467165                                                                               \n",
       "pos.dna.01.S17.001.fa                                          -2.714916                                                                               \n",
       "MEEN.93.S11.001.fa                                             -3.596317                                                                               \n",
       "DETH.41.S27.001.fa                                             -0.362683                                                                               \n",
       "TECA.45.S15.001.fa                                             -2.342227                                                                               \n",
       "BODA.52.S9.001.fa                                              -0.727855                                                                               \n",
       "CHMA.61.S3.001.fa                                               0.590895                                                                               \n",
       "ASTER.02.S34.001.fa                                             1.110924                                                                               \n",
       "BIRO.74.S19.001.fa                                              1.265571                                                                               \n",
       "CAAU.07.S2.001.fa                                               0.328233                                                                               \n",
       "DIGMA.71.S21.001.fa                                             0.604328                                                                               \n",
       "FAFE.16.S10.001.fa                                             -0.797606                                                                               \n",
       "OPSE.53.S5.001.fa                                               0.762726                                                                               \n",
       "LAES.91.S29.001.fa                                             -2.009986                                                                               \n",
       "DOEC.81.S30.001.fa                                             -2.681997                                                                               \n",
       "EKJO.64.S32.001.fa                                             -3.517006                                                                               \n",
       "DIAL.15.S23.001.fa                                             -1.053558                                                                               \n",
       "GIFI.21.S22.001.fa                                             -2.909653                                                                               \n",
       "TAOT.65.S24.001.fa                                             -2.264877                                                                               \n",
       "DEFI.14.S6.001.fa                                              -3.090619                                                                               \n",
       "MUNA.98.S28.001.fa                                             -4.406562                                                                               \n",
       "ROAL.75.S33.001.fa                                             -2.515834                                                                               \n",
       "MERI.28.S20.001.fa                                             -3.217738                                                                               \n",
       "BUMA.05.S4.001.fa                                              -2.212961                                                                               \n",
       "MOEM.48.S16.001.fa                                             -0.922624                                                                               \n",
       "PAMA.46.S18.001.fa                                             -1.335513                                                                               \n",
       "VUGE.37.S14.001.fa                                             -2.197253                                                                               \n",
       "TOSA.76.S25.001.fa                                             -3.086498                                                                               \n",
       "DEIS.70.S26.001.fa                                             -1.246573                                                                               \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_oleiciplenus  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                          \n",
       "CAPMA.83.S12.001.fa                                            -3.675290                                                                          \n",
       "FOGI.63.S13.001.fa                                             -1.696893                                                                          \n",
       "COCL.11.S7.001.fa                                              -0.489334                                                                          \n",
       "WACH.38.S1.001.fa                                              -2.093441                                                                          \n",
       "PRPI.42.S31.001.fa                                             -0.051808                                                                          \n",
       "RODE.80.S8.001.fa                                              -2.030491                                                                          \n",
       "pos.dna.01.S17.001.fa                                          -2.155300                                                                          \n",
       "MEEN.93.S11.001.fa                                             -1.389233                                                                          \n",
       "DETH.41.S27.001.fa                                             -2.735072                                                                          \n",
       "TECA.45.S15.001.fa                                             -1.738692                                                                          \n",
       "BODA.52.S9.001.fa                                              -1.865448                                                                          \n",
       "CHMA.61.S3.001.fa                                               0.945615                                                                          \n",
       "ASTER.02.S34.001.fa                                            -0.424424                                                                          \n",
       "BIRO.74.S19.001.fa                                              1.526685                                                                          \n",
       "CAAU.07.S2.001.fa                                               0.564564                                                                          \n",
       "DIGMA.71.S21.001.fa                                            -0.262493                                                                          \n",
       "FAFE.16.S10.001.fa                                             -1.121834                                                                          \n",
       "OPSE.53.S5.001.fa                                              -0.606017                                                                          \n",
       "LAES.91.S29.001.fa                                              0.495566                                                                          \n",
       "DOEC.81.S30.001.fa                                             -0.129765                                                                          \n",
       "EKJO.64.S32.001.fa                                             -1.907569                                                                          \n",
       "DIAL.15.S23.001.fa                                              0.031977                                                                          \n",
       "GIFI.21.S22.001.fa                                             -2.573181                                                                          \n",
       "TAOT.65.S24.001.fa                                             -2.552559                                                                          \n",
       "DEFI.14.S6.001.fa                                              -1.534714                                                                          \n",
       "MUNA.98.S28.001.fa                                             -1.715021                                                                          \n",
       "ROAL.75.S33.001.fa                                             -1.401189                                                                          \n",
       "MERI.28.S20.001.fa                                             -3.448850                                                                          \n",
       "BUMA.05.S4.001.fa                                              -1.443736                                                                          \n",
       "MOEM.48.S16.001.fa                                              0.683690                                                                          \n",
       "PAMA.46.S18.001.fa                                             -2.977741                                                                          \n",
       "VUGE.37.S14.001.fa                                             -1.232786                                                                          \n",
       "TOSA.76.S25.001.fa                                             -2.393351                                                                          \n",
       "DEIS.70.S26.001.fa                                              0.276297                                                                          \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__CAG-81;s__CAG-81_sp900066785  \n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                               \n",
       "CAPMA.83.S12.001.fa                                            -3.829440                                                               \n",
       "FOGI.63.S13.001.fa                                             -2.913825                                                               \n",
       "COCL.11.S7.001.fa                                              -3.856630                                                               \n",
       "WACH.38.S1.001.fa                                              -1.995148                                                               \n",
       "PRPI.42.S31.001.fa                                             -1.882509                                                               \n",
       "RODE.80.S8.001.fa                                              -1.091414                                                               \n",
       "pos.dna.01.S17.001.fa                                          -1.703315                                                               \n",
       "MEEN.93.S11.001.fa                                             -1.218305                                                               \n",
       "DETH.41.S27.001.fa                                              0.163711                                                               \n",
       "TECA.45.S15.001.fa                                             -0.965607                                                               \n",
       "BODA.52.S9.001.fa                                               0.195402                                                               \n",
       "CHMA.61.S3.001.fa                                              -0.637804                                                               \n",
       "ASTER.02.S34.001.fa                                             0.490222                                                               \n",
       "BIRO.74.S19.001.fa                                             -1.350001                                                               \n",
       "CAAU.07.S2.001.fa                                              -0.054863                                                               \n",
       "DIGMA.71.S21.001.fa                                             0.456438                                                               \n",
       "FAFE.16.S10.001.fa                                              0.994868                                                               \n",
       "OPSE.53.S5.001.fa                                              -1.502694                                                               \n",
       "LAES.91.S29.001.fa                                             -1.028367                                                               \n",
       "DOEC.81.S30.001.fa                                             -0.393402                                                               \n",
       "EKJO.64.S32.001.fa                                             -1.843030                                                               \n",
       "DIAL.15.S23.001.fa                                             -0.745484                                                               \n",
       "GIFI.21.S22.001.fa                                             -3.132796                                                               \n",
       "TAOT.65.S24.001.fa                                             -1.208824                                                               \n",
       "DEFI.14.S6.001.fa                                              -1.908701                                                               \n",
       "MUNA.98.S28.001.fa                                             -0.154564                                                               \n",
       "ROAL.75.S33.001.fa                                             -2.638173                                                               \n",
       "MERI.28.S20.001.fa                                             -2.688044                                                               \n",
       "BUMA.05.S4.001.fa                                              -1.600636                                                               \n",
       "MOEM.48.S16.001.fa                                              0.319989                                                               \n",
       "PAMA.46.S18.001.fa                                             -2.327154                                                               \n",
       "VUGE.37.S14.001.fa                                             -0.935806                                                               \n",
       "TOSA.76.S25.001.fa                                             -2.393351                                                               \n",
       "DEIS.70.S26.001.fa                                             -3.350707                                                               \n",
       "\n",
       "[35 rows x 181 columns]"
      ]
     },
     "execution_count": 63,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_taxa_clr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [],
   "source": [
    "from skbio.stats.composition import ancom"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [],
   "source": [
    "mask_controls = [not _.startswith(\"pos\") for _ in df_taxa.index]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa.index = sample_ids\n",
    "\n",
    "df_taxa_nc = df_taxa.copy()\n",
    "\n",
    "df_taxa_nc[\"sample_id\"] = sample_ids\n",
    "\n",
    "df_taxa_nc = df_taxa_nc.loc[mask_controls].copy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(33, 182)"
      ]
     },
     "execution_count": 67,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_taxa_nc.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_nc = pd.merge(df_taxa_nc, df_metadata_file[[\"RDV\", \"sample_id\"]], on=\"sample_id\", how=\"inner\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [],
   "source": [
    "ancom_results = ancom(df_taxa_nc.drop(columns=[\"RDV\", \"sample_id\"]) + 1, df_taxa_nc[\"RDV\"])[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 70,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ancom_results[\"Reject null hypothesis\"].sum()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 119,
   "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>#OTU ID</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Phocaeicola;s__Phocaeicola_vulgatus</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_fragilis</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli_D</th>\n",
       "      <th>k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_flexneri</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_bouchesdurhonensis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Tannerellaceae;g__Parabacteroides;s__Parabacteroides_distasonis_A</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_uniformis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_faecis</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_thetaiotaomicron</th>\n",
       "      <th>...</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Porphyromonadaceae;g__Porphyromonas;s__Porphyromonas_uenonis_A</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_sp001553265</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_oris</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp900066825</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Faecalimonas;s__Faecalimonas_sp900556835</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Hungatella_A;s__Hungatella_A_sp003478355</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp000177015</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Paraprevotella;s__Paraprevotella_sp003477995</th>\n",
       "      <th>k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_oleiciplenus</th>\n",
       "      <th>k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__CAG-81;s__CAG-81_sp900066785</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>pos.dna.02.S35.001.fa</th>\n",
       "      <td>12.564805</td>\n",
       "      <td>12.596202</td>\n",
       "      <td>10.432703</td>\n",
       "      <td>10.382420</td>\n",
       "      <td>10.717464</td>\n",
       "      <td>-0.529343</td>\n",
       "      <td>-0.683494</td>\n",
       "      <td>5.178716</td>\n",
       "      <td>1.714401</td>\n",
       "      <td>2.002083</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "      <td>-2.475253</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CAPMA.83.S12.001.fa</th>\n",
       "      <td>1.502473</td>\n",
       "      <td>0.108900</td>\n",
       "      <td>8.026002</td>\n",
       "      <td>6.006419</td>\n",
       "      <td>4.688253</td>\n",
       "      <td>-3.136293</td>\n",
       "      <td>0.814951</td>\n",
       "      <td>1.198817</td>\n",
       "      <td>3.197282</td>\n",
       "      <td>2.198838</td>\n",
       "      <td>...</td>\n",
       "      <td>3.138940</td>\n",
       "      <td>4.841704</td>\n",
       "      <td>3.443421</td>\n",
       "      <td>-2.730828</td>\n",
       "      <td>-2.124692</td>\n",
       "      <td>-0.929852</td>\n",
       "      <td>-4.522588</td>\n",
       "      <td>-2.443146</td>\n",
       "      <td>-3.675290</td>\n",
       "      <td>-3.829440</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>FOGI.63.S13.001.fa</th>\n",
       "      <td>3.594296</td>\n",
       "      <td>-0.611240</td>\n",
       "      <td>0.322862</td>\n",
       "      <td>-0.505445</td>\n",
       "      <td>-2.583584</td>\n",
       "      <td>-0.913521</td>\n",
       "      <td>-0.274768</td>\n",
       "      <td>2.714890</td>\n",
       "      <td>0.887994</td>\n",
       "      <td>2.404863</td>\n",
       "      <td>...</td>\n",
       "      <td>2.407628</td>\n",
       "      <td>1.353132</td>\n",
       "      <td>-0.121505</td>\n",
       "      <td>-2.117494</td>\n",
       "      <td>-3.004797</td>\n",
       "      <td>-1.820056</td>\n",
       "      <td>-3.537014</td>\n",
       "      <td>-2.626143</td>\n",
       "      <td>-1.696893</td>\n",
       "      <td>-2.913825</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>COCL.11.S7.001.fa</th>\n",
       "      <td>1.536088</td>\n",
       "      <td>1.074962</td>\n",
       "      <td>4.663100</td>\n",
       "      <td>3.188711</td>\n",
       "      <td>1.774582</td>\n",
       "      <td>-2.081677</td>\n",
       "      <td>-2.133863</td>\n",
       "      <td>1.997295</td>\n",
       "      <td>2.537797</td>\n",
       "      <td>4.726931</td>\n",
       "      <td>...</td>\n",
       "      <td>2.444156</td>\n",
       "      <td>4.645267</td>\n",
       "      <td>1.707124</td>\n",
       "      <td>-2.901118</td>\n",
       "      <td>-3.214776</td>\n",
       "      <td>-2.330573</td>\n",
       "      <td>-4.213305</td>\n",
       "      <td>-3.594266</td>\n",
       "      <td>-0.489334</td>\n",
       "      <td>-3.856630</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>WACH.38.S1.001.fa</th>\n",
       "      <td>4.240597</td>\n",
       "      <td>-0.421227</td>\n",
       "      <td>2.999156</td>\n",
       "      <td>1.341303</td>\n",
       "      <td>0.023436</td>\n",
       "      <td>-1.574201</td>\n",
       "      <td>-1.208627</td>\n",
       "      <td>2.685130</td>\n",
       "      <td>0.050759</td>\n",
       "      <td>1.814724</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.117196</td>\n",
       "      <td>3.141126</td>\n",
       "      <td>1.670049</td>\n",
       "      <td>-2.016654</td>\n",
       "      <td>-1.179972</td>\n",
       "      <td>-1.957565</td>\n",
       "      <td>-4.815676</td>\n",
       "      <td>-1.953475</td>\n",
       "      <td>-2.093441</td>\n",
       "      <td>-1.995148</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PRPI.42.S31.001.fa</th>\n",
       "      <td>3.169797</td>\n",
       "      <td>-1.438640</td>\n",
       "      <td>7.172793</td>\n",
       "      <td>4.708191</td>\n",
       "      <td>3.824696</td>\n",
       "      <td>-0.532085</td>\n",
       "      <td>0.522833</td>\n",
       "      <td>2.376027</td>\n",
       "      <td>1.803052</td>\n",
       "      <td>1.172110</td>\n",
       "      <td>...</td>\n",
       "      <td>2.369289</td>\n",
       "      <td>-5.137255</td>\n",
       "      <td>-6.075525</td>\n",
       "      <td>0.473331</td>\n",
       "      <td>-1.156355</td>\n",
       "      <td>-0.917108</td>\n",
       "      <td>-1.554945</td>\n",
       "      <td>-1.227844</td>\n",
       "      <td>-0.051808</td>\n",
       "      <td>-1.882509</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>RODE.80.S8.001.fa</th>\n",
       "      <td>5.918568</td>\n",
       "      <td>-1.343871</td>\n",
       "      <td>4.211689</td>\n",
       "      <td>4.242758</td>\n",
       "      <td>4.215763</td>\n",
       "      <td>-0.235021</td>\n",
       "      <td>1.030012</td>\n",
       "      <td>2.403992</td>\n",
       "      <td>1.124451</td>\n",
       "      <td>2.924393</td>\n",
       "      <td>...</td>\n",
       "      <td>3.704593</td>\n",
       "      <td>1.527456</td>\n",
       "      <td>-1.924759</td>\n",
       "      <td>-0.191832</td>\n",
       "      <td>-1.928109</td>\n",
       "      <td>-1.326253</td>\n",
       "      <td>-1.318927</td>\n",
       "      <td>-1.467165</td>\n",
       "      <td>-2.030491</td>\n",
       "      <td>-1.091414</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>pos.dna.01.S17.001.fa</th>\n",
       "      <td>10.758325</td>\n",
       "      <td>10.466699</td>\n",
       "      <td>8.427199</td>\n",
       "      <td>8.411471</td>\n",
       "      <td>8.758813</td>\n",
       "      <td>-1.616304</td>\n",
       "      <td>-0.923157</td>\n",
       "      <td>3.486100</td>\n",
       "      <td>2.323658</td>\n",
       "      <td>1.673341</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.923157</td>\n",
       "      <td>-1.010168</td>\n",
       "      <td>-1.798626</td>\n",
       "      <td>-0.667223</td>\n",
       "      <td>-0.635475</td>\n",
       "      <td>-1.010168</td>\n",
       "      <td>-2.309451</td>\n",
       "      <td>-2.714916</td>\n",
       "      <td>-2.155300</td>\n",
       "      <td>-1.703315</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MEEN.93.S11.001.fa</th>\n",
       "      <td>4.110142</td>\n",
       "      <td>6.091029</td>\n",
       "      <td>4.004449</td>\n",
       "      <td>3.807656</td>\n",
       "      <td>2.422215</td>\n",
       "      <td>-0.461099</td>\n",
       "      <td>0.197382</td>\n",
       "      <td>3.913378</td>\n",
       "      <td>1.410490</td>\n",
       "      <td>4.790671</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.942135</td>\n",
       "      <td>-2.301677</td>\n",
       "      <td>-2.368428</td>\n",
       "      <td>0.770305</td>\n",
       "      <td>1.697052</td>\n",
       "      <td>0.457184</td>\n",
       "      <td>-0.315932</td>\n",
       "      <td>-3.596317</td>\n",
       "      <td>-1.389233</td>\n",
       "      <td>-1.218305</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DETH.41.S27.001.fa</th>\n",
       "      <td>3.183976</td>\n",
       "      <td>-3.761992</td>\n",
       "      <td>-0.928224</td>\n",
       "      <td>-2.920100</td>\n",
       "      <td>-3.462587</td>\n",
       "      <td>0.307518</td>\n",
       "      <td>-3.401551</td>\n",
       "      <td>1.633435</td>\n",
       "      <td>0.197562</td>\n",
       "      <td>-0.076657</td>\n",
       "      <td>...</td>\n",
       "      <td>1.070765</td>\n",
       "      <td>0.845063</td>\n",
       "      <td>0.288315</td>\n",
       "      <td>1.436088</td>\n",
       "      <td>-1.135138</td>\n",
       "      <td>-0.484491</td>\n",
       "      <td>-1.364955</td>\n",
       "      <td>-0.362683</td>\n",
       "      <td>-2.735072</td>\n",
       "      <td>0.163711</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TECA.45.S15.001.fa</th>\n",
       "      <td>3.134289</td>\n",
       "      <td>4.631672</td>\n",
       "      <td>2.388486</td>\n",
       "      <td>3.961840</td>\n",
       "      <td>-0.605582</td>\n",
       "      <td>-0.851349</td>\n",
       "      <td>0.515884</td>\n",
       "      <td>3.626052</td>\n",
       "      <td>-0.435956</td>\n",
       "      <td>-0.125998</td>\n",
       "      <td>...</td>\n",
       "      <td>1.255953</td>\n",
       "      <td>-2.019154</td>\n",
       "      <td>-2.957723</td>\n",
       "      <td>0.560763</td>\n",
       "      <td>-1.981214</td>\n",
       "      <td>0.301699</td>\n",
       "      <td>-1.118817</td>\n",
       "      <td>-2.342227</td>\n",
       "      <td>-1.738692</td>\n",
       "      <td>-0.965607</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BODA.52.S9.001.fa</th>\n",
       "      <td>4.091882</td>\n",
       "      <td>1.095134</td>\n",
       "      <td>0.574600</td>\n",
       "      <td>-0.936630</td>\n",
       "      <td>-1.537898</td>\n",
       "      <td>0.710902</td>\n",
       "      <td>-0.204385</td>\n",
       "      <td>2.611502</td>\n",
       "      <td>1.548782</td>\n",
       "      <td>-0.011895</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.812368</td>\n",
       "      <td>-2.907754</td>\n",
       "      <td>-2.278180</td>\n",
       "      <td>0.806878</td>\n",
       "      <td>-1.435304</td>\n",
       "      <td>0.612192</td>\n",
       "      <td>-0.467510</td>\n",
       "      <td>-0.727855</td>\n",
       "      <td>-1.865448</td>\n",
       "      <td>0.195402</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CHMA.61.S3.001.fa</th>\n",
       "      <td>6.751274</td>\n",
       "      <td>1.924646</td>\n",
       "      <td>3.074376</td>\n",
       "      <td>2.536282</td>\n",
       "      <td>1.136575</td>\n",
       "      <td>1.722341</td>\n",
       "      <td>1.558442</td>\n",
       "      <td>6.117136</td>\n",
       "      <td>3.203964</td>\n",
       "      <td>5.417908</td>\n",
       "      <td>...</td>\n",
       "      <td>-6.522518</td>\n",
       "      <td>-5.223235</td>\n",
       "      <td>-4.676691</td>\n",
       "      <td>0.755190</td>\n",
       "      <td>0.602765</td>\n",
       "      <td>0.508339</td>\n",
       "      <td>-0.216243</td>\n",
       "      <td>0.590895</td>\n",
       "      <td>0.945615</td>\n",
       "      <td>-0.637804</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ASTER.02.S34.001.fa</th>\n",
       "      <td>7.511852</td>\n",
       "      <td>-0.131282</td>\n",
       "      <td>-3.570534</td>\n",
       "      <td>-3.458830</td>\n",
       "      <td>-3.434732</td>\n",
       "      <td>0.951255</td>\n",
       "      <td>1.994439</td>\n",
       "      <td>3.549830</td>\n",
       "      <td>2.011697</td>\n",
       "      <td>4.680516</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.781043</td>\n",
       "      <td>-1.652275</td>\n",
       "      <td>-2.837292</td>\n",
       "      <td>0.171671</td>\n",
       "      <td>-0.867344</td>\n",
       "      <td>-0.210790</td>\n",
       "      <td>-1.297471</td>\n",
       "      <td>1.110924</td>\n",
       "      <td>-0.424424</td>\n",
       "      <td>0.490222</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BIRO.74.S19.001.fa</th>\n",
       "      <td>6.603156</td>\n",
       "      <td>4.111350</td>\n",
       "      <td>-1.495615</td>\n",
       "      <td>-1.615759</td>\n",
       "      <td>-2.541918</td>\n",
       "      <td>-0.057356</td>\n",
       "      <td>-0.108232</td>\n",
       "      <td>4.973692</td>\n",
       "      <td>1.036538</td>\n",
       "      <td>4.116220</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.852802</td>\n",
       "      <td>-4.475508</td>\n",
       "      <td>-3.103602</td>\n",
       "      <td>1.275340</td>\n",
       "      <td>-0.624752</td>\n",
       "      <td>-0.903965</td>\n",
       "      <td>1.473348</td>\n",
       "      <td>1.265571</td>\n",
       "      <td>1.526685</td>\n",
       "      <td>-1.350001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CAAU.07.S2.001.fa</th>\n",
       "      <td>5.567414</td>\n",
       "      <td>-1.544217</td>\n",
       "      <td>-4.652267</td>\n",
       "      <td>-5.025247</td>\n",
       "      <td>-6.750757</td>\n",
       "      <td>-0.619196</td>\n",
       "      <td>1.828631</td>\n",
       "      <td>5.706051</td>\n",
       "      <td>1.065846</td>\n",
       "      <td>3.684255</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.491175</td>\n",
       "      <td>-0.625905</td>\n",
       "      <td>-4.495425</td>\n",
       "      <td>0.780260</td>\n",
       "      <td>0.033700</td>\n",
       "      <td>-1.256430</td>\n",
       "      <td>-2.549514</td>\n",
       "      <td>0.328233</td>\n",
       "      <td>0.564564</td>\n",
       "      <td>-0.054863</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DIGMA.71.S21.001.fa</th>\n",
       "      <td>5.374936</td>\n",
       "      <td>0.454759</td>\n",
       "      <td>-1.006752</td>\n",
       "      <td>0.006649</td>\n",
       "      <td>-1.588852</td>\n",
       "      <td>-0.308543</td>\n",
       "      <td>-0.375151</td>\n",
       "      <td>6.540831</td>\n",
       "      <td>0.876363</td>\n",
       "      <td>1.213984</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.997182</td>\n",
       "      <td>-5.422697</td>\n",
       "      <td>-2.208633</td>\n",
       "      <td>4.121940</td>\n",
       "      <td>-0.361749</td>\n",
       "      <td>1.308846</td>\n",
       "      <td>-0.674639</td>\n",
       "      <td>0.604328</td>\n",
       "      <td>-0.262493</td>\n",
       "      <td>0.456438</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>FAFE.16.S10.001.fa</th>\n",
       "      <td>5.584302</td>\n",
       "      <td>-0.067135</td>\n",
       "      <td>-1.416951</td>\n",
       "      <td>-2.181225</td>\n",
       "      <td>-2.503694</td>\n",
       "      <td>1.577985</td>\n",
       "      <td>0.749284</td>\n",
       "      <td>5.106302</td>\n",
       "      <td>6.007668</td>\n",
       "      <td>3.278542</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.479546</td>\n",
       "      <td>-5.285171</td>\n",
       "      <td>-5.385255</td>\n",
       "      <td>1.517647</td>\n",
       "      <td>-0.069201</td>\n",
       "      <td>1.602478</td>\n",
       "      <td>0.442297</td>\n",
       "      <td>-0.797606</td>\n",
       "      <td>-1.121834</td>\n",
       "      <td>0.994868</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>OPSE.53.S5.001.fa</th>\n",
       "      <td>7.120223</td>\n",
       "      <td>1.635535</td>\n",
       "      <td>1.917310</td>\n",
       "      <td>-0.146372</td>\n",
       "      <td>-1.667841</td>\n",
       "      <td>1.399707</td>\n",
       "      <td>3.893911</td>\n",
       "      <td>4.110070</td>\n",
       "      <td>1.131665</td>\n",
       "      <td>4.004854</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.781932</td>\n",
       "      <td>-4.310000</td>\n",
       "      <td>-4.449762</td>\n",
       "      <td>-0.532751</td>\n",
       "      <td>-0.095620</td>\n",
       "      <td>-1.119344</td>\n",
       "      <td>-2.203747</td>\n",
       "      <td>0.762726</td>\n",
       "      <td>-0.606017</td>\n",
       "      <td>-1.502694</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>LAES.91.S29.001.fa</th>\n",
       "      <td>3.183280</td>\n",
       "      <td>-0.423308</td>\n",
       "      <td>5.162297</td>\n",
       "      <td>3.618158</td>\n",
       "      <td>2.340659</td>\n",
       "      <td>-0.533104</td>\n",
       "      <td>-2.047337</td>\n",
       "      <td>4.861458</td>\n",
       "      <td>-0.527922</td>\n",
       "      <td>1.678388</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.659404</td>\n",
       "      <td>-3.815749</td>\n",
       "      <td>-2.204774</td>\n",
       "      <td>1.430708</td>\n",
       "      <td>-2.543399</td>\n",
       "      <td>0.343855</td>\n",
       "      <td>0.111238</td>\n",
       "      <td>-2.009986</td>\n",
       "      <td>0.495566</td>\n",
       "      <td>-1.028367</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DOEC.81.S30.001.fa</th>\n",
       "      <td>4.092462</td>\n",
       "      <td>-0.679445</td>\n",
       "      <td>3.202142</td>\n",
       "      <td>3.084641</td>\n",
       "      <td>3.229331</td>\n",
       "      <td>-0.325333</td>\n",
       "      <td>1.312362</td>\n",
       "      <td>3.325615</td>\n",
       "      <td>4.108773</td>\n",
       "      <td>1.343689</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.364446</td>\n",
       "      <td>-2.452479</td>\n",
       "      <td>-3.483327</td>\n",
       "      <td>0.449270</td>\n",
       "      <td>-1.440480</td>\n",
       "      <td>-0.737190</td>\n",
       "      <td>-1.007263</td>\n",
       "      <td>-2.681997</td>\n",
       "      <td>-0.129765</td>\n",
       "      <td>-0.393402</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EKJO.64.S32.001.fa</th>\n",
       "      <td>2.758383</td>\n",
       "      <td>1.899094</td>\n",
       "      <td>3.268204</td>\n",
       "      <td>1.459727</td>\n",
       "      <td>1.149572</td>\n",
       "      <td>2.914325</td>\n",
       "      <td>-1.524576</td>\n",
       "      <td>1.952459</td>\n",
       "      <td>1.117723</td>\n",
       "      <td>-0.472484</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.004701</td>\n",
       "      <td>2.080180</td>\n",
       "      <td>5.461192</td>\n",
       "      <td>-0.952057</td>\n",
       "      <td>-2.130712</td>\n",
       "      <td>-0.723798</td>\n",
       "      <td>-2.536177</td>\n",
       "      <td>-3.517006</td>\n",
       "      <td>-1.907569</td>\n",
       "      <td>-1.843030</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DIAL.15.S23.001.fa</th>\n",
       "      <td>4.740969</td>\n",
       "      <td>4.392127</td>\n",
       "      <td>1.668718</td>\n",
       "      <td>0.472244</td>\n",
       "      <td>-1.063350</td>\n",
       "      <td>-0.037714</td>\n",
       "      <td>1.751842</td>\n",
       "      <td>3.866595</td>\n",
       "      <td>-0.562924</td>\n",
       "      <td>2.009991</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.826564</td>\n",
       "      <td>-1.847991</td>\n",
       "      <td>-2.486585</td>\n",
       "      <td>0.842740</td>\n",
       "      <td>-1.756216</td>\n",
       "      <td>1.252917</td>\n",
       "      <td>-0.041654</td>\n",
       "      <td>-1.053558</td>\n",
       "      <td>0.031977</td>\n",
       "      <td>-0.745484</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>GIFI.21.S22.001.fa</th>\n",
       "      <td>3.184368</td>\n",
       "      <td>2.202335</td>\n",
       "      <td>-0.188357</td>\n",
       "      <td>-0.344703</td>\n",
       "      <td>-0.019281</td>\n",
       "      <td>2.197504</td>\n",
       "      <td>-1.151795</td>\n",
       "      <td>2.699819</td>\n",
       "      <td>1.567684</td>\n",
       "      <td>0.284930</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.117893</td>\n",
       "      <td>3.923379</td>\n",
       "      <td>5.297858</td>\n",
       "      <td>0.153738</td>\n",
       "      <td>-2.727331</td>\n",
       "      <td>-1.341037</td>\n",
       "      <td>-3.132796</td>\n",
       "      <td>-2.909653</td>\n",
       "      <td>-2.573181</td>\n",
       "      <td>-3.132796</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TAOT.65.S24.001.fa</th>\n",
       "      <td>2.689188</td>\n",
       "      <td>1.152850</td>\n",
       "      <td>0.797345</td>\n",
       "      <td>-1.012114</td>\n",
       "      <td>-2.398408</td>\n",
       "      <td>-0.606649</td>\n",
       "      <td>-1.571730</td>\n",
       "      <td>3.131588</td>\n",
       "      <td>-0.432295</td>\n",
       "      <td>-0.878582</td>\n",
       "      <td>...</td>\n",
       "      <td>-2.552559</td>\n",
       "      <td>3.106923</td>\n",
       "      <td>5.611904</td>\n",
       "      <td>0.568337</td>\n",
       "      <td>-2.958024</td>\n",
       "      <td>-1.399879</td>\n",
       "      <td>-2.734880</td>\n",
       "      <td>-2.264877</td>\n",
       "      <td>-2.552559</td>\n",
       "      <td>-1.208824</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DEFI.14.S6.001.fa</th>\n",
       "      <td>3.475730</td>\n",
       "      <td>5.737301</td>\n",
       "      <td>3.139103</td>\n",
       "      <td>5.188512</td>\n",
       "      <td>0.212743</td>\n",
       "      <td>-1.604646</td>\n",
       "      <td>0.610018</td>\n",
       "      <td>3.058043</td>\n",
       "      <td>-0.740189</td>\n",
       "      <td>2.086746</td>\n",
       "      <td>...</td>\n",
       "      <td>-3.038194</td>\n",
       "      <td>1.291514</td>\n",
       "      <td>-0.052186</td>\n",
       "      <td>-0.119842</td>\n",
       "      <td>1.011532</td>\n",
       "      <td>-1.295560</td>\n",
       "      <td>-2.430593</td>\n",
       "      <td>-3.090619</td>\n",
       "      <td>-1.534714</td>\n",
       "      <td>-1.908701</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MUNA.98.S28.001.fa</th>\n",
       "      <td>2.218994</td>\n",
       "      <td>-1.491165</td>\n",
       "      <td>-1.683421</td>\n",
       "      <td>-2.960675</td>\n",
       "      <td>-3.069716</td>\n",
       "      <td>-2.292634</td>\n",
       "      <td>-0.164856</td>\n",
       "      <td>4.237870</td>\n",
       "      <td>-0.495943</td>\n",
       "      <td>1.740706</td>\n",
       "      <td>...</td>\n",
       "      <td>0.302494</td>\n",
       "      <td>0.470776</td>\n",
       "      <td>-0.653966</td>\n",
       "      <td>-0.514026</td>\n",
       "      <td>-0.816727</td>\n",
       "      <td>-0.916989</td>\n",
       "      <td>-1.873143</td>\n",
       "      <td>-4.406562</td>\n",
       "      <td>-1.715021</td>\n",
       "      <td>-0.154564</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>ROAL.75.S33.001.fa</th>\n",
       "      <td>2.617716</td>\n",
       "      <td>0.620977</td>\n",
       "      <td>1.251969</td>\n",
       "      <td>0.998512</td>\n",
       "      <td>1.236305</td>\n",
       "      <td>-0.857894</td>\n",
       "      <td>-0.566466</td>\n",
       "      <td>2.787953</td>\n",
       "      <td>-0.119041</td>\n",
       "      <td>0.706430</td>\n",
       "      <td>...</td>\n",
       "      <td>2.554822</td>\n",
       "      <td>1.559750</td>\n",
       "      <td>0.630927</td>\n",
       "      <td>0.524060</td>\n",
       "      <td>-2.359089</td>\n",
       "      <td>-2.246501</td>\n",
       "      <td>-3.651487</td>\n",
       "      <td>-2.515834</td>\n",
       "      <td>-1.401189</td>\n",
       "      <td>-2.638173</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MERI.28.S20.001.fa</th>\n",
       "      <td>1.645391</td>\n",
       "      <td>-2.516686</td>\n",
       "      <td>6.880540</td>\n",
       "      <td>5.298956</td>\n",
       "      <td>5.544301</td>\n",
       "      <td>-2.042753</td>\n",
       "      <td>0.192413</td>\n",
       "      <td>2.897733</td>\n",
       "      <td>0.579711</td>\n",
       "      <td>-1.315868</td>\n",
       "      <td>...</td>\n",
       "      <td>0.439290</td>\n",
       "      <td>2.206092</td>\n",
       "      <td>-0.342471</td>\n",
       "      <td>-0.432335</td>\n",
       "      <td>-3.599673</td>\n",
       "      <td>-0.984146</td>\n",
       "      <td>-0.108111</td>\n",
       "      <td>-3.217738</td>\n",
       "      <td>-3.448850</td>\n",
       "      <td>-2.688044</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BUMA.05.S4.001.fa</th>\n",
       "      <td>5.184688</td>\n",
       "      <td>0.726147</td>\n",
       "      <td>-1.573127</td>\n",
       "      <td>-1.091536</td>\n",
       "      <td>-2.000538</td>\n",
       "      <td>-1.530629</td>\n",
       "      <td>0.112467</td>\n",
       "      <td>4.452070</td>\n",
       "      <td>-0.874077</td>\n",
       "      <td>1.570038</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.666852</td>\n",
       "      <td>-1.504957</td>\n",
       "      <td>-2.146903</td>\n",
       "      <td>-0.352957</td>\n",
       "      <td>-1.530629</td>\n",
       "      <td>-0.746264</td>\n",
       "      <td>-1.368205</td>\n",
       "      <td>-2.212961</td>\n",
       "      <td>-1.443736</td>\n",
       "      <td>-1.600636</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>MOEM.48.S16.001.fa</th>\n",
       "      <td>5.631256</td>\n",
       "      <td>1.947917</td>\n",
       "      <td>-1.263681</td>\n",
       "      <td>-1.753768</td>\n",
       "      <td>-2.266186</td>\n",
       "      <td>0.196565</td>\n",
       "      <td>2.102101</td>\n",
       "      <td>4.604422</td>\n",
       "      <td>2.247472</td>\n",
       "      <td>3.520333</td>\n",
       "      <td>...</td>\n",
       "      <td>-4.459434</td>\n",
       "      <td>-5.226689</td>\n",
       "      <td>-7.098492</td>\n",
       "      <td>0.976780</td>\n",
       "      <td>2.587275</td>\n",
       "      <td>1.486921</td>\n",
       "      <td>-1.131064</td>\n",
       "      <td>-0.922624</td>\n",
       "      <td>0.683690</td>\n",
       "      <td>0.319989</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PAMA.46.S18.001.fa</th>\n",
       "      <td>2.826640</td>\n",
       "      <td>1.425925</td>\n",
       "      <td>0.920859</td>\n",
       "      <td>-0.381243</td>\n",
       "      <td>-0.867528</td>\n",
       "      <td>-1.068199</td>\n",
       "      <td>-2.690059</td>\n",
       "      <td>2.445372</td>\n",
       "      <td>-0.587451</td>\n",
       "      <td>-0.174381</td>\n",
       "      <td>...</td>\n",
       "      <td>-1.335513</td>\n",
       "      <td>1.219461</td>\n",
       "      <td>-0.602835</td>\n",
       "      <td>-0.996740</td>\n",
       "      <td>-0.908771</td>\n",
       "      <td>-0.037698</td>\n",
       "      <td>-2.690059</td>\n",
       "      <td>-1.335513</td>\n",
       "      <td>-2.977741</td>\n",
       "      <td>-2.327154</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>VUGE.37.S14.001.fa</th>\n",
       "      <td>4.700142</td>\n",
       "      <td>-1.786927</td>\n",
       "      <td>4.158713</td>\n",
       "      <td>1.681301</td>\n",
       "      <td>0.816839</td>\n",
       "      <td>-1.146786</td>\n",
       "      <td>0.389677</td>\n",
       "      <td>2.790468</td>\n",
       "      <td>2.500457</td>\n",
       "      <td>0.155391</td>\n",
       "      <td>...</td>\n",
       "      <td>0.423770</td>\n",
       "      <td>-1.154333</td>\n",
       "      <td>-3.232563</td>\n",
       "      <td>0.427640</td>\n",
       "      <td>-0.313027</td>\n",
       "      <td>0.374880</td>\n",
       "      <td>-0.862795</td>\n",
       "      <td>-2.197253</td>\n",
       "      <td>-1.232786</td>\n",
       "      <td>-0.935806</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TOSA.76.S25.001.fa</th>\n",
       "      <td>2.213781</td>\n",
       "      <td>6.897451</td>\n",
       "      <td>-1.166429</td>\n",
       "      <td>-0.946432</td>\n",
       "      <td>-1.813533</td>\n",
       "      <td>-1.793730</td>\n",
       "      <td>1.110166</td>\n",
       "      <td>2.397261</td>\n",
       "      <td>2.802775</td>\n",
       "      <td>0.494040</td>\n",
       "      <td>...</td>\n",
       "      <td>0.618325</td>\n",
       "      <td>1.556518</td>\n",
       "      <td>1.321092</td>\n",
       "      <td>0.577063</td>\n",
       "      <td>-2.275568</td>\n",
       "      <td>-0.003606</td>\n",
       "      <td>-2.781117</td>\n",
       "      <td>-3.086498</td>\n",
       "      <td>-2.393351</td>\n",
       "      <td>-2.393351</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>DEIS.70.S26.001.fa</th>\n",
       "      <td>2.781823</td>\n",
       "      <td>1.181892</td>\n",
       "      <td>1.358823</td>\n",
       "      <td>0.169754</td>\n",
       "      <td>-0.047490</td>\n",
       "      <td>4.855276</td>\n",
       "      <td>-1.869103</td>\n",
       "      <td>2.909639</td>\n",
       "      <td>-1.153483</td>\n",
       "      <td>-1.433785</td>\n",
       "      <td>...</td>\n",
       "      <td>-0.406268</td>\n",
       "      <td>3.801718</td>\n",
       "      <td>10.340782</td>\n",
       "      <td>1.103640</td>\n",
       "      <td>-3.573851</td>\n",
       "      <td>0.557308</td>\n",
       "      <td>-3.350707</td>\n",
       "      <td>-1.246573</td>\n",
       "      <td>0.276297</td>\n",
       "      <td>-3.350707</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>35 rows × 181 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Phocaeicola;s__Phocaeicola_vulgatus  \\\n",
       "pos.dna.02.S35.001.fa                                          12.564805                                                                      \n",
       "CAPMA.83.S12.001.fa                                             1.502473                                                                      \n",
       "FOGI.63.S13.001.fa                                              3.594296                                                                      \n",
       "COCL.11.S7.001.fa                                               1.536088                                                                      \n",
       "WACH.38.S1.001.fa                                               4.240597                                                                      \n",
       "PRPI.42.S31.001.fa                                              3.169797                                                                      \n",
       "RODE.80.S8.001.fa                                               5.918568                                                                      \n",
       "pos.dna.01.S17.001.fa                                          10.758325                                                                      \n",
       "MEEN.93.S11.001.fa                                              4.110142                                                                      \n",
       "DETH.41.S27.001.fa                                              3.183976                                                                      \n",
       "TECA.45.S15.001.fa                                              3.134289                                                                      \n",
       "BODA.52.S9.001.fa                                               4.091882                                                                      \n",
       "CHMA.61.S3.001.fa                                               6.751274                                                                      \n",
       "ASTER.02.S34.001.fa                                             7.511852                                                                      \n",
       "BIRO.74.S19.001.fa                                              6.603156                                                                      \n",
       "CAAU.07.S2.001.fa                                               5.567414                                                                      \n",
       "DIGMA.71.S21.001.fa                                             5.374936                                                                      \n",
       "FAFE.16.S10.001.fa                                              5.584302                                                                      \n",
       "OPSE.53.S5.001.fa                                               7.120223                                                                      \n",
       "LAES.91.S29.001.fa                                              3.183280                                                                      \n",
       "DOEC.81.S30.001.fa                                              4.092462                                                                      \n",
       "EKJO.64.S32.001.fa                                              2.758383                                                                      \n",
       "DIAL.15.S23.001.fa                                              4.740969                                                                      \n",
       "GIFI.21.S22.001.fa                                              3.184368                                                                      \n",
       "TAOT.65.S24.001.fa                                              2.689188                                                                      \n",
       "DEFI.14.S6.001.fa                                               3.475730                                                                      \n",
       "MUNA.98.S28.001.fa                                              2.218994                                                                      \n",
       "ROAL.75.S33.001.fa                                              2.617716                                                                      \n",
       "MERI.28.S20.001.fa                                              1.645391                                                                      \n",
       "BUMA.05.S4.001.fa                                               5.184688                                                                      \n",
       "MOEM.48.S16.001.fa                                              5.631256                                                                      \n",
       "PAMA.46.S18.001.fa                                              2.826640                                                                      \n",
       "VUGE.37.S14.001.fa                                              4.700142                                                                      \n",
       "TOSA.76.S25.001.fa                                              2.213781                                                                      \n",
       "DEIS.70.S26.001.fa                                              2.781823                                                                      \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_fragilis  \\\n",
       "pos.dna.02.S35.001.fa                                          12.596202                                                                      \n",
       "CAPMA.83.S12.001.fa                                             0.108900                                                                      \n",
       "FOGI.63.S13.001.fa                                             -0.611240                                                                      \n",
       "COCL.11.S7.001.fa                                               1.074962                                                                      \n",
       "WACH.38.S1.001.fa                                              -0.421227                                                                      \n",
       "PRPI.42.S31.001.fa                                             -1.438640                                                                      \n",
       "RODE.80.S8.001.fa                                              -1.343871                                                                      \n",
       "pos.dna.01.S17.001.fa                                          10.466699                                                                      \n",
       "MEEN.93.S11.001.fa                                              6.091029                                                                      \n",
       "DETH.41.S27.001.fa                                             -3.761992                                                                      \n",
       "TECA.45.S15.001.fa                                              4.631672                                                                      \n",
       "BODA.52.S9.001.fa                                               1.095134                                                                      \n",
       "CHMA.61.S3.001.fa                                               1.924646                                                                      \n",
       "ASTER.02.S34.001.fa                                            -0.131282                                                                      \n",
       "BIRO.74.S19.001.fa                                              4.111350                                                                      \n",
       "CAAU.07.S2.001.fa                                              -1.544217                                                                      \n",
       "DIGMA.71.S21.001.fa                                             0.454759                                                                      \n",
       "FAFE.16.S10.001.fa                                             -0.067135                                                                      \n",
       "OPSE.53.S5.001.fa                                               1.635535                                                                      \n",
       "LAES.91.S29.001.fa                                             -0.423308                                                                      \n",
       "DOEC.81.S30.001.fa                                             -0.679445                                                                      \n",
       "EKJO.64.S32.001.fa                                              1.899094                                                                      \n",
       "DIAL.15.S23.001.fa                                              4.392127                                                                      \n",
       "GIFI.21.S22.001.fa                                              2.202335                                                                      \n",
       "TAOT.65.S24.001.fa                                              1.152850                                                                      \n",
       "DEFI.14.S6.001.fa                                               5.737301                                                                      \n",
       "MUNA.98.S28.001.fa                                             -1.491165                                                                      \n",
       "ROAL.75.S33.001.fa                                              0.620977                                                                      \n",
       "MERI.28.S20.001.fa                                             -2.516686                                                                      \n",
       "BUMA.05.S4.001.fa                                               0.726147                                                                      \n",
       "MOEM.48.S16.001.fa                                              1.947917                                                                      \n",
       "PAMA.46.S18.001.fa                                              1.425925                                                                      \n",
       "VUGE.37.S14.001.fa                                             -1.786927                                                                      \n",
       "TOSA.76.S25.001.fa                                              6.897451                                                                      \n",
       "DEIS.70.S26.001.fa                                              1.181892                                                                      \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli  \\\n",
       "pos.dna.02.S35.001.fa                                          10.432703                                                                                   \n",
       "CAPMA.83.S12.001.fa                                             8.026002                                                                                   \n",
       "FOGI.63.S13.001.fa                                              0.322862                                                                                   \n",
       "COCL.11.S7.001.fa                                               4.663100                                                                                   \n",
       "WACH.38.S1.001.fa                                               2.999156                                                                                   \n",
       "PRPI.42.S31.001.fa                                              7.172793                                                                                   \n",
       "RODE.80.S8.001.fa                                               4.211689                                                                                   \n",
       "pos.dna.01.S17.001.fa                                           8.427199                                                                                   \n",
       "MEEN.93.S11.001.fa                                              4.004449                                                                                   \n",
       "DETH.41.S27.001.fa                                             -0.928224                                                                                   \n",
       "TECA.45.S15.001.fa                                              2.388486                                                                                   \n",
       "BODA.52.S9.001.fa                                               0.574600                                                                                   \n",
       "CHMA.61.S3.001.fa                                               3.074376                                                                                   \n",
       "ASTER.02.S34.001.fa                                            -3.570534                                                                                   \n",
       "BIRO.74.S19.001.fa                                             -1.495615                                                                                   \n",
       "CAAU.07.S2.001.fa                                              -4.652267                                                                                   \n",
       "DIGMA.71.S21.001.fa                                            -1.006752                                                                                   \n",
       "FAFE.16.S10.001.fa                                             -1.416951                                                                                   \n",
       "OPSE.53.S5.001.fa                                               1.917310                                                                                   \n",
       "LAES.91.S29.001.fa                                              5.162297                                                                                   \n",
       "DOEC.81.S30.001.fa                                              3.202142                                                                                   \n",
       "EKJO.64.S32.001.fa                                              3.268204                                                                                   \n",
       "DIAL.15.S23.001.fa                                              1.668718                                                                                   \n",
       "GIFI.21.S22.001.fa                                             -0.188357                                                                                   \n",
       "TAOT.65.S24.001.fa                                              0.797345                                                                                   \n",
       "DEFI.14.S6.001.fa                                               3.139103                                                                                   \n",
       "MUNA.98.S28.001.fa                                             -1.683421                                                                                   \n",
       "ROAL.75.S33.001.fa                                              1.251969                                                                                   \n",
       "MERI.28.S20.001.fa                                              6.880540                                                                                   \n",
       "BUMA.05.S4.001.fa                                              -1.573127                                                                                   \n",
       "MOEM.48.S16.001.fa                                             -1.263681                                                                                   \n",
       "PAMA.46.S18.001.fa                                              0.920859                                                                                   \n",
       "VUGE.37.S14.001.fa                                              4.158713                                                                                   \n",
       "TOSA.76.S25.001.fa                                             -1.166429                                                                                   \n",
       "DEIS.70.S26.001.fa                                              1.358823                                                                                   \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_coli_D  \\\n",
       "pos.dna.02.S35.001.fa                                          10.382420                                                                                     \n",
       "CAPMA.83.S12.001.fa                                             6.006419                                                                                     \n",
       "FOGI.63.S13.001.fa                                             -0.505445                                                                                     \n",
       "COCL.11.S7.001.fa                                               3.188711                                                                                     \n",
       "WACH.38.S1.001.fa                                               1.341303                                                                                     \n",
       "PRPI.42.S31.001.fa                                              4.708191                                                                                     \n",
       "RODE.80.S8.001.fa                                               4.242758                                                                                     \n",
       "pos.dna.01.S17.001.fa                                           8.411471                                                                                     \n",
       "MEEN.93.S11.001.fa                                              3.807656                                                                                     \n",
       "DETH.41.S27.001.fa                                             -2.920100                                                                                     \n",
       "TECA.45.S15.001.fa                                              3.961840                                                                                     \n",
       "BODA.52.S9.001.fa                                              -0.936630                                                                                     \n",
       "CHMA.61.S3.001.fa                                               2.536282                                                                                     \n",
       "ASTER.02.S34.001.fa                                            -3.458830                                                                                     \n",
       "BIRO.74.S19.001.fa                                             -1.615759                                                                                     \n",
       "CAAU.07.S2.001.fa                                              -5.025247                                                                                     \n",
       "DIGMA.71.S21.001.fa                                             0.006649                                                                                     \n",
       "FAFE.16.S10.001.fa                                             -2.181225                                                                                     \n",
       "OPSE.53.S5.001.fa                                              -0.146372                                                                                     \n",
       "LAES.91.S29.001.fa                                              3.618158                                                                                     \n",
       "DOEC.81.S30.001.fa                                              3.084641                                                                                     \n",
       "EKJO.64.S32.001.fa                                              1.459727                                                                                     \n",
       "DIAL.15.S23.001.fa                                              0.472244                                                                                     \n",
       "GIFI.21.S22.001.fa                                             -0.344703                                                                                     \n",
       "TAOT.65.S24.001.fa                                             -1.012114                                                                                     \n",
       "DEFI.14.S6.001.fa                                               5.188512                                                                                     \n",
       "MUNA.98.S28.001.fa                                             -2.960675                                                                                     \n",
       "ROAL.75.S33.001.fa                                              0.998512                                                                                     \n",
       "MERI.28.S20.001.fa                                              5.298956                                                                                     \n",
       "BUMA.05.S4.001.fa                                              -1.091536                                                                                     \n",
       "MOEM.48.S16.001.fa                                             -1.753768                                                                                     \n",
       "PAMA.46.S18.001.fa                                             -0.381243                                                                                     \n",
       "VUGE.37.S14.001.fa                                              1.681301                                                                                     \n",
       "TOSA.76.S25.001.fa                                             -0.946432                                                                                     \n",
       "DEIS.70.S26.001.fa                                              0.169754                                                                                     \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Proteobacteria;c__Gammaproteobacteria;o__Enterobacterales;f__Enterobacteriaceae;g__Escherichia;s__Escherichia_flexneri  \\\n",
       "pos.dna.02.S35.001.fa                                          10.717464                                                                                       \n",
       "CAPMA.83.S12.001.fa                                             4.688253                                                                                       \n",
       "FOGI.63.S13.001.fa                                             -2.583584                                                                                       \n",
       "COCL.11.S7.001.fa                                               1.774582                                                                                       \n",
       "WACH.38.S1.001.fa                                               0.023436                                                                                       \n",
       "PRPI.42.S31.001.fa                                              3.824696                                                                                       \n",
       "RODE.80.S8.001.fa                                               4.215763                                                                                       \n",
       "pos.dna.01.S17.001.fa                                           8.758813                                                                                       \n",
       "MEEN.93.S11.001.fa                                              2.422215                                                                                       \n",
       "DETH.41.S27.001.fa                                             -3.462587                                                                                       \n",
       "TECA.45.S15.001.fa                                             -0.605582                                                                                       \n",
       "BODA.52.S9.001.fa                                              -1.537898                                                                                       \n",
       "CHMA.61.S3.001.fa                                               1.136575                                                                                       \n",
       "ASTER.02.S34.001.fa                                            -3.434732                                                                                       \n",
       "BIRO.74.S19.001.fa                                             -2.541918                                                                                       \n",
       "CAAU.07.S2.001.fa                                              -6.750757                                                                                       \n",
       "DIGMA.71.S21.001.fa                                            -1.588852                                                                                       \n",
       "FAFE.16.S10.001.fa                                             -2.503694                                                                                       \n",
       "OPSE.53.S5.001.fa                                              -1.667841                                                                                       \n",
       "LAES.91.S29.001.fa                                              2.340659                                                                                       \n",
       "DOEC.81.S30.001.fa                                              3.229331                                                                                       \n",
       "EKJO.64.S32.001.fa                                              1.149572                                                                                       \n",
       "DIAL.15.S23.001.fa                                             -1.063350                                                                                       \n",
       "GIFI.21.S22.001.fa                                             -0.019281                                                                                       \n",
       "TAOT.65.S24.001.fa                                             -2.398408                                                                                       \n",
       "DEFI.14.S6.001.fa                                               0.212743                                                                                       \n",
       "MUNA.98.S28.001.fa                                             -3.069716                                                                                       \n",
       "ROAL.75.S33.001.fa                                              1.236305                                                                                       \n",
       "MERI.28.S20.001.fa                                              5.544301                                                                                       \n",
       "BUMA.05.S4.001.fa                                              -2.000538                                                                                       \n",
       "MOEM.48.S16.001.fa                                             -2.266186                                                                                       \n",
       "PAMA.46.S18.001.fa                                             -0.867528                                                                                       \n",
       "VUGE.37.S14.001.fa                                              0.816839                                                                                       \n",
       "TOSA.76.S25.001.fa                                             -1.813533                                                                                       \n",
       "DEIS.70.S26.001.fa                                             -0.047490                                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_bouchesdurhonensis  \\\n",
       "pos.dna.02.S35.001.fa                                          -0.529343                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -3.136293                                                                                \n",
       "FOGI.63.S13.001.fa                                             -0.913521                                                                                \n",
       "COCL.11.S7.001.fa                                              -2.081677                                                                                \n",
       "WACH.38.S1.001.fa                                              -1.574201                                                                                \n",
       "PRPI.42.S31.001.fa                                             -0.532085                                                                                \n",
       "RODE.80.S8.001.fa                                              -0.235021                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -1.616304                                                                                \n",
       "MEEN.93.S11.001.fa                                             -0.461099                                                                                \n",
       "DETH.41.S27.001.fa                                              0.307518                                                                                \n",
       "TECA.45.S15.001.fa                                             -0.851349                                                                                \n",
       "BODA.52.S9.001.fa                                               0.710902                                                                                \n",
       "CHMA.61.S3.001.fa                                               1.722341                                                                                \n",
       "ASTER.02.S34.001.fa                                             0.951255                                                                                \n",
       "BIRO.74.S19.001.fa                                             -0.057356                                                                                \n",
       "CAAU.07.S2.001.fa                                              -0.619196                                                                                \n",
       "DIGMA.71.S21.001.fa                                            -0.308543                                                                                \n",
       "FAFE.16.S10.001.fa                                              1.577985                                                                                \n",
       "OPSE.53.S5.001.fa                                               1.399707                                                                                \n",
       "LAES.91.S29.001.fa                                             -0.533104                                                                                \n",
       "DOEC.81.S30.001.fa                                             -0.325333                                                                                \n",
       "EKJO.64.S32.001.fa                                              2.914325                                                                                \n",
       "DIAL.15.S23.001.fa                                             -0.037714                                                                                \n",
       "GIFI.21.S22.001.fa                                              2.197504                                                                                \n",
       "TAOT.65.S24.001.fa                                             -0.606649                                                                                \n",
       "DEFI.14.S6.001.fa                                              -1.604646                                                                                \n",
       "MUNA.98.S28.001.fa                                             -2.292634                                                                                \n",
       "ROAL.75.S33.001.fa                                             -0.857894                                                                                \n",
       "MERI.28.S20.001.fa                                             -2.042753                                                                                \n",
       "BUMA.05.S4.001.fa                                              -1.530629                                                                                \n",
       "MOEM.48.S16.001.fa                                              0.196565                                                                                \n",
       "PAMA.46.S18.001.fa                                             -1.068199                                                                                \n",
       "VUGE.37.S14.001.fa                                             -1.146786                                                                                \n",
       "TOSA.76.S25.001.fa                                             -1.793730                                                                                \n",
       "DEIS.70.S26.001.fa                                              4.855276                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Tannerellaceae;g__Parabacteroides;s__Parabacteroides_distasonis_A  \\\n",
       "pos.dna.02.S35.001.fa                                          -0.683494                                                                                  \n",
       "CAPMA.83.S12.001.fa                                             0.814951                                                                                  \n",
       "FOGI.63.S13.001.fa                                             -0.274768                                                                                  \n",
       "COCL.11.S7.001.fa                                              -2.133863                                                                                  \n",
       "WACH.38.S1.001.fa                                              -1.208627                                                                                  \n",
       "PRPI.42.S31.001.fa                                              0.522833                                                                                  \n",
       "RODE.80.S8.001.fa                                               1.030012                                                                                  \n",
       "pos.dna.01.S17.001.fa                                          -0.923157                                                                                  \n",
       "MEEN.93.S11.001.fa                                              0.197382                                                                                  \n",
       "DETH.41.S27.001.fa                                             -3.401551                                                                                  \n",
       "TECA.45.S15.001.fa                                              0.515884                                                                                  \n",
       "BODA.52.S9.001.fa                                              -0.204385                                                                                  \n",
       "CHMA.61.S3.001.fa                                               1.558442                                                                                  \n",
       "ASTER.02.S34.001.fa                                             1.994439                                                                                  \n",
       "BIRO.74.S19.001.fa                                             -0.108232                                                                                  \n",
       "CAAU.07.S2.001.fa                                               1.828631                                                                                  \n",
       "DIGMA.71.S21.001.fa                                            -0.375151                                                                                  \n",
       "FAFE.16.S10.001.fa                                              0.749284                                                                                  \n",
       "OPSE.53.S5.001.fa                                               3.893911                                                                                  \n",
       "LAES.91.S29.001.fa                                             -2.047337                                                                                  \n",
       "DOEC.81.S30.001.fa                                              1.312362                                                                                  \n",
       "EKJO.64.S32.001.fa                                             -1.524576                                                                                  \n",
       "DIAL.15.S23.001.fa                                              1.751842                                                                                  \n",
       "GIFI.21.S22.001.fa                                             -1.151795                                                                                  \n",
       "TAOT.65.S24.001.fa                                             -1.571730                                                                                  \n",
       "DEFI.14.S6.001.fa                                               0.610018                                                                                  \n",
       "MUNA.98.S28.001.fa                                             -0.164856                                                                                  \n",
       "ROAL.75.S33.001.fa                                             -0.566466                                                                                  \n",
       "MERI.28.S20.001.fa                                              0.192413                                                                                  \n",
       "BUMA.05.S4.001.fa                                               0.112467                                                                                  \n",
       "MOEM.48.S16.001.fa                                              2.102101                                                                                  \n",
       "PAMA.46.S18.001.fa                                             -2.690059                                                                                  \n",
       "VUGE.37.S14.001.fa                                              0.389677                                                                                  \n",
       "TOSA.76.S25.001.fa                                              1.110166                                                                                  \n",
       "DEIS.70.S26.001.fa                                             -1.869103                                                                                  \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_uniformis  \\\n",
       "pos.dna.02.S35.001.fa                                           5.178716                                                                       \n",
       "CAPMA.83.S12.001.fa                                             1.198817                                                                       \n",
       "FOGI.63.S13.001.fa                                              2.714890                                                                       \n",
       "COCL.11.S7.001.fa                                               1.997295                                                                       \n",
       "WACH.38.S1.001.fa                                               2.685130                                                                       \n",
       "PRPI.42.S31.001.fa                                              2.376027                                                                       \n",
       "RODE.80.S8.001.fa                                               2.403992                                                                       \n",
       "pos.dna.01.S17.001.fa                                           3.486100                                                                       \n",
       "MEEN.93.S11.001.fa                                              3.913378                                                                       \n",
       "DETH.41.S27.001.fa                                              1.633435                                                                       \n",
       "TECA.45.S15.001.fa                                              3.626052                                                                       \n",
       "BODA.52.S9.001.fa                                               2.611502                                                                       \n",
       "CHMA.61.S3.001.fa                                               6.117136                                                                       \n",
       "ASTER.02.S34.001.fa                                             3.549830                                                                       \n",
       "BIRO.74.S19.001.fa                                              4.973692                                                                       \n",
       "CAAU.07.S2.001.fa                                               5.706051                                                                       \n",
       "DIGMA.71.S21.001.fa                                             6.540831                                                                       \n",
       "FAFE.16.S10.001.fa                                              5.106302                                                                       \n",
       "OPSE.53.S5.001.fa                                               4.110070                                                                       \n",
       "LAES.91.S29.001.fa                                              4.861458                                                                       \n",
       "DOEC.81.S30.001.fa                                              3.325615                                                                       \n",
       "EKJO.64.S32.001.fa                                              1.952459                                                                       \n",
       "DIAL.15.S23.001.fa                                              3.866595                                                                       \n",
       "GIFI.21.S22.001.fa                                              2.699819                                                                       \n",
       "TAOT.65.S24.001.fa                                              3.131588                                                                       \n",
       "DEFI.14.S6.001.fa                                               3.058043                                                                       \n",
       "MUNA.98.S28.001.fa                                              4.237870                                                                       \n",
       "ROAL.75.S33.001.fa                                              2.787953                                                                       \n",
       "MERI.28.S20.001.fa                                              2.897733                                                                       \n",
       "BUMA.05.S4.001.fa                                               4.452070                                                                       \n",
       "MOEM.48.S16.001.fa                                              4.604422                                                                       \n",
       "PAMA.46.S18.001.fa                                              2.445372                                                                       \n",
       "VUGE.37.S14.001.fa                                              2.790468                                                                       \n",
       "TOSA.76.S25.001.fa                                              2.397261                                                                       \n",
       "DEIS.70.S26.001.fa                                              2.909639                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_faecis  \\\n",
       "pos.dna.02.S35.001.fa                                           1.714401                                                                    \n",
       "CAPMA.83.S12.001.fa                                             3.197282                                                                    \n",
       "FOGI.63.S13.001.fa                                              0.887994                                                                    \n",
       "COCL.11.S7.001.fa                                               2.537797                                                                    \n",
       "WACH.38.S1.001.fa                                               0.050759                                                                    \n",
       "PRPI.42.S31.001.fa                                              1.803052                                                                    \n",
       "RODE.80.S8.001.fa                                               1.124451                                                                    \n",
       "pos.dna.01.S17.001.fa                                           2.323658                                                                    \n",
       "MEEN.93.S11.001.fa                                              1.410490                                                                    \n",
       "DETH.41.S27.001.fa                                              0.197562                                                                    \n",
       "TECA.45.S15.001.fa                                             -0.435956                                                                    \n",
       "BODA.52.S9.001.fa                                               1.548782                                                                    \n",
       "CHMA.61.S3.001.fa                                               3.203964                                                                    \n",
       "ASTER.02.S34.001.fa                                             2.011697                                                                    \n",
       "BIRO.74.S19.001.fa                                              1.036538                                                                    \n",
       "CAAU.07.S2.001.fa                                               1.065846                                                                    \n",
       "DIGMA.71.S21.001.fa                                             0.876363                                                                    \n",
       "FAFE.16.S10.001.fa                                              6.007668                                                                    \n",
       "OPSE.53.S5.001.fa                                               1.131665                                                                    \n",
       "LAES.91.S29.001.fa                                             -0.527922                                                                    \n",
       "DOEC.81.S30.001.fa                                              4.108773                                                                    \n",
       "EKJO.64.S32.001.fa                                              1.117723                                                                    \n",
       "DIAL.15.S23.001.fa                                             -0.562924                                                                    \n",
       "GIFI.21.S22.001.fa                                              1.567684                                                                    \n",
       "TAOT.65.S24.001.fa                                             -0.432295                                                                    \n",
       "DEFI.14.S6.001.fa                                              -0.740189                                                                    \n",
       "MUNA.98.S28.001.fa                                             -0.495943                                                                    \n",
       "ROAL.75.S33.001.fa                                             -0.119041                                                                    \n",
       "MERI.28.S20.001.fa                                              0.579711                                                                    \n",
       "BUMA.05.S4.001.fa                                              -0.874077                                                                    \n",
       "MOEM.48.S16.001.fa                                              2.247472                                                                    \n",
       "PAMA.46.S18.001.fa                                             -0.587451                                                                    \n",
       "VUGE.37.S14.001.fa                                              2.500457                                                                    \n",
       "TOSA.76.S25.001.fa                                              2.802775                                                                    \n",
       "DEIS.70.S26.001.fa                                             -1.153483                                                                    \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_thetaiotaomicron  \\\n",
       "pos.dna.02.S35.001.fa                                           2.002083                                                                              \n",
       "CAPMA.83.S12.001.fa                                             2.198838                                                                              \n",
       "FOGI.63.S13.001.fa                                              2.404863                                                                              \n",
       "COCL.11.S7.001.fa                                               4.726931                                                                              \n",
       "WACH.38.S1.001.fa                                               1.814724                                                                              \n",
       "PRPI.42.S31.001.fa                                              1.172110                                                                              \n",
       "RODE.80.S8.001.fa                                               2.924393                                                                              \n",
       "pos.dna.01.S17.001.fa                                           1.673341                                                                              \n",
       "MEEN.93.S11.001.fa                                              4.790671                                                                              \n",
       "DETH.41.S27.001.fa                                             -0.076657                                                                              \n",
       "TECA.45.S15.001.fa                                             -0.125998                                                                              \n",
       "BODA.52.S9.001.fa                                              -0.011895                                                                              \n",
       "CHMA.61.S3.001.fa                                               5.417908                                                                              \n",
       "ASTER.02.S34.001.fa                                             4.680516                                                                              \n",
       "BIRO.74.S19.001.fa                                              4.116220                                                                              \n",
       "CAAU.07.S2.001.fa                                               3.684255                                                                              \n",
       "DIGMA.71.S21.001.fa                                             1.213984                                                                              \n",
       "FAFE.16.S10.001.fa                                              3.278542                                                                              \n",
       "OPSE.53.S5.001.fa                                               4.004854                                                                              \n",
       "LAES.91.S29.001.fa                                              1.678388                                                                              \n",
       "DOEC.81.S30.001.fa                                              1.343689                                                                              \n",
       "EKJO.64.S32.001.fa                                             -0.472484                                                                              \n",
       "DIAL.15.S23.001.fa                                              2.009991                                                                              \n",
       "GIFI.21.S22.001.fa                                              0.284930                                                                              \n",
       "TAOT.65.S24.001.fa                                             -0.878582                                                                              \n",
       "DEFI.14.S6.001.fa                                               2.086746                                                                              \n",
       "MUNA.98.S28.001.fa                                              1.740706                                                                              \n",
       "ROAL.75.S33.001.fa                                              0.706430                                                                              \n",
       "MERI.28.S20.001.fa                                             -1.315868                                                                              \n",
       "BUMA.05.S4.001.fa                                               1.570038                                                                              \n",
       "MOEM.48.S16.001.fa                                              3.520333                                                                              \n",
       "PAMA.46.S18.001.fa                                             -0.174381                                                                              \n",
       "VUGE.37.S14.001.fa                                              0.155391                                                                              \n",
       "TOSA.76.S25.001.fa                                              0.494040                                                                              \n",
       "DEIS.70.S26.001.fa                                             -1.433785                                                                              \n",
       "\n",
       "#OTU ID                ...  \\\n",
       "pos.dna.02.S35.001.fa  ...   \n",
       "CAPMA.83.S12.001.fa    ...   \n",
       "FOGI.63.S13.001.fa     ...   \n",
       "COCL.11.S7.001.fa      ...   \n",
       "WACH.38.S1.001.fa      ...   \n",
       "PRPI.42.S31.001.fa     ...   \n",
       "RODE.80.S8.001.fa      ...   \n",
       "pos.dna.01.S17.001.fa  ...   \n",
       "MEEN.93.S11.001.fa     ...   \n",
       "DETH.41.S27.001.fa     ...   \n",
       "TECA.45.S15.001.fa     ...   \n",
       "BODA.52.S9.001.fa      ...   \n",
       "CHMA.61.S3.001.fa      ...   \n",
       "ASTER.02.S34.001.fa    ...   \n",
       "BIRO.74.S19.001.fa     ...   \n",
       "CAAU.07.S2.001.fa      ...   \n",
       "DIGMA.71.S21.001.fa    ...   \n",
       "FAFE.16.S10.001.fa     ...   \n",
       "OPSE.53.S5.001.fa      ...   \n",
       "LAES.91.S29.001.fa     ...   \n",
       "DOEC.81.S30.001.fa     ...   \n",
       "EKJO.64.S32.001.fa     ...   \n",
       "DIAL.15.S23.001.fa     ...   \n",
       "GIFI.21.S22.001.fa     ...   \n",
       "TAOT.65.S24.001.fa     ...   \n",
       "DEFI.14.S6.001.fa      ...   \n",
       "MUNA.98.S28.001.fa     ...   \n",
       "ROAL.75.S33.001.fa     ...   \n",
       "MERI.28.S20.001.fa     ...   \n",
       "BUMA.05.S4.001.fa      ...   \n",
       "MOEM.48.S16.001.fa     ...   \n",
       "PAMA.46.S18.001.fa     ...   \n",
       "VUGE.37.S14.001.fa     ...   \n",
       "TOSA.76.S25.001.fa     ...   \n",
       "DEIS.70.S26.001.fa     ...   \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Porphyromonadaceae;g__Porphyromonas;s__Porphyromonas_uenonis_A  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                               \n",
       "CAPMA.83.S12.001.fa                                             3.138940                                                                               \n",
       "FOGI.63.S13.001.fa                                              2.407628                                                                               \n",
       "COCL.11.S7.001.fa                                               2.444156                                                                               \n",
       "WACH.38.S1.001.fa                                              -2.117196                                                                               \n",
       "PRPI.42.S31.001.fa                                              2.369289                                                                               \n",
       "RODE.80.S8.001.fa                                               3.704593                                                                               \n",
       "pos.dna.01.S17.001.fa                                          -0.923157                                                                               \n",
       "MEEN.93.S11.001.fa                                             -2.942135                                                                               \n",
       "DETH.41.S27.001.fa                                              1.070765                                                                               \n",
       "TECA.45.S15.001.fa                                              1.255953                                                                               \n",
       "BODA.52.S9.001.fa                                              -1.812368                                                                               \n",
       "CHMA.61.S3.001.fa                                              -6.522518                                                                               \n",
       "ASTER.02.S34.001.fa                                            -0.781043                                                                               \n",
       "BIRO.74.S19.001.fa                                             -4.852802                                                                               \n",
       "CAAU.07.S2.001.fa                                              -0.491175                                                                               \n",
       "DIGMA.71.S21.001.fa                                            -3.997182                                                                               \n",
       "FAFE.16.S10.001.fa                                             -4.479546                                                                               \n",
       "OPSE.53.S5.001.fa                                              -3.781932                                                                               \n",
       "LAES.91.S29.001.fa                                             -0.659404                                                                               \n",
       "DOEC.81.S30.001.fa                                             -2.364446                                                                               \n",
       "EKJO.64.S32.001.fa                                             -1.004701                                                                               \n",
       "DIAL.15.S23.001.fa                                             -2.826564                                                                               \n",
       "GIFI.21.S22.001.fa                                             -1.117893                                                                               \n",
       "TAOT.65.S24.001.fa                                             -2.552559                                                                               \n",
       "DEFI.14.S6.001.fa                                              -3.038194                                                                               \n",
       "MUNA.98.S28.001.fa                                              0.302494                                                                               \n",
       "ROAL.75.S33.001.fa                                              2.554822                                                                               \n",
       "MERI.28.S20.001.fa                                              0.439290                                                                               \n",
       "BUMA.05.S4.001.fa                                              -0.666852                                                                               \n",
       "MOEM.48.S16.001.fa                                             -4.459434                                                                               \n",
       "PAMA.46.S18.001.fa                                             -1.335513                                                                               \n",
       "VUGE.37.S14.001.fa                                              0.423770                                                                               \n",
       "TOSA.76.S25.001.fa                                              0.618325                                                                               \n",
       "DEIS.70.S26.001.fa                                             -0.406268                                                                               \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_sp001553265  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                       \n",
       "CAPMA.83.S12.001.fa                                             4.841704                                                                       \n",
       "FOGI.63.S13.001.fa                                              1.353132                                                                       \n",
       "COCL.11.S7.001.fa                                               4.645267                                                                       \n",
       "WACH.38.S1.001.fa                                               3.141126                                                                       \n",
       "PRPI.42.S31.001.fa                                             -5.137255                                                                       \n",
       "RODE.80.S8.001.fa                                               1.527456                                                                       \n",
       "pos.dna.01.S17.001.fa                                          -1.010168                                                                       \n",
       "MEEN.93.S11.001.fa                                             -2.301677                                                                       \n",
       "DETH.41.S27.001.fa                                              0.845063                                                                       \n",
       "TECA.45.S15.001.fa                                             -2.019154                                                                       \n",
       "BODA.52.S9.001.fa                                              -2.907754                                                                       \n",
       "CHMA.61.S3.001.fa                                              -5.223235                                                                       \n",
       "ASTER.02.S34.001.fa                                            -1.652275                                                                       \n",
       "BIRO.74.S19.001.fa                                             -4.475508                                                                       \n",
       "CAAU.07.S2.001.fa                                              -0.625905                                                                       \n",
       "DIGMA.71.S21.001.fa                                            -5.422697                                                                       \n",
       "FAFE.16.S10.001.fa                                             -5.285171                                                                       \n",
       "OPSE.53.S5.001.fa                                              -4.310000                                                                       \n",
       "LAES.91.S29.001.fa                                             -3.815749                                                                       \n",
       "DOEC.81.S30.001.fa                                             -2.452479                                                                       \n",
       "EKJO.64.S32.001.fa                                              2.080180                                                                       \n",
       "DIAL.15.S23.001.fa                                             -1.847991                                                                       \n",
       "GIFI.21.S22.001.fa                                              3.923379                                                                       \n",
       "TAOT.65.S24.001.fa                                              3.106923                                                                       \n",
       "DEFI.14.S6.001.fa                                               1.291514                                                                       \n",
       "MUNA.98.S28.001.fa                                              0.470776                                                                       \n",
       "ROAL.75.S33.001.fa                                              1.559750                                                                       \n",
       "MERI.28.S20.001.fa                                              2.206092                                                                       \n",
       "BUMA.05.S4.001.fa                                              -1.504957                                                                       \n",
       "MOEM.48.S16.001.fa                                             -5.226689                                                                       \n",
       "PAMA.46.S18.001.fa                                              1.219461                                                                       \n",
       "VUGE.37.S14.001.fa                                             -1.154333                                                                       \n",
       "TOSA.76.S25.001.fa                                              1.556518                                                                       \n",
       "DEIS.70.S26.001.fa                                              3.801718                                                                       \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Prevotella;s__Prevotella_oris  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                \n",
       "CAPMA.83.S12.001.fa                                             3.443421                                                                \n",
       "FOGI.63.S13.001.fa                                             -0.121505                                                                \n",
       "COCL.11.S7.001.fa                                               1.707124                                                                \n",
       "WACH.38.S1.001.fa                                               1.670049                                                                \n",
       "PRPI.42.S31.001.fa                                             -6.075525                                                                \n",
       "RODE.80.S8.001.fa                                              -1.924759                                                                \n",
       "pos.dna.01.S17.001.fa                                          -1.798626                                                                \n",
       "MEEN.93.S11.001.fa                                             -2.368428                                                                \n",
       "DETH.41.S27.001.fa                                              0.288315                                                                \n",
       "TECA.45.S15.001.fa                                             -2.957723                                                                \n",
       "BODA.52.S9.001.fa                                              -2.278180                                                                \n",
       "CHMA.61.S3.001.fa                                              -4.676691                                                                \n",
       "ASTER.02.S34.001.fa                                            -2.837292                                                                \n",
       "BIRO.74.S19.001.fa                                             -3.103602                                                                \n",
       "CAAU.07.S2.001.fa                                              -4.495425                                                                \n",
       "DIGMA.71.S21.001.fa                                            -2.208633                                                                \n",
       "FAFE.16.S10.001.fa                                             -5.385255                                                                \n",
       "OPSE.53.S5.001.fa                                              -4.449762                                                                \n",
       "LAES.91.S29.001.fa                                             -2.204774                                                                \n",
       "DOEC.81.S30.001.fa                                             -3.483327                                                                \n",
       "EKJO.64.S32.001.fa                                              5.461192                                                                \n",
       "DIAL.15.S23.001.fa                                             -2.486585                                                                \n",
       "GIFI.21.S22.001.fa                                              5.297858                                                                \n",
       "TAOT.65.S24.001.fa                                              5.611904                                                                \n",
       "DEFI.14.S6.001.fa                                              -0.052186                                                                \n",
       "MUNA.98.S28.001.fa                                             -0.653966                                                                \n",
       "ROAL.75.S33.001.fa                                              0.630927                                                                \n",
       "MERI.28.S20.001.fa                                             -0.342471                                                                \n",
       "BUMA.05.S4.001.fa                                              -2.146903                                                                \n",
       "MOEM.48.S16.001.fa                                             -7.098492                                                                \n",
       "PAMA.46.S18.001.fa                                             -0.602835                                                                \n",
       "VUGE.37.S14.001.fa                                             -3.232563                                                                \n",
       "TOSA.76.S25.001.fa                                              1.321092                                                                \n",
       "DEIS.70.S26.001.fa                                             10.340782                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp900066825  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -2.730828                                                                                \n",
       "FOGI.63.S13.001.fa                                             -2.117494                                                                                \n",
       "COCL.11.S7.001.fa                                              -2.901118                                                                                \n",
       "WACH.38.S1.001.fa                                              -2.016654                                                                                \n",
       "PRPI.42.S31.001.fa                                              0.473331                                                                                \n",
       "RODE.80.S8.001.fa                                              -0.191832                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -0.667223                                                                                \n",
       "MEEN.93.S11.001.fa                                              0.770305                                                                                \n",
       "DETH.41.S27.001.fa                                              1.436088                                                                                \n",
       "TECA.45.S15.001.fa                                              0.560763                                                                                \n",
       "BODA.52.S9.001.fa                                               0.806878                                                                                \n",
       "CHMA.61.S3.001.fa                                               0.755190                                                                                \n",
       "ASTER.02.S34.001.fa                                             0.171671                                                                                \n",
       "BIRO.74.S19.001.fa                                              1.275340                                                                                \n",
       "CAAU.07.S2.001.fa                                               0.780260                                                                                \n",
       "DIGMA.71.S21.001.fa                                             4.121940                                                                                \n",
       "FAFE.16.S10.001.fa                                              1.517647                                                                                \n",
       "OPSE.53.S5.001.fa                                              -0.532751                                                                                \n",
       "LAES.91.S29.001.fa                                              1.430708                                                                                \n",
       "DOEC.81.S30.001.fa                                              0.449270                                                                                \n",
       "EKJO.64.S32.001.fa                                             -0.952057                                                                                \n",
       "DIAL.15.S23.001.fa                                              0.842740                                                                                \n",
       "GIFI.21.S22.001.fa                                              0.153738                                                                                \n",
       "TAOT.65.S24.001.fa                                              0.568337                                                                                \n",
       "DEFI.14.S6.001.fa                                              -0.119842                                                                                \n",
       "MUNA.98.S28.001.fa                                             -0.514026                                                                                \n",
       "ROAL.75.S33.001.fa                                              0.524060                                                                                \n",
       "MERI.28.S20.001.fa                                             -0.432335                                                                                \n",
       "BUMA.05.S4.001.fa                                              -0.352957                                                                                \n",
       "MOEM.48.S16.001.fa                                              0.976780                                                                                \n",
       "PAMA.46.S18.001.fa                                             -0.996740                                                                                \n",
       "VUGE.37.S14.001.fa                                              0.427640                                                                                \n",
       "TOSA.76.S25.001.fa                                              0.577063                                                                                \n",
       "DEIS.70.S26.001.fa                                              1.103640                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Faecalimonas;s__Faecalimonas_sp900556835  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                            \n",
       "CAPMA.83.S12.001.fa                                            -2.124692                                                                            \n",
       "FOGI.63.S13.001.fa                                             -3.004797                                                                            \n",
       "COCL.11.S7.001.fa                                              -3.214776                                                                            \n",
       "WACH.38.S1.001.fa                                              -1.179972                                                                            \n",
       "PRPI.42.S31.001.fa                                             -1.156355                                                                            \n",
       "RODE.80.S8.001.fa                                              -1.928109                                                                            \n",
       "pos.dna.01.S17.001.fa                                          -0.635475                                                                            \n",
       "MEEN.93.S11.001.fa                                              1.697052                                                                            \n",
       "DETH.41.S27.001.fa                                             -1.135138                                                                            \n",
       "TECA.45.S15.001.fa                                             -1.981214                                                                            \n",
       "BODA.52.S9.001.fa                                              -1.435304                                                                            \n",
       "CHMA.61.S3.001.fa                                               0.602765                                                                            \n",
       "ASTER.02.S34.001.fa                                            -0.867344                                                                            \n",
       "BIRO.74.S19.001.fa                                             -0.624752                                                                            \n",
       "CAAU.07.S2.001.fa                                               0.033700                                                                            \n",
       "DIGMA.71.S21.001.fa                                            -0.361749                                                                            \n",
       "FAFE.16.S10.001.fa                                             -0.069201                                                                            \n",
       "OPSE.53.S5.001.fa                                              -0.095620                                                                            \n",
       "LAES.91.S29.001.fa                                             -2.543399                                                                            \n",
       "DOEC.81.S30.001.fa                                             -1.440480                                                                            \n",
       "EKJO.64.S32.001.fa                                             -2.130712                                                                            \n",
       "DIAL.15.S23.001.fa                                             -1.756216                                                                            \n",
       "GIFI.21.S22.001.fa                                             -2.727331                                                                            \n",
       "TAOT.65.S24.001.fa                                             -2.958024                                                                            \n",
       "DEFI.14.S6.001.fa                                               1.011532                                                                            \n",
       "MUNA.98.S28.001.fa                                             -0.816727                                                                            \n",
       "ROAL.75.S33.001.fa                                             -2.359089                                                                            \n",
       "MERI.28.S20.001.fa                                             -3.599673                                                                            \n",
       "BUMA.05.S4.001.fa                                              -1.530629                                                                            \n",
       "MOEM.48.S16.001.fa                                              2.587275                                                                            \n",
       "PAMA.46.S18.001.fa                                             -0.908771                                                                            \n",
       "VUGE.37.S14.001.fa                                             -0.313027                                                                            \n",
       "TOSA.76.S25.001.fa                                             -2.275568                                                                            \n",
       "DEIS.70.S26.001.fa                                             -3.573851                                                                            \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__Hungatella_A;s__Hungatella_A_sp003478355  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                            \n",
       "CAPMA.83.S12.001.fa                                            -0.929852                                                                            \n",
       "FOGI.63.S13.001.fa                                             -1.820056                                                                            \n",
       "COCL.11.S7.001.fa                                              -2.330573                                                                            \n",
       "WACH.38.S1.001.fa                                              -1.957565                                                                            \n",
       "PRPI.42.S31.001.fa                                             -0.917108                                                                            \n",
       "RODE.80.S8.001.fa                                              -1.326253                                                                            \n",
       "pos.dna.01.S17.001.fa                                          -1.010168                                                                            \n",
       "MEEN.93.S11.001.fa                                              0.457184                                                                            \n",
       "DETH.41.S27.001.fa                                             -0.484491                                                                            \n",
       "TECA.45.S15.001.fa                                              0.301699                                                                            \n",
       "BODA.52.S9.001.fa                                               0.612192                                                                            \n",
       "CHMA.61.S3.001.fa                                               0.508339                                                                            \n",
       "ASTER.02.S34.001.fa                                            -0.210790                                                                            \n",
       "BIRO.74.S19.001.fa                                             -0.903965                                                                            \n",
       "CAAU.07.S2.001.fa                                              -1.256430                                                                            \n",
       "DIGMA.71.S21.001.fa                                             1.308846                                                                            \n",
       "FAFE.16.S10.001.fa                                              1.602478                                                                            \n",
       "OPSE.53.S5.001.fa                                              -1.119344                                                                            \n",
       "LAES.91.S29.001.fa                                              0.343855                                                                            \n",
       "DOEC.81.S30.001.fa                                             -0.737190                                                                            \n",
       "EKJO.64.S32.001.fa                                             -0.723798                                                                            \n",
       "DIAL.15.S23.001.fa                                              1.252917                                                                            \n",
       "GIFI.21.S22.001.fa                                             -1.341037                                                                            \n",
       "TAOT.65.S24.001.fa                                             -1.399879                                                                            \n",
       "DEFI.14.S6.001.fa                                              -1.295560                                                                            \n",
       "MUNA.98.S28.001.fa                                             -0.916989                                                                            \n",
       "ROAL.75.S33.001.fa                                             -2.246501                                                                            \n",
       "MERI.28.S20.001.fa                                             -0.984146                                                                            \n",
       "BUMA.05.S4.001.fa                                              -0.746264                                                                            \n",
       "MOEM.48.S16.001.fa                                              1.486921                                                                            \n",
       "PAMA.46.S18.001.fa                                             -0.037698                                                                            \n",
       "VUGE.37.S14.001.fa                                              0.374880                                                                            \n",
       "TOSA.76.S25.001.fa                                             -0.003606                                                                            \n",
       "DEIS.70.S26.001.fa                                              0.557308                                                                            \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Oscillospirales;f__Oscillospiraceae;g__Lawsonibacter;s__Lawsonibacter_sp000177015  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                                \n",
       "CAPMA.83.S12.001.fa                                            -4.522588                                                                                \n",
       "FOGI.63.S13.001.fa                                             -3.537014                                                                                \n",
       "COCL.11.S7.001.fa                                              -4.213305                                                                                \n",
       "WACH.38.S1.001.fa                                              -4.815676                                                                                \n",
       "PRPI.42.S31.001.fa                                             -1.554945                                                                                \n",
       "RODE.80.S8.001.fa                                              -1.318927                                                                                \n",
       "pos.dna.01.S17.001.fa                                          -2.309451                                                                                \n",
       "MEEN.93.S11.001.fa                                             -0.315932                                                                                \n",
       "DETH.41.S27.001.fa                                             -1.364955                                                                                \n",
       "TECA.45.S15.001.fa                                             -1.118817                                                                                \n",
       "BODA.52.S9.001.fa                                              -0.467510                                                                                \n",
       "CHMA.61.S3.001.fa                                              -0.216243                                                                                \n",
       "ASTER.02.S34.001.fa                                            -1.297471                                                                                \n",
       "BIRO.74.S19.001.fa                                              1.473348                                                                                \n",
       "CAAU.07.S2.001.fa                                              -2.549514                                                                                \n",
       "DIGMA.71.S21.001.fa                                            -0.674639                                                                                \n",
       "FAFE.16.S10.001.fa                                              0.442297                                                                                \n",
       "OPSE.53.S5.001.fa                                              -2.203747                                                                                \n",
       "LAES.91.S29.001.fa                                              0.111238                                                                                \n",
       "DOEC.81.S30.001.fa                                             -1.007263                                                                                \n",
       "EKJO.64.S32.001.fa                                             -2.536177                                                                                \n",
       "DIAL.15.S23.001.fa                                             -0.041654                                                                                \n",
       "GIFI.21.S22.001.fa                                             -3.132796                                                                                \n",
       "TAOT.65.S24.001.fa                                             -2.734880                                                                                \n",
       "DEFI.14.S6.001.fa                                              -2.430593                                                                                \n",
       "MUNA.98.S28.001.fa                                             -1.873143                                                                                \n",
       "ROAL.75.S33.001.fa                                             -3.651487                                                                                \n",
       "MERI.28.S20.001.fa                                             -0.108111                                                                                \n",
       "BUMA.05.S4.001.fa                                              -1.368205                                                                                \n",
       "MOEM.48.S16.001.fa                                             -1.131064                                                                                \n",
       "PAMA.46.S18.001.fa                                             -2.690059                                                                                \n",
       "VUGE.37.S14.001.fa                                             -0.862795                                                                                \n",
       "TOSA.76.S25.001.fa                                             -2.781117                                                                                \n",
       "DEIS.70.S26.001.fa                                             -3.350707                                                                                \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Paraprevotella;s__Paraprevotella_sp003477995  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                               \n",
       "CAPMA.83.S12.001.fa                                            -2.443146                                                                               \n",
       "FOGI.63.S13.001.fa                                             -2.626143                                                                               \n",
       "COCL.11.S7.001.fa                                              -3.594266                                                                               \n",
       "WACH.38.S1.001.fa                                              -1.953475                                                                               \n",
       "PRPI.42.S31.001.fa                                             -1.227844                                                                               \n",
       "RODE.80.S8.001.fa                                              -1.467165                                                                               \n",
       "pos.dna.01.S17.001.fa                                          -2.714916                                                                               \n",
       "MEEN.93.S11.001.fa                                             -3.596317                                                                               \n",
       "DETH.41.S27.001.fa                                             -0.362683                                                                               \n",
       "TECA.45.S15.001.fa                                             -2.342227                                                                               \n",
       "BODA.52.S9.001.fa                                              -0.727855                                                                               \n",
       "CHMA.61.S3.001.fa                                               0.590895                                                                               \n",
       "ASTER.02.S34.001.fa                                             1.110924                                                                               \n",
       "BIRO.74.S19.001.fa                                              1.265571                                                                               \n",
       "CAAU.07.S2.001.fa                                               0.328233                                                                               \n",
       "DIGMA.71.S21.001.fa                                             0.604328                                                                               \n",
       "FAFE.16.S10.001.fa                                             -0.797606                                                                               \n",
       "OPSE.53.S5.001.fa                                               0.762726                                                                               \n",
       "LAES.91.S29.001.fa                                             -2.009986                                                                               \n",
       "DOEC.81.S30.001.fa                                             -2.681997                                                                               \n",
       "EKJO.64.S32.001.fa                                             -3.517006                                                                               \n",
       "DIAL.15.S23.001.fa                                             -1.053558                                                                               \n",
       "GIFI.21.S22.001.fa                                             -2.909653                                                                               \n",
       "TAOT.65.S24.001.fa                                             -2.264877                                                                               \n",
       "DEFI.14.S6.001.fa                                              -3.090619                                                                               \n",
       "MUNA.98.S28.001.fa                                             -4.406562                                                                               \n",
       "ROAL.75.S33.001.fa                                             -2.515834                                                                               \n",
       "MERI.28.S20.001.fa                                             -3.217738                                                                               \n",
       "BUMA.05.S4.001.fa                                              -2.212961                                                                               \n",
       "MOEM.48.S16.001.fa                                             -0.922624                                                                               \n",
       "PAMA.46.S18.001.fa                                             -1.335513                                                                               \n",
       "VUGE.37.S14.001.fa                                             -2.197253                                                                               \n",
       "TOSA.76.S25.001.fa                                             -3.086498                                                                               \n",
       "DEIS.70.S26.001.fa                                             -1.246573                                                                               \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Bacteroidota;c__Bacteroidia;o__Bacteroidales;f__Bacteroidaceae;g__Bacteroides;s__Bacteroides_oleiciplenus  \\\n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                                          \n",
       "CAPMA.83.S12.001.fa                                            -3.675290                                                                          \n",
       "FOGI.63.S13.001.fa                                             -1.696893                                                                          \n",
       "COCL.11.S7.001.fa                                              -0.489334                                                                          \n",
       "WACH.38.S1.001.fa                                              -2.093441                                                                          \n",
       "PRPI.42.S31.001.fa                                             -0.051808                                                                          \n",
       "RODE.80.S8.001.fa                                              -2.030491                                                                          \n",
       "pos.dna.01.S17.001.fa                                          -2.155300                                                                          \n",
       "MEEN.93.S11.001.fa                                             -1.389233                                                                          \n",
       "DETH.41.S27.001.fa                                             -2.735072                                                                          \n",
       "TECA.45.S15.001.fa                                             -1.738692                                                                          \n",
       "BODA.52.S9.001.fa                                              -1.865448                                                                          \n",
       "CHMA.61.S3.001.fa                                               0.945615                                                                          \n",
       "ASTER.02.S34.001.fa                                            -0.424424                                                                          \n",
       "BIRO.74.S19.001.fa                                              1.526685                                                                          \n",
       "CAAU.07.S2.001.fa                                               0.564564                                                                          \n",
       "DIGMA.71.S21.001.fa                                            -0.262493                                                                          \n",
       "FAFE.16.S10.001.fa                                             -1.121834                                                                          \n",
       "OPSE.53.S5.001.fa                                              -0.606017                                                                          \n",
       "LAES.91.S29.001.fa                                              0.495566                                                                          \n",
       "DOEC.81.S30.001.fa                                             -0.129765                                                                          \n",
       "EKJO.64.S32.001.fa                                             -1.907569                                                                          \n",
       "DIAL.15.S23.001.fa                                              0.031977                                                                          \n",
       "GIFI.21.S22.001.fa                                             -2.573181                                                                          \n",
       "TAOT.65.S24.001.fa                                             -2.552559                                                                          \n",
       "DEFI.14.S6.001.fa                                              -1.534714                                                                          \n",
       "MUNA.98.S28.001.fa                                             -1.715021                                                                          \n",
       "ROAL.75.S33.001.fa                                             -1.401189                                                                          \n",
       "MERI.28.S20.001.fa                                             -3.448850                                                                          \n",
       "BUMA.05.S4.001.fa                                              -1.443736                                                                          \n",
       "MOEM.48.S16.001.fa                                              0.683690                                                                          \n",
       "PAMA.46.S18.001.fa                                             -2.977741                                                                          \n",
       "VUGE.37.S14.001.fa                                             -1.232786                                                                          \n",
       "TOSA.76.S25.001.fa                                             -2.393351                                                                          \n",
       "DEIS.70.S26.001.fa                                              0.276297                                                                          \n",
       "\n",
       "#OTU ID                k__Bacteria;p__Firmicutes_A;c__Clostridia;o__Lachnospirales;f__Lachnospiraceae;g__CAG-81;s__CAG-81_sp900066785  \n",
       "pos.dna.02.S35.001.fa                                          -2.475253                                                               \n",
       "CAPMA.83.S12.001.fa                                            -3.829440                                                               \n",
       "FOGI.63.S13.001.fa                                             -2.913825                                                               \n",
       "COCL.11.S7.001.fa                                              -3.856630                                                               \n",
       "WACH.38.S1.001.fa                                              -1.995148                                                               \n",
       "PRPI.42.S31.001.fa                                             -1.882509                                                               \n",
       "RODE.80.S8.001.fa                                              -1.091414                                                               \n",
       "pos.dna.01.S17.001.fa                                          -1.703315                                                               \n",
       "MEEN.93.S11.001.fa                                             -1.218305                                                               \n",
       "DETH.41.S27.001.fa                                              0.163711                                                               \n",
       "TECA.45.S15.001.fa                                             -0.965607                                                               \n",
       "BODA.52.S9.001.fa                                               0.195402                                                               \n",
       "CHMA.61.S3.001.fa                                              -0.637804                                                               \n",
       "ASTER.02.S34.001.fa                                             0.490222                                                               \n",
       "BIRO.74.S19.001.fa                                             -1.350001                                                               \n",
       "CAAU.07.S2.001.fa                                              -0.054863                                                               \n",
       "DIGMA.71.S21.001.fa                                             0.456438                                                               \n",
       "FAFE.16.S10.001.fa                                              0.994868                                                               \n",
       "OPSE.53.S5.001.fa                                              -1.502694                                                               \n",
       "LAES.91.S29.001.fa                                             -1.028367                                                               \n",
       "DOEC.81.S30.001.fa                                             -0.393402                                                               \n",
       "EKJO.64.S32.001.fa                                             -1.843030                                                               \n",
       "DIAL.15.S23.001.fa                                             -0.745484                                                               \n",
       "GIFI.21.S22.001.fa                                             -3.132796                                                               \n",
       "TAOT.65.S24.001.fa                                             -1.208824                                                               \n",
       "DEFI.14.S6.001.fa                                              -1.908701                                                               \n",
       "MUNA.98.S28.001.fa                                             -0.154564                                                               \n",
       "ROAL.75.S33.001.fa                                             -2.638173                                                               \n",
       "MERI.28.S20.001.fa                                             -2.688044                                                               \n",
       "BUMA.05.S4.001.fa                                              -1.600636                                                               \n",
       "MOEM.48.S16.001.fa                                              0.319989                                                               \n",
       "PAMA.46.S18.001.fa                                             -2.327154                                                               \n",
       "VUGE.37.S14.001.fa                                             -0.935806                                                               \n",
       "TOSA.76.S25.001.fa                                             -2.393351                                                               \n",
       "DEIS.70.S26.001.fa                                             -3.350707                                                               \n",
       "\n",
       "[35 rows x 181 columns]"
      ]
     },
     "execution_count": 119,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_taxa_clr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "df_taxa_clr.index.name = \"sample_id\"\n",
    "\n",
    "df_metadata_file.loc[df_metadata_file[\"sample_id\"] == \"MEIR_28\", \"sample_id\"] = \"MERI_28\"\n",
    "\n",
    "df_merged = pd.merge(df_taxa_clr, df_metadata_file, on=\"sample_id\", how=\"left\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# MERI 28 is misspelled\n",
    "set(sample_ids).difference(set(df_metadata_file[\"sample_id\"]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "[col for col in df_merged.columns if \"RDV\" in col]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# hospitalization_day\n",
    "# outcome severe\n",
    "# CONC_N1\n",
    "\n",
    "# RDV"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Veillonella"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'df_merged' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-80-44d86050273b>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mdf_merged\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"RDV\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m: name 'df_merged' is not defined"
     ]
    }
   ],
   "source": [
    "df_merged[\"RDV\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [],
   "source": [
    "for col in genera:\n",
    "    sns.boxplot(x=\"RDV\", y=col, data=df_merged)\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[]"
      ]
     },
     "execution_count": 75,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "species = \"Enterococcus_faecalis\"\n",
    "\n",
    "[col for col in df_taxa_nc.columns if species in col]"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:.conda-type_1]",
   "language": "python",
   "name": "conda-env-.conda-type_1-py"
  },
  "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.8.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}