-
Notifications
You must be signed in to change notification settings - Fork 4
/
definition.c
executable file
·60 lines (42 loc) · 1.4 KB
/
definition.c
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
#include "utils.h"
int definition_main(int argc, char *argv[]){
if ( argc != 2 ) {
fprintf(stderr, "\nUsage: emapper-utils definition <emapper>\n\n");
return 1;
}
gzFile fp;
fp = strcmp(argv[1], "-")? gzopen(argv[1], "r") : gzdopen(fileno(stdin), "r");
if(fp){
kstream_t *ks;
ks = ks_init(fp);
int k;
kstring_t kt = {0, 0, 0};
kstring_t kv = {0, 0, 0};
printf("#gene\tdefinition\n");
while( ks_getuntil( ks, '\n', &kt, 0) >= 0){
if( kt.s[0] == '#' ) continue;
ks_tokaux_t aux;
char *p;
kv.l = 0;
for (p = kstrtok(kt.s, "\t", &aux), k = 0; p; p = kstrtok(0, 0, &aux), ++k){
if(k == 0){
kt.s[aux.p -p] = '\0';
}else if(k == 21){
p[aux.p - p] = '\0';
kputsn(p, aux.p -p , &kv);
break;
}
}
if(kv.l == 0) continue;
printf("%s\t%s\n", kt.s, kv.s);
}
free(kt.s);
free(kv.s);
ks_destroy(ks);
gzclose(fp);
}else{
fprintf(stderr, "[ERR]: can't open file %s\n", argv[1]);
exit(1);
}
return 0;
}