statgen / METAL

Meta-analysis of genomewide association scans
Other
39 stars 12 forks source link

What does the warning: "Bad alleles for marker '', expecting 'a/c' found 'a/g' mean? #23

Open ghost opened 2 years ago

ghost commented 2 years ago

Hi, I am getting a lot of warnings with the above text.

I found the source code origin of the warning:

 if (allele1[marker] == "")
                allele1[marker] = tokens[firstColumn],
                allele2[marker] = tokens[secondColumn],
                original_flipped[marker] = flip ? 'Y' : 'N';
            else if (tokens[firstColumn] != allele1[marker] ||
                     tokens[secondColumn] != allele2[marker]) {
                if (++badAlleles > maxWarnings) continue;
                printf("## WARNING: Bad alleles for marker '%s', expecting '%s/%s' found '%s/%s'\n",
                       (const char *) tokens[markerColumn],
                       (const char *) allele1[marker],
                       (const char *) allele2[marker],
                       (const char *) tokens[firstColumn],
                       (const char *) tokens[secondColumn]);

But I can't understand what comparison is being made. It happens when I process the very GWAS file, so it is not comparing alleles with the other GWAS results. It reads as if it is being compared to a reference build but I dont see that in the source code.

Any further info would be appreciated.

cwnag-c commented 1 year ago

I alsohad the same problem, did you solve this problem?

dbw198211 commented 1 year ago

I met this problem too, how to solve this problem?

Stacy33333 commented 4 months ago

Did anyone solve this problem? I also met this error.