Class
Gaps
Efficient storage of gaps for a sequence.
Gaps wrap a SequenceConcept and allows to (1) insert gaps into the sequence and (2) select an infix of the gapped sequence (clipping). The gaps are not inserted into the underlying sequence (source) but stored separately. Using the clipping is optional and meant for selecting parts of the alignment as a part of the result of a local alignment algorithm.
In the figure above, the source sequence has seven characters, the gapped sequence has four gaps and thus consists of eleven characters. The gapped sequence is clipped to start at position 0 in the gapped sequence and to end at position 8 in the gapped sequence (the positions given as half-open intervals [begin, end)).
The figure shows the three coordinate systems that are used with Gaps objects. The source position is the position in the underlying sequence. The unclipped view position is the position in the gapped sequence without gaps. The view position is the position in the gapped sequence but including the clipping: All (clipped) view positions have the clipping begin position subtracted from them.
Illustration of Gaps object and positions with clipping.
Gaps<TSequence, TSpec>
Include Headers
seqan/align.h
Parameters
TSequence
The type of the underlying sequence.
TSpec
Specialization tag.
Implements
Specializations
AnchorGapsStores gaps anchors of the first characters behind gaps.
ArrayGapsStores length of gap- and non-gapped runs in an array.
Metafunctions
GetValueType for reading values.
IsSequenceDetermines whether a container stores its elements in sequential order.
IteratorType of iterator objects that are used to traverse the container.
PositionType of an object that represents a position in a container.
ReferenceReference type.
SizeType of an object that is suitable to hold size information.
SourceReturn underlying sequence of Gaps/Alignments.
ValueType of the items in the container or behind an iterator.
Member Functions
GapsConstructor
Functions
assignSourceAssign the source of a Gaps object, copying data.
assignValueAssigns value to item. (ContainerConcept)
backThe last item in container. (ContainerConcept)
bandedChainAlignmentComputes the best global pairwise alignment between two sequences given a non-empty seed chain.
beginThe begin of a container. (ContainerConcept)
beginPositionBegin position of object in host. (ContainerConcept)
clearClippingClear clipping from Gaps object.
clearGapsClear gaps and clipping from Gaps object.
clippedBeginPositionReturn the begin position of the clipping in the unclipped gapped sequence.
clippedEndPositionReturn the end position of the clipping in the unclipped gapped sequence.
endThe end of a container. (ContainerConcept)
endPositionEnd position of object in host. (ContainerConcept)
frontThe first item in container. (ContainerConcept)
getValueAccess to the value. (ContainerConcept)
globalAlignmentComputes the best global pairwise alignment.
insertGapInsert one gap into a gapped sequence.
insertGapsInsert multiple gaps into a gapped sequence.
isGapQuery whether a given clipped view position is a gap.
lengthThe number of items/characters. (ContainerConcept)
localAlignmentComputes the best pairwise local alignment using the Smith-Waterman algorithm.
moveValueAssigns value to item. (ContainerConcept)
nextLocalAlignmentCompute next suboptimal local alignment.
positionPosition of an iterator. (ContainerConcept)
removeGapRemove one gap from a gapped sequence.
removeGapsRemove multiple gaps from a gapped sequence.
setBeginPositionSet the begin position of the clipped gapped sequence, given a source position.
setClippedBeginPositionSets the begin position of the clipping.
setClippedEndPositionSets the end position of the clipping.
setEndPositionSet the end position of the clipped gapped sequence, given a source position.
setSourceSet the source of a Gaps object, do not copy if possible.
sourceReturn underlying object.
splitAlignmentCompute split alignments.
toSourcePositionTransforms view to source position, if the view position is a gap, the original position of the next non-gap entry is returned.
toViewPositionTransforms source to view position.
unclippedLengthReturn length of the gapped sequence without clipping.
SeqAn - Sequence Analysis Library - www.seqan.de
 

Page built @2013/07/11 09:12:16