Epinano workflow, how to get per-read results

[adalisan]

I have a question about the epinano workflow. As I understand it, the code uses Epinano_variants.py script to compute reference misalignment frequencies at positions in order to compute per-site statistics.
There is a second step in Epinano library (implemented in Epinano_differr.R) which computes different statistics including per_read outliers and generates some plots. MOP2 doesn't use this script, am I correct?
Is there any replacement for this second step in MOP2? Was this step removed/not implemented on purpose?
Thank you for your work in this useful package. :)

[lucacozzuto]

@ADelgadoT