MoldMitochondrialGeneticCode.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: The Bio++ Development Group
//
// SPDX-License-Identifier: CECILL-2.1
 
#ifndef BPP_SEQ_GENETICCODE_MOLDMITOCHONDRIALGENETICCODE_H
#define BPP_SEQ_GENETICCODE_MOLDMITOCHONDRIALGENETICCODE_H
 
 
#include "../Alphabet/NucleicAlphabet.h"
#include "GeneticCode.h"
 
namespace bpp
{
class MoldMitochondrialGeneticCode :
  public virtual GeneticCode
{
public:
  MoldMitochondrialGeneticCode(std::shared_ptr<const NucleicAlphabet> alphabet) :
    GeneticCode(alphabet)
  {
    init_();
  }
 
  virtual ~MoldMitochondrialGeneticCode() {}
 
  virtual MoldMitochondrialGeneticCode* clone() const
  {
    return new MoldMitochondrialGeneticCode(*this);
  }
 
public:
  size_t getNumberOfStopCodons() const { return 2; }
 
  std::vector<int> getStopCodonsAsInt() const
  {
    return std::vector<int>({48, 50});
  }
 
  std::vector<std::string> getStopCodonsAsChar() const
  {
    return std::vector<std::string>({"TAA", "TAG"});
  }
 
  bool isStop(int state) const
  {
    // Test:
    codonAlphabet_->intToChar(state); // throw exception if invalid state!
    return state == 48 || state == 50;
  }
 
  bool isStop(const std::string& state) const
  {
    int i = codonAlphabet_->charToInt(state);
    return i == 48 || i == 50;
  }
 
  bool isAltStart(int state) const
  {
    // Test:
    codonAlphabet_->intToChar(state); // throw exception if invalid state!
    return state == 12 || state == 13 || state == 15 || state == 30 || state == 46 || state == 60 || state == 62;
  }
 
  bool isAltStart(const std::string& state) const
  {
    int i = codonAlphabet_->charToInt(state);
    return i == 12 || i == 13 || i == 15 || i == 30 || i == 46 || i == 60 || i == 62;
  }
 
private:
  void init_();
};
} // end of namespace bpp.
#endif // BPP_SEQ_GENETICCODE_MOLDMITOCHONDRIALGENETICCODE_H