void WindowsStorage::add(junction window, Sequence sequence, int status, Germline *germline) { if (! hasWindow(window)) { // First time we see that window: init status_by_window[window].resize(STATS_SIZE); seqs_by_window[window].init(nb_bins, max_value_bins, &scorer); seqs_by_window[window].setMaxNbReadsStored(getMaximalNbReadsPerWindow()); } seqs_by_window[window].add(sequence); status_by_window[window][status]++; germline_by_window[window] = germline; }
size_t WindowsStorage::getNbReads(junction window) { assert(hasWindow(window)); return seqs_by_window[window].getNbInserted(); }
float WindowsStorage::getAverageLength(junction window) { assert(hasWindow(window)); return seqs_by_window[window].getAverageLength(); }