This may reflect: the inherent difficulties of using NGS bioinformatics approaches to achieve accurate variant annotation in repetitive sequences, such as polyAla stretches; the high GC content of the FOXE1 gene (~73%), which could lead to a low sequencing depth and coverage (i.e., 51.7%) [13]; and/or the inclusion of CH patients regardless of their thyroid phenotype [6,7,13]. This evidence concerns the gene FOXE1 and thyroid gland disorder.