6 #define HEX_STR "0123456789abcdef"
8 KSEQ_INIT(gzFile, gzread)
10 static void md5_one(const char *fn)
12 MD5_CTX md5_one, md5_all;
16 unsigned char unordered[16], digest[16];
18 for (l = 0; l < 16; ++l) unordered[l] = 0;
19 fp = strcmp(fn, "-")? gzopen(fn, "r") : gzdopen(fileno(stdin), "r");
21 fprintf(stderr, "md5fa: %s: No such file or directory\n", fn);
27 while ((l = kseq_read(seq)) >= 0) {
28 for (i = k = 0; i < seq->seq.l; ++i) {
29 if (islower(seq->seq.s[i])) seq->seq.s[k++] = toupper(seq->seq.s[i]);
30 else if (isupper(seq->seq.s[i])) seq->seq.s[k++] = seq->seq.s[i];
33 MD5Update(&md5_one, (unsigned char*)seq->seq.s, k);
34 MD5Final(digest, &md5_one);
35 for (l = 0; l < 16; ++l) {
36 printf("%c%c", HEX_STR[digest[l]>>4&0xf], HEX_STR[digest[l]&0xf]);
37 unordered[l] ^= digest[l];
39 printf(" %s %s\n", fn, seq->name.s);
40 MD5Update(&md5_all, (unsigned char*)seq->seq.s, k);
42 MD5Final(digest, &md5_all);
44 for (l = 0; l < 16; ++l)
45 printf("%c%c", HEX_STR[digest[l]>>4&0xf], HEX_STR[digest[l]&0xf]);
46 printf(" %s >ordered\n", fn);
47 for (l = 0; l < 16; ++l)
48 printf("%c%c", HEX_STR[unordered[l]>>4&0xf], HEX_STR[unordered[l]&0xf]);
49 printf(" %s >unordered\n", fn);
52 int main(int argc, char *argv[])
55 if (argc == 1) md5_one("-");
56 else for (i = 1; i < argc; ++i) md5_one(argv[i]);