SEARCH · Search PubMed

Results for “Genome alignment”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66Linked to original sources

PubMed1992

An insect picornavirus may have genome organization similar to that of caliciviruses.

Computer-assisted analysis of the amino acid sequence of the product encoded by the sequenced 3' portion of the cricket paralysis virus (CrPV), an insect picornavirus, genome showed that this protein is homologous not to the RNA-directed RNA polymerases, as originally suggested, but to the capsid proteins of mammalian picornaviruses. Alignment of the CrPV protein sequence with those of picornavirus and calicivirus capsid proteins demonstrated that the sequenced portion of the insect picornavirus genome encodes the C-terminal part of VP3 and the entire VP1. Thus CrPV seems to have a genome organization distinct from that of other picornaviruses but closely resembling that of caliciviruses, with the capsid proteins encoded in the 3' part of the genome. On the other hand, the tentative phylogenetic trees generated from the VP3 alignment revealed grouping of CrPV with hepatitis A virus, a true picornavirus, not with caliciviruses. Thus CrPV may be a picornavirus with a calicivirus-like genome organization. Different options for CrPV genome expression are discussed.

Amino Acid Sequence↗

PubMed2026

Reducing haystacks to needles - ViralClust: A Nextflow pipeline to cluster viral sequences.

BACKGROUND: The rapid accumulation of viral genome sequences presents major challenges for downstream analysis tools, including tools for multiple sequence alignments, phylogeny, and genome/alignment visualization, due to computational constraints and sampling biases caused by outbreak-driven over-representation. Selecting representative genomes through clustering offers a principled alternative to random subsampling, yet choosing appropriate clustering strategies remains non-trivial and context-dependent. RESULTS: Here, we present ViralClust, a modular Nextflow pipeline for bias-aware representative selection from large viral genome datasets. ViralClust integrates five distinct clustering algorithms (CD-HIT-EST, SUMACLUST, VSEARCH, MMSeqs2, and HDBSCAN) within a unified workflow, enabling direct comparison of clustering outcomes and flexible adaptation to diverse biological questions, considering a balanced phylogenetic distribution of the selected sequences. We evaluated ViralClust on six RNA and DNA virus datasets ranging from 632 to 156,586 sequences and spanning genome lengths from 890 to 197,185 nucleotides. Across all datasets, clustering reduced dataset size by ~95 % or more while preserving genetic diversity across species, genera, and families, and effectively mitigating biases introduced by outbreaks, partial genomes, and sequence orientation artifacts. CONCLUSIONS: By supporting whole-genome clustering and scalable representative selection, ViralClust enables efficient and reproducible downstream analyses that would otherwise be computationally infeasible. Rather than offering a prescriptive, guided analysis engine, our framework functions as a flexible comparative collection of complementary strategies, allowing users to empirically evaluate trade-offs and choose the ideal method tailored to their specific analytical endpoints.

Bioinformatics↗

PubMed2026

Chromosome-Level Reference Genome of the Desert Night Lizard Xantusia vigilis.

We present a reference-quality genome assembly for the desert night lizard (Xantusia vigilis). The night lizards (Xantusiidae) are a family of small-bodied lizards found in North America (Xantusia), Central America (Lepidophyma), and Cuba (Cricosaura). The night lizard family has an independent evolutionary history of at least 80 million years from its sister taxa within Scincoidea. The Xantusiids have several unique ecological, behavioral and evolutionary characteristics. For instance, the family contains the only squamate species that form diploid, unisexual, parthenogenic lineages. In addition, most night lizards are viviparous and form stable kin groups that are maintained over multiple years, an unusual life history strategy among lizards. Combining PacBio long-read sequencing, Hi-C, and RNAseq data we developed a reference-quality genome for the desert night lizard, X. vigilis. We assembled a complete mitochondrion and ~ 2.2 Gb nuclear genome, with 20 scaffolds that correlate in size to the X. vigilis karyotype. In addition, we found that X. vigilis chromosome 1 aligns with gene content of both of macrochromosome 1 and microchromosome 9 from a genome assembly of a species in the sister family Cordylidae (Hemicordylus capensis).

Xantusia↗

PubMed1999

Transformation distances: a family of dissimilarity measures based on movements of segments.

MOTIVATION: Evolution acts in several ways on DNA: either by mutating a base, or by inserting, deleting or copying a segment of the sequence (Ruddle, 1997; Russell, 1994; Li and Grauer, 1991). Classical alignment methods deal with point mutations (Waterman, 1995), genome-level mutations are studied using genome rearrangement distances (Bafna and Pevzner, 1993, 1995; Kececioglu and Sankoff, 1994; Kececioglu and Ravi, 1995). The latter distances generally operate, not on the sequences, but on an ordered list of genes. To our knowledge, no measure of distance attempts to compare sequences using a general set of segment-based operations. RESULTS: Here we define a new family of distances, called transformation distances, which quantify the dissimilarity between two sequences in terms of segment-based events. We focus on the case where segment-copy, -reverse-copy and -insertion are allowed in our set of operations. Those events are weighted by their description length, but other sets of weights are possible when biological information is available. The transformation distance from sequence S to sequence T is then the Minimum Description Length among all possible scripts that build T knowing S with segment-based operations. The underlying idea is related to Kolmogorov complexity theory. We present an algorithm which, given two sequences S and T, computes exactly and efficiently the transformation distance from S to T. Unlike alignment methods, the method we propose does not necessarily respect the order of the residues within the compared sequences and is therefore able to account for duplications and translocations that cannot be properly described by sequence alignment. A biological application on Tnt1 tobacco retrotransposon is presented. AVAILABILITY: The algorithm and the graphical interface can be downloaded at http://www.lifl.fr/ approximately varre/TD

Algorithms↗

PubMed2004

DbW: automatic update of a functional family-specific multiple alignment.

MOTIVATION: Recent advances in gene sequencing have provided complete sequence information for a number of genomes and as a result the amount of data in the sequence databases is growing at an exponential rate. We introduce here a new program, DbW, to automate the update of a functional family-specific multiple alignment that tries to include relevant sequences. The program is based on the use of different sources of information: sequences and annotations in databases. RESULTS: The advantages of DbW are demonstrated using the 20 families of aminoacyl-tRNA synthetases, where DbW detects a maximum of homologous sequences in the Swiss-Prot and SPTREMBL databases. The global specificity of DbW in this test is 98.4% (1.6% of the sequences included in the alignment did not belong to the family according to their function), and the global sensitivity of DbW is estimated to be 95.2%. Thus, DbW provides a reliable basis for the many applications that rely on accurate multiple alignments, e.g. functional residue identification, 2D/3D structure prediction or homology modeling. AVAILABILITY: The DbW software is available for download at ftp://ftp-igbmc.u-strasbg.fr/pub/DbW/DbW.tar and online at http://titus.u-strasbg.fr/DbW CONTACT: prigent@igbmc.u-strasbg.fr.

Algorithms↗

PubMed2026

Secondary metabolite profiling of rare Micromonospora spp. from cold desert of NW Himalayas via multi-omics analysis.

INTRODUCTION: The genus Micromonospora is a prolific producer of specialized metabolites with pharmacological and agronomic relevance. Natural products derived from the genus Micromonospora have a distinctive chemical diversity and enormous therapeutic potential, thus represent a potential source for drugs and drug leads. OBJECTIVE: To explore the biosynthetic potential of four Micromonospora strains isolated from cold desert of NW Himalayas through genome mining and to correlate predicted biosynthetic gene clusters with chemical features detected by untargeted LC-HRMS metabolomics. METHOD: High-quality genomes were annotated for BGCs and matched against untargeted LC-HRMS features (peak picking, alignment, and annotation to chemical classes). Each isolate was grown in triplicate, and fermented broth was pooled for further metabolomic studies. RESULTS: By integrating genomic and metabolomic approaches, specialized biosynthetic gene clusters and strain-based putative metabolite classes were identified. LRS1 showed elevated xanthines (RiPP/siderophore), LRS3 had phenolic glycosides (hybrid PKS/NRPS), LRS4 showed 70-fold hydroxycinnamate enrichment (Type II PKS), and LRS5 displayed p-benzoquinone enrichment (Type III PKS). The metabolite profile of each strain aligned with its predicted biosynthetic gene cluster composition. CONCLUSION: Under a single growth regime, each Micromonospora strain exhibits a distinct metabolomic profile. This metabologenomics workflow can be further explored to isolate specialized metabolites with potential therapeutic and agricultural value.

Micromonospora↗

PubMed2004

Development and large scale benchmark testing of the PROSPECTOR_3 threading algorithm.

This article describes the PROSPECTOR_3 threading algorithm, which combines various scoring functions designed to match structurally related target/template pairs. Each variant described was found to have a Z-score above which most identified templates have good structural (threading) alignments, Z(struct) (Z(good)). 'Easy' targets with accurate threading alignments are identified as single templates with Z > Z(good) or two templates, each with Z > Z(struct), having a good consensus structure in mutually aligned regions. 'Medium' targets have a pair of templates lacking a consensus structure, or a single template for which Z(struct) < Z < Z(good). PROSPECTOR_3 was applied to a comprehensive Protein Data Bank (PDB) benchmark composed of 1491 single domain proteins, 41-200 residues long and no more than 30% identical to any threading template. Of the proteins, 878 were found to be easy targets, with 761 having a root mean square deviation (RMSD) from native of less than 6.5 A. The average contact prediction accuracy was 46%, and on average 17.6 residue continuous fragments were predicted with RMSD values of 2.0 A. There were 606 medium targets identified, 87% (31%) of which had good structural (threading) alignments. On average, 9.1 residue, continuous fragments with RMSD of 2.5 A were predicted. Combining easy and medium sets, 63% (91%) of the targets had good threading (structural) alignments compared to native; the average target/template sequence identity was 22%. Only nine targets lacked matched templates. Moreover, PROSPECTOR_3 consistently outperforms PSIBLAST. Similar results were predicted for open reading frames (ORFS) < or =200 residues in the M. genitalium, E. coli and S. cerevisiae genomes. Thus, progress has been made in identification of weakly homologous/analogous proteins, with very high alignment coverage, both in a comprehensive PDB benchmark as well as in genomes.

Algorithms↗

PubMed2002

Molecular phylogenetics and mitochondrial genomic evolution in the chamaeleonidae (Reptilia, Squamata).

A phylogenetic hypothesis for the lizard family Chamaeleonidae is generated from 1503 aligned base positions (883 parsimony-informative) of mitochondrial DNA for specimens representing 59 species (57 ingroup and two outgroup). Sequences are reported for a genomic segment encoding eight transfer RNAs, NADH dehydrogenase component 2 (ND2), and portions of NADH dehydrogenase component 1 (ND1) and cytochrome c oxidase subunit 1 (COI). Newly reported genomic rearrangements and duplications support the hypothesis that mitochondrial gene order and content are destabilized by phylogenetic loss of a functional origin for light-strand replication between the genes encoding tRNA(Asn) and tRNA(Cys). A novel gene order characterizes all sampled Brookesia except B. nasus. Brookesia nasus, the apparent sister taxon of a clade formed by all other Brookesia, has the ancestral gene order but contains a large tandem duplication. An apparently noncoding 220 base pair insertion between the genes encoding ND2 and tRNA(Trp) is reported for Bradypodion tavetanum. Phylogenetic analysis identifies nine clades whose ancestral lineages diverged early in chamaeleonid evolutionary history: (1) Brookesia (possibly excluding B. nasus), (2) Chamaeleo subgenus Chamaeleo (excluding C. namaquensis), (3) Chamaeleo subgenus Trioceros, (4) viviparous Bradypodion, (5) oviparous Bradypodion, (6) genus Furcifer (except F. balteatus), and (7-9) three distinct clades of Calumma. Chamaeleo namaquensis, Brookesia nasus, Furcifer balteatus, Rhampholeon brevicaudatus, and R. spectrum represent ancient lineages dating to approximately the same time. Multiple independent losses and a possible secondary gain of horns are inferred for Trioceros. Viviparity has at least two separate origins in chameleons, one in Bradypodion and

Animals↗

PubMed2003

Comparative map between chicken chromosome 15 and human chromosomal region 12q24 and 22q11-q12.

The physical and comparative map of GGA15 was improved by the construction of 9 BAC contigs around loci previously mapped on GGA15 by linkage analysis. In total, 240 BAC clones were isolated, covering 30-35% of GGA15, and 120 STS were developed (104 STS derived from BAC end sequences and 18 STS derived within genes). Seventeen chicken orthologues of human genes located on human Chr 22q11-q12 were directly mapped within BAC contigs of GGA15. Furthermore, the partial sequences of the chicken BAC clones were compared with sequences present in the EMBL/GenBank databases and revealed matches to 26 genes, ESTs, and genomic clones located on HSA22q11-q12 and HSA12q24. These results provide a better alignment of GGA15 with the corresponding regions in human and mouse, and improve our knowledge of the evolution and dynamics of the vertebrate genome.

Animals↗

PubMed2003

Putative proteins related to group II intron reverse transcriptase/maturases are encoded by nuclear genes in higher plants.

The Arabidopsis thaliana nuclear genome sequence revealed several open reading frames encoding proteins related to group II intron-encoded reverse transcriptase/maturases. Here, we show via sequence alignments that at least four such open reading frames are conserved in the nuclear genomes of A.thaliana and Oryza sativa (rice) and that they encode putative proteins belonging to two different classes (nMat-1 and nMat-2), neither of which is associated with a group II intron RNA structure. The two nMat-1 proteins have reverse transcriptase, maturase and DNA endonuclease domains characteristic of canonical group II intron-encoded proteins, while the two nMat-2 proteins have reverse transcriptase and maturase domains linked to a novel C-terminal domain. Although some nMat proteins have mutations expected to inactivate intron mobility functions, all could potentially retain the RNA splicing function. These nuclear maturase-like proteins may be imported into organelles to function in group II intron splicing and/or they may have assumed other cellular functions. Nuclear-encoded maturases could regulate organellar gene expression and may reflect a step in the evolution of mobile group II introns into spliceosomal introns.

Amino Acid Sequence↗

PubMed2025

GiantHunter: accurate detection of giant virus in metagenomic data using reinforcement-learning and Monte Carlo tree search.

MOTIVATION: Nucleocytoplasmic large DNA viruses (NCLDVs) are notable for their large genomes and extensive gene repertoires, which contribute to their widespread environmental presence and critical roles in processes such as host metabolic reprogramming and nutrient cycling. Metagenomic sequencing has emerged as a powerful tool for uncovering novel NCLDVs in environmental samples. However, identifying NCLDV sequences in metagenomic data remains challenging due to their high genomic diversity, limited reference genomes, and shared regions with other microbes. Existing alignment-based and machine learning methods struggle with achieving optimal trade-offs between sensitivity and precision. RESULTS: In this work, we present GiantHunter, a reinforcement learning-based tool for identifying NCLDVs from metagenomic data. By employing a Monte Carlo tree search strategy, GiantHunter dynamically selects representative non-NCLDV sequences as the negative training data, enabling the model to establish a robust decision boundary. Benchmarking on rigorously designed experiments shows that GiantHunter achieves high precision while maintaining competitive sensitivity, improving the F1-score by 10% and reducing computational cost by 90% compared to the second-best method. To demonstrate its real-world utility, we applied GiantHunter to 60 metagenomic datasets collected from six cities along the Yangtze River, located both upstream and downstream of the Three Gorges Dam. The results reveal significant differences in NCLDV diversity correlated with proximity to the dam, likely influenced by reduced flow velocity caused by the dam. These findings highlight GiantHunter's potential to advance our understanding of NCLDVs and their ecological roles in diverse environments. AVAILABILITY AND IMPLEMENTATION: The source code of GiantHunter is available via: https://github.com/FuchuanQu/GiantHunter.

Metagenomics↗

PubMed2004

ATRX, a member of the SNF2 family of helicase/ATPases, is required for chromosome alignment and meiotic spindle organization in metaphase II stage mouse oocytes.

ATRX is a centromeric heterochromatin binding protein belonging to the SNF2 family of helicase/ATPases with chromatin remodeling activity. Mutations in the human ATRX gene result in X-linked alpha-thalassaemia with mental retardation (ATRX) syndrome and correlate with changes in methylation of repetitive DNA sequences. We show here that ATRX also functions to regulate key stages of meiosis in mouse oocytes. At the germinal vesicle (GV) stage, ATRX was found associated with the perinucleolar heterochromatin rim in transcriptionally quiescent oocytes. Phosphorylation of ATRX during meiotic maturation is dependent upon calcium calmodulin kinase (CamKII) activity. Meiotic resumption also coincides with deacetylation of histone H4 at lysine 5 (H4K5 Ac) while ATRX and histone H3 methylated on lysine 9 (H3K9) remained bound to the centromeres and interstitial regions of condensing chromosomes, respectively. Inhibition of histone deacetylases (HDACs) with trichostatin A (TSA) disrupted ATRX binding to the centromeres of hyperacetylated chromosomes resulting in abnormal chromosome alignments at metaphase II (MII). Similarly, while selective ablation of ATRX by antibody microinjection and RNA interference (RNAi) had no effect on the progression of meiosis, it had severe consequences for the alignment of chromosomes on the metaphase II spindle. These results suggest that genome-wide epigenetic modifications such as global histone deacetylation are essential for the binding of ATRX to centromeric heterochromatin. Moreover, centromeric ATRX is required for correct chromosome alignment and organization of a bipolar meiotic metaphase II spindle.

Animals↗

PubMed2002

Single nucleotide polymorphism mapping using genome-wide unique sequences.

As more and more genomic DNAs are sequenced to characterize human genetic variations, the demand for a very fast and accurate method to genomically position these DNA sequences is high. We have developed a new mapping method that does not require sequence alignment. In this method, we first identified DNA fragments of 15 bp in length that are unique in the human genome and then used them to position single nucleotide polymorphism (SNP) sequences. By use of four desktop personal computers with AMD K7 (1 GHz) processors, our new method mapped more than 1.6 million SNP sequences in 20 hr and achieved a very good agreement with mapping results from alignment-based methods.

Chromosome Mapping↗

PubMed2002

The tmRNA Website: invasion by an intron.

tmRNA (also known as 10Sa RNA or SsrA) plays a central role in an unusual mode of translation, whereby a stalled ribosome switches from a problematic mRNA to a short reading frame within tmRNA during translation of a single polypeptide chain. Research on the mechanism, structure and biology of tmRNA is served by the tmRNA Website, a collection of sequences for tmRNA and the encoded proteolysis-inducing peptide tags, alignments, careful documentation and other information; the URL is http://www.indiana.edu/~tmrna. Four pseudoknots are usually present in each tmRNA, so the database is rich with information on pseudoknot variability. Since last year it has doubled (227 tmRNA sequences as of September 2001), a sequence alignment for the tmRNA cofactor SmpB has been included, and genomic data for Clostridium botulinum has revealed a group I (subgroup IA3) intron interrupting the tmRNA T-loop.

Base Sequence↗

PubMed2025

Complete genome sequence and genomic characterization of the probiotic Limosilactobacillus reuteri PSC102.

BACKGROUND: Gut microbiota are potential sources of probiotics and play an essential role in maintaining intestinal health. Limosilactobacillus reuteri PSC102 (L. reuteri PSC102), which was isolated from the feces of healthy pigs, exhibited health-beneficial properties. AIM: We aimed to conduct a whole-genome sequencing analysis of L. reuteri PSC102 to determine its molecular characteristics as a probiotic strain. METHODS: Limosilactobacillus reuteri PSC102 cells were cultured in De Man-Rogosa-Sharpe medium, followed by DNA extraction for genomic analysis using the PacBio-Illumina sequencing platform. The EzBioCloud software was used to perform gene assembly, and the genes were interpreted by the National Center for Biotechnology Information (NCBI) and the Glimmer program. Core and pan-genomic analyses were performed to assess the extent of functional conservation in the genomic sequence. Moreover, the NCBI database and the Basic Local Alignment Search Tool software were used to identify antimicrobial resistance genes and virulence factors. RESULTS: Limosilactobacillus reuteri PSC102 consists of a single circular chromosome with 2,048,626 bp, a guanine- cytosine of 38.9%, 18 rRNA genes, and 69 tRNA genes. Among the 1,846 protein-coding sequences, genes associated with probiotic characteristics were identified, including genes involved in host-microbe interactions, stress tolerance, biogenesis, and defense mechanisms. Furthermore, the genome of L. reuteri PSC102 comprises 2,446 pan-genome and 1,222 core-genome orthologous gene clusters. A total of 74 unique genes were identified in L. reuteri PSC102 genome. These genes mostly encode proteins potentially involved in the transport and metabolism of amino acids and carbohydrates. Moreover, antibacterial resistance genes and virulence factors were absent in L. reuteri PSC102. CONCLUSION: The results of the molecular insight into L. reuteri PSC102 corroborates its use as a probiotic in humans and other animals.

Limosilactobacillus reuteri↗

PubMed2000

Comparison between Pyrococcus horikoshii and Pyrococcus abyssi genome sequences reveals linkage of restriction-modification genes with large genome polymorphisms.

Recent work suggests that restriction-modification gene complexes are mobile genetic elements that insert themselves into the genome and cause various genome rearrangements. In the present work, the complete genome sequences of Pyrococcus horikoshii and Pyrococcus abyssi, two species in a genus of hyperthermophilic archaeon (archaebacterium), were compared to detect large genome polymorphisms linked with restriction-modification gene homologs. Sequence alignments, GC content analysis, and codon usage analysis demonstrated the diversity of these homologs and revealed a possible case of relatively recent acquisition (horizontal transfer). In two cases out of the six large polymorphisms identified, there was insertion of a DNA segment with a modification gene homolog, accompanied by target deletion (simple substitution). In two other cases, homologous DNA segments carrying a modification gene homolog were present at different locations in the two genomes (transposition). In both cases, substitution (insertion/deletion) in one of the two loci was accompanied by inversion of adjacent chromosomal segment. In the fifth case, substitution by a DNA segment carrying type I restriction, modification, and specificity gene homologs was likewise accompanied by adjacent inversion. In the last case, two homologous DNA segments, were found at different loci in the two genomes (transposition), but only one of them had insertion of a modification homolog and an unknown ORF. The possible relationship of these polymorphisms to attack by restriction enzymes on the chromosome will be discussed.

Chromosome Inversion↗

PubMed2026

Whole genome sequencing reveals a specific microbiota in subglottic stenosis C. acnes may contribute to inflammation.

PURPOSE: Subglottic stenosis (SGS) progressively reduces the airway below the vocal folds. The cause is not known and there is a recurrent need of surgical treatment. Including all phenotypes, SGS affects 1/400 000/yr, with a female dominance. Previous studies have revealed a possible role of the Mycobacterium complex in SGS development. Our hypothesis is that microbiota is associated with the inflammation in SGS, if true it might affect the prevailing treatment options. METHODS: This prospective cross-sectional study included biopsies from 34 patients with subglottic stenosis, collected between 2020 and 2023. Nucleic acids were extracted from the tissue samples and analysed using whole genome sequencing. Microbial composition was characterized using taxonomic profiling of sequencing data. Species with sufficient read counts were selected for further validation using sequence alignment methods to ensure accuracy of identification. RESULTS: Using the most comprehensive form of genomic testing currently in clinical use, we present curated and stable data on the presence of Cutibacterium acnes in 28 out of the 34 cases. CONCLUSION: Cutibacterium acnes may serve as a driver of the inflammation characterizing SGS and should be considered in therapeutically oriented future studies.

Cutibacterium acnes↗

PubMed2004

The BAX PCR assay for screening Listeria monocytogenes targets a partial putative gene lmo2234.

The BAX PCR for screening Listeria monocytogenes is a commercial PCR assay for specifically targeting L. monocytogenes, a foodborne pathogen that can contaminate a variety of foods and cause a potentially fatal disease, listeriosis, among high-risk populations. The high specificity (> 98%) of this PCR assay is achieved by targeting a species-specific genomic region (approximately 400 bp) presumably found only in L. monocytogenes. In this study, the identity of the BAX PCR-targeted genomic region was determined by using PCR cloning, DNA sequencing, and basic local alignment search tool (BLAST) analysis of the amplicon sequences of an L. monocytogenes serotype 1/2a strain. BLAST analysis identified the BAX PCR amplicon (GenBank accession no. AY364605) as a 423-bp genomic region between nucleotides 224,409 and 224,831 in the genome of L. monocytogenes (serotype 1/2a strain EGD-e), including a 145-bp noncoding region and a 278-bp partial coding sequence of a putative gene, lmo2234. The translated amino acid sequence (92 amino acids) of this partial coding region is highly conserved between L. monocytogenes and Listeria innocua (93% homology). Reverse-position-specific BLAST analysis identified a conserved domain in Lmo2234 that was similar (95.3% aligned, E value = 9E-18) to the consensus amino acid sequence of sugar phosphate isomerases/epimerases (National Center for Biotechnology Information conserved domain database accession no. COG 1082.1, IolE), indicating that Lmo2234 might be involved in bacterial carbohydrate transport and metabolism.

Amino Acid Sequence↗