In [ ]:

    

%matplotlib inline


    

In [ ]:

    

import os
import sys

from bokeh.plotting import *

pyrnafold_path = os.path.abspath(os.path.join('..'))
if pyrnafold_path not in sys.path:
    sys.path.append(pyrnafold_path)
    
from pyrnafold.pyrnafold import trange_df
from pyrnafold.parse import parse_fasta
from pyrnafold.plotting import heatmap_T_range, sequence_heatmap


    

In [ ]:

    

output_notebook()


    

In [ ]:

    

df = trange_df('../data/hHSR', trange=range(35,45))

fasta = parse_fasta('../data/hHSR.fa')
for rec in fasta:
    id, seq = rec
    
df['seq'] = [seq[x-1] for x in df['pos']]
df


    

In [ ]:

    

temps = set(df['Temp'])
heatmap_T_range(df, temps, title='Human HSR1')


    

In [ ]:

    

sequence_heatmap(df, temp=42)


    

In [ ]: