Nanopore Sequencing of SARS-CoV-2: Comparison of Short and Long PCR-tiling Amplicon Protocols ============================================================================================= * Broňa Brejová * Kristína Boršová * Viktória Hodorová * Viktória Čabanová * Askar Gafurov * Dominika Fričová * Martina Neboháčová * Tomáš Vinař * Boris Klempa * Jozef Nosek ## Abstract Surveillance of the SARS-CoV-2 variants including the quickly spreading mutants by rapid and near real-time sequencing of the viral genome provides an important tool for effective health policy decision making in the ongoing COVID-19 pandemic. Here we evaluated PCR-tiling of short (∼400-bp) and long (∼2 and ∼2.5-kb) amplicons combined with nanopore sequencing on a MinION device for analysis of the SARS-CoV-2 genome sequences. Analysis of several sequencing runs demonstrated that using the long amplicon schemes outperforms the original protocol based on the 400-bp amplicons. It also illustrated common artefacts and problems associated with this approach, such as uneven genome coverage, variable fraction of discarded sequencing reads, as well as the reads derived from the viral sub-genomic RNAs and/or human and bacterial contamination. Keywords * SARS-CoV-2 * mutants * PCR-tiling * nanopore sequencing * multiplexing ## Introduction Massive spreading of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) within the human population began in December 2019 in Wuhan, Hubei Province, China (1-3). In the following weeks, the virus has been quickly transmitted all over the globe. As of April 2021, it infected more than 147 million humans and caused over 3.1 million deaths ([https://arcg.is/0fHmTX](https://arcg.is/0fHmTX); 4). The 29,903-nt long genomic RNA sequence of the SARS-CoV-2 strain Wuhan-Hu-1 (Genbank/RefSeq acc.nos. MN908947 / NC_045512; 1) and related isolates (2,3) were determined early in 2020 and facilitated rapid development of molecular diagnostics as well as the analysis of additional isolates from other geographical regions of the world. Currently, more than 1 million SARS-CoV-2 genome sequences are available in the GISAID repository ([http://www.gisaid.org](http://www.gisaid.org)), thus representing an unprecedented resource for the scientific community and public health officials. The virus sequences were determined using a range of experimental approaches based on metagenomics, sequence capture or enrichment, amplicon pools by deploying short (*e*.*g*., Illumina) or long-read (*e*.*g*., Pacific Biosciences, Oxford Nanopore Technologies) sequencing platforms. Of these, the nanopore sequencing becomes increasingly popular as in addition to sequencing of viral genomic RNA it also permits transcriptome mapping, characterization of sub-genomic RNA molecules, and identification of modified nucleotides in the viral genome (5-7). Rapid, cost-effective, and near real-time genome sequencing of the SARS-CoV-2 variants combined with epidemiological data provides an important resource not only for understanding the virus transmission, its genetic alterations and evolution, but also for making the policy decisions in combating the pandemic (8). Monitoring sequence diversification plays an essential role in continual refinement of molecular diagnostics (*e*.*g*., redesigning the primers for nucleic acid amplification techniques (9) or development of screening tools for variants of concerns (VoC) and those evading the immune response (10, 11)). This underscores the importance of genomic epidemiology, although the elucidation of direct links between particular mutation(s) and the virus spreading or clinical implications still represents a challenging task (12-22). Nanopore sequencing of tiled PCR-generated amplicon pools represents a powerful tool for investigating viral genomes. The protocol has been developed by the Artic Network ([https://artic.network/](https://artic.network/)) for sequencing of Ebola, Zika, and Chikungunya genomes (23, 24). In January 2020, the original protocol was promptly adjusted for rapid sequence determination of SARS-CoV-2 RNA prepared directly from clinical samples such as nasopharyngeal or oropharyngeal swabs. Additional studies described its modifications including alternative primer schemes and different amplicon sizes or different sequencing chemistries (25-36). Its further improvements resulted in simplification of the sequencing library preparation, shortened hands-on time, and increased sample multiplexing (up to 96) that decreased the reagent costs to about £10 per sample, making this approach affordable for epidemiologic surveillance of the pandemics (36). Importantly, the rigorous comparison of nanopore sequencing with Illumina short reads technology demonstrated that in spite of relatively high error rates in individual nanopore reads, the highly accurate consensus single nucleotide variant (SNV) calling with >99% sensitivity and >99% precision can be achieved with a minimum of about 60-fold coverage (37). In this study, we compare the performance of several PCR-tiling based protocols which were evaluated as part of our efforts to sequence isolates of SARS-CoV-2 from Slovakia collected between March 2020 and March 2021. Using the generated sequence data, we investigate the nature of common problems and artefacts associated with this approach. We compared the sequencing results obtained from the libraries containing multiplexed barcoded SARS-CoV-2 samples made of ∼400-bp, ∼2-kb, and ∼2.5-kb long overlapping amplicon pools as well as the combination of short and long amplicons. Our results show that sequencing of long amplicons clearly outperforms the original protocol based on shorter amplicons in terms of lower coverage variation and overall quality of the final sequence consensus. We also compare the performance of MinION runs with the standard (FLO-MIN106, FLO-MIN106D) and flongle (FLO-FLG001) flow cells differing by nominal pore counts, *i*.*e*. 2048 (split into four sets of 512 each) and 126, respectively. ## Results and Discussion ### Comparison of PCR-tiling amplification protocols The PCR-tiling amplification combined with nanopore sequencing was employed for genome sequence analysis of 152 SARS-CoV-2 isolates from Slovakia (**S1 Table**). The genome sequences were obtained using primer schemes generating either ∼400-bp (Artic Network version V3, [https://github.com/artic-network/artic-ncov2019](https://github.com/artic-network/artic-ncov2019)), ∼2-kb (35), or ∼2.5-kb long amplicons (27). In some experiments, we loaded the same sequencing libraries to both the standard and Flongle flow cells. This allowed us to compare and evaluate the MinION runs and to analyze the problems affecting the performance of nanopore sequencing. To compare primer sets for short and long amplicons and sequencing devices, three different batches (UKBA-2, UKBA-3 and UKBA-4 in **Table 1**) consisting of 10-12 multiplexed samples were sequenced using multiple strategies. **Fig 1** shows the comparison of the fraction of samples in a sequencing run successfully sequenced at various cut-offs measuring the total amount of sequencing normalized by the number of samples in the run. In both batches UKBA-2 and UKBA-3, 2-kb amplicons clearly outperform 400-bp amplicons. Sequencing of a mixture of longer and shorter amplicon pools provided comparable results to sequencing longer amplicons alone, perhaps because the mixture was enriched in the long amplicons. Finally, the Flongle and standard flow cells are similarly successful at comparable sequencing volumes. However, there are two disadvantages to using the Flongle flow cells. First, the Flongle cannot be washed and reused, its entire capacity is used for a single experiment. Second, since there is a large variance in the amount of data produced by a single Flongle flow cell (in our experiments, the number of active pores in Flongles ranged between 18 and 67 pores and produced between 110 and 830 Mbp - see Table 1), the capacity may be insufficient to completely recover sequences of 10 or more multiplexed samples. We consider as an important advantage that the runs using the standard flow cells can be terminated when sufficient data is collected, and thus these flow cells can be reused in further experiments after washing with nuclease containing buffer (*i*.*e*., EXP-WSH003 or EXP-WSH004). Moreover, the standard flow cells allow simultaneous sequencing of a greater number of barcoded samples with a longer run. View this table: [Table 1.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/T1) Table 1. Overview of the MinION sequencing runs. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F1) Fig 1. The fraction of successfully sequenced multiplexed samples over time. A sample is considered as successfully sequenced if the resulting sequence produced by the Artic pipeline has fewer than 500-bp (a) or 3-kb (b) marked as missing bases. Each run is represented by several time points, each point showing the percentage of successfully sequenced barcodes (y-axis) upon reaching a specified amount of sequenced data per barcode (x-axis). Note that batch UKBA-2 included samples with low product concentrations after PCR amplification. As a result, three samples (barcodes 02, 06 and 11) could not be completed reliably even after combining data from all six sequencing runs. Batches UKBA-3 and UKBA-4 contained only samples with Cq values from RT-qPCR below 26. **S1 Figure** shows the amount of missing sequence in individual samples plotted against possible explanatory variables, namely the Cq values, amplicon concentration, and RNA sample storage time prior to amplification. Although the expected trends are in some cases observable, they are not followed universally. Using the Artic pipeline for further analysis, sequencing reads must first pass a series of filters to ensure no barcode bleeding and to remove possible contamination. The number of reads passing these filters and used for the identification of variants in the final step of the pipeline varied between runs. In our experiments their fraction comprise between 14 and 55% (**Fig 2A**). Majority of failed reads are due to the low quality or incompleteness (groups A-C) and comprise between 41-78% of all reads. While there are no clear differences between short and long amplicon protocols, with 2-kb amplicons these low-quality reads are apparently more prevalent on the Flongle runs compared to the standard flow cells. ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F2) Fig 2. Reasons for discarding reads in the Artic pipeline. The sequencing reads must pass through a series of filters to ensure correct sample assignment and the read quality. A: reads without barcode identification. B: reads with only one barcode (Artic pipeline requires barcodes on both ends to ensure that the whole read was sequenced and to decrease the probability of barcode bleeding). C: low-quality reads (base caller quality less than 7). D: reads do not align to the SARS-CoV-2 reference. E: reads are too short (likely due to fragmentation) or F: too long (i.e. chimeric reads); lengths between 1500 and 3000 for 2-kb amplicons, between 350 and 619 for 400-bp amplicons. The reads passing all filters are included in group G. Panel (A): Summary per run. Panel (B): Detailed per-barcode analysis for UKBA-2 samples, 2-kb amplicons, standard flow cell. Interestingly, in some runs, up to 6% of reads that pass the base quality filters do not map to the target reference genome. In particular, four samples in batch UKBA-2 of 2-kb amplicon run (barcodes 02, 07, 08 and 11) have a very high fraction of non-target reads (**Fig 2B**). The majority (82-96%) of these reads map to the human genome, and a smaller fraction (0.3-9%) map to bacterial genomes, including the species colonizing human oral cavity and respiratory tract (*e*.*g*., *Actinomyces graevenitzii, Haemophilus parainfluenzae, Leptotrichia spp*., *Prevotella spp*., *Pseudomonas aeruginosa, Rothia mucilaginosa, Streptococcus pneumoniae, S. mitis, S. parasanguinis, S. salivarius, Tannerella forsythia, Veillonella parvula*). All four samples showed a lower viral load (*i*.*e*., Cq value > 30) in RT-qPCR assays, and the amplification in the PCR-tiling protocol resulted in lower product yield. Human and bacterial reads represent artefacts apparently resulting from a non-specific amplification of contaminating nucleic acids present in clinical samples. We have also observed that some amplicons originate from sub-genomic RNAs that co-purify with the SARS-CoV-2 genomic RNA. It has been demonstrated that the amount of sub-genomic RNAs correlates with the disease severity. As these molecules are strongly repressed in asymptomatic patients (38), their proportion in the sequencing data can serve as a molecular marker. The most abundant reads are derived from the N mRNA (39). The sub-genomic RNAs are generated in the process of the virus replication/transcription (5) and start with a leader sequence originating from the untranslated 5’ end of the viral genome, followed by a downstream sequence containing a particular open reading frame. The leftmost primer in both 400-bp and 2-kb primer sets investigated in this study is contained within the leader sequence. This facilitates amplification of sub-genomic RNAs with appropriate right primers (**Fig 3**). ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F3) Fig 3. Sub-genomic RNAs (black) and amplicons of primer pool 1 from the 2-kb primer set (blue), visualized by the UCSC genome browser (40). **Table 2** lists the fraction of selected sub-genomic RNAs among reads that could be aligned to the SARS-CoV-2 genome. These fractions are relatively low, with the remaining sub-genomic RNAs being even more rare. However, the fractions vary among the samples. In UKBA-2 run with 2-kb amplicons, the highest fraction of 14.3% was observed for the gene N mRNA in barcode 07 and the fraction of 7.5% was observed for the ORF3a mRNA in barcode 11. Some of these sub-genomic amplicons are discarded from the analysis as too short, while others lead to lower coverage in areas not covered by the sub-genomic RNA (**Fig 4**). View this table: [Table 2.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/T2) Table 2. Sub-genomic RNA fractions (in %) among different MinION runs with the standard flow cells (A) and in different barcodes of batch UKBA-2, 2-kb amplicons (B). Only reads that align to the SARS-CoV-2 genome and can be demultiplexed were considered. Only genes with the highest numbers of sub-genomic RNA reads are shown. ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F4) Fig 4. Coverage along the genome in two MinION runs for batch UKBA-2. In both runs, an initial portion of the run containing on average 40-Mbp of sequencing data per barcode was used. Coverage values higher than 1000 were clipped at this value and are shown in blue. Coverage below 20 (default Artic cutoff) is shown in red. Medians of 10-bp windows are shown for smoothing. The very starts and ends of the genome are not covered by amplicons, and are thus displayed in red. Shaded area in the left column corresponds to amplicon 13. Some barcodes have a visible dip in the coverage at the left end of this amplicon; this difference in coverage is caused by reads originating from sub-genomic RNAs corresponding to the gene S. Similar plots for additional runs are shown in **S2 Fig**. From these pilot experiments, we conclude that even though 400-bp amplicons have a lower percentage of discarded reads (**Fig 2**), they produce fewer finished sequences at a comparable overall amount of sequence data (**Fig 1**). The reason is a very uneven coverage of individual amplicons (**Fig 4**). This is observed in both sets of primers, but for the 400-bp amplicons we see a much lower coverage in the worst covered regions (**Fig 5**). Additional sequencing runs (UKBA-6, UKBA-10, UKBA-11, and UKBA-12) were performed with long 2-kb amplicons on standard MinION flow cells with similar results (**Fig 2A; S2 Fig**). ![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F5) Fig 5. Coverage distribution in different sequencing runs. For each barcode, coverage by reads passing the Artic filter was computed along the genome and the distribution of the coverage values was summarized as a violin plot (blue), cropped at coverage 1000. Orange dots represent median coverage and green dots 10% percentile (corresponding to a approx. 3-kb portion of the genome with the lowest coverage). In all runs, an initial portion containing on average 40-Mbp of sequencing data per barcode was used. To investigate if a different primer scheme for generating long amplicons can solve the problem with uneven coverage (in particular, see amplicon 13 in **Fig 4** which partially covers the S gene region important for identification of the SARS-CoV-2 Variants of Concerns), we also tested the 2.5-kb primer panel (27). Except for the leftmost primer, the primer positions in this panel differ from those of the 2-kb scheme. We have performed two sequencing runs with 2.5-kb primer set (UKBA-19, UKBA-21). In the first experiment, we have noticed an almost complete drop of coverage in the last amplicon derived from the 3’ end of the genome; for the second experiment, we have replaced the primers for the right-most amplicon with the right-most primers from the 2-kb panel, which mitigated the issue. Comparing the coverage of individual amplicons between the 2-kb and 2.5-kb schemes (**Fig 6**), the coverage in the 2.5-kb scheme indeed appears to be more even. **Fig 5** illustrates that our modification of the 2.5-kb scheme leads to a particularly small difference between the median coverage and coverage of the lowest 10% of the genome, which may result in fewer regions with insufficient coverage. However, we have also noticed a higher percentage of failed reads, with only 24% (UKBA-19) and 16% (UKBA-21) reads passing all filters and being usable for variant identification. Further analysis revealed a notable increase in single-barcode reads (group B) and shorter than expected reads (group E), pointing to difficulties in amplifying and sequencing longer fragments. More experiments are required to determine whether the 2.5-kb scheme results in more fully-assembled genomes over the 2-kb scheme. ![Fig 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/05/13/2021.05.12.21256693/F6.medium.gif) [Fig 6.](http://medrxiv.org/content/early/2021/05/13/2021.05.12.21256693/F6) Fig 6. Average coverage along the genome for seven runs with 2-kb amplicons (batches UKBA-2,3,4,6,10,11,12) and two runs with 2.5-kb amplicons (UKBA-19 with the original primer set and UKBA-21 with the last primer pair replaced by its counterpart from the 2-kb scheme). Each line depicts the average coverage over all samples in a run at the time point when 40-Mbp per sample was sequenced on average. Medians of 50-bp windows are shown for smoothing. Note a drop-out in the amplicon 13 (2-kb scheme) which covers a 3’ end of orf1b and about a third of the S gene including the region associated with mutations in Variants of Concern such as B.1.1.7. ## Conclusions In this paper, we have compared three versions of PCR-tiling protocol for sequencing SARS-CoV-2 genomes from clinical samples on MinION platform. Our results have shown that even though the protocol based on short 400-bp amplicons generally produces more usable data, the coverage of individual amplicons varies widely which may result in difficulties in recovering individual mutations in under-represented amplicons. In comparison, longer amplicons tend to produce close-to-finished genomes more quickly, generally requiring smaller amounts of raw data produced per barcode sequenced. However, protocols based on long amplicons produce a higher percentage of reads that are unsuitable for further analysis with the Artic pipeline, likely due to a combination of fragmentation of synthesized molecules and prematurely aborted molecules during sequencing. The longer amplicon protocols are also less suitable for applications, where original RNA molecules in clinical samples may already be fragmented. Generally, the Flongle flow cells performed worse in sequencing multiplexed libraries containing barcoded samples than regular MinION flow cells, which have an added advantage of ability to adjust the length of the run based on the library and individual sample quality. Interestingly, PCR-tiling protocols were able to also pick up sub-genomic RNA transcripts, and the proportion of these transcripts varied between samples. Since increased levels of sub-genomic transcripts are correlated with severe cases of COVID-19, these protocols could be optimized to detect the levels of sub-genomic transcripts more accurately and used as a biomarker for disease severity. It is evident that effective epidemiologic surveillance of the pandemic is strongly dependent on systematic sequencing of SARS-CoV-2 isolates. The MinION platform from Oxford Nanopore Technologies is one of the most powerful and versatile means for acquisition of viral sequences. Yet, as demonstrated in this study, the pros and cons of a particular protocol must be taken into account to ensure that the sequencing results will be of a highest quality, which is an essential prerequisite for their utility in fighting the pandemic. ## Materials and Methods ### Collection of samples and RNA preparation Oropharyngeal swabs of patients with suspected COVID-19, collected between March 30, 2020 and March 19, 2021, were preserved in 2-3 ml of viral transport medium and delivered to the laboratory of the Biomedical Research Centre of the Slovak Academy of Sciences in Bratislava, Slovakia in frame of the routine RT-qPCR diagnostics for SARS-CoV-2. Initially (UKBA-2 samples), 100 µl of the swab medium was used for the RNA extraction using the Zymo Research Quick-RNA™ Viral 96 Kit (Zymo Research, Irvin, California, USA). Resulting RNA was eluted to 20 µl of nuclease free water. For all other specimens, the Biomek i5 Automated Workstation (Beckman Coulter, Indianapolis, Indiana, USA) was employed using the RNAdvanced Viral kit (Beckman Coulter, Indianapolis, Indiana, USA). In this case, RNA was extracted from 200 µl of swab medium and eluted to 40 µl of nuclease free water. ### Real-time quantitative PCR (RT-qPCR) In frame of the routine RT-qPCR diagnostics, presence of SARS-CoV-2 RNA was detected by vDETECT COVID-19 RT-qPCR kit, rTEST COVID-19 RT-qPCR kit or rTEST COVID-19 RT-qPCR ALLPLEX kit (MultiplexDX, Bratislava Slovakia) targeting RNA-dependent RNA polymerase (RdRp) and Envelope (E) genes. The RT-qPCR assays were carried on QuantStudio(tm) 5 Real-Time PCR System (Applied Biosystem, Foster City, California, USA). ### Library preparation and DNA sequencing The sequencing libraries were constructed using a ligation kit (SQK-LSK109) essentially as described in a PCR-tiling of COVID-19 virus protocol (PTC\_9096_v109_revF_06Feb2020; Oxford Nanopore Technologies) with minor modifications. Briefly, RNA samples extracted from swabs positive for the presence of SARS-CoV-2 in RT-qPCR assay (quantification cycle (Cq) values 13.46 - 32.03; **S1 Table**) were converted into cDNA using a SuperScript IV reverse transcriptase (Thermo Fisher Scientific) or LunaScript® RT SuperMix Kit (New England Biolabs). For each sample, the overlapping amplicons were generated using a Q5® Hot Start High-Fidelity DNA polymerase (New England Biolabs) and the primer pools spanning the SARS-CoV-2 genome sequence (*i*.*e*., 400-bp Artic nCoV-2019 V3 panel ([https://github.com/artic-network/artic-ncov2019](https://github.com/artic-network/artic-ncov2019)) purchased from Integrated DNA Technologies (IDT, cat.no. 10006788) and the 2-kb (35) and 2.5-kb schemes (27), custom synthesized by Microsynth). The same cycling program was used for all amplicon types (*i*.*e*., 30 sec initial denaturation at 98°C, followed by 25 to 35 cycles of 15 sec at 98°C (denaturation) and 5 min at 65°C (combined annealing and polymerization), and cooling to 4°C). The amplifications were performed in two separate reactions and the overlapping amplicons were pooled, purified using an equal volume of AMPure XP magnetic beads (Beckman Coulter) and quantified using a Qubit 3.0 spectrophotometer and dsDNA Broad Range Assay Kit (Thermo Fisher Scientific). About 50-75 ng (400-bp amplicons) and 250-300 ng (2 and 2.5-kb amplicons) of each SARS-CoV-2 isolate were treated with NEBNext Ultra II End repair / dA-tailing Module (New England Biolabs). The samples were then barcoded using EXP-NBD104 (barcodes 1-12) or EXP-NBD114 (barcodes 13-24) kits (Oxford Nanopore Technologies) and NEBNext Ultra II Ligation Master Mix (New England Biolabs). Barcoded samples were pooled and purified using 0.6 volume of AMPure XP magnetic beads. The AMII sequencing adapter (Oxford Nanopore Technologies) was ligated to about 75 ng (400-bp amplicons) or 300 ng (2 and 2.5-kb amplicons) of barcoded pools using Quick T4 DNA ligase (New England Biolabs) and the sequencing libraries were purified using 0.6 volume of AMPure XP magnetic beads. About 20 ng (400-bp amplicons) and 90 ng (2 and 2.5-kb amplicons) of the libraries were loaded on an R9.4.1 flow cell (FLO-MIN106). The sequencing was performed using a MinION Mk-1b device (Oxford Nanopore Technologies). For sequencing on the Flongle flow cells (FLO-FLG001), the library preparation was the same, except that one third to one half of the library was loaded compared to the amount used for the standard flow cell. ### Data processing Nanopore sequencing data were base called and demultiplexed using Guppy v.3.4.4. Variant analysis was performed using Artic analysis pipeline v.1.1.3. ([https://artic.network/ncov-2019/ncov2019-bioinformatics-sop.html](https://artic.network/ncov-2019/ncov2019-bioinformatics-sop.html)) using recommended settings. Minimum and maximum read lengths in the Artic guppyplex filter were set to 350 and 619 for the 400-bp amplicons and 1500 and 3000 for both the 2 and 2.5-kb amplicons, respectively. For batch UKBA-2, the final sequences were produced by first combining sequencing reads from both standard and Flongle runs with the same primer set and running the Artic pipeline. Subsequently the results for the two primer sets were combined so that regions sufficiently covered by at least one amplicon set were considered as finished. The same process was used in batch UKBA-3, but there was only data from standard flow cells available. Subsequent batches were based on 2 or 2.5-kb amplicons sequenced on a standard flow cell. To compare different primer sets and sequencing devices, reads were also demultiplexed at the less strict default Guppy settings and aligned to various reference genomes by minimap2 v. 2.13-r852-dirty (41). Reference genomes include the SARS-CoV-2 genome MN908947.3 (1), the human genome version hg19 downloaded from the UCSC genome browser (40), and the database for bacterial species typing included in the Japsa software (42). To detect sub-genomic RNAs, reads were aligned to transcripts downloaded from the UCSC genome browser by minimap2, and classified as sub-genomic, if the alignment to a sub-genomic RNA has at least 5 matches more than the best alignment to the reference genome. Read coverage was computed using genomecov tool from BEDTools (43) with options -bga -split. To compare the results for various sequencing data volumes, reads were ordered by the sequencing finish time and the initial portion with the desired total length was selected and used for the analysis in the Artic pipeline. To compare batches with a different number of samples, the cutoffs were expressed as the average amount per barcode. ### Data availability The SARS-CoV-2 genome sequences determined in this study were deposited in the GISAID database ([http://www.gisaid.org](http://www.gisaid.org)). The accession IDs are shown in **S1 Table**. The sequencing reads mapping to the reference genome were deposited to ENA under study PRJEB44303. ## Supporting information Supporting Information [[supplements/256693_file02.pdf]](pending:yes) ## Data Availability The SARS-CoV-2 genome sequences determined in this study were deposited in the GISAID database ([http://www.gisaid.org](http://www.gisaid.org)). The accession IDs are shown in S1 Table. The sequencing reads mapping to the reference genome were deposited to ENA under study PRJEB44303. ## Ethics statement All clinical specimens used within this study were previously collected for the purpose of primary diagnosis of SARS-CoV-2 and were transferred to the Biomedical Research Center of the Slovak Academy of Sciences, Bratislava, Slovakia while made unidentifiable for the researchers performing this study. The study has been approved by the Ethics committee of Biomedical Research Center of the Slovak Academy of Sciences, Bratislava, Slovakia (Ethics committee statement No. EK/BmV-02/2020). ## Author contributions Conceived and designed the experiments: BB, JN, TV. Performed the experiments: KB, VC, DF, VH, JN. Analyzed the data: BB, AG, TV. Wrote the first draft of manuscript: BB, JN, TV. Edited and approved the manuscript: BB, KB, VC, DF, AG, VH, BK, MN, JN, TV. Obtained funding: BB, BK, MN, JN, TV. ## Competing interests The authors declare that they have no conflict of interest. ## Acknowledgements The authors wish to thank Lubomir Tomaska (Comenius University in Bratislava) for critical reading of the manuscript and discussions. The project was supported by grants from the Slovak Research and Development Agency (APVV-18-0239 to JN, PP-COVID-20-0017 to BK), the Scientific Grant Agency (VEGA 1/0463/20 to BB, VEGA 1/0458/18 to TV, VEGA 1/0027/19 to JN, VEGA 1/0136/20 to MN), and the European Union’s Horizon 2020 research and innovation program (EVA-GLOBAL project #871029 to BK and PANGAIA project #872539). The research was also supported in part by the Operation Program of Integrated Infrastructure (OPII) projects ITMS2014: 313011ATL7 and ITMS2014+: 313021×329 (Advancing University Capacity and Competence in Research, Development and Innovation), co-financed by the European Regional Development Fund. * Received May 12, 2021. * Revision received May 12, 2021. * Accepted May 13, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG, et al. A new coronavirus associated with human respiratory disease in China. Nature 2020;579(7798): 265–269. doi: 10.1038/s41586-020-2008-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2008-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32015508&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 2. 2.Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020;382(8): 727–733. doi: 10.1056/NEJMoa2001017 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31978945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 3. 3.Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020;579(7798): 270–273. doi: 10.1038/s41586-020-2012-7 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2012-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32015507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 4. 4.Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020;20(5): 533–534. doi: 10.1016/S1473-3099(20)30120-1 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32087114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 5. 5.Kim D, Lee JY, Yang JS, Kim JW, Kim VN, Chang H. The architecture of SARS-CoV-2 transcriptome. Cell 2020;181(4): 914–921.e10. doi: 10.1016/j.cell.2020.04.011 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.04.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32330414&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 6. 6.Taiaroa G, Rawlinson D, Featherstone L, Pitt M, Caly L, Druce J, et al. Direct RNA sequencing and early evolution of SARS-CoV-2. bioRxiv 2020.03.05.976167 [Preprint]. 2020 [posted 2020 Apr 03; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.03.05.976167](https://doi.org/10.1101/2020.03.05.976167) 7. 7.Miladi M, Fuchs J, Maier W, Weigang S, Pedrosa NDi, Weiss L, et al. The landscape of SARS-CoV-2 RNA modifications. bioRxiv 2020.07.18.204362 [Preprint]. 2020 [posted 2020 Jul 8; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.07.18.204362](https://doi.org/10.1101/2020.07.18.204362) 8. 8.Oude Munnink BB, Nieuwenhuijse DF, Stein M, O’Toole Á, Haverkate M, Mollers M, et al. Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands. Nat Med. 2020;26(9): 1405–1410. doi: 10.1038/s41591-020-0997-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0997-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32678356&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 9. 9.Nayar G, Seabolt E, Kunitomi M, Agarwal A, Beck KL, Mukherjee V, et al. Analysis and forecasting of global RT-PCR primers for SARS-CoV-2. bioRxiv 2020.12.26.424429 [Preprint]. 2021 [posted 2021 Jan 07; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.12.26.424429](https://doi.org/10.1101/2020.12.26.424429) 10. 10.Bal A, Destras G, Gaymard A, Stefic K, Marlet J, Eymieux S, et al. COVID-Diagnosis HCL Study Group. Two-step strategy for the identification of SARS-CoV-2 variant of concern 202012/01 and other variants with spike deletion H69-V70, France, August to December 2020. Euro Surveill. 2021;26(3): 2100008. doi: 10.2807/1560-7917.ES.2021.26.3.2100008 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2021.26.3.2100008&link_type=DOI) 11. 11.KováČová V, Boršová K, Paul ED, Radvánszka M, Hajdu R, Čabanová V, et al. A novel, multiplexed RT-qPCR assay to distinguish lineage B.1.1.7 from the remaining SARS-CoV-2 lineages. medRxiv 2021.02.09.21251168 [Preprint]. 2021 [posted 2021 Mar 3; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2021.02.09.21251168](https://doi.org/10.1101/2021.02.09.21251168) 12. 12.Choi B, Choudhary MC, Regan J, Sparks JA, Padera RF, Qiu X, et al. Persistence and evolution of SARS-CoV-2 in an immunocompromised host. N Engl J Med. 2020;383(23): 2291–2293. doi: 10.1056/NEJMc2031364 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2031364&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33176080&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 13. 13.Eskier D, Karakülah G, Suner A, Oktay Y. RdRp mutations are associated with SARS-CoV-2 genome evolution. Peer J. 2020;8: e9587. doi: 10.7717/peerj.9587 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7717/peerj.9587&link_type=DOI) 14. 14.Hodcroft EB, Zuber M, Nadeau S, Crawford KHD, Bloom JD, Veesler D, et al. Emergence and spread of a SARS-CoV-2 variant through Europe in the summer of 2020. medRxiv 2020.10.25.20219063 [Preprint]. 2021 [posted 2021 Mar 24; cited 2021 Apr 14]. Available from: [https://www.medrxiv.org/content/10.1101/2020.10.25.20219063v3](https://www.medrxiv.org/content/10.1101/2020.10.25.20219063v3) 15. 15.Hou YJ, Chiba S, Halfmann P, Ehre C, Kuroda M, Dinnon KH 3rd, et al. SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo. Science 2020; 370(6523): 1464–1468. doi: 10.1126/science.abe8499 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNzAvNjUyMy8xNDY0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDUvMTMvMjAyMS4wNS4xMi4yMTI1NjY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. 16.Kemp SA, Collier DA, Datir RP, Ferreira Iatm, Gayed S, Jahun A, et al. SARS-CoV-2 evolution during treatment of chronic infection. Nature 2021; 592(7853): 277–282. doi: 10.1038/s41586-021-03291-y [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03291-y&link_type=DOI) 17. 17.Kemp SA, Harvey WT, Datir RP, Collier DA, Ferreira Iatm, Carabelli AM, Robertson DL, et al. Recurrent emergence and transmission of a SARS-CoV-2 Spike deletion ΔH69/V70. bioRxiv 2020.12.14.422555 [Preprint]. 2021 [posted 2021 Mar 08; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.12.14.422555](https://doi.org/10.1101/2020.12.14.422555) 18. 18.Lemieux JE, Siddle KJ, Shaw BM, Loreth C, Schaffner SF, Gladden-Young A, et al. Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events. Science 2021;371(6529): eabe3261. doi: 10.1126/science.abe3261 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNzEvNjUyOS9lYWJlMzI2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA1LzEzLzIwMjEuMDUuMTIuMjEyNTY2OTMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 19. 19.Plante JA, Liu Y, Liu J, Xia H, Johnson BA, Lokugamage KG, et al. Spike mutation D614G alters SARS-CoV-2 fitness. Nature 2021;592(7852): 116–121. doi: 10.1038/s41586-020-2895-3 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2895-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33106671&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 20. 20.Thomson EC, Rosen LE, Shepherd JG, Spreafico R, da Silva Filipe A, Wojcechowskyj JA, et al. Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity. Cell 2021;184(5): 1171–1187.e20. doi: 10.1016/j.cell.2021.01.037 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2021.01.037&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33621484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 21. 21.Velazquez-Salinas L, Zarate S, Eberl S, Gladue DP, Novella I, Borca MV. Positive selection of ORF1ab, ORF3a, and ORF8 genes drives the early evolutionary trends of SARS-CoV-2 during the 2020 COVID-19 pandemic. Front Microbiol. 2020;11: 550674. doi: 10.3389/fmicb.2020.550674 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fmicb.2020.550674&link_type=DOI) 22. 22.Weisblum Y, Schmidt F, Zhang F, DaSilva J, Poston D, Lorenzi JC, et al. Escape from neutralizing antibodies by SARS-CoV-2 spike protein variants. Elife 2020;9: e61312. doi: 10.7554/eLife.61312 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.61312&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33112236&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 23. 23.Quick J, Loman NJ, Duraffour S, Simpson JT, Severi E, Cowley L, et al. Real-time, portable genome sequencing for Ebola surveillance. Nature 2016;530(7589): 228–232. doi: 10.1038/nature16996 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature16996&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26840485&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 24. 24.Quick J, Grubaugh ND, Pullan ST, Claro IM, Smith AD, Gangavarapu K, et al. Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples. Nat Protoc. 2017;12(6): 1261–1276. doi: 10.1038/nprot.2017.066 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nprot.2017.066&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28538739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 25. 25.Baker DJ, Kay GL, Aydin A, Le-Viet T, Kay GL, Rudder S, et al. CoronaHiT: high-throughput sequencing of SARS-CoV-2 genomes. Genome Med. 2021;13: 21. doi: 10.1186/s13073-021-00839-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-021-00839-5&link_type=DOI) 26. 26.Charre C, Ginevra C, Sabatier M, Regue H, Destras G, Brun S, et al. Evaluation of NGS-based approaches for SARS-CoV-2 whole genome characterisation. Virus Evol. 2020;6(2): veaa075. doi: 10.1093/ve/veaa075 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/veaa075&link_type=DOI) 27. 27.Eden J-S, Sim E. SARS-CoV-2 Genome Sequencing Using Long Pooled Amplicons on Illumina Platforms. 2020. doi: 10.17504/protocols.io.befyjbpw [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.17504/protocols.io.befyjbpw&link_type=DOI) 28. 28.Freed NE, Vlková M, Faisal MB, Silander OK. Rapid and inexpensive whole-genome sequencing of SARS-CoV-2 using 1200 bp tiled amplicons and Oxford Nanopore rapid barcoding. Biol Methods Protoc. 2020;5(1): bpaa014. doi: 10.1093/biomethods/bpaa014 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomethods/bpaa014&link_type=DOI) 29. 29.Gohl DM, Garbe J, Grady P, Daniel J, Watson RHB, Auch B, et al. A rapid, cost-effective tailed amplicon method for sequencing SARS-CoV-2. BMC Genomics 2020; 21(1): 863. doi: 10.1186/s12864-020-07283-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12864-020-07283-6&link_type=DOI) 30. 30.González-Recio O, Gutiérrez-Rivas M, Peiró-Pastor R, Aguilera-Sepúlveda P, Cano-Gómez C, Jiménez-Clavero MÁ, et al. Sequencing of SARS-CoV-2 genome using different Nanopore chemistries. Appl Microbiol Biotechnol. 2021 Apr 1;1–10. doi: 10.1007/s00253-021-11250-w [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00253-021-11250-w&link_type=DOI) 31. 31.Itokawa K, Sekizuka T, Hashino M, Tanaka R, Kuroda M. Disentangling primer interactions improves SARS-CoV-2 genome sequencing by multiplex tiling PCR. PLoS One 2020; 15(9): e0239403. doi: 10.1371/journal.pone.0239403 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0239403&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32946527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 32. 32.Moore SC, Penrice-Randal R, Alruwaili M, Dong X, Pullan ST, Carter D, et al. Amplicon based MinION sequencing of SARS-CoV-2 and metagenomic characterisation of nasopharyngeal swabs from patients with COVID-19. medRxiv 2020.03.05.20032011 [Preprint]. 2020 [posted 2020 Mar 08; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.03.05.20032011](https://doi.org/10.1101/2020.03.05.20032011) 33. 33.Nasir JA, Kozak RA, Aftanas P, Raphenya AR, Smith KM, Maguire F, et al. A Comparison of Whole Genome Sequencing of SARS-CoV-2 Using amplicon-based sequencing, random hexamers, and bait capture. Viruses 2020;12(8): 895. doi: 10.3390/v12080895 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/v12080895&link_type=DOI) 34. 34.Paden CR, Tao Y, Queen K, Zhang J, Li Y, Uehara A, et al. Rapid, sensitive, full-genome sequencing of severe acute respiratory syndrome coronavirus 2. Emerg Infect Dis. 2020;26(10): 2401–2405. doi: 10.3201/eid2610.201800 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2610.201800&link_type=DOI) 35. 35.Resende PC, Motta FC, Roy S, Appolinario L, Fabri A, Xavier J, et al. SARS-CoV-2 genomes recovered by long amplicon tiling multiplex approach using nanopore sequencing and applicable to other sequencing platforms. bioRxiv 2020.04.30.069039 [Preprint]. 2020 [posted 2020 May 01; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.04.30.069039](https://doi.org/10.1101/2020.04.30.069039) 36. 36.Tyson JR, James P, Stoddart D, Sparks N, Wickenhagen A, Hall G, et al. Improvements to the ARTIC multiplex PCR method for SARS-CoV-2 genome sequencing using nanopore. bioRxiv 2020.09.04.283077 [Preprint]. 2020 [posted 2020 Sep 04; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.09.04.283077](https://doi.org/10.1101/2020.09.04.283077) 37. 37.Bull RA, Adikari TN, Ferguson JM, Hammond JM, Stevanovski I, Beukers AG, et al. Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis. Nat Commun. 2020;11(1): 6272. doi: 10.1038/s41467-020-20075-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-20075-6&link_type=DOI) 38. 38.Wong CH, Ngan CY, Goldfeder RL, Idol J, Kuhlberg C, Maurya R, et al. Subgenomic RNAs as molecular indicators of asymptomatic SARS-CoV-2 infection. bioRxiv 2021.02.06.430041 [Preprint] 2021 [posted 2021 Feb 06; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2021.02.06.430041](https://doi.org/10.1101/2021.02.06.430041) 39. 39.Parker MD, Lindsey BB, Leary S, Gaudieri S, Chopra A, Wyles M, et al. periscope: sub-genomic RNA identification in SARS-CoV-2 ARTIC Network nanopore sequencing data. bioRxiv 2020.07.01.181867 [Preprint]. 2020 [posted 2020 Nov 06; cited 2021 Apr 14]. Available from: [https://doi.org/10.1101/2020.07.01.181867](https://doi.org/10.1101/2020.07.01.181867) 40. 40.Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6): 996–1006. doi: 10.1101/gr.229102 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjEyLzYvOTk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDUvMTMvMjAyMS4wNS4xMi4yMTI1NjY5My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 41. 41.Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 2018;34: 3094–3100. doi:10.1093/bioinformatics/bty191 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty191&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29750242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) 42. 42.Cao MD, Ganesamoorthy D, Elliott AG, Zhang H, Cooper MA, Coin LJM. Streaming algorithms for identification of pathogens and antibiotic resistance potential from real-time MinION(TM) sequencing. Gigascience 2016;5(1): 32. doi: 10.1186/s13742-016-0137-2 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13742-016-0137-2&link_type=DOI) 43. 43.Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 2010;26(6): 841–842. doi: 10.1093/bioinformatics/btq033 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20110278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F05%2F13%2F2021.05.12.21256693.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275243500019&link_type=ISI)