Number of Mismatches and Length of Longest Match Correlate with Alignment Score in Swalign Built-in Function in MATLAB

Article Information

Wenfa Ng*

Department of Biomedical Engineering, National University of Singapore, Lower Kent Ridge, Singapore

*Corresponding author: Wenfa Ng, Department of Biomedical Engineering, National University of Singapore, Lower Kent Ridge, Singapore

Received: 07 October 2021; Accepted: 21 October 2021; Published: 25 October 2021

Citation: Wenfa Ng. Number of Mismatches and Length of Longest Match Correlate with Alignment Score in Swalign Built-in Function in MATLAB. Journal of Bioinformatics and Systems Biology 4 (2021): 117-121.

Share at Facebook

Abstract

Understanding how one sequence relates to another at the nucleotide or amino acid level allows the derivation of new knowledge regarding the provenance of particular sequence as well as the determination of consensus sequence motifs that informs biological conservation at the sequence level. To this end, local or multiple sequence alignments tools in bioinformatics have been developed to automatically profile two or more nucleotide or amino acid sequence in search of matches in stretches of nucleotides or amino acid sequence that yield an alignment. While alignment score is a common metric for assessing alignment quality, relative difference between alignment scores does not readily correlate with concrete measures such as number of mismatches and length of longest match in alignment. Thus, using swalign local sequence alignment function in MATLAB on 200 alignments between RNA-seq sequence read and reference Escherichia coli K-12 MG1655 genome sequence in the sense and antisense direction, this work sought to shed some light on how alignment score from swalign correlates with number of mismatches and length of longest match. Results revealed that number of mismatches negatively correlate with alignment score; thereby, validating theoretical predictions that larger number of mismatches would result in a poorer alignment and lower alignment score. However, dependence of alignment score on other factors such as length of longest match and gap penalty from opening an alignment gap prevents linear relationship to be obtained between number of mismatches and alignment score. On the other hand, length of longest match was found to positively correlate with alignment score as predicted from theoretical understanding. But, data obtained revealed that clusters of data points gather at two regions of the scatter plot involving short matches and low alignment score, as well as long matches and high alignment score. Such clustering and sparseness of

Keywords

Swalign; Local sequence alignment tool; Number of mismatches; Length of longest match; Escherichia coli K-12 MG1655; Alignment score; Correlation; RNA Sequencing; Sequence read; Reference genome

Swalign articles; Local sequence alignment tool articles; Number of mismatches articles; Length of longest match articles; Escherichia coli K-12 MG1655 articles; Alignment score articles; Correlation articles; RNA Sequencing articles; Sequence read articles; Reference genome articles

Swalign articles Swalign Research articles Swalign review articles Swalign PubMed articles Swalign PubMed Central articles Swalign 2023 articles Swalign 2024 articles Swalign Scopus articles Swalign impact factor journals Swalign Scopus journals Swalign PubMed journals Swalign medical journals Swalign free journals Swalign best journals Swalign top journals Swalign free medical journals Swalign famous journals Swalign Google Scholar indexed journals Local sequence alignment tool articles Local sequence alignment tool Research articles Local sequence alignment tool review articles Local sequence alignment tool PubMed articles Local sequence alignment tool PubMed Central articles Local sequence alignment tool 2023 articles Local sequence alignment tool 2024 articles Local sequence alignment tool Scopus articles Local sequence alignment tool impact factor journals Local sequence alignment tool Scopus journals Local sequence alignment tool PubMed journals Local sequence alignment tool medical journals Local sequence alignment tool free journals Local sequence alignment tool best journals Local sequence alignment tool top journals Local sequence alignment tool free medical journals Local sequence alignment tool famous journals Local sequence alignment tool Google Scholar indexed journals Number of mismatches articles Number of mismatches Research articles Number of mismatches review articles Number of mismatches PubMed articles Number of mismatches PubMed Central articles Number of mismatches 2023 articles Number of mismatches 2024 articles Number of mismatches Scopus articles Number of mismatches impact factor journals Number of mismatches Scopus journals Number of mismatches PubMed journals Number of mismatches medical journals Number of mismatches free journals Number of mismatches best journals Number of mismatches top journals Number of mismatches free medical journals Number of mismatches famous journals Number of mismatches Google Scholar indexed journals Length of longest match articles Length of longest match Research articles Length of longest match review articles Length of longest match PubMed articles Length of longest match PubMed Central articles Length of longest match 2023 articles Length of longest match 2024 articles Length of longest match Scopus articles Length of longest match impact factor journals Length of longest match Scopus journals Length of longest match PubMed journals Length of longest match medical journals Length of longest match free journals Length of longest match best journals Length of longest match top journals Length of longest match free medical journals Length of longest match famous journals Length of longest match Google Scholar indexed journals Escherichia coli K-12 MG1655 articles Escherichia coli K-12 MG1655 Research articles Escherichia coli K-12 MG1655 review articles Escherichia coli K-12 MG1655 PubMed articles Escherichia coli K-12 MG1655 PubMed Central articles Escherichia coli K-12 MG1655 2023 articles Escherichia coli K-12 MG1655 2024 articles Escherichia coli K-12 MG1655 Scopus articles Escherichia coli K-12 MG1655 impact factor journals Escherichia coli K-12 MG1655 Scopus journals Escherichia coli K-12 MG1655 PubMed journals Escherichia coli K-12 MG1655 medical journals Escherichia coli K-12 MG1655 free journals Escherichia coli K-12 MG1655 best journals Escherichia coli K-12 MG1655 top journals Escherichia coli K-12 MG1655 free medical journals Escherichia coli K-12 MG1655 famous journals Escherichia coli K-12 MG1655 Google Scholar indexed journals Alignment score articles Alignment score Research articles Alignment score review articles Alignment score PubMed articles Alignment score PubMed Central articles Alignment score 2023 articles Alignment score 2024 articles Alignment score Scopus articles Alignment score impact factor journals Alignment score Scopus journals Alignment score PubMed journals Alignment score medical journals Alignment score free journals Alignment score best journals Alignment score top journals Alignment score free medical journals Alignment score famous journals Alignment score Google Scholar indexed journals Correlation articles Correlation Research articles Correlation review articles Correlation PubMed articles Correlation PubMed Central articles Correlation 2023 articles Correlation 2024 articles Correlation Scopus articles Correlation impact factor journals Correlation Scopus journals Correlation PubMed journals Correlation medical journals Correlation free journals Correlation best journals Correlation top journals Correlation free medical journals Correlation famous journals Correlation Google Scholar indexed journals RNA Sequencing articles RNA Sequencing Research articles RNA Sequencing review articles RNA Sequencing PubMed articles RNA Sequencing PubMed Central articles RNA Sequencing 2023 articles RNA Sequencing 2024 articles RNA Sequencing Scopus articles RNA Sequencing impact factor journals RNA Sequencing Scopus journals RNA Sequencing PubMed journals RNA Sequencing medical journals RNA Sequencing free journals RNA Sequencing best journals RNA Sequencing top journals RNA Sequencing free medical journals RNA Sequencing famous journals RNA Sequencing Google Scholar indexed journals Sequence read articles Sequence read Research articles Sequence read review articles Sequence read PubMed articles Sequence read PubMed Central articles Sequence read 2023 articles Sequence read 2024 articles Sequence read Scopus articles Sequence read impact factor journals Sequence read Scopus journals Sequence read PubMed journals Sequence read medical journals Sequence read free journals Sequence read best journals Sequence read top journals Sequence read free medical journals Sequence read famous journals Sequence read Google Scholar indexed journals Reference genome articles Reference genome Research articles Reference genome review articles Reference genome PubMed articles Reference genome PubMed Central articles Reference genome 2023 articles Reference genome 2024 articles Reference genome Scopus articles Reference genome impact factor journals Reference genome Scopus journals Reference genome PubMed journals Reference genome medical journals Reference genome free journals Reference genome best journals Reference genome top journals Reference genome free medical journals Reference genome famous journals Reference genome Google Scholar indexed journals

Article Details

1. Introduction

Sequence alignment is a foundational tool in bioinformatics for understanding provenance of different genes and determining consensus sequence of particular motifs such as the -35 and -10 box of promoters. The method is also used in determining how sequence reads from transcriptome sequencing matches to their corresponding genes in the reference genome. To help in this endeavour, multiple sequence alignment tools using different algorithms have been published. Typically, sequence alignment tools would present the result of the alignment in both graphical format and an alignment score. However, alignment score remains a nebulous concept without detailed knowledge of the sequence alignment algorithm. Specifically, how the alignment score relates to the number of mismatches and length of longest match are key questions that determine how an alignment score is to be interpreted. Theoretically, alignment score should negatively correlate with number of mismatches; thus, a high alignment score should correlate with a sequence alignment with few mismatches. On the other hand, alignment score typically would be higher if the length of the longest match in sequence alignment is long; thereby, yielding a positive correlation between alignment score and length of longest match. Overall, in understanding the significance of an alignment score, several parameters that feed into the calculation of alignment score would need to be understood. Besides number of mismatches and length of longest match in the alignment, penalty for opening a gap in the alignment also affects the alignment score.  

Thus, the objective of this work is to understand how number of mismatch and length of longest match affects the alignment score of swalign in-built sequence alignment tool in MATLAB. Specifically, the type of alignment investigated in this work concerns how well a sequence read from an RNA-sequencing experiment matches to the reference genome from which RNA transcripts derive. Thus, sequence alignments encountered in this work are those pertaining to short sequence reads (~50 bases) to a reference genome. Hence, limits exist in which longer matches (hundreds of base-pairs) could not be used to understand how length of longest match influence alignment score. In addition, there is also a limit on the number of mismatch possible in the alignment, which meant that the data points obtained could not fully represent the full range of alignment score possible. Such limitations represent the extend in which results derived from this work should be understood.

2. Materials and Methods

Sequence reads from RNA-sequencing data of Escherichia coli was aligned to the genome sequence of E. coli K-12 MG1655 using the swalign built-in function in MATLAB. Swalign is a local sequence alignment tool developed by Smith-Waterman [1]. Both the genome sequence and its reverse complement was used in the alignment; thereby, yielding two sequence alignment for each sequence read. The obtained alignment was used in further analysing the number of mismatches and length of longest match in the sequence alignment. Specifically, number of mismatches in the alignment was calculated and length of longest match was determined by a MATLAB programme. A total of 100 sequence reads was analysed in this work. Number of mismatches was determined by calculating the total number of mismatches in the alignment. On the other hand, length of longest match was calculated by first determining the length of each match within the alignment, and subsequently calculating the match that has the longest length of consecutive nucleotides. A zip file of the MATLAB programme could be found in the Supplementary information of this preprint.

3. Results and Discussion

Analysis of the results revealed that number of mismatches was negatively correlated with alignment score of swalign (Figure 1). Specifically, it could be observed that the higher the number of mismatches, the lower the alignment score. Given that the data points covered the range of mismatches evenly, a negative correlation could be established between alignment score and number of mismatches, which meant that alignment score of swalign function is diagnostic of how well the sequence read align to the genome sequence. The result obtained is also consistent with the theoretical understanding that a greater number of mismatches would result in a lower alignment score. However, even though a negative correlation exists between alignment score and number of mismatches, a direct quantitative relationship between the two variables could not be obtained as alignment score also depends on many other factors, one of which is length of longest match in alignment. 

For a sequence read of defined length, the greater the number of mismatches, the probability of having a long match would be correspondingly reduced, except for cases where the mismatches are clustered together at a particular section of the sequence read. This is seen in Figure 2 where clusters of data points gather at the short length of match and low alignment score region, as well as the long length of match and high alignment score region. Data available reveals that a positive correlation exists between length of longest match and alignment score, which is consistent with theoretical understanding of how well instances of long matches influence the alignment score.

fortune-biomass-feedstock

Figure 1: Negative correlation could be seen between number of mismatches in alignment and score of the alignment for sequence reads aligned to a reference genome using swalign function in MATLAB.

fortune-biomass-feedstock

Figure 2: Scatter plot of length of longest match and alignment score for alignment of 100 sequence reads to both the genome sequence of E. coli K-12 MG1655 and its reverse complement using swalign function in MATLAB.

However, the data obtained is sparse in the region of intermediate length of longest match and intermediate alignment score, which together with the clustering of data points prevents the elucidation of a quantitative positive correlation between length of longest match and alignment score. More importantly, as mentioned above, alignment score also depends on number of mismatches and the penalty imposed on opening a gap in the alignment; thus, a linear fit could not be obtained between length of longest match and alignment score.

4. Conclusions

Overall, the data obtained from aligning 100 sequence reads to the sense and antisense strand of the genome sequence of E. coli K-12 MG1655 revealed that number of mismatches in the alignment negatively correlate with alignment score readable from the swalign function which is based on the Smith and Waterman algorithm. However, dependencies of alignment score on other parameters prevent the elucidation of a linear correlation between number of mismatch and alignment score. On the other hand, length of longest match was found to cluster at two regions of the scatter plot relating length of longest match with alignment score, which also precludes the elucidation of a linear relationship between the variables through linear regression. While length of longest match positively correlates with alignment score, dependence of alignment score on other parameters such as gap penalty of opening an alignment gap and number of mismatches also hampers a linear relationship to be obtained between length of longest match and alignment score. 

Collectively, elucidation of the dependencies of alignment score on number of mismatches and length of longest match helped verify that swalign function could detect and quantify the relative effect of different number of mismatches and length of longest match on alignment quality. However, given that alignment score also depends on other variables difficult to quantify and the inability to obtain a linear relationship between alignment score and number of mismatches and length of longest match respectively, alignment score remains a nebulous quantity to the end user where its application space is limited to a qualitative understanding of alignment quality of different sequence reads to a reference genome sequence. Finally, while local sequence alignment tools are expected to evolve to better align two sequences with less computational power, utility of alignment score for obtaining a quantitative understanding of quality of alignment remain difficult as the variable is inherently dependent on a multitude of factors.    

Supplementary Information

Source code of the MATLAB programme used for this work is in the programme file that is part of the zip file 

Conflicts of Interest

The author declares no conflicts of interest.

Funding

No funding was used in this work.

References

  1. Smith T F, Waterman M S. Identification of common molecular subsequences. Journal of Molecular Biology 147 (1981): 195-197.

© 2016-2024, Copyrights Fortune Journals. All Rights Reserved