content: '''
  Copyleft Oct 14, 2016 Arya Iranmehr, PhD Student, Bafna Lab, UC San Diego,  Email: airanmehr@gmail.com
  '''
  import numpy as np;
  
  np.set_printoptions(linewidth=200, precision=5, suppress=True)
  import pandas as pd;
  
  pd.options.display.max_rows = 20;
  pd.options.display.expand_frame_repr = False
  import seaborn as sns
  import matplotlib as mpl
  import os;
  
  home = os.path.expanduser('~') + '/'
  import Utils.Util as utl
  import Scripts.TimeSeriesPaper.RealData.Utils as rutl
  
  a = rutl.loadAllScores().groupby(level='h', axis=1).apply(rutl.HstatisticAll)
  df = pd.read_pickle(utl.outpath + 'real/scores.df')
  i = df.lrd.sort_values().index[-1]
  df.loc[i]
  
  cd = pd.read_pickle(utl.outpath + 'real/CD.F59.df')
  
  import Utils.Plots as pplt
  import pylab as plt
  
  names = rutl.loadSNPIDs()
  sns.set_style("white", {"grid.color": "0.9", 'axes.linewidth': .5, "grid.linewidth": "9.99"})
  mpl.rc('font', **{'family': 'serif', 'serif': ['Computer Modern']});
  mpl.rc('text', usetex=True)
  reload(pplt)
  f, ax = plt.subplots(1, 2, sharey=True, dpi=300, figsize=(4, 2))
  i = a[0.5].sort_values().index[-1]
  sns.set_context("notebook", font_scale=1, rc={"lines.linewidth": 1.2})
  
  pplt.plotSiteReal(cd.loc[i], ax=ax[0], legend=True)
  ax[0].set_title('{}:{:.0f} ({})'.format(i[0], i[1], names.loc[i]), fontsize=8)
  
  i = df.lrdiff.sort_values().index[-1]
  pplt.plotSiteReal(cd.loc[i], ax=ax[1])
  sns.set_context("notebook", font_scale=1, rc={"lines.linewidth": 1.2})
  
  ax[1].set_title('{}:{:.0f} ({})'.format(i[0], i[1], names.loc[i]), fontsize=8)
  plt.gcf().subplots_adjust(bottom=0.2)
  pplt.savefig('topSNPs', 300)
  plt.show()
  