SuperlocusIterator.hpp

Go to the documentation of this file.

// Copyright 2010 Complete Genomics, Inc.
//
// Licensed under the Apache License, Version 2.0 (the "License"); you
// may not use this file except in compliance with the License. You
// may obtain a copy of the License at
//
// http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
// implied. See the License for the specific language governing
// permissions and limitations under the License.

#ifndef CGATOOLS_VARIANTS_SUPERLOCUSITERATOR_HPP_
#define CGATOOLS_VARIANTS_SUPERLOCUSITERATOR_HPP_ 1


#include "cgatools/core.hpp"
#include "cgatools/variants/Superlocus.hpp"
#include "cgatools/variants/VariantFileIterator.hpp"

#include <deque>

namespace cgatools { namespace variants {

    class SuperlocusQueueIterator;

    class SuperlocusIterator : private boost::noncopyable
    {
    public:
        SuperlocusIterator(uint32_t extend3Mers = 4, uint32_t extendBases = 0);

        void setVariantFile(VariantFileIterator& iter);

        void setVariantFiles(VariantFileIterator& iterA, VariantFileIterator& iterB);

        void setVariantFiles(const std::vector<VariantFileIterator*> iters);

        void skipToVariant(const reference::Range& range, const std::string& alleleSeq);

        void seekFirst();

        bool eof() const
        {
            return eof_;
        }

        const Superlocus& operator*() const
        {
            return sl_;
        }

        const Superlocus* operator->() const
        {
            return &sl_;
        }

        reference::Location getPrecedingRefStart() const
        {
            return precedingRefStart_;
        }

        SuperlocusIterator& operator++()
        {
            next();
            return *this;
        }

    private:
        void next();
        bool findVariant(
            const reference::Location& loc, bool trackRefStart,
            reference::Range& range, reference::Range& searchRange);
        void retireLoci();
        void extendVariant();
        void extendVariant(reference::Range& slRange, reference::Range& searchRange);

        reference::Range extendSearchRange(
            const reference::Range& range, const Locus& locus) const;
        uint32_t extendLeftBySuffixMatching(const Locus& locus) const;
        uint32_t extendRightByPrefixMatching(const Locus& locus) const;
        uint32_t extendLeftBySuffixMatching(
            const reference::Location& loc, const std::string& sequence) const;
        uint32_t extendRightByPrefixMatching(
            const reference::Location& loc, const std::string& sequence) const;
        size_t findContig(const reference::Location& loc) const;

        Superlocus sl_;
        std::vector<VariantFileIterator*> iters_;
        std::vector< std::deque<Locus> >& queues_;
        std::vector< boost::shared_ptr<SuperlocusQueueIterator> > qIters_;
        bool started_;
        bool eof_;
        uint32_t extend3Mers_;
        uint32_t extendBases_;
        uint32_t padBases_;
        reference::Location queueCutoff_;
        reference::Location slRetired_;
        const reference::CrrFile* crr_;
        std::vector<reference::Range> contigs_;
        // If the superlocus is adjacent to a reference call, this
        // location tracks the beginning of that call; if this location
        // is equal or greater than the superlocus range start, then
        // there's a non-reference call to the left of the superlocus.
        reference::Location precedingRefStart_;
    };

} } // cgatools::variants

#endif // CGATOOLS_VARIANTS_SUPERLOCUSITERATOR_HPP_

---