File: read.h

package info (click to toggle)
maq 0.7.1-3
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 1,360 kB
  • ctags: 993
  • sloc: cpp: 5,025; ansic: 3,329; sh: 3,282; perl: 2,547; makefile: 27
file content (48 lines) | stat: -rw-r--r-- 1,123 bytes parent folder | download | duplicates (5)
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
#ifndef LH3_READ_H
#define LH3_READ_H

#include <zlib.h>
#include "dword.hh"
#include "const.h"
#include "seq.h"

#define MA_NONAME

typedef struct __longreads_t
{
	int size_l, size_r;
	int n_reads;
	bit8_t **seq;
	char **name;
} longreads_t;

struct __match_data_t;
struct __match_aux_t;

#ifdef __cplusplus
extern "C" {
#endif
	longreads_t *ma_load_reads(gzFile fp_l, int size_l, gzFile fp_r, int size_r);
	void delete_longreads(longreads_t *lr);
	void ma_init_match_data(struct __match_data_t *d, gzFile fp_l, int *size_l, gzFile fp_r, int *size_r, int is_trim,
							const char *adapter, gzFile hits_fp);
	void maq_methy_modify(struct __match_data_t *d, struct __match_aux_t *o);
	static inline int ma_load_1read(gzFile fp, seq_t *seq, char *name)
	{
		int len;
		if (gzread(fp, &len, sizeof(int)) == 0) return -1;
		gzread(fp, name, sizeof(char) * len);
		gzread(fp, &len, sizeof(int));
		if (seq->m < len + 1) {
			seq->m = len + 1;
			seq->s = (unsigned char*)realloc(seq->s, sizeof(char) * seq->m);
		}
		seq->l = len;
		gzread(fp, seq->s, sizeof(char) * len);
		return len;
	}
#ifdef __cplusplus
}
#endif

#endif