Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Jun 13;38(12):3291-3293.
doi: 10.1093/bioinformatics/btac311.

ReadItAndKeep: rapid decontamination of SARS-CoV-2 sequencing reads

Martin Hunt  1   2 , Jeremy Swann  2 , Bede Constantinides  2 , Philip W Fowler  2 , Zamin Iqbal  1
Affiliations

ReadItAndKeep: rapid decontamination of SARS-CoV-2 sequencing reads

Martin Hunt et al. Bioinformatics. .

Abstract

Summary: Viral sequence data from clinical samples frequently contain contaminating human reads, which must be removed prior to sharing for legal and ethical reasons. To enable host read removal for SARS-CoV-2 sequencing data on low-specification laptops, we developed ReadItAndKeep, a fast lightweight tool for Illumina and nanopore data that only keeps reads matching the SARS-CoV-2 genome. Peak RAM usage is typically below 10 MB, and runtime less than 1 min. We show that by excluding the polyA tail from the viral reference, ReadItAndKeep prevents bleed-through of human reads, whereas mapping to the human genome lets some reads escape. We believe our test approach (including all possible reads from the human genome, human samples from each of the 26 populations in the 1000 genomes data and a diverse set of SARS-CoV-2 genomes) will also be useful for others.

Availability and implementation: ReadItAndKeep is implemented in C++, released under the MIT license, and available from https://github.com/GenomePathogenAnalysisService/read-it-and-keep.

Supplementary information: Supplementary data are available at Bioinformatics online.

PubMed Disclaimer

References

    1. Bush S.J. et al. (2020) Evaluation of methods for detecting human reads in microbial sequencing datasets. Microb. Genomics, 6, e000393. - PMC - PubMed
    1. Byrska-Bishop M. et al. (2021) High coverage whole genome sequencing of the expanded 1000 genomes project cohort including 602 trios. bioRxiv. - PMC - PubMed
    1. Jain M. et al. (2018) Nanopore sequencing and assembly of a human genome with ultra-long reads. Nat. Biotechnol., 36, 338–345. - PMC - PubMed
    1. Langmead B., Salzberg S.L. (2012) Fast gapped-read alignment with Bowtie 2. Nat. Methods, 9, 357–359. - PMC - PubMed
    1. Li H. (2013) Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997 [q-bio]. arXiv: 1303.3997.

Publication types