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/aln__writer_8cpp_source.html below:

NCBI C++ ToolKit: src/objtools/writers/aln_writer.cpp Source File

64  unsigned int

uFlags) :

82  unsigned int

uFlags) :

133  "Empty bioseq handle"

);

147  "Empty sequence string"

);

165  const auto

num_rows = denseg.

GetDim

();

166  const auto

num_segs = denseg.

GetNumseg

();

169  for

(

int row

=0;

row

<num_rows; ++

row

)

179  "Unable to fetch Bioseq"

);

191  for

(

int

seg=0; seg<num_segs; ++seg)

193  const auto

start = denseg.

GetStarts

()[seg*num_rows +

row

];

202  string

defline =

">"

+ best_id;

256  "Unable to determine product length"

);

262  "Unable to determine product length"

);

267  return

(length >= 0) ? length : -length;

279  string

prev_genomic_id;

280  string

prev_product_id;

284

exon->IsSetGenomic_id() ?

285

exon->GetGenomic_id() :

289

exon->IsSetProduct_id() ?

290

exon->GetProduct_id() :

295

exon->IsSetGenomic_strand() ?

296

exon->GetGenomic_strand() :

297

default_genomic_strand;

300

exon->IsSetProduct_strand() ?

301

exon->GetProduct_strand() :

302

default_product_strand;

304  const auto

genomic_start = exon->GetGenomic_start();

305  const auto

genomic_end = exon->GetGenomic_end();

307  if

(genomic_end < genomic_start) {

310  "Bad genomic location: end < start"

);

312  const int

genomic_length = genomic_end - genomic_start;

314  const int

product_start = exon->GetProduct_start().AsSeqPos();

315  const int

product_end = exon->GetProduct_end().AsSeqPos();

318  if

(product_end < product_start) {

321  "Bad product location: end < start"

);

324  const int

product_length = product_end - product_start;

333  "Unable to resolve genomic sequence"

);

337  GetSeqString

(genomic_bsh, genomic_range, genomic_strand, genomic_seq);

346  "Unable to resolve product sequence"

);

350  GetSeqString

(product_bsh, product_range, product_strand, product_seq);

352  if

(exon->IsSetParts()) {

353  AddGaps

(product_type, exon->GetParts(), genomic_seq, product_seq);

356  if

(product_length != genomic_length) {

359  "Lengths of genomic and product sequences don't match"

);

374  string

& genomic_seq,

375  string

& product_seq)

378  if

(exon_chunks.empty()) {

382  string

genomic_string;

383  string

product_string;

385  const unsigned int

res_width =

393  unsigned int

interval_width = 0;

395  auto

chunk_type = exon_chunk->Which();

396  switch

(chunk_type) {

400  switch

(chunk_type) {

405

interval_width = exon_chunk->GetMatch();

408

interval_width = exon_chunk->GetMismatch();

411

interval_width = exon_chunk->GetDiag();

414

genomic_string.append(genomic_seq, genomic_pos, interval_width);

415

product_string.append(product_seq, product_pos, (interval_width + (res_width-1))/res_width);

416

genomic_pos += interval_width;

417

product_pos += interval_width/res_width;

421

interval_width = exon_chunk->GetGenomic_ins();

422

genomic_string.append(genomic_seq, genomic_pos, interval_width);

423

product_string.append(interval_width/res_width,

'-'

);

424

genomic_pos += interval_width;

428

interval_width = exon_chunk->GetProduct_ins();

429

genomic_string.append(interval_width,

'-'

);

430

product_string.append(product_seq, product_pos, interval_width/res_width);

431

product_pos += interval_width/res_width;

437

genomic_seq = genomic_string;

438

product_seq = product_string;

451  if

(start >= seq_plus.size()) {

454  "Bad location: impossible start"

);

457  return

seq_plus.substr(start,

len

);

488  const auto

num_segs = sparse_align.

GetNumseg

();

509  for

(

int

seg=0; seg<num_segs; ++seg) {

511  const auto len

= sparse_align.

GetLens

()[seg];

538  for

(

int

seg=0; seg<num_segs; ++seg) {

540  const auto len

= sparse_align.

GetLens

()[seg];

553  if

(defline.back() ==

'|'

&& defline.size() > 1)

555  const auto

length = defline.size();

556  m_Os

<< defline.substr(0,length-1) <<

"\n"

;

559  m_Os

<< defline <<

"\n"

;

562  while

(pos < seqdata.size()) {

568  "Processing terminated by user"

);

User-defined methods of the data storage class.

User-defined methods of the data storage class.

User-defined methods of the data storage class.

User-defined methods of the data storage class.

User-defined methods of the data storage class.

User-defined methods of the data storage class.

User-defined methods of the data storage class.

unsigned int s_ProductLength(const CProduct_pos &start, const CProduct_pos &end)

void ProcessSeqId(const CSeq_id &id, CBioseq_Handle &bsh, CRange< TSeqPos > &range)

bool WriteSparseAlign(const CSparse_align &sparse_aln)

void WriteContiguous(const string &defline, const string &seqdata)

bool WriteSplicedExons(const list< CRef< CSpliced_exon >> &exons, CSpliced_seg::TProduct_type product_type, CRef< CSeq_id > default_genomic_id, ENa_strand default_genomic_strand, CRef< CSeq_id > default_product_id, ENa_strand default_product_strand)

CAlnWriter(CScope &scope, CNcbiOstream &ostr, unsigned int uFlags)

void GetSeqString(CBioseq_Handle bsh, const CRange< TSeqPos > &range, ENa_strand strand, string &seq)

bool WriteAlignSparseSeg(const CSparse_seg &sparse_seg)

bool WriteAlign(const CSeq_align &align, const string &name="", const string &descr="") override

Write a raw Seq-align to the internal output stream.

string GetSegString(const string &seq_plus, CSeqUtil::ECoding coding, ENa_strand strand, int start, size_t len)

void AddGaps(CSpliced_seg::TProduct_type product_type, const CSpliced_exon::TParts &exon_chunks, string &genomic_seq, string &product_seq)

bool WriteAlignSplicedSeg(const CSpliced_seg &spliced_seg)

string GetBestId(const CSeq_id &)

bool WriteAlignDenseSeg(const CDense_seg &denseg)

const CSeq_id & GetSeq_id(TDim row) const

static SIZE_TYPE ReverseComplement(const string &src, TCoding src_coding, TSeqPos pos, TSeqPos length, string &dst)

Defines and provides stubs for a general interface to a variety of file formatters.

bool GetBestId(CSeq_id_Handle idh, CScope &scope, string &best_id) const

unique_ptr< CGenbankIdResolve > mpIdResolve

#define NCBI_THROW(exception_class, err_code, message)

Generic macro to throw an exception, given the exception class, error code and message string.

const string AsFastaString(void) const

virtual void Assign(const CSerialObject &source, ESerialRecursionMode how=eRecursive)

Optimized implementation of CSerialObject::Assign, which is not so efficient.

static CSeq_id_Handle GetHandle(const CSeq_id &id)

Normal way of getting a handle, works for any seq-id.

@ eFasta

Tagged ID in NCBI's traditional FASTA style.

CBioseq_Handle GetBioseqHandle(const CSeq_id &id)

Get bioseq handle by seq-id.

bool IsNucleotide(void) const

TSeqPos GetBioseqLength(void) const

CSeqVector GetSeqVector(EVectorCoding coding, ENa_strand strand=eNa_strand_plus) const

Get sequence: Iupacna or Iupacaa if use_iupac_coding is true.

@ eCoding_Iupac

Set coding to printable coding (Iupacna or Iupacaa)

void GetSeqData(TSeqPos start, TSeqPos stop, string &buffer) const

Fill the buffer string with the sequence data for the interval [start, stop).

CRef< C > Ref(C *object)

Helper functions to get CRef<> and CConstRef<> objects.

void Reset(void)

Reset reference object.

TThisType & SetToOpen(position_type toOpen)

#define END_NCBI_SCOPE

End previously defined NCBI scope.

#define BEGIN_NCBI_SCOPE

Define ncbi namespace.

IO_PREFIX::ostream CNcbiOstream

Portable alias for ostream.

static bool IsBlank(const CTempString str, SIZE_TYPE pos=0)

Check if a string is blank (has no text).

static const char label[]

void SetFrom(TFrom value)

Assign a value to From data member.

TTo GetTo(void) const

Get the To member data.

TFrom GetFrom(void) const

Get the From member data.

bool CanGetDim(void) const

Check if it is safe to call GetDim method.

const TDenseg & GetDenseg(void) const

Get the variant data.

E_Choice Which(void) const

Which variant is currently selected.

const TGenomic_id & GetGenomic_id(void) const

Get the Genomic_id member data.

E_Choice Which(void) const

Which variant is currently selected.

bool IsSetStrands(void) const

Check if a value has been assigned to Strands data member.

bool IsSetProduct_strand(void) const

should be 'plus' or 'minus' Check if a value has been assigned to Product_strand data member.

const TStarts & GetStarts(void) const

Get the Starts member data.

const TProduct_id & GetProduct_id(void) const

Get the Product_id member data.

const TLens & GetLens(void) const

Get the Lens member data.

bool CanGetNumseg(void) const

Check if it is safe to call GetNumseg method.

const TFirst_id & GetFirst_id(void) const

Get the First_id member data.

TProduct_type GetProduct_type(void) const

Get the Product_type member data.

TGenomic_strand GetGenomic_strand(void) const

Get the Genomic_strand member data.

bool CanGetIds(void) const

Check if it is safe to call GetIds method.

const TSpliced & GetSpliced(void) const

Get the variant data.

TDim GetDim(void) const

Get the Dim member data.

bool IsSetGenomic_strand(void) const

Check if a value has been assigned to Genomic_strand data member.

const TLens & GetLens(void) const

Get the Lens member data.

const TExons & GetExons(void) const

Get the Exons member data.

bool IsSetSecond_strands(void) const

Check if a value has been assigned to Second_strands data member.

const TFirst_starts & GetFirst_starts(void) const

Get the First_starts member data.

TNumseg GetNumseg(void) const

Get the Numseg member data.

bool IsSetExons(void) const

set of segments involved each segment corresponds to one exon exons are always in biological order Ch...

const TSecond_strands & GetSecond_strands(void) const

Get the Second_strands member data.

list< CRef< CSpliced_exon_chunk > > TParts

const TSecond_id & GetSecond_id(void) const

Get the Second_id member data.

bool CanGetStarts(void) const

Check if it is safe to call GetStarts method.

TNumseg GetNumseg(void) const

Get the Numseg member data.

TProduct_strand GetProduct_strand(void) const

Get the Product_strand member data.

const TSparse & GetSparse(void) const

Get the variant data.

const TRows & GetRows(void) const

Get the Rows member data.

const TStrands & GetStrands(void) const

Get the Strands member data.

const TSegs & GetSegs(void) const

Get the Segs member data.

bool CanGetLens(void) const

Check if it is safe to call GetLens method.

bool IsSetGenomic_id(void) const

Check if a value has been assigned to Genomic_id data member.

@ e_Product_ins

insertion in product sequence (i.e. gap in the genomic sequence)

@ e_Diag

both sequences are represented, there is sufficient similarity between product and genomic sequences....

@ e_Genomic_ins

insertion in genomic sequence (i.e. gap in the product sequence)

@ e_Match

both sequences represented, product and genomic sequences match

@ e_Mismatch

both sequences represented, product and genomic sequences do not match

@ e_not_set

No variant selected.

@ eProduct_type_transcript

ENa_strand

strand of nucleic acid

#define row(bind, expected)


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