annotate bwa-0.6.2/main.c @ 0:dd1186b11b3b draft

Uploaded BWA
author ashvark
date Fri, 18 Jul 2014 07:55:14 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
1 #include <stdio.h>
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
2 #include <string.h>
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
3 #include "main.h"
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
4 #include "utils.h"
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
5
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
6 #ifndef PACKAGE_VERSION
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
7 #define PACKAGE_VERSION "0.6.2-r126"
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
8 #endif
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
9
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
10 static int usage()
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
11 {
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
12 fprintf(stderr, "\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
13 fprintf(stderr, "Program: bwa (alignment via Burrows-Wheeler transformation)\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
14 fprintf(stderr, "Version: %s\n", PACKAGE_VERSION);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
15 fprintf(stderr, "Contact: Heng Li <lh3@sanger.ac.uk>\n\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
16 fprintf(stderr, "Usage: bwa <command> [options]\n\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
17 fprintf(stderr, "Command: index index sequences in the FASTA format\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
18 fprintf(stderr, " aln gapped/ungapped alignment\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
19 fprintf(stderr, " samse generate alignment (single ended)\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
20 fprintf(stderr, " sampe generate alignment (paired ended)\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
21 fprintf(stderr, " bwasw BWA-SW for long queries\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
22 fprintf(stderr, " fastmap identify super-maximal exact matches\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
23 fprintf(stderr, "\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
24 fprintf(stderr, " fa2pac convert FASTA to PAC format\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
25 fprintf(stderr, " pac2bwt generate BWT from PAC\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
26 fprintf(stderr, " pac2bwtgen alternative algorithm for generating BWT\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
27 fprintf(stderr, " bwtupdate update .bwt to the new format\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
28 fprintf(stderr, " bwt2sa generate SA from BWT and Occ\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
29 fprintf(stderr, " pac2cspac convert PAC to color-space PAC\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
30 fprintf(stderr, " stdsw standard SW/NW alignment\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
31 fprintf(stderr, "\n");
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
32 return 1;
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
33 }
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
34
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
35 void bwa_print_sam_PG()
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
36 {
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
37 printf("@PG\tID:bwa\tPN:bwa\tVN:%s\n", PACKAGE_VERSION);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
38 }
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
39
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
40 int main(int argc, char *argv[])
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
41 {
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
42 int i, ret;
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
43 double t_real;
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
44 t_real = realtime();
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
45 if (argc < 2) return usage();
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
46 if (strcmp(argv[1], "fa2pac") == 0) ret = bwa_fa2pac(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
47 else if (strcmp(argv[1], "pac2bwt") == 0) ret = bwa_pac2bwt(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
48 else if (strcmp(argv[1], "pac2bwtgen") == 0) ret = bwt_bwtgen_main(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
49 else if (strcmp(argv[1], "bwtupdate") == 0) ret = bwa_bwtupdate(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
50 else if (strcmp(argv[1], "bwt2sa") == 0) ret = bwa_bwt2sa(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
51 else if (strcmp(argv[1], "index") == 0) ret = bwa_index(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
52 else if (strcmp(argv[1], "aln") == 0) ret = bwa_aln(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
53 else if (strcmp(argv[1], "sw") == 0) ret = bwa_stdsw(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
54 else if (strcmp(argv[1], "samse") == 0) ret = bwa_sai2sam_se(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
55 else if (strcmp(argv[1], "sampe") == 0) ret = bwa_sai2sam_pe(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
56 else if (strcmp(argv[1], "pac2cspac") == 0) ret = bwa_pac2cspac(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
57 else if (strcmp(argv[1], "stdsw") == 0) ret = bwa_stdsw(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
58 else if (strcmp(argv[1], "bwtsw2") == 0) ret = bwa_bwtsw2(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
59 else if (strcmp(argv[1], "dbwtsw") == 0) ret = bwa_bwtsw2(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
60 else if (strcmp(argv[1], "bwasw") == 0) ret = bwa_bwtsw2(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
61 else if (strcmp(argv[1], "fastmap") == 0) ret = main_fastmap(argc-1, argv+1);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
62 else {
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
63 fprintf(stderr, "[main] unrecognized command '%s'\n", argv[1]);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
64 return 1;
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
65 }
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
66 err_fflush(stdout);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
67 err_fclose(stdout);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
68 if (ret == 0) {
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
69 fprintf(stderr, "[%s] Version: %s\n", __func__, PACKAGE_VERSION);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
70 fprintf(stderr, "[%s] CMD:", __func__);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
71 for (i = 0; i < argc; ++i)
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
72 fprintf(stderr, " %s", argv[i]);
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
73 fprintf(stderr, "\n[%s] Real time: %.3f sec; CPU: %.3f sec\n", __func__, realtime() - t_real, cputime());
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
74 }
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
75 return 0;
dd1186b11b3b Uploaded BWA
ashvark
parents:
diff changeset
76 }