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
|
#ifndef PAIREDENDHIT_H_
#define PAIREDENDHIT_H_
#include<iostream>
#include "SingleHit.h"
class PairedEndHit : public SingleHit {
public:
PairedEndHit() : SingleHit() {
insertL = 0;
}
PairedEndHit(int sid, int pos, int insertL, double conprb = 0.0) : SingleHit(sid, pos, conprb) {
this->insertL = insertL;
}
int getInsertL() const { return insertL; }
bool read(std::istream&);
void write(std::ostream&);
private:
int insertL; // insert length
};
bool PairedEndHit::read(std::istream& in) {
conprb = 0.0;
return ((bool)(in>>sid>>pos>>insertL));
}
void PairedEndHit::write(std::ostream& out) {
out<<" "<<sid<<" "<<pos<<" "<<insertL;
}
#endif /* PAIREDENDHIT_H_ */
|
