void check_interval(const Interval& actual, const Interval& truth) { BOOST_CHECK_EQUAL(actual.chr(), truth.chr()); BOOST_CHECK_EQUAL(actual.start(), truth.start()); BOOST_CHECK_EQUAL(actual.stop(), truth.stop()); }
string ReferenceMap::get_sequence(const Interval& interval, const bool reverse_strand) const { const auto& seq = this->at(interval.chr()); auto result = seq.substr(interval.start()-1, interval.size()); return reverse_strand ? utils::complement(result) : result; }