read_bases.h

Go to the documentation of this file.

#ifndef __gamgee__read_bases__
#define __gamgee__read_bases__

#include "htslib/sam.h"

#include <memory>
#include <map>

namespace gamgee {

enum class Base { A = 1, C = 2, G = 4, T = 8, N = 15 };

class ReadBases {
public:
  explicit ReadBases(const std::shared_ptr<bam1_t>& sam_record);
  ReadBases(const ReadBases& other);
  ReadBases(ReadBases&& other) noexcept;
  ReadBases& operator=(const ReadBases& other);
  ReadBases& operator=(ReadBases&& other) noexcept;
  ~ReadBases() = default; 

  Base operator[](const uint32_t index) const;   
  uint32_t size() const { return m_num_bases; }; 
  std::string to_string() const;

private:
  std::shared_ptr<bam1_t> m_sam_record; 
  uint8_t* m_bases;                     
  uint32_t m_num_bases;                 

  static const std::map<Base, const char*> base_to_string_map; 
};

}

#endif /* __gamgee__read_bases__ */