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
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
#ifndef RRESOLVER_CONTIGS_H
#define RRESOLVER_CONTIGS_H 1
#include "Common/Options.h"
#include "Common/StringUtil.h"
#include "Common/IOUtil.h"
#include "Common/Sequence.h"
#include "Common/ConstString.h"
#include "Common/ContigPath.h"
#include "Common/ContigNode.h"
#include "Common/ContigProperties.h"
#include "Common/Histogram.h"
#include "DataLayer/FastaReader.h"
#include "Graph/ContigGraph.h"
#include "Graph/ContigGraphAlgorithms.h"
#include "Graph/DirectedGraph.h"
#include "Graph/GraphIO.h"
#include "Graph/GraphUtil.h"
#include <string>
#include <vector>
#include <set>
typedef ContigGraph<DirectedGraph<ContigProperties, Distance> > Graph;
typedef Graph::vertex_descriptor vertex_descriptor;
typedef Graph::edge_descriptor edge_descriptor;
typedef Graph::adjacency_iterator adjacency_iterator;
typedef Graph::in_edge_iterator in_edge_iterator;
typedef Graph::out_edge_iterator out_edge_iterator;
class ContigSequence : public const_string {
public:
ContigSequence(const std::string& s): const_string(s) {
strSize = const_string::size();
}
size_t size() const { return strSize; }
private:
size_t strSize;
};
typedef std::vector<std::pair<std::string, std::string>> ContigSequencesInfo;
typedef std::vector<ContigSequence> ContigSequences;
typedef std::vector<ContigNode> ContigPath;
typedef std::vector<ContigPath> ContigPaths;
typedef std::vector<std::pair<ContigNode, int>> ImaginaryContigPath; // Each element is a pair of contig and distance to the next
typedef std::set<ImaginaryContigPath> ImaginaryContigPaths;
extern Graph g_contigGraph;
extern std::vector<std::string> g_contigComments;
extern ContigSequences g_contigSequences;
int distanceBetween(const ContigNode& node1, const ContigNode& node2);
const ContigSequence& getContigSequence(const ContigNode &node);
int getContigSize(const ContigNode& node);
const std::string& getContigComment(const ContigNode &node);
Sequence getPathSequence(const ContigPath &path);
Sequence getPathSequence(const ImaginaryContigPath &path);
long getContigNumOfKmers(const ContigNode& node);
double getContigBaseCoverage(const ContigNode& node);
void loadContigGraph(const std::string &contigGraphPath);
void loadContigs(const std::string &contigsPath);
void storeContigGraph(const std::string &contigGraphPath, const std::string &program, const std::string& commandLine);
void storeContigs(const std::string &contigsPath);
unsigned num_vertices_removed(const Graph& contigGraph);
void assembleContigs();
#endif
|
