Keep in mind that the GATK variant callers are designed to be as sensitive as possible and will therefore include many false positives, so you need to apply some filtering steps after calling to remove those false positives, as described in the GATK Best Practices. It's essentially impossible to answer your question without knowing more about how you did the variant calling in both cases, and what kind of filtering and evaluation you did on the results.
It's also important to understand that QUAL scores are calculated differently by different variant callers, so it's tricky to compare them directly. You'll get more insights from evaluating your results relative to known callsets or truth sets.
Of interest: Variant callers reporting different read depth on same alignment
Thanks! vcftools has a feature called --diff
It creates a neat outfile with the following contents;
I'm trying to find some consistent themes as to the variant calling discrepancies.
When examining the tail end of data found in the INFO column, you'll notice a difference between GATK and Freebayes:
Can anyone decipher this information?
Hello,
Have a look at the header of your vcf files. All these entrys should be described under FORMAT.
fin swimmer
How big is the difference of the number of variants between them?
One reason can be that freebayes describes multiple variants that are close together as one haplotype if they can be asigned to one allele. Whereas GATK maybe report every change seperately.
Fin swimmer