A RetroSearch Logo

Home - News ( United States | United Kingdom | Italy | Germany ) - Football scores

Search Query:

Showing content from http://www.ncbi.nlm.nih.gov/IEB/ToolBox/CPP_DOC/doxyhtml/traceback_8cpp_source.html below:

NCBI C++ ToolKit: src/algo/cobalt/traceback.cpp Source File

53  if

(m_Script.empty()) {

58  if

(last_op.

op_type

== op_type) {

129  for

(

int i

= tback_range.

GetFrom

();

i

<= tback_range.

GetTo

();

i

++) {

130  switch

(tback[

i

]) {

160

TScriptOps::iterator itr(

m_Script

.begin());

162  if

(link_start + itr->num_ops > tback_range.

GetFrom

())

165

link_start += itr->num_ops;

173  int

curr_tback = tback_range.

GetFrom

();

175

curr_tback <= tback_range.

GetTo

()) {

181  int

num_ops =

min

(tback_range.

GetTo

() - curr_tback + 1,

182

link_start + itr->num_ops - curr_tback);

184

new_script.

AddOps

(itr->op_type, num_ops);

185

curr_tback += num_ops;

186

link_start += itr->num_ops;

199  bool

go_past_seq1_gap)

202  TOffset

q = start_offsets.first;

203  TOffset

s = start_offsets.second;

204

TScriptOps::iterator itr(

m_Script

.begin());

216  if

(s + itr->num_ops > seq2_target) {

217  if

(go_past_seq1_gap ==

false

) {

218

q--; s--; curr_tback--;

222

curr_tback += itr->num_ops;

229

curr_tback += itr->num_ops;

236  if

(s + itr->num_ops > seq2_target) {

237

curr_tback += seq2_target - s;

238

q += seq2_target - s;

244

curr_tback += itr->num_ops;

253

curr_tback += itr->num_ops;

263

new_offsets.first = q;

264

new_offsets.second = s;

265

new_tback = curr_tback;

273  bool

go_past_seq2_gap)

276  TOffset

q = start_offsets.first;

277  TOffset

s = start_offsets.second;

278

TScriptOps::iterator itr(

m_Script

.begin());

291

curr_tback += itr->num_ops;

298  if

(q + itr->num_ops > seq1_target) {

299

curr_tback += seq1_target - q;

300

s += seq1_target - q;

306

curr_tback += itr->num_ops;

314  if

(q + itr->num_ops > seq1_target) {

315  if

(go_past_seq2_gap ==

false

) {

316

q--; s--; curr_tback--;

320

curr_tback += itr->num_ops;

327

curr_tback += itr->num_ops;

337

new_offsets.first = q;

338

new_offsets.second = s;

339

new_tback = curr_tback;

346  int

**seq2_pssm,

int

gap_open,

int

gap_extend)

348  TOffset

q = start_offsets.first;

349  TOffset

s = start_offsets.second;

356

TScriptOps::iterator itr(

m_Script

.begin());

358  if

(link_start + itr->num_ops > tback_range.

GetFrom

())

361  switch

(itr->op_type) {

376

link_start += itr->num_ops;

383  switch

(itr->op_type) {

385

s += tback_range.

GetFrom

() - link_start;

388

s += tback_range.

GetFrom

() - link_start;

389

q += tback_range.

GetFrom

() - link_start;

392

q += tback_range.

GetFrom

() - link_start;

401

curr_tback <= tback_range.

GetTo

()) {

407  int

num_ops =

min

(tback_range.

GetTo

() - curr_tback + 1,

408

link_start + itr->num_ops - curr_tback);

411  switch

(itr->op_type) {

413

score += gap_open + num_ops * gap_extend;

417  for

(

int i

= 0;

i

< num_ops;

i

++) {

419

score += seq2_pssm[s +

i

][c];

425

score += gap_open + num_ops * gap_extend;

432

curr_tback += num_ops;

433

link_start += itr->num_ops;

446

vector<TOffsetPair> region_list;

447  int

q = start_offsets.first;

448  int

s = start_offsets.second;

450

TScriptOps::iterator itr(

m_Script

.begin());

452  switch

(itr->op_type) {

int TOffset

Basic data type for offsets into a sequence.

pair< TOffset, TOffset > TOffsetPair

Basic type specifying a range on a sequence.

Interface for the traceback from blast hits.

void VerifyScript(TRange seq1_range, TRange seq2_range)

Validate that the alignment described by the CEditScript has the same size for each sequence as the i...

int GetScore(TRange tback_range, TOffsetPair start_offsets, CSequence &seq1, int **seq2_pssm, int gap_open, int gap_extend)

Compute the score associated with (a portion of) an alignment Assumes that seq1 is a sequence and tha...

void FindOffsetFromSeq2(TOffsetPair start_offsets, TOffsetPair &new_offsets, TOffset seq2_target, TOffset &new_tback, bool go_past_seq1_gap)

Given a subject offset, find the corresponding query offset.

TScriptOps m_Script

The list of edit operations in the current edit script.

CEditScript()

Generate empty traceback.

void FindOffsetFromSeq1(TOffsetPair start_offsets, TOffsetPair &new_offsets, TOffset seq1_target, TOffset &new_tback, bool go_past_seq2_gap)

Given a query offset, find the corresponding subject offset.

CEditScript MakeEditScript(TRange tback_range)

Return an edit script corresponding to a subset of the complete traceback available.

vector< TOffsetPair > ListMatchRegions(TOffsetPair start_offsets)

Compile a list of regions in the current edit script that contain substitutions.

void AddOps(EGapAlignOpType op_type, int num_ops)

Add a new edit operation to the current list.

Class for representing protein sequences.

unsigned char GetLetter(int pos) const

Access the sequence letter at a specified position.

EGapAlignOpType

Operation types within the edit script.

@ eGapAlignIns

Insertion: a gap in subject.

@ eGapAlignSub

Substitution.

@ eGapAlignDel

Deletion: a gap in query.

vector< ETranscriptSymbol > TTranscript

int TSignedSeqPos

Type for signed sequence position.

int32_t Int4

4-byte (32-bit) signed integer

position_type GetLength(void) const

#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.

const TStarts & GetStarts(void) const

Get the Starts member data.

const TLens & GetLens(void) const

Get the Lens member data.

TLen GetLen(void) const

Get the Len member data.

vector< TSignedSeqPos > TStarts

TDim GetDim(void) const

Get the Dim member data.

TDim GetDim(void) const

Get the Dim member data.

TNumseg GetNumseg(void) const

Get the Numseg member data.

unsigned int

A callback function used to compare two keys in a database.

Runlength-encoded representation of a traceback operation.

int num_ops

number of such operations

EGapAlignOpType op_type

type of operation

Edit script: linked list of correspondencies between two sequences.

Int4 * num

Array of number of operations.

Int4 size

Size of above arrays.

EGapAlignOpType * op_type

Array of type of operation.

Interface for CEditScript class.


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