For cfDNA sequencing library preparation, the targeted DNA covering 0.69 Mb was enriched by designed biotinylated probes with the NimbleGen SeqCap EZ Choice (Roche) according to the manufacturer’s standard protocol. The cfDNA was end repaired, A-tailed and adaptor ligated with KAPA Hyper Prep Kit (KAPA Biosystems). Then, the 8-bp barcode was added to each sample through PCR and followed by hybridization. Then the captured DNA was amplified by PCR and pooled together for sequencing.
The DNA libraries were subjected to 150-bp paired-end sequencing on NovaSeq platform (Illumina). All sequenced reads, after low quality reads had been filtered, were aligned on GRCh37 human reference genome with BWA-MEM (version 0.7.17-r1188). Duplicated reads were deleted. Variants were called with GATK 3.7 from sequencing data of gDNA. The sequencing depth was determined for each base from sequencing data of cfDNA.