Viral RNA extracted from original swab samples and virus culture using the RNeasy minikit (Qiagen) was reverse transcribed to cDNA using a LunaScript RT SuperMix kit (New England Biolabs). Primer pools targeting SARS-CoV-2 were designed using the PrimalScheme tool (66 (link)), and PCR was conducted using PhusionFlash PCR master mix (Themo Fisher). Sequencing libraries were prepared using an NEBNext Ultra II FS DNA library kit (New England Biolabs) according to the manufacturer’s instructions and were sequenced using an Illumina Miseq with a v3 sequencing kit. Raw sequence reads were trimmed, and low quality (quality score of <30) and short (<25 nt) sequences were removed using Trimmomatic (67 (link)). The trimmed sequence reads were assembled to the reference sequence (NC_045512.2) using the BWA-MEM (68 ) algorithm implemented in SAMTools version 1.8 (69 (link)). Mutation frequencies in virus populations were estimated based on minority variant calling with LoFreq (70 (link)).
Free full text: Click here