def plot_matrix_muts(reactivity_type, max_val = 0.025):
matrix_data = cotrans_matrix_data(reactivity_type, max_val = max_val)
plt.matshow(matrix_data)
ax = plt.gca()
ax.grid(color='grey',linestyle='-',linewidth='0.5')
ax.xaxis.set_ticks_position('bottom')
plt.xlabel("Nucleotide (nt)")
plt.ylabel("RNA Length (nt)")
cbar = plt.colorbar(fraction=0.046, pad=0.04)
cbar.set_label(reactivity_type)
return plt
def plot_cotrans_untreated():
matrix_data = cotrans_matrix_data('untreated_counts', max_val = 5000)
plt.matshow(matrix_data,cmap='jet')
ax = plt.gca()
ax.grid(color='grey',linestyle='-',linewidth='0.5')
ax.xaxis.set_ticks_position('bottom')
plt.xlabel("Nucleotide (nt)")
plt.ylabel("RNA Length (nt)")
cbar = plt.colorbar(fraction=0.046, pad=0.04)
cbar.set_label('Counts')
plt.title('Untreated Channel Counts')
return plt