I've sampled multiple populations of a species of frog (a typical outbreeding species) and for every population I'm getting observed heterozygosities (e.g 0.66) much larger than expected heterozygosities (e.g. 0.38) and correspondingly very negative Fis values (e.g. -0.65).
I can think of no biologically realistic explanations for this degree of heterozygosity excess. Could this be an indication that the RAD data is bad--that there's something wrong with the library/sequencing? Is there something I can do to evaluate this possibility and gain trust in my data?
I performed de novo assembly in dDocent followed by the recommended filtering steps.
Thanks for any insights.