sequence_impl.hh

Go to the documentation of this file.

//------------------------------------------------------------------------------
// This file is part of the OpenStructure project <www.openstructure.org>
//
// Copyright (C) 2008-2020 by the OpenStructure authors
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 3.0 of the License, or (at your option)
// any later version.
// This library is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License for more
// details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this library; if not, write to the Free Software Foundation, Inc.,
// 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//------------------------------------------------------------------------------
#ifndef OST_SEQUENCE_HH
#define OST_SEQUENCE_HH
 
/*
  Author: Marco Biasini
 */
#include <ost/message.hh> 
#include <list>
 
#include <boost/shared_ptr.hpp>
#include <ost/generic_property.hh>
#include <ost/config.hh>
#if(OST_INFO_ENABLED)
#include <ost/info/info_fw.hh>
#endif
#include <ost/mol/residue_prop.hh>
#include <ost/mol/entity_view.hh>
#include <ost/mol/residue_view.hh>
 
#include <ost/seq/module_config.hh>
#include <ost/seq/invalid_sequence.hh>
namespace ost { namespace seq { namespace impl {
 
class SequenceImpl;
 
typedef boost::shared_ptr<SequenceImpl> SequenceImplPtr;
 
class DLLEXPORT_OST_SEQ SequenceImpl : public GenericPropContainerImpl {
public:
  static SequenceImplPtr FromString(const String& seq_name,
                                    const String& sequence_string,
                                    const String& role="UNKNOWN");
 
  int GetResidueIndex(int pos) const;
 
  int GetPos(int residue_index) const;
 
  int GetFirstNonGap() const;
  
  int GetLastNonGap() const;  
  
  const String& GetName() const;
 
  void SetName(const String& name);
  
  void Cut(int start, int n);
 
  void ShiftRegion(int start, int end, int amount);
  const String& GetString() const {
    return seq_string_;
  }
  void SetString(const String& seq);
 
  void Normalise();
 
  void Replace(const String& str,int start, int end);
 
  String GetGaplessString() const;
 
  int GetOffset() const;
  
  int GetIndex(const String& substr) const;
  void SetOffset(int offset);
 
  int GetLength() const;
  SequenceImpl(const String& seq_name, const String& sequence_string, 
               const String& role);
 
  char GetOneLetterCode(int position) const;
 
  void SetOneLetterCode(int position, char new_char);
  mol::ResidueView GetResidue(int position) const;
 
  mol::EntityView GetAttachedView() const;
 
  void AttachView(const mol::EntityView& view);
 
  SequenceImplPtr Copy() const;
  void AttachView(const mol::EntityView& view, const String& chain_name);
 
  bool HasAttachedView() const;
  
  void Append(char olc);
  
  char& operator[](size_t index)
  {
    return seq_string_[index];
  }
  char operator[](size_t index) const
  {
    return seq_string_[index];
  }
  const String& GetRole() const
  {
    return seq_role_;
  }
  
  void SetRole(const String& role)
  {
    seq_role_=role;
  }
private:
 
  void ShiftsFromSequence();
 
 
  int GetPosNoBounds(int index) const;
 
  static bool IsSequenceStringSane(const String& seq_string);
 
  typedef struct {
    int start;
    int shift;
  } Shift;
  String              seq_name_;
  String              seq_string_;
  String              seq_role_;
  std::list<Shift>    shifts_;
  bool                editing_;
  int                 offset_;
  mol::EntityView     attached_view_;
};
 
typedef std::vector<SequenceImplPtr>  SequenceList;
 
#if(OST_INFO_ENABLED)
SequenceImplPtr DLLEXPORT_OST_SEQ SequenceImplFromInfo(const info::InfoGroup& group);
 
void DLLEXPORT_OST_SEQ SequenceImplToInfo(const SequenceImplPtr& sequence,
                                     info::InfoGroup& group);
#endif
 
 
}}} //ns
#endif