sequence_handle.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_SEQ_SEQUENCE_HANDLE_HH
#define OST_SEQ_SEQUENCE_HANDLE_HH

/*
  Author: Marco Biasini
 */

#include <ost/base.hh>
#include <ost/generic_property.hh>
#include <ost/config.hh>
#if(OST_INFO_ENABLED)
#include <ost/info/info_fw.hh>
#endif
#include <ost/seq/module_config.hh>

#include <ost/mol/residue_view.hh>
#include <ost/mol/entity_view.hh>

#include <ost/seq/impl/sequence_impl_fw.hh>

namespace ost { namespace seq {

class SequenceHandle;
class AlignmentHandle;

class DLLEXPORT_OST_SEQ ConstSequenceHandle : 
  public ConstGenericPropContainer<ConstSequenceHandle> {
public:
  friend class ConstGenericPropContainer<ConstSequenceHandle>;
  friend class AlignmentHandle;
  friend class ConstSequenceList;
  friend class SequenceList;
  
  typedef String::const_iterator iterator;
  
  ConstSequenceHandle();
  
  int GetResidueIndex(int pos) const;

  int GetPos(int residue_index) const;

  int GetFirstNonGap() const;
  
  int GetLastNonGap() const;  
  
  const String& GetName() const;  
  
  const String& GetString() const;
  
  String GetGaplessString() const;

  int GetOffset() const;  
  
  int GetLength() const;
  
  int GetIndex(const String& substr) const;

  char GetOneLetterCode(int position) const;

  mol::ResidueView GetResidue(int position) const;

  mol::EntityView GetAttachedView() const;
  
  SequenceHandle Copy() const;
  
  bool HasAttachedView() const;
  
  const String& GetRole() const;

  bool operator==(const ConstSequenceHandle& rhs) const;
  bool operator!=(const ConstSequenceHandle& rhs) const;  
  
  char operator[](int index) const;
  
  iterator begin() const { return this->GetString().begin(); }
  iterator end() const { return this->GetString().end(); }
  
  bool IsValid() const;
protected:
  GenericPropContainerImpl* GpImpl();

  const GenericPropContainerImpl* GpImpl() const;
public:
  ConstSequenceHandle(const impl::SequenceImplPtr& impl);
  impl::SequenceImplPtr& Impl() const;  
private:
  void CheckValidity() const;
  mutable impl::SequenceImplPtr impl_;  
};

//
class DLLEXPORT_OST_SEQ SequenceHandle : 
  public GenericPropContainer<SequenceHandle> {
public:
  typedef String::const_iterator iterator;
  friend class GenericPropContainer<SequenceHandle>;  
  friend class SequenceList;

  friend class AlignmentHandle;
  friend class ConstSequenceList;
  
  int GetResidueIndex(int pos) const;

  int GetPos(int residue_index) const;

  int GetFirstNonGap() const;
  
  int GetLastNonGap() const;  
  
  const String& GetName() const;  
  
  const String& GetString() const;
  
  String GetGaplessString() const;

  int GetOffset() const;  
  
  int GetLength() const;

  char GetOneLetterCode(int position) const;

  int GetIndex(const String& substr) const;
  
  mol::ResidueView GetResidue(int position) const;

  mol::EntityView GetAttachedView() const;

  void Normalise();

  SequenceHandle Copy() const;
  
  bool HasAttachedView() const;
  
  bool operator==(const SequenceHandle& rhs) const;
  bool operator!=(const SequenceHandle& rhs) const;  
  
  bool IsValid() const;
    
  SequenceHandle();
  
  void SetName(const String& name);

  void SetString(const String& seq);

  void SetOffset(int offset);

  void AttachView(const mol::EntityView& view);
  
  void SetOneLetterCode(int position, char new_char);
  
  char operator[](size_t index) const;
    
  iterator begin() const { return this->GetString().begin(); }
  iterator end() const { return this->GetString().end(); }
  
  void Append(char olc);

  operator ConstSequenceHandle() const;
  void AttachView(const mol::EntityView& view, const String& chain_name);
   
  SequenceHandle(const impl::SequenceImplPtr& impl);  
  
  const String& GetRole() const;
  
  void SetRole(const String& role) const;
  
  impl::SequenceImplPtr& Impl() const;  

  GenericPropContainerImpl* GpImpl();

  const GenericPropContainerImpl* GpImpl() const;
private:
  void CheckValidity() const;  
  mutable impl::SequenceImplPtr impl_;
};

SequenceHandle DLLEXPORT_OST_SEQ CreateSequence(const String& name, 
                                                const String& seq, 
                                                const String& role="UNKNOWN");
#if(OST_INFO_ENABLED)                                                
void DLLEXPORT_OST_SEQ SequenceToInfo(const ConstSequenceHandle& sequence,
                                      info::InfoGroup& group);

SequenceHandle DLLEXPORT_OST_SEQ SequenceFromInfo(info::InfoGroup& group);
#endif


DLLEXPORT_OST_SEQ std::ostream& operator<<(std::ostream& os, 
                                           const ConstSequenceHandle& sequence);

bool DLLEXPORT_OST_SEQ Match(const ConstSequenceHandle& s1,
                             const ConstSequenceHandle& s2);

bool DLLEXPORT_OST_SEQ Match(const String& s1,
                             const String& s2);
}}

#endif