Augustus 3.4.0
|
global multiple sequence alignment More...
#include <alignment.hh>
Public Member Functions | |
AlignmentRow (string seqID, int chrPos, Strand strand, string rowbuffer) | |
int | chrStart () const |
int | chrEnd () const |
int | aliEnd () const |
int | getSeqLen () const |
int | getCumFragLen () const |
void | setCumFragLen (int len) |
int | gapLenAfterFrag (size_t i) const |
void | addFragment (int chrPos, int aliPos, int len) |
void | addFragment (fragment &f) |
string | getSignature () const |
void | pack () |
int | getAliPos (int chrPos, vector< fragment >::const_iterator from) |
int | getAliPos (int chrPos, vector< fragment >::const_iterator *from) |
int | getAliPos (int chrPos) |
int | getChrPos (int aliPos, vector< fragment >::const_iterator from) |
int | getChrPos (int aliPos) |
Public Attributes | |
string | seqID |
Strand | strand |
vector< fragment > | frags |
Friends | |
ostream & | operator<< (ostream &strm, const AlignmentRow &row) |
void | appendRow (AlignmentRow **r1, const AlignmentRow *r2, int aliLen1, string sigstr) |
global multiple sequence alignment
alignment does not contain the sequence itself, but the info, where gaps are.
Generation of exon candidates.
int AlignmentRow::getAliPos | ( | int | chrPos, |
vector< fragment >::const_iterator | from | ||
) |
convert from chromosomal to alignment position start search from the fragment 'from' on, i.e. assume that aliPos is not to the left of fragment *from return -1, if chrPos is outside the range of these fragments return -2, if position is otherwise not mappable: no fragment contains the chrPos, i.e. chrPos is in a gap
|
friend |
append row r2 to r1, thereby shifting all alignment coordinates of r2 by aliLen1 if signature string (chr and strand) is not empty, treat everything else as missing