featureinfo.cc

Go to the documentation of this file.

#include "featureinfo.h"
#include "featuredata.h"
#include <cassert>
#include <sstream>

std::string FeatureInfo::getPosInfo( TR_DNA pos )
{
  if ( pos < m_data->startPos() || pos > m_data->endPos() )
    return std::string();
  
  return myInfoPos(pos);
}

std::string FeatureInfo::getRangeInfo( TR_DNA begin, TR_DNA end )
{
  if ( end < m_data->startPos() || begin > m_data->endPos() )
    return std::string();
  
  return myInfoRange(begin, end);
}

std::string FeatureInfo::myInfoPos( TR_DNA pos )
{
  return std::string();
}

std::string FeatureInfo::myInfoRange( TR_DNA begin, TR_DNA end )
{
  return std::string();
}

//DnpData
std::string DnpInfo::myInfoPos( TR_DNA pos ) 
{
  assert( pos == dnpData()->startPos() );
  assert( pos == dnpData()->endPos() );
  
  ostringstream os;
  
  os<<"DNP ID: "<<dnpData()->get_dnpID();
  
  return os.str();
}

std::string DnpInfo::myInfoRange( TR_DNA begin, TR_DNA end ) 
{
  //Sanity check
  assert( dnpData()->startPos() == dnpData()->endPos() );
  ostringstream os;
  os.str("");
  if ( dnpData()->startPos() >= begin && dnpData()->startPos() <= end) {
    os<<"DNP ID: "<<dnpData()->get_dnpID()<<" at pos "<<dnpData()->startPos();

  }
  
  return os.str();
}

DnpData * DnpInfo::dnpData()
{
    return dynamic_cast<DnpData *>( m_data );
}

//DnaStrData
std::string DnaStrInfo::myInfoPos( TR_DNA pos ) 
{
  assert( pos >= 0 );
  
  assert( pos < dnaStrData()->dnaVector.stlVector().size() );

  ostringstream os;
  
  os<<"Base: "<<dnaStrData()->dnaVector.stlVector()[ pos ];
  
  return os.str();
}

DnaStrData * DnaStrInfo::dnaStrData()
{
    return dynamic_cast<DnaStrData *>( m_data );
}

//QualityData
std::string QualityInfo::myInfoPos( TR_DNA pos ) 
{
  assert( pos >= 0 );
  
  assert( pos < qualityData()->qualityVector.stlVector().size() );

  ostringstream os;
  
  os<<"Qual: "<<qualityData()->qualityVector.stlVector()[ pos ];
  
  return os.str();
}

QualityData * QualityInfo::qualityData()
{
    return dynamic_cast<QualityData *>( m_data );
}

---