SeqPanther: Sequence manipulation and mutation statistics toolset

Abstract

Pathogen genomes harbor critical information necessary to support genomic investigations that inform public health interventions such as treatment, control, and eradication. To extract this information, their sequences are analysed to identify structural variations such as single nucleotide polymorphisms (SNPs) and insertions and deletions (indels) that may be associated with phenotypes of interest. Typically, this involves generating a consensus sequence from raw reads, aligning it to a reference and identifying positions where variations occur. Several pipelines exist to map raw reads and assemble whole genomes for downstream analysis. However, there is no easy-to-use, freely available bioinformatics quality control (QC) tool to explore mappings for both positional codons and nucleotide distributions in mapped short reads of microbial genomes. To address this problem, we have developed a fast and accurate tool to summarise read counts associated with codons, nucleotides, and indels in mapped next-generation sequencing (NGS) short reads. The tool, developed in Python, also provides a visualization of the genome sequencing depth and coverage. Furthermore, the tool can be run in single or batch mode, where several genomes need to be analysed. Our tool produces a text-based report that enables quick review or can be imported into any analytical tool for upstream analysis. Additionally, the tool also provides functionality to modify the consensus sequences by adding, masking, or restoring to wild-type mutations specified by the user. Availability: SeqPanther is available at https://github.com/codemeleon/seqPanther, along with the necessary documentation for installation and usage.