From the VCF, Plink was used to generate .ped and .map files. (http://pngu.mgh.harvard.edu/purcell/plink/) [58 (link)]. To detect and characterize the stretches of heterozygosity the package “detectRUNS” in R was then used. (https://github.com/bioinformatics-ptp/detectRUNS/tree/master/detectRUNS). We used the function slidingRuns.run with the following parameters: WindowSize=10, threshold=0.05, RoHet=True, minDensity=1/100, rest as default.
Free full text: Click here