本文整理汇总了C++中Genome::getReadAt方法的典型用法代码示例。如果您正苦于以下问题：C++ Genome::getReadAt方法的具体用法？C++ Genome::getReadAt怎么用？C++ Genome::getReadAt使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Genome的用法示例。
在下文中一共展示了Genome::getReadAt方法的2个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: addUpstreamRead
void ReadContainer::addUpstreamRead ( const Genome& genome, const chr_num_t chr, const chr_pos_t pos ) {
assume(pos <= 0, "3' ends need to be given as negative numbers!");
int ipos = pos - 1; // fix difference in one-based offsets
auto upstr = genome.getReadAt(chr, ipos); // 3' ends are linked at negative indices
if (upstr != nullptr) {
// check if this is a multistrand splice junction
if (genome.multistrand != nullptr &&
(upstr->chromosome != chromosome || (upstr->flags & MULTISTRAND) == MULTISTRAND)) {
upstr->flags |= MULTISTRAND;
flags |= MULTISTRAND;
genome.multistrand->insert(shared_from_this());
genome.multistrand->insert(upstr);
}
int link = findLink(upstr, false); // search link upstr <= this
int backLink = upstr->findLink(shared_from_this()); // search link upstr => this
if (link == -1) { // upstr <= this unknown yet
fivePrimeRead->push_back(upstr);
}
if (backLink >= 0) { // upstr => this known
upstr->threePrimeRefs->at(backLink)++;
}
else { // upstr => this new
upstr->threePrimeRead->push_back(shared_from_this());
upstr->threePrimeRefs->push_back(1);
}
}
}示例2: addDownstreamRead
void ReadContainer::addDownstreamRead ( const Genome& genome, const chr_num_t chr, const chr_pos_t pos ) {
auto dnstr = genome.getReadAt(chr, pos);
if (dnstr != nullptr) {
// check if this is a multistrand splice junction
if (dnstr->chromosome != chromosome && !(dnstr->flags & MULTISTRAND) ) {
flags |= MULTISTRAND;
genome.multistrand->insert(shared_from_this()); // adding one multistrand seed suffices
}
int link = findLink(dnstr); // search link this => dnstr
int backLink = dnstr->findLink(shared_from_this(), false); // search link this <= dnstr
if (link == -1) { // this => dnstr unknown yet
threePrimeRead->push_back(dnstr);
threePrimeRefs->push_back(1);
}
else { // this => dnstr known
threePrimeRefs->at(link)++;
}
if (backLink == -1) { // this <= dnstr unknown yet
dnstr->fivePrimeRead->push_back(shared_from_this());
}
}
}