Mercurial > repos > vipints > rdiff

annotate rDiff/src/tools/read_utils/get_reads_direct.h @ 0:0f80a5141704

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
version 0.3 uploaded
author vipints
date Thu, 14 Feb 2013 23:38:36 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
1 /* written by Jonas Behr, Regina Bohnert and Gunnar Raetsch, FML Tuebingen, Germany, 2010 */
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
2
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
3 #ifndef __GET_READS_DIRECT_H__
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
4 #define __GET_READS_DIRECT_H__
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
5
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
6 #include <vector>
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
7 using std::vector;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
8 #include "read.h"
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
9 #include <stdlib.h>
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
10 #include "sam.h"
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
11
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
12 //static int g_flag_on = 0, g_flag_off = 0;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
13 static int left_flag_mask = strtol((char*) "0x40", 0, 0);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
14 static int right_flag_mask = strtol((char*) "0x80", 0, 0);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
15 static int reverse_flag_mask = strtol((char*) "0x10", 0, 0);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
16
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
17 static int subsample = 1000;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
18 //static int collapse = 0;
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
19
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
20 int get_reads_from_bam(char* filename, char* region, vector<CRead*>* reads, char strand, int lsubsample);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
21 void parse_cigar(bam1_t* b, CRead* read, bam_header_t* header);
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
22
0f80a5141704 version 0.3 uploaded
vipints
parents:
diff changeset
23 #endif