-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathAlignmentSeq.H
57 lines (49 loc) · 1.52 KB
/
AlignmentSeq.H
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
/***********************************************************************
AlignmentSeq.H
BOOM : Bioinformatics Object Oriented Modules
Copyright (C)2012 William H. Majoros ([email protected]).
This is OPEN SOURCE SOFTWARE governed by the Gnu General Public
License (GPL) version 3, as described at www.opensource.org.
***********************************************************************/
#ifndef INCL_AlignmentSeq_H
#define INCL_AlignmentSeq_H
#include "String.H"
#include "Sequence.H"
#include "BitSet.H"
// a sequence (i.e., "track") in an alignment
namespace BOOM {
class AlignmentSeq
{
public:
AlignmentSeq(const BOOM::String &name,int ID,const BitSet &gapSymbols);
const BOOM::String &getName() const;
int getID() const;
void setID(int id) {ID=id;}
int getBegin() const; // relative to source genome, not to alignment!
int getLength() const;
void addHSP(int begin,const BOOM::Sequence &seq);
void extendToLength(int length,Symbol gapSymbol);
const BOOM::Sequence &getSeq() const;
BOOM::Sequence &getSeq();
inline Symbol &operator[](int);
void deleteSequence();
void setTarget();
bool isTarget() const;
void append(const AlignmentSeq &);
void append(const AlignmentSeq &,int begin,int length);
String &getAnnoTrack();
private:
BOOM::String name;
int begin; // relative to source genome, not to alignment!
BOOM::Sequence seq;
bool isTargetTrack;
int ID;
BitSet gapSymbols;
String annoTrack;
};
Symbol &AlignmentSeq::operator[](int index)
{
return seq[index];
}
}
#endif