Boltzmann sampling plot script needs updating

[maxieds]

I need to update this script using the details in the Slack discussion threads on this topic.

@ceheitsch
Is there a good way to compute the full set of possible maximal helix classes given just a base sequence with no secondary pairing data?

Here is what it produces at current:

[maxieds]

A compendia of the discussions and suggestions internally by @ceheitsch and @hbanoscervantes:

1. You might find some of the python scripts useful (changed to public version): https://github.com/gtDMMB/ipam-wbio-scripts
2. And Forrest has code to index helix classes consistently across different samples: https://github.gatech.edu/fhurley6/HelixEnumeration
3. And here's the python version: https://github.gatech.edu/fhurley6/HelixAnalysis/blob/master/Enumerate.py
4. Q: Figure 1 of the following paper: https://www.cell.com/biophysj/fulltext/S0006-3495(17)30565-9
   
   A: Here is a brief explanation from a paper:
   • Boltzmann sample of 1,000 structures through GTboltzmann (and here for docs).
   • Next, it assigns every helix present in the sample to an equivalence class known as a helix class. Every time a helix is assigned to a helix class, the frequency of the helix class is augmented by 1.
   • See this paper