Early transmission patterns of coronavirus disease 2019 (COVID-19) in travellers from Wuhan to Thailand, January 2020 - Amazon S3
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Rapid communication Early transmission patterns of coronavirus disease 2019 (COVID-19) in travellers from Wuhan to Thailand, January 2020 Pilailuk Okada¹ , Rome Buathong² , Siripaporn Phuygun¹ , Thanutsapa Thanadachakul¹ , Sittiporn Parnmen¹ , Warawan Wongboot¹ , Sunthareeya Waicharoen¹ , Supaporn Wacharapluesadee³ , Sumonmal Uttayamakul² , Apichart Vachiraphan2 , Malinee Chittaganpitch¹ , Nanthawan Mekha¹ , Noppavan Janejai¹ , Sopon Iamsirithaworn² , Raphael TC Lee⁴ , Sebastian Maurer- Stroh4,5 1. Department of Medical Sciences, Ministry of Public Health, Thailand 2. Department of Disease Control, Ministry of Public Health, Thailand 3. Thai Red Cross Emerging Infectious Diseases - Health Science Centre, Chulalongkorn University, Thailand 4. Bioinformatics Institute, Agency for Science Technology and Research, Singapore 5. Department of Biological Sciences, National University of Singapore, Singapore Correspondence: Pilailuk Okada (pilailuk.o@dmsc.mail.go.th) Citation style for this article: Okada Pilailuk , Buathong Rome , Phuygun Siripaporn , Thanadachakul Thanutsapa , Parnmen Sittiporn , Wongboot Warawan , Waicharoen Sunthareeya , Wacharapluesadee Supaporn , Uttayamakul Sumonmal , Vachiraphan Apichart , Chittaganpitch Malinee , Mekha Nanthawan , Janejai Noppavan , Iamsirithaworn Sopon , Lee Raphael TC , Maurer-Stroh Sebastian . Early transmission patterns of coronavirus disease 2019 (COVID-19) in travellers from Wuhan to Thailand, January 2020. Euro Surveill. 2020;25(8):pii=2000097. https://doi.org/10.2807/1560-7917.ES.2020.25.8.2000097 Article submitted on 06 Feb 2020 / accepted on 27 Feb 2020 / published on 27 Feb 2020 We report two cases of coronavirus disease 2019 arrival gate was 38.6 °C by thermoscanner, and con- (COVID-19) in travellers from Wuhan, China to firmed with a tympanic thermometer. After being inter- Thailand. Both were independent introductions on viewed by quarantine officers, she was transferred to separate flights, discovered with thermoscanners Bamrasnaradura Infectious Disease Institute (BIDI) and confirmed with RT-PCR and genome sequencing. Hospital, Nonthaburi, for isolation and laboratory inves- Both cases do not seem directly linked to the Huanan tigations. She reported a runny nose and sore throat Seafood Market in Hubei but the viral genomes are but no dyspnea or diarrhoea. Upon admission, her vital identical to four other sequences from Wuhan, sug- signs were normal except for elevated blood pressure. gesting early spread within the city already in the first Her physical examination was unremarkable including week of January. inconspicuous lung sounds. Her complete blood count suggested a viral infection from relatively decreased In late December 2019, an outbreak with an initially neutrophil (48%; norm: 35–75%) to lymphocyte (40%; undiagnosed pneumonia was reported in the city of norm: 20–59%) ratio [8]. The chest X-ray (CXR) results Wuhan, Hubei Province, China, and linked to the Huanan on 8 January were compatible with pneumonia with Seafood Market [1,2]. The causative pathogen was mild thickening lung marking at both lower lungs pos- identified as a novel betacoronavirus within the severe sibly because of crowded lung rather than infiltration. acute respiratory syndrome (SARS) coronavirus (CoV) It also showed borderline cardiomegaly. Repeat CXR family, recently termed SARS-CoV-2 [3-7]. In response after 7 days showed no remarkable changes. to the outbreak, several countries including Thailand, established thermal screening at the airport for travel- In the interview, the patient explicitly stated that she lers from Wuhan since 3 January. On 8 January and 13 had not visited the Huanan Seafood Market but she January, suspected cases of infection with SARS-CoV-2 was living near another local market, which she visited were identified at Bangkok Suvarnabhumi airport. We daily until illness onset. This market was at ca 7.5 km report the investigation, basic clinical characteristics distance from the Huanan Seafood Market. The patient and viral genomes derived from these cases. also reported that she had not purchased live animals and or visited stalls with live animals. She did not visit Case 1 Jinyintai Hospital or other hospitals in Wuhan. However, A woman in her early 60s from Wuhan developed a she visited local dispensaries in Wuhan to obtain medi- fever with chills, sore throat and headache on 5 January cation. She also reported no contact with persons with 2020. She went to a local health facility in China and respiratory symptoms. received undisclosed medication. On 8 January 2020, she took a direct ca 4 h flight to Thailand from Wuhan, Clinical specimens collected on admission included with five family members, as part of a tour group of 16 the upper respiratory tract secretions and sputum. (including the case). Her measured temperature at the These specimens tested positive on 12 January for the www.eurosurveillance.org 1
Figure 1 Phylogenetic trees of Thai sequences in context of all coronavirus families (A) and structural mapping of mutations in the spike glycoprotein between SARS CoV (PDB:6CG [12]) and the current SARS-CoV-2 using YASARA [20] (B) A. B. 100 Bovine respiratory coronavirus AH187/NC 012948/1-30969 100 Bovine coronavirus/NC 003045/1-31028 100 Bovine respiratory coronavirus bovine/US/OH-440-TC/1996/NC 012949/1-30953 100 Human enteric coronavirus strain 4408/NC 012950/1-30953 100 Human coronavirus OC43/NC 005147/1-30738 100 Porcine hemagglutinating encephalomyelitis virus/NC 007732/1-30480 Equine coronavirus/NC 010327/1-30992 100 Rabbit coronavirus HKU14/NC 017083/1-31100 Rat coronavirus Parker/NC 012936/1-31250 100 Murine hepatitis virus strain A59/NC 001846/1-31357 100 100 Murine hepatitis virus strain JHM/NC 006852/1-31526 94 Human coronavirus HKU1/NC 006577/1-29926 Bat coronavirus HKU9-1/NC 009021/1-29114 100 Human betacoronavirus 2c EMC 2012/JX869059/1-30118 100 Bat coronavirus HKU5-1/NC 009020/1-30482 100 Bat coronavirus HKU4-1/NC 009019/1-30286 100 Bat coronavirus BtCoV/133/2005/NC 008315/1-30307 100 100 Bat coronavirus BM48-31/BGR/2008/NC 014470/1-29276 SARS coronavirus/NC 004718/1-29751 100 100 MG772933.1 Bat SARS-like coronavirus isolate bat-SL-CoVZC45/1-29802 100 MG772934.1 Bat SARS-like coronavirus isolate bat-SL-CoVZXC21/1-29732 100 BetaCoV/Nonthaburi/61/2020|EPI ISL 403962/1-29848 100 BetaCoV/Nonthaburi/74/2020|EPI ISL 403963/1-29859 Feline infectious peritonitis virus/NC 002306/1-29355 49 Human coronavirus NL63/NC 005831/1-27553 100 Human coronavirus 229E/NC 002645/1-27317 Bat coronavirus HKU2/NC 009988/1-27165 52 100 Scotophilus bat coronavirus 512/NC 009657/1-28203 81 Porcine epidemic diarrhea virus/NC 003436/1-28033 85 Bat coronavirus HKU8/NC 010438/1-28773 Bat coronavirus 1B/NC 010436/1-28476 100 100 Bat coronavirus 1A/NC 010437/1-28326 100 Turkey coronavirus/NC 010800/1-27657 100 Avian infectious bronchitis virus/NC 001451/1-27608 Beluga Whale coronavirus SW1/NC 010646/1-31686 44 Night-heron coronavirus HKU19/NC 016994/1-26077 Wigeon coronavirus HKU20/NC 016995/1-26227 Common-moorhen coronavirus HKU21/NC 016996/1-26223 100 Thrush coronavirus HKU12-600/NC 011549/1-26396 49 White-eye coronavirus HKU16/NC 016991/1-26041 100 100 Munia coronavirus HKU13-3514/NC 011550/1-26552 100 Magpie-robin coronavirus HKU18/NC 016993/1-26689 100 Sparrow coronavirus HKU17/NC 016992/1-26083 100 Porcine coronavirus HKU15/NC 016990/1-25437 CoV: coronavirus; MERS: Middle East respiratory syndrome coronavirus; SARS: severe acute respiratory syndrome. Panel A: blue: SARS-CoV-2; red: SARS; purple: MERS; green: common cold. Panel B: cyan: ACE2 human host receptor; grey: CoV spike glycoprotein trimer (PDB:6ACG); red: mutations between SARS-CoV vs SARS-CoV-2. The phylogenetic tree was created from whole genome alignment with MAFFT using the neighbour-joining method with maximum composite likelihood (MCL) model, uniform site rates and 500 bootstrap tests using MEGA X. CoV family by using a conventional nested RT-PCR [9]. thermoscanner detected a fever of 38.0 °C that was Genomic sequencing analysis included Sanger and confirmed with a tympanic thermometer. The patient next generation sequencing were performed by the reported a sore throat, that her fever onset date was 13 Emerging Infectious Diseases Health Science Center, January and that she had a cough for an undetermined the Thai Red Cross Society and the Thai National period. The patient was hospitalised at BIDI. Upon Institute of Health, Department of Medical Sciences admission, she reported mild tachypnoea, and her and the sequence shared via the Global Initiative CXR was compatible with pneumonia. Similar to case on Sharing All Influenza Data (GISAID) EpiCoV data- 1, the first CXR taken on 13 January, showed thickening base (EPI_ISL_403962). The sequencing protocol and interstitial lung marking at both lower lung fields and details are provided in the Supplementary Materials. both perihilar regions because of interstitial infiltra- The patient recovered after testing negative for SARS- tion or crowded lung marking, mild cardiomegaly and CoV-2, and returned to China without signs and symp- dilated aorta. Follow-up CXR on 17 January additionally toms on 18 January 2020. showed recent hazy with reticular opacities at left mid- dle lung field. The patient was not in severe condition Case 2 but stable. A woman in her mid-70s from Wuhan landed at Suvarnabhumi airport on 13 January. She travelled In the patient interview, she reported that she did not to Thailand with three family members as part of visit the Huanan Seafood Market or other markets. a tour group of 20 (including the case). An airport She also reported no contact with pigs, camels, other 2 www.eurosurveillance.org
Figure 2 Within-outbreak SARS-CoV-2 sequence divergence and clusters, China and Thailand, January 2020 2019-12- 2019-12- 2019-12- 2019-12- 2019-12- 2019-12- 2019-12- 2019-12- 2019-12-2019-12-2019-12-2019-12- 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 2020-01 24 30 30 30 30 30 30 30 30 30 30 30 -01 -08 -13 -14 -15 -15 -15 -16 -17 -17 -18 WUHAN- WUHAN- WUHAN-WUHAN-WUHAN- WUHAN-WUHAN-WUHAN-WUHAN-WUHAN-WUHAN-WUHAN-WUHAN-NONTH- NONTH- GUANG- GUANG- GUANG- GUANG- ZHEJI- GUANG- ZHEJI- GUANG- /IPB- /HB- /IPB- /IPB- /IPB- /IVDC-H /IVDC-H /WIV02 /WIV04 /WIV05 /WIV06 /WIV07 /IVDC-H ABURI/ ABURI/7DONG/2DONG/2 DONG/2 DONG/2 ANG/W DONG/2 ANG/W DONG/2 CAMS-W CDC-HB- CAMS- CAMS- CAMS-W B-01 B-05 B-04 61 4 0SF012 0SF013 0SF014 0SF025 Z-01 0SF028 Z-02 0SF040 H-01 01 WH-02 WH-03 H-04 2019-12-24 WUHAN/IPBCAMS-WH-01 2019-12-30 WUHAN/HBCDC-HB-01 4 2019-12-30 WUHAN/IPBCAMS-WH-02 9 7 2019-12-30 WUHAN/IPBCAMS-WH-03 4 2 7 2019-12-30 WUHAN/IPBCAMS-WH-04 3 1 6 1 2019-12-30 WUHAN/IVDC-HB-01 3 1 6 1 0 2019-12-30 WUHAN/IVDC-HB-05 5 3 8 3 2 2 2019-12-30 WUHAN/WIV02 5 3 8 3 2 2 4 2019-12-30 WUHAN/WIV04 3 1 6 1 0 0 2 2 2019-12-30 WUHAN/WIV05 5 3 8 3 2 2 4 4 2 2019-12-30 WUHAN/WIV06 3 1 6 1 0 0 2 2 0 2 2019-12-30 WUHAN/WIV07 5 3 8 3 2 2 4 4 2 4 2 2020-01-01 WUHAN/IVDC-HB-04 6 4 9 4 3 3 5 5 3 5 3 5 2020-01-08 NONTHABURI/61 3 1 6 1 0 0 2 2 0 2 0 2 3 2020-01-13 NONTHABURI/74 3 1 6 1 0 0 2 2 0 2 0 2 3 0 2020-01-14 GUANGDONG/20SF012 6 4 9 4 3 3 5 5 3 5 3 5 6 3 3 2020-01-15 GUANGDONG/20SF013 6 4 9 4 3 3 5 5 3 5 3 5 6 3 3 0 2020-01-15 GUANGDONG/20SF014 4 2 7 2 1 1 3 3 1 3 1 3 4 1 1 4 4 2020-01-15 GUANGDONG/20SF025 6 4 9 4 3 3 5 5 3 5 3 5 6 3 3 0 0 4 2020-01-16 ZHEJIANG/WZ-01 5 3 8 3 2 2 4 3 2 4 2 4 5 2 2 5 5 3 5 2020-01-17 GUANGDONG/20SF028 4 2 7 2 1 1 3 3 1 3 1 3 4 1 1 4 4 2 4 3 2020-01-17 ZHEJIANG/WZ-02 3 1 6 1 0 0 2 2 0 2 0 2 3 0 0 3 3 1 3 2 1 2020-01-18 GUANGDONG/20SF040 4 2 7 2 1 1 3 3 1 3 1 3 4 1 1 4 4 2 4 3 0 1 Number of pair-wise nt differences across whole genomes colour-coded from zero (green) to nine (red). Blue: Thai sample names. Orange: samples with sequences identical to each other and the Thai sequences. We gratefully acknowledge the authors, the originating and submitting laboratories for their sequence and metadata shared through GISAID, on which this research is based (as listed in Supplementary Table 1). mammals (or areas with dead birds), or any consump- at the protein level (Figure 1B). This surface protein is tion of raw or undercooked foods. She stated that critical for ACE2 host receptor interaction and is also she was not in contact with persons with respiratory a target of the immune response [12,13]. Given several symptoms. mutations in the binding interface, it may differ in host cell binding efficiency compared with SARS-CoV which A conventional nested RT-PCR test of this patient was could result in differences in virulence and transmis- positive for the CoV family [9]. Subsequent genome sion potential [14,15]. sequencing was again compatible with the SARS- CoV-2 and shared via the GISAID EpiCoV database The genomes of the two separate cases of coronavi- (EPI_ISL_403963). A nasopharyngeal swab also tested rus disease 2019 (COVID-19) are identical over the full positive by RT-PCR for adenovirus. The patient was no length of close to 30 kb and are furthermore identical longer febrile as of 17 January, and after testing nega- to five other sequences (four from Wuhan and one from tive for the CoV family by conventional nested RT-PCR, Zhejiang); together these sequences form the largest she was discharged and she returned to China. cluster of identical cases within the early outbreak, comprising a core of at least indirectly linked cases Genome sequence analysis (Figure 2). Within-outbreak sequence divergence is Comparing the two genome sequences with a non- generally low with 0–9 nt differences over the whole redundant selection of representatives from all known genome and mutations unique to individual strains are CoV families by alignment [10] and phylogenetic tree possibly related to quality differences of the samples (Figure 1A) [11] shows that they belong to the SARS and noise of the methods used for sequencing. family of betacoronaviruses and while related to SARS- CoV (80% genome identity), they were most closely Follow-up of contacts related to SARS-like bat CoV from China (88% identity) Case 1 travelled in a tour group with 15 others and 40 as closest known sequence at the time of emergence. close contacts were identified: 15 members of the tour Structural mapping of mutations in the spike glyco- group, 15 people sitting within two rows in front and protein between SARS CoV and the two cases of the back of the seat of case 1 on the same airplane, nine SARS-CoV-2 reported here shows only 76% identity crew members and one port health officer. They were www.eurosurveillance.org 3
Conflict of interest monitored for COVID-19 with RT-PCR tests on days 1, 12 and 14 ; all tested negative. None declared. Case 2 travelled in a tour group with 19 others and 44 close contacts were identified: 19 members of the tour Authors’ contributions group, 15 within two rows front and back of the seat of PO, SPh, TT, SPa, WW, SWai, MC, NM and NJ contributed case 2, nine crew members and one port health officer. work in the Department of Medical Sciences; RB, AV, SU and All were monitored for COVID-19 with RT-PCR tests on SI in the Department of Disease Control; and SWac in the Health Science Centre. RTCL and SMS provided the sequence days 1, 12 and 14 and all tested negative. analysis and figures. PO and SMS drafted the manuscript. All authors saw and agreed to the final version. Discussion and conclusion According to the cases’ history, the two imported COVID-19 cases described here are not directly linked, References yet their genomes are identical. Further, according to 1. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical information provided by the cases, they have no direct features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497-506. https://doi. link to the Huanan Seafood Market but their genomes org/10.1016/S0140-6736(20)30183-5 PMID: 31986264 are identical to four sequences from Wuhan collected 2. Mahase E. China coronavirus: what do we know so far? BMJ. on 30 December 2019 and sequenced by three differ- 2020;368:m308. https://doi.org/10.1136/bmj.m308 PMID: 31980434 ent laboratories, indicating potential wider distribution 3. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A Novel in the city. Although within-outbreak sequence diver- Coronavirus from Patients with Pneumonia in China, 2019. N gence is low for this virus, identical genomes of up to Engl J Med. 2020;NEJMoa2001017. https://doi.org/10.1056/ NEJMoa2001017 PMID: 31978945 30 kb are rare and a strong sign of recent transmis- 4. Du Toit A. Outbreak of a novel coronavirus. Nat Rev Microbiol. sion linkage, albeit with unknown number of genera- 2020;18(3):123. https://doi.org/10.1038/s41579-020-0332-0 PMID: 31988490 tions within the transmission chain. The missing link to 5. Wang C, Horby PW, Hayden FG, Gao GF. A novel the market, the time of travel and onset of symptoms coronavirus outbreak of global health concern. Lancet. in the two COVID-19 cases together with the incuba- 2020;395(10223):470-3. https://doi.org/10.1016/S0140- 6736(20)30185-9 PMID: 31986257 tion period of mean 6.4 days (range: 5.6–7.7) [16] and 6. Munster VJ, Koopmans M, van Doremalen N, van Riel D, de the time of closest genome neighbours obtained from Wit E. A Novel Coronavirus Emerging in China - Key Questions sequences in Wuhan, suggest that transmission within for Impact Assessment. N Engl J Med. 2020;NEJMp2000929. https://doi.org/10.1056/NEJMp2000929 PMID: 31978293 Wuhan beyond the Huanan Seafood Market is likely to 7. Gorbalenya AE, Baker SC, Baric RS, de Groot RJ, Drosten C, have occurred in the first week of January or earlier. Gulyaeva AA, et al. Severe acute respiratory syndrome-related coronavirus: The species and its viruses – a statement of the People travelling out of the city since then may have Coronavirus Study Group. bioRxiv. 2020;937862: (Preprint). spread the virus further before travel restrictions were Available from: https://doi.org/10.1101/2020.02.07.937862 enforced on 23 January [17]. 8. Naess A, Nilssen SS, Mo R, Eide GE, Sjursen H. Role of neutrophil to lymphocyte and monocyte to lymphocyte ratios in the diagnosis of bacterial infection in patients with fever. Thailand implemented measures for screening Infection. 2017;45(3):299-307. https://doi.org/10.1007/s15010- 016-0972-1 PMID: 27995553 patients travelling from Wuhan since 3 January 2020 9. Department of Medical Sciences, Ministry of Public Health, at Suvarnabhumi Airport, Don Mueang, Phuket and Thailand. Diagnostic detection of Novel coronavirus 2019 Chiang Mai airports, and stepped up surveillance at by real time RT-PCR. 23 Jan 2020. [Accessed 26 Feb 2020]. Available from: https://www.who.int/docs/default-source/ public and private hospitals. The two cases described coronaviruse/conventional-rt-pcr-followed-by-sequencing-for- here who tested positive for SARS-CoV-2 were the first detection-of-ncov-rirl-nat-inst-health-t.pdf?sfvrsn=42271c6d_4 confirmed exported cases from China, suggesting 10. Nakamura T, Yamada KD, Tomii K, Katoh K. Parallelization of MAFFT for large-scale multiple sequence alignments. early international spread. Therefore, while one can- Bioinformatics. 2018;34(14):2490-2. https://doi.org/10.1093/ not exclude the possibility that asymptomatic cases in bioinformatics/bty121 PMID: 29506019 11. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: their incubation period would be missed and become Molecular Evolutionary Genetics Analysis across Computing infectious later, screening of incoming travellers from Platforms. Mol Biol Evol. 2018;35(6):1547-9. https://doi. org/10.1093/molbev/msy096 PMID: 29722887 an affected area proved successful at least in these 12. Song W, Gui M, Wang X, Xiang Y. Cryo-EM structure of the two instances, when COVID-19 cases were detected, SARS coronavirus spike glycoprotein in complex with its isolated, observed and only discharged once they host cell receptor ACE2. PLoS Pathog. 2018;14(8):e1007236. https://doi.org/10.1371/journal.ppat.1007236 PMID: 30102747 tested negative for SARS-CoV-2. Rapid response includ- 13. Du L, He Y, Zhou Y, Liu S, Zheng B-J, Jiang S. The spike ing genome sequencing and sharing via GISAID [18,19] protein of SARS-CoV--a target for vaccine and therapeutic (https://www.gisaid.org/) enabled fast dissemination development. Nat Rev Microbiol. 2009;7(3):226-36. https:// doi.org/10.1038/nrmicro2090 PMID: 19198616 of results and provided a glimpse of early transmission 14. Zhou P, Yang X-L, Wang X-G, Hu B, Zhang L, Zhang W, et al. patterns in this outbreak. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020. https://doi.org/10.1038/ s41586-020-2012-7 PMID: 32015507 15. Wan Y, Shang J, Graham R, Baric RS, Li F. Receptor Acknowledgements recognition by novel coronavirus from Wuhan: An analysis based on decade-long structural studies of SARS. J Virol. The authors wish to thank all the participants in this study. 2020;JVI.00127-20. https://doi.org/10.1128/JVI.00127-20 We would like to thank GISAID President Peter Bogner for PMID: 31996437 help in sharing with GISAID and Dr Opart Karnkawinpong 16. Backer JA, Klinkenberg D, Wallinga J. Incubation period (Director General of Department of Medical Sciences) for his of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China, 20-28 January 2020. Euro support and guidance during the investigation. 4 www.eurosurveillance.org
Surveill. 2020;25(5). https://doi.org/10.2807/1560-7917. ES.2020.25.5.2000062 PMID: 32046819 17. Crossley G. Wuhan lockdown ‘unprecedented’, shows commitment to contain virus: WHO representative in China. Reuters. 23 Jan 2020. Available from: https://www.reuters. com/article/us-china-health-who-idUSKBN1ZM1G9 18. Shu Y, McCauley J. GISAID: Global initiative on sharing all influenza data - from vision to reality. Euro Surveill. 2017;22(13):22. https://doi.org/10.2807/1560-7917. ES.2017.22.13.30494 PMID: 28382917 19. Elbe S, Buckland-Merrett G. Data, disease and diplomacy: GISAID’s innovative contribution to global health. Glob Chall. 2017;1(1):33-46. https://doi.org/10.1002/gch2.1018 PMID: 31565258 20. Krieger E, Vriend G. YASARA View - molecular graphics for all devices - from smartphones to workstations. Bioinformatics. 2014;30(20):2981-2. https://doi.org/10.1093/bioinformatics/ btu426 PMID: 24996895 License, supplementary material and copyright This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY 4.0) Licence. You may share and adapt the material, but must give appropriate credit to the source, provide a link to the licence and indicate if changes were made. Any supplementary material referenced in the article can be found in the online version. This article is copyright of the authors or their affiliated in- stitutions, 2020. www.eurosurveillance.org 5
You can also read