sam_pair_iterator.h

Go to the documentation of this file.

#ifndef __gamgee__sam_pair_iterator__
#define __gamgee__sam_pair_iterator__

#include "sam.h"

#include "htslib/sam.h"

#include <fstream>
#include <queue>
#include <memory>

namespace gamgee {

class SamPairIterator {
  public:

    SamPairIterator();

    SamPairIterator(samFile * sam_file_ptr, const std::shared_ptr<bam_hdr_t>& sam_header_ptr);

    SamPairIterator(SamPairIterator&& other);
    
    bool operator!=(const SamPairIterator& rhs);

    std::pair<Sam,Sam> operator*();

    std::pair<Sam,Sam> operator++();

  private:
    using SamPtrQueue = std::queue<std::shared_ptr<bam1_t>>;

    SamPtrQueue m_supp_alignments;                     
    samFile * m_sam_file_ptr;                          
    const std::shared_ptr<bam_hdr_t> m_sam_header_ptr; 
    std::shared_ptr<bam1_t> m_sam_record_ptr1;         
    std::shared_ptr<bam1_t> m_sam_record_ptr2;         
    std::pair<Sam,Sam> m_sam_records;                  

    std::pair<Sam,Sam> fetch_next_pair();              
    bool read_sam(std::shared_ptr<bam1_t>& record_ptr);                 
    Sam make_sam(std::shared_ptr<bam1_t>& record_ptr);                  
    Sam next_primary_alignment(std::shared_ptr<bam1_t>& record_ptr);
    std::pair<Sam,Sam> next_supplementary_alignment();
};

}  // end namespace gamgee

#endif