Abstract:Many algorithms for sequence analysis rely on word matching or word statistics. Often, these approaches can be improved if binary patterns representing match and don't-care positions are used as a filter, such that only those positions of words are considered that correspond to the match positions of the patterns. The performance of these approaches, however, depends on the underlying patterns. Herein, we show that the overlap complexity of a pattern set that was introduced by Ilie and Ilie is closely related to the variance of the number of matches between two evolutionarily related sequences with respect to this pattern set. We propose a modified hill-climbing algorithm to optimize pattern sets for database searching, read mapping and alignment-free sequence comparison of nucleic-acid sequences; our implementation of this algorithm is called rasbhari. Depending on the application at hand, rasbhari can either minimize the overlap complexity of pattern sets, maximize their sensitivity in database searching or minimize the variance of the number of pattern-based matches in alignment-free sequence comparison. We show that, for database searching, rasbhari generates pattern sets with slightly higher sensitivity than existing approaches. In our Spaced Words approach to alignment-free sequence comparison, pattern sets calculated with rasbhari led to more accurate estimates of phylogenetic distances than the randomly generated pattern sets that we previously used. Finally, we used rasbhari to generate patterns for short read classification with CLARK-S. Here too, the sensitivity of the results could be improved, compared to the default patterns of the program. We integrated rasbhari into Spaced Words; the source code of rasbhari is freely available at <a class="link-external link-http" href="http://rasbhari.gobics.de/" rel="external noopener nofollow">this http URL</a>

Pebblescout is an easy-to-use tool for fast sequence search in petabase-scale nucleotide resources

Indexing and searching petabase-scale nucleotide resources

Indexing All Life's Known Biological Sequences

Reference Sequence Browser: An R application with a User-Friendly GUI to rapidly query sequence databases

BnaSNPDB: an Interactive Web Portal for the Efficient Retrieval and Analysis of SNPs among 1,007 Rapeseed Accessions

Sequenceserver: A Modern Graphical User Interface for Custom BLAST Databases

Lapidary: Identifying and reporting amino acid sequences in metagenomes using sequence reads and Diamond

NucleoSeeker - Precision filtering of RNA databases to curate high-quality datasets

Indexing and real-time user-friendly queries in terabyte-sized complex genomic datasets with kmindex and ORA

Gencube: Efficient retrieval, download, and unification of genomic data from leading biodiversity databases

Proksee: in-depth characterization and visualization of bacterial genomes

Ingap: an Integrated Next-Generation Genome Analysis Pipeline

piBUSS: a parallel BEAST/BEAGLE utility for sequence simulation under complex evolutionary scenarios

iSeq: An integrated tool to fetch public sequencing data

Seqminer2: an efficient tool to query and retrieve genotypes for statistical genetics analyses from biobank scale sequence dataset

RasBhari: optimizing spaced seeds for database searching, read mapping and alignment-free sequence comparison

mbctools: A User-Friendly Metabarcoding and Cross-Platform Pipeline for Analyzing Multiple Amplicon Sequencing Data across a Large Diversity of Organisms

Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor

Mantis: A Fast, Small, and Exact Large-Scale Sequence-Search Index

BusyBee Web: metagenomic data analysis by bootstrapped supervised binning and annotation

SeqCAT: Sequence Conversion and Analysis Toolbox