call and intensity sanity check

[dkj]

ensure called base is the one with highest corrected intensity when applying (and calibrating?) qualities - an assert?

[srl147]

The bases are called by RTA, which treats the intensities as floating point numbers so the channel with the maximum intensity is unique. The cif/dif files contain integers, in practice there are a significant number of bases where multiple channels match the maximum intensity. In src/calibration_pu.c uncomment

//#define CHECK_BASECALL

if you want to do this check. This finds the largest intensity and then checks that called base corresponds to a channel with this intensity

[dkj]

We need to turn this into a assert which runs all the time (to catch problems arising from when we've altered other bits of the pipeline e.g. composing read from a subset of cycles).

[dkj]

And we need this check in predictor_pu i.e. it's most important to confirm base called is base with max intensity where we actually set the qualities....