->m_intron_params->m_intronlen.MinLen();
82 return m_data->m_intergenic_params->m_intergeniclen.MinLen();
87 return m_data->m_intron_params->m_intronlen.MaxLen();
92 doublep = exp(
m_data->m_intron_params->m_intronlen.ClosingScore(
l));
99 if(
m_data->m_range.Empty() )
112 const intGC_RANGE_SIZE = 200000;
113 TSignedSeqRangegc_range(middle-GC_RANGE_SIZE/2, middle+GC_RANGE_SIZE/2);
115gc_range +=
m_data->m_range;
120 if(c ==
enC|| c ==
enG) {
126 m_data->m_gccontent =
static_cast<int>(
m_data->m_gccontent*100.0/length + 0.5);
140 m_data->m_intron_params->SetSeqLen(
m_data->m_range.GetLength() );
141 m_data->GetHMMParameter(
m_data->m_intergenic_params);
142 m_data->m_intergenic_params->SetSeqLen(
m_data->m_range.GetLength() );
153 return m_data->m_gccontent;
162 return Run( cls, leftwall, rightwall,
false,
false,
163mpp,
BadScore(), notbridgeable_gaps_len, ggapinfo, pcsf_slice
168 boolleftwall,
boolrightwall,
boolleftanchor,
boolrightanchor,
doublempp,
doubleconsensuspenalty,
178*
m_data->m_intron_params,
180cls, initial_fshifts, mpp, *
this)
182 m_data->m_ss->Init(
m_data->m_seq, leftwall, rightwall, consensuspenalty, notbridgeable_gaps_len, ggapinfo, pcsf_slice);
184*
m_data->m_intron_params,
185*
m_data->m_intergenic_params,
187leftanchor, rightanchor) );
188 return m_data->m_parse->Path()->Score();
195 return m_data->m_parse->GetGenes();
203 if(!genes.empty() && !genes.back().RightComplete()) {
204 TSignedSeqPospartial_start = genes.back().Limits().GetFrom();
207 if(!genes.empty()) {
208right = genes.back().Limits().GetTo();
211right = partial_start-100;
233 m_data->m_parse->PrintInfo();
TSignedSeqPos MapEditedToOrig(TSignedSeqPos edited_pos) const
int GetMinIntergenicLen() const
double GetChanceOfIntronLongerThan(int l) const
double Run(bool leftwall=true, bool rightwall=true, double mpp=10, SPhyloCSFSlice *pcsf_slice=nullptr)
list< CGeneModel > GetGenes() const
unique_ptr< SGnomonEngineImplData > m_data
CGnomonEngine(CConstRef< CHMMParameters > hmm_params, const CResidueVec &sequence, TSignedSeqRange range=TSignedSeqRange::GetWhole(), SPhyloCSFSlice *pcsf_slice=nullptr)
TSignedSeqPos PartialModelStepBack(list< CGeneModel > &genes) const
const CResidueVec & GetSeq() const
int GetMaxIntronLen() const
int GetMinIntronLen() const
void ResetRange(TSignedSeqRange range)
HMM model parameters just create it and pass to a Gnomon engine.
vector< TResidue > CResidueVec
list< CGeneModel > TGeneModelList
vector< CInDelInfo > TInDels
int TSignedSeqPos
Type for signed sequence position.
#define NCBI_THROW(exception_class, err_code, message)
Generic macro to throw an exception, given the exception class, error code and message string.
CRange< TSignedSeqPos > TSignedSeqRange
#define END_NCBI_SCOPE
End previously defined NCBI scope.
#define END_SCOPE(ns)
End the previously defined scope.
#define BEGIN_NCBI_SCOPE
Define ncbi namespace.
#define BEGIN_SCOPE(ns)
Define a new scope.
TTo GetTo(void) const
Get the To member data.
TFrom GetFrom(void) const
Get the From member data.
Defines classes: CDirEntry, CFile, CDir, CSymLink, CMemoryFile, CFileUtil, CFileLock,...
static SLJIT_INLINE sljit_ins l(sljit_gpr r, sljit_s32 d, sljit_gpr x, sljit_gpr b)
SGnomonEngineImplData(CConstRef< CHMMParameters > hmm_params, CResidueVec &&sequence, TSignedSeqRange range, SPhyloCSFSlice *pcsf_slice)
const value_slice::CValueConvert< value_slice::SRunTimeCP, FROM > Convert(const FROM &value)
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4