DNA-seq of muscle
After DNA sample extracted, a bioruptor (Covaris S220) was used to fragment sheep genomic DNA to a mean size of approximately 500 bp. After  fragmentation, DNA fragments were blunt-ended, added dA at the 3-ends, and then added adapters. The procedure was carried out in strict accordance with the Illumina manufacturer instructions. Adapters-added DNA were given a bisulfite conversion step, described in a previous study. Bisulfite-treated DNA was PCR amplified and then pair-end sequenced with 125 readlength of each end using Illumina high throughput sequencing system (Hiseq 2500)
A-B1