The application relates to the field of
bioinformatics, and relates to a method for improving
gene sequence alignment based on BWT transformation by using a multiplication
algorithm. First, a
suffix array is directly calculated by using the multiplication
algorithm through a reference long sequence, and then data structures required by BWT transformation are obtained. The BWT matrix does not need to be stored in space, N*N two-dimensional space is saved, and time is greatly improved. Seed sequences are subjected to pretreatment operation, the seed sequences are sorted, the same suffixes are gathered together, the intermediate result of the previous
seed sequence alignment is conveniently saved, and the number of alignment times is reduced. After the sa interval is found, a search is conducted by using a BFS (Breadth-First Search) method, the BFS can guarantee that the optimal sequence
score is found as soon as possible, the optimal sa is not missed, then a heap
data structure is used for storage, the top of the heap is the highest
score each time, the highest
score item is obtained each time to conduct the next layer search, and the
optimal matching result is obtained until the matching is finished.