Session 14: converting some Perl to Python

[pvanheus]

Our starting point for this session is @kamzune's Perl script calculate_nuc.pl. The challenge: convert this to Python (with some BioPython as well I expect).

[pvanheus]

The original code can be found here: https://github.com/kamzune/nucleotide-frequencies

[pvanheus]

The resultant code is in https://github.com/SANBI-SA/studyGroupWork/tree/master/session14 .

We started with calculate_nuc.pl and did a pretty straightforward conversion to calculate_nuc.py.orig and then switched to using BioPython (installed via conda install -y biopython because my Anaconda install has bioconda channels) for calculate_nuc.py.bioperl and finally added a nicer command line interface with click to get calculate_nuc.py. Try it with:

python calculate_nuc.py seq.fasta