定義：タンデムリピート 記号の集合を $\Sigma$（例えば $\{A,C,G,T\}$）、長さ1以上の記号列の集合を $\Sigma^{+}$ と表記する。記号列 $p\in\Sigma^{+}$ が $k(\geq1)$ 回連続して出現する記号列を、正確なタンデムリピート(tandem repeat)と呼び、$p^k$ と表記する。ある記号列を正確なタンデムリピートとアラインメントした時に、ある閾値以上の配列類似度を持つ時、タンデムリピートと呼ぶ。

Ref: Kerasy Documentation: Tandem Repeats

演習問題1(必須)

　記号列 $q$ が正確なタンデムリピート $p^k$ を部分列として含むとする。$p$ と $k$ を検出するアルゴリズムを考えよ。

　例題として $10$ 個の配列を与えるので、各配列から $p$ と $k$ を決定せよ。

解答

In [1]:

import time
import matplotlib.pyplot as plt

In [2]:

from kerasy.Bio.tandem import find_tandem
from kerasy.Bio.string import StringSearch
from kerasy.utils import readMultiSeq

In [3]:

sequences = readMultiSeq("exact.fasta")

In [4]:

len_sequences = [len(seq) for seq in sequences]

Results¶

In [8]:

fig, (ax_sais, ax_dp) = plt.subplots(1,2,figsize=(12,6))

ax_sais.scatter(len_sequences, pt_sais, color="r")
ax_sais.set_xlabel("the length of sequence"), ax_sais.set_ylabel("Processing time."), ax_sais.set_title("Suffix Array", fontsize=18)

ax_dp.scatter(len_sequences, pt_dp, color="b")
ax_dp.set_xlabel("the length of sequence"), ax_dp.set_ylabel("Processing time."), ax_dp.set_title("Dynamic Programming", fontsize=18)

plt.tight_layout()
plt.show()

Suffix Array	Dynamic Programming
$$O(n)$$	$$O(n^2)$$

find tandem repeat DP¶

In [5]:

pt_dp = []
for seq in sequences:
    s = time.time()
    score, tandems = find_tandem(seq, method="DP")
    pt_dp.append(time.time()-s)
    print(f"score : {score:<4}")
    print(f"tandem:\n{tandems[0]}")
    print()

score : 33  
tandem:
TTG

score : 259 
tandem:
TTTTGCC

score : 493 
tandem:
ATCAGTTTAAATTACCA

score : 3127
tandem:
TCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGC

score : 1001
tandem:
ACTGGTATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTG

score : 2369
tandem:
CCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCG

score : 2869
tandem:
CGTTCTGATTTTACTTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCACGGAAGATTATGAGCACTTCGCGG

score : 2223
tandem:
TAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGC

score : 8517
tandem:
CATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGAC

score : 7007
tandem:
ATTAGTACCCCATCGGTTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTGCTACAGTCTTACTCCCCCTATAGCTGGGTGACT

find tandem repeat using Suffix Array¶

In [6]:

pt_sais = []
for seq in sequences:
    s = time.time()
    score, tandems = find_tandem(seq, method="SAIS")
    pt_sais.append(time.time()-s)
    print(f"score : {score:<4}")
    print(f"tandem:\n{tandems[0]}")
    print()

score : 33  
tandem:
TGT

score : 259 
tandem:
TTTGCCT

score : 493 
tandem:
TCAGTTTAAATTACCAA

score : 3127
tandem:
TCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGC

score : 1001
tandem:
GTATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTGACTG

score : 2369
tandem:
CCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCG

score : 2869
tandem:
GTTCTGATTTTACTTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCACGGAAGATTATGAGCACTTCGCGGC

score : 2223
tandem:
TAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGC

score : 8517
tandem:
CATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGAC

score : 7007
tandem:
ATTAGTACCCCATCGGTTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTGCTACAGTCTTACTCCCCCTATAGCTGGGTGACT

display tandem repeat¶

In [7]:

for seq in sequences:
    db = StringSearch(seq, verbose=0)
    tandem = db.find_tandem()
    db.where(tandem, width=50)
    print("\n")

Model: Suffix Array
Number of matches: 11
Query: TGT
===============================================================
S: [00] GCCCCTGCACTCTCGAGTTGCGCAAAGTCTTGTTGTTGTTGTTGTTGTTG [50]
 : [00]                               <-><-><-><-><-><-><- [50]

S: [50] TTGTTGTTGTTGTGGGAATTGGC [72]
 : [50] ><-><-><-><->           [72]
===============================================================


Model: Suffix Array
Number of matches: 37
Query: TTTGCCT
=================================================================
S: [000] ATTACTATTCGCCGATGGTTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTT [050]
 : [000]                     <-----><-----><-----><-----><- [050]

S: [050] TGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTT [100]
 : [050] ----><-----><-----><-----><-----><-----><-----><-- [100]

S: [100] GCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTG [150]
 : [100] ---><-----><-----><-----><-----><-----><-----><--- [150]

S: [150] CCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGC [200]
 : [150] --><-----><-----><-----><-----><-----><-----><---- [200]

S: [200] CTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCCTTTTGCC [250]
 : [200] -><-----><-----><-----><-----><-----><-----><----- [250]

S: [250] TTTTGCCTTTTGCCTTTTGCCTTTTGCCTCAGTGTTTGCCCATGTAATG [298]
 : [250] ><-----><-----><-----><----->                     [298]
=================================================================


Model: Suffix Array
Number of matches: 29
Query: TCAGTTTAAATTACCAA
=================================================================
S: [000] GATAGACGGATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTT [050]
 : [000]           <---------------><---------------><----- [050]

S: [050] TAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGT [100]
 : [050] ----------><---------------><---------------><---- [100]

S: [100] TTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAG [150]
 : [100] -----------><---------------><---------------><--- [150]

S: [150] TTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCA [200]
 : [150] ------------><---------------><---------------><-- [200]

S: [200] GTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAATC [250]
 : [200] -------------><---------------><---------------><- [250]

S: [250] AGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAAT [300]
 : [250] --------------><---------------><--------------->< [300]

S: [300] CAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCAA [350]
 : [300] ---------------><---------------><---------------> [350]

S: [350] TCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACCA [400]
 : [350] <---------------><---------------><--------------- [400]

S: [400] ATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTACC [450]
 : [400] ><---------------><---------------><-------------- [450]

S: [450] AATCAGTTTAAATTACCAATCAGTTTAAATTACCAATCAGTTTAAATTAC [500]
 : [450] -><---------------><---------------><------------- [500]

S: [500] CAACTATAGGAGCAAGGTTGGTTTCTAACTCCT [532]
 : [500] -->                               [532]
=================================================================


Model: Suffix Array
Number of matches: 53
Query: TCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGC
===================================================================
S: [0000] ACCACTACGGTTTGATGAAAGTAGTCAAGGAGGTCGACGTGTAGGGCGTG [0050]
 : [0000]                                                    [0050]

S: [0050] TCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCT [0100]
 : [0050] <------------------------------------------------- [0100]

S: [0100] CCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTG [0150]
 : [0100] --------><---------------------------------------- [0150]

S: [0150] ACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCC [0200]
 : [0150] -----------------><------------------------------- [0200]

S: [0200] CGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAG [0250]
 : [0200] --------------------------><---------------------- [0250]

S: [0250] GCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTA [0300]
 : [0250] -----------------------------------><------------- [0300]

S: [0300] GGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCAC [0350]
 : [0300] --------------------------------------------><---- [0350]

S: [0350] CCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATC [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] ATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTC [0450]
 : [0400] ---><--------------------------------------------- [0450]

S: [0450] AGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGG [0500]
 : [0450] ------------><------------------------------------ [0500]

S: [0500] CCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGC [0550]
 : [0500] ---------------------><--------------------------- [0550]

S: [0550] AGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAAC [0600]
 : [0550] ------------------------------><------------------ [0600]

S: [0600] GAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGA [0650]
 : [0600] ---------------------------------------><--------- [0650]

S: [0650] CTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCT [0700]
 : [0650] ------------------------------------------------>< [0700]

S: [0700] CCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTC [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] CATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGA [0800]
 : [0750] -------><----------------------------------------- [0800]

S: [0800] CGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCC [0850]
 : [0800] ----------------><-------------------------------- [0850]

S: [0850] GTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGG [0900]
 : [0850] -------------------------><----------------------- [0900]

S: [0900] CAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAG [0950]
 : [0900] ----------------------------------><-------------- [0950]

S: [0950] GAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACC [1000]
 : [0950] -------------------------------------------><----- [1000]

S: [1000] CCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCA [1050]
 : [1000] -------------------------------------------------- [1050]

S: [1050] TGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCA [1100]
 : [1050] --><---------------------------------------------- [1100]

S: [1100] GCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGC [1150]
 : [1100] -----------><------------------------------------- [1150]

S: [1150] CTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCA [1200]
 : [1150] --------------------><---------------------------- [1200]

S: [1200] GCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACG [1250]
 : [1200] -----------------------------><------------------- [1250]

S: [1250] AAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGAC [1300]
 : [1250] --------------------------------------><---------- [1300]

S: [1300] TTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTC [1350]
 : [1300] -----------------------------------------------><- [1350]

S: [1350] CACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCC [1400]
 : [1350] -------------------------------------------------- [1400]

S: [1400] ATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGAC [1450]
 : [1400] ------><------------------------------------------ [1450]

S: [1450] GTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCG [1500]
 : [1450] ---------------><--------------------------------- [1500]

S: [1500] TGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGC [1550]
 : [1500] ------------------------><------------------------ [1550]

S: [1550] AGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGG [1600]
 : [1550] ---------------------------------><--------------- [1600]

S: [1600] AACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCC [1650]
 : [1600] ------------------------------------------><------ [1650]

S: [1650] CGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCAT [1700]
 : [1650] -------------------------------------------------- [1700]

S: [1700] GCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAG [1750]
 : [1700] -><----------------------------------------------- [1750]

S: [1750] CTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCC [1800]
 : [1750] ----------><-------------------------------------- [1800]

S: [1800] TGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAG [1850]
 : [1800] -------------------><----------------------------- [1850]

S: [1850] CCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGA [1900]
 : [1850] ----------------------------><-------------------- [1900]

S: [1900] AGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACT [1950]
 : [1900] -------------------------------------><----------- [1950]

S: [1950] TAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCC [2000]
 : [1950] ----------------------------------------------><-- [2000]

S: [2000] ACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCA [2050]
 : [2000] -------------------------------------------------- [2050]

S: [2050] TCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACG [2100]
 : [2050] -----><------------------------------------------- [2100]

S: [2100] TCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGT [2150]
 : [2100] --------------><---------------------------------- [2150]

S: [2150] GGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCA [2200]
 : [2150] -----------------------><------------------------- [2200]

S: [2200] GCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGA [2250]
 : [2200] --------------------------------><---------------- [2250]

S: [2250] ACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCC [2300]
 : [2250] -----------------------------------------><------- [2300]

S: [2300] GACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATG [2350]
 : [2300] -------------------------------------------------- [2350]

S: [2350] CTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGC [2400]
 : [2350] ><------------------------------------------------ [2400]

S: [2400] TCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCT [2450]
 : [2400] ---------><--------------------------------------- [2450]

S: [2450] GACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGC [2500]
 : [2450] ------------------><------------------------------ [2500]

S: [2500] CCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAA [2550]
 : [2500] ---------------------------><--------------------- [2550]

S: [2550] GGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTT [2600]
 : [2550] ------------------------------------><------------ [2600]

S: [2600] AGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCA [2650]
 : [2600] ---------------------------------------------><--- [2650]

S: [2650] CCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCAT [2700]
 : [2650] -------------------------------------------------- [2700]

S: [2700] CATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGT [2750]
 : [2700] ----><-------------------------------------------- [2750]

S: [2750] CAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTG [2800]
 : [2750] -------------><----------------------------------- [2800]

S: [2800] GCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAG [2850]
 : [2800] ----------------------><-------------------------- [2850]

S: [2850] CAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAA [2900]
 : [2850] -------------------------------><----------------- [2900]

S: [2900] CGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGCTCCACCCCG [2950]
 : [2900] ----------------------------------------><-------- [2950]

S: [2950] ACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCTCCATCATGC [3000]
 : [2950] -------------------------------------------------> [3000]

S: [3000] TCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTGACGTCAGCT [3050]
 : [3000] <------------------------------------------------- [3050]

S: [3050] CCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCCCGTGGCCTG [3100]
 : [3050] --------><---------------------------------------- [3100]

S: [3100] ACGTCAGCTCCATCATGCTCCACCCCGACTTAGGAACGAAGGCAGCAGCC [3150]
 : [3100] -----------------><------------------------------- [3150]

S: [3150] CGTGGCCTGACGTCAGCTCCATCATGCATTGGCGACTTCGTTATTGCTAC [3200]
 : [3150] -------------------------->                        [3200]

S: [3200] CACCTGTCTTACAGAGG [3216]
 : [3200]                   [3216]
===================================================================


Model: Suffix Array
Number of matches: 11
Query: GTATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTGACTG
===================================================================
S: [0000] TATTCAATCAATCCGTCGTTTTAACAACACTCAGAAGCCAACTGGTATAT [0050]
 : [0000]                                             <----- [0050]

S: [0050] GGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAG [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] TTTTACTCATCTTGCAATTGGGACGTCCCTGACTGGTATATGGAGAGTAC [0150]
 : [0100] ----------------------------------><-------------- [0150]

S: [0150] CACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCA [0200]
 : [0150] -------------------------------------------------- [0200]

S: [0200] TCTTGCAATTGGGACGTCCCTGACTGGTATATGGAGAGTACCACGCCACA [0250]
 : [0200] -------------------------><----------------------- [0250]

S: [0250] AAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAAT [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] TGGGACGTCCCTGACTGGTATATGGAGAGTACCACGCCACAAAACGTTCT [0350]
 : [0300] ----------------><-------------------------------- [0350]

S: [0350] CGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTC [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] CCTGACTGGTATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGG [0450]
 : [0400] -------><----------------------------------------- [0450]

S: [0450] GCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTGACTGG [0500]
 : [0450] ------------------------------------------------>< [0500]

S: [0500] TATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTA [0550]
 : [0500] -------------------------------------------------- [0550]

S: [0550] ACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTGACTGGTATATGGAG [0600]
 : [0550] ---------------------------------------><--------- [0600]

S: [0600] AGTACCACGCCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTT [0650]
 : [0600] -------------------------------------------------- [0650]

S: [0650] ACTCATCTTGCAATTGGGACGTCCCTGACTGGTATATGGAGAGTACCACG [0700]
 : [0650] ------------------------------><------------------ [0700]

S: [0700] CCACAAAACGTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTT [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] GCAATTGGGACGTCCCTGACTGGTATATGGAGAGTACCACGCCACAAAAC [0800]
 : [0750] ---------------------><--------------------------- [0800]

S: [0800] GTTCTCGTGAACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGG [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] ACGTCCCTGACTGGTATATGGAGAGTACCACGCCACAAAACGTTCTCGTG [0900]
 : [0850] ------------><------------------------------------ [0900]

S: [0900] AACGGGCCGCATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTG [0950]
 : [0900] -------------------------------------------------- [0950]

S: [0950] ACTGGTATATGGAGAGTACCACGCCACAAAACGTTCTCGTGAACGGGCCG [1000]
 : [0950] ---><--------------------------------------------- [1000]

S: [1000] CATTAACAAGTTTTACTCATCTTGCAATTGGGACGTCCCTGACTGCATAG [1050]
 : [1000] -------------------------------------------->      [1050]

S: [1050] AACAGACATACATCCGTATAA [1070]
 : [1050]                       [1070]
===================================================================


Model: Suffix Array
Number of matches: 23
Query: CCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCG
===================================================================
S: [0000] GTATAAAACTGAAGTAGTGAATCCAGTGATCCCGGGTTAAGTCATCATTC [0050]
 : [0000]                               <------------------- [0050]

S: [0050] AACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTT [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] CTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCA [0150]
 : [0100] --------------------------------><---------------- [0150]

S: [0150] TTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAA [0200]
 : [0150] -------------------------------------------------- [0200]

S: [0200] GTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCA [0250]
 : [0200] -----------------------------------><------------- [0250]

S: [0250] TCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACG [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] AAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAG [0350]
 : [0300] --------------------------------------><---------- [0350]

S: [0350] TCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAG [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] ACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTT [0450]
 : [0400] -----------------------------------------><------- [0450]

S: [0450] AAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACT [0500]
 : [0450] -------------------------------------------------- [0500]

S: [0500] AAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGG [0550]
 : [0500] --------------------------------------------><---- [0550]

S: [0550] GTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGC [0600]
 : [0550] -------------------------------------------------- [0600]

S: [0600] ACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCC [0650]
 : [0600] -----------------------------------------------><- [0650]

S: [0650] CGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCT [0700]
 : [0650] -------------------------------------------------- [0700]

S: [0700] CGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCC [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] GCCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGAT [0800]
 : [0750] ><------------------------------------------------ [0800]

S: [0800] CCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCG [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] GCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTG [0900]
 : [0850] ---><--------------------------------------------- [0900]

S: [0900] GATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAG [0950]
 : [0900] -------------------------------------------------- [0950]

S: [0950] TCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTAT [1000]
 : [0950] ------><------------------------------------------ [1000]

S: [1000] TTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGC [1050]
 : [1000] -------------------------------------------------- [1050]

S: [1050] GAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCACGG [1100]
 : [1050] ---------><--------------------------------------- [1100]

S: [1100] TATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGGAAC [1150]
 : [1100] -------------------------------------------------- [1150]

S: [1150] TGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAATACGCA [1200]
 : [1150] ------------><------------------------------------ [1200]

S: [1200] CGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATACTGG [1250]
 : [1200] -------------------------------------------------- [1250]

S: [1250] AACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAATAC [1300]
 : [1250] ---------------><--------------------------------- [1300]

S: [1300] GCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGATAC [1350]
 : [1300] -------------------------------------------------- [1350]

S: [1350] TGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAGTAA [1400]
 : [1350] ------------------><------------------------------ [1400]

S: [1400] TACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACAAGA [1450]
 : [1400] -------------------------------------------------- [1450]

S: [1450] TACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCCGAG [1500]
 : [1450] ---------------------><--------------------------- [1500]

S: [1500] TAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTTACA [1550]
 : [1500] -------------------------------------------------- [1550]

S: [1550] AGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAACCC [1600]
 : [1550] ------------------------><------------------------ [1600]

S: [1600] GAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCTTTT [1650]
 : [1600] -------------------------------------------------- [1650]

S: [1650] ACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATTCAA [1700]
 : [1650] ---------------------------><--------------------- [1700]

S: [1700] CCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGTTCT [1750]
 : [1700] -------------------------------------------------- [1750]

S: [1750] TTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATCATT [1800]
 : [1750] ------------------------------><------------------ [1800]

S: [1800] CAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAAAGT [1850]
 : [1800] -------------------------------------------------- [1850]

S: [1850] TCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTCATC [1900]
 : [1850] ---------------------------------><--------------- [1900]

S: [1900] ATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGACGAA [1950]
 : [1900] -------------------------------------------------- [1950]

S: [1950] AGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAAGTC [2000]
 : [1950] ------------------------------------><------------ [2000]

S: [2000] ATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAAGAC [2050]
 : [2000] -------------------------------------------------- [2050]

S: [2050] GAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGTTAA [2100]
 : [2050] ---------------------------------------><--------- [2100]

S: [2100] GTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCACTAA [2150]
 : [2100] -------------------------------------------------- [2150]

S: [2150] GACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCGGGT [2200]
 : [2150] ------------------------------------------><------ [2200]

S: [2200] TAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCGCAC [2250]
 : [2200] -------------------------------------------------- [2250]

S: [2250] TAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGCCCG [2300]
 : [2250] ---------------------------------------------><--- [2300]

S: [2300] GGTTAAGTCATCATTCAACCCGAGTAATACGCACGGTATTTGGATCCTCG [2350]
 : [2300] -------------------------------------------------- [2350]

S: [2350] CACTAAGACGAAAGTTCTTTTACAAGATACTGGAACTGCGAGTCGGCCGT [2400]
 : [2350] ------------------------------------------------>  [2400]

S: [2400] AGATTCAATTCTGTCTGAG [2418]
 : [2400]                     [2418]
===================================================================


Model: Suffix Array
Number of matches: 19
Query: GTTCTGATTTTACTTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCACGGAAGATTATGAGCACTTCGCGGC
===================================================================
S: [0000] TACTAGCGACATTACGCGTACGTTCTGATTTTACTTGCCATGCAACCACG [0050]
 : [0000]                      <---------------------------- [0050]

S: [0050] CACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGATT [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] TTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCACGG [0150]
 : [0100] -------------------------------------------------- [0150]

S: [0150] AAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCAACCAC [0200]
 : [0150] ---------------------><--------------------------- [0200]

S: [0200] GCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGAT [0250]
 : [0200] -------------------------------------------------- [0250]

S: [0250] TTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCACG [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] GAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCAACCA [0350]
 : [0300] ----------------------><-------------------------- [0350]

S: [0350] CGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAGA [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] TTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCAC [0450]
 : [0400] -------------------------------------------------- [0450]

S: [0450] GGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCAACC [0500]
 : [0450] -----------------------><------------------------- [0500]

S: [0500] ACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTAG [0550]
 : [0500] -------------------------------------------------- [0550]

S: [0550] ATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATCA [0600]
 : [0550] -------------------------------------------------- [0600]

S: [0600] CGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCAAC [0650]
 : [0600] ------------------------><------------------------ [0650]

S: [0650] CACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTTA [0700]
 : [0650] -------------------------------------------------- [0700]

S: [0700] GATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAATC [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] ACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCAA [0800]
 : [0750] -------------------------><----------------------- [0800]

S: [0800] CCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGTT [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] AGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAAT [0900]
 : [0850] -------------------------------------------------- [0900]

S: [0900] CACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGCA [0950]
 : [0900] --------------------------><---------------------- [0950]

S: [0950] ACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAGT [1000]
 : [0950] -------------------------------------------------- [1000]

S: [1000] TAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAAA [1050]
 : [1000] -------------------------------------------------- [1050]

S: [1050] TCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATGC [1100]
 : [1050] ---------------------------><--------------------- [1100]

S: [1100] AACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAAG [1150]
 : [1100] -------------------------------------------------- [1150]

S: [1150] TTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGAA [1200]
 : [1150] -------------------------------------------------- [1200]

S: [1200] ATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCATG [1250]
 : [1200] ----------------------------><-------------------- [1250]

S: [1250] CAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTAA [1300]
 : [1250] -------------------------------------------------- [1300]

S: [1300] GTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATGA [1350]
 : [1300] -------------------------------------------------- [1350]

S: [1350] AATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCAT [1400]
 : [1350] -----------------------------><------------------- [1400]

S: [1400] GCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGTA [1450]
 : [1400] -------------------------------------------------- [1450]

S: [1450] AGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTATG [1500]
 : [1450] -------------------------------------------------- [1500]

S: [1500] AAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCCA [1550]
 : [1500] ------------------------------><------------------ [1550]

S: [1550] TGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAGT [1600]
 : [1550] -------------------------------------------------- [1600]

S: [1600] AAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTAT [1650]
 : [1600] -------------------------------------------------- [1650]

S: [1650] GAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGCC [1700]
 : [1650] -------------------------------><----------------- [1700]

S: [1700] ATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGAG [1750]
 : [1700] -------------------------------------------------- [1750]

S: [1750] TAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATTA [1800]
 : [1750] -------------------------------------------------- [1800]

S: [1800] TGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTGC [1850]
 : [1800] --------------------------------><---------------- [1850]

S: [1850] CATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAGA [1900]
 : [1850] -------------------------------------------------- [1900]

S: [1900] GTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACATT [1950]
 : [1900] -------------------------------------------------- [1950]

S: [1950] ATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTTG [2000]
 : [1950] ---------------------------------><--------------- [2000]

S: [2000] CCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTAG [2050]
 : [2000] -------------------------------------------------- [2050]

S: [2050] AGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACAT [2100]
 : [2050] -------------------------------------------------- [2100]

S: [2100] TATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACTT [2150]
 : [2100] ----------------------------------><-------------- [2150]

S: [2150] GCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCTA [2200]
 : [2150] -------------------------------------------------- [2200]

S: [2200] GAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGACA [2250]
 : [2200] -------------------------------------------------- [2250]

S: [2250] TTATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTACT [2300]
 : [2250] -----------------------------------><------------- [2300]

S: [2300] TGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCCT [2350]
 : [2300] -------------------------------------------------- [2350]

S: [2350] AGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGAC [2400]
 : [2350] -------------------------------------------------- [2400]

S: [2400] ATTATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTAC [2450]
 : [2400] ------------------------------------><------------ [2450]

S: [2450] TTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGCC [2500]
 : [2450] -------------------------------------------------- [2500]

S: [2500] TAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACGA [2550]
 : [2500] -------------------------------------------------- [2550]

S: [2550] CATTATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTTA [2600]
 : [2550] -------------------------------------><----------- [2600]

S: [2600] CTTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTGC [2650]
 : [2600] -------------------------------------------------- [2650]

S: [2650] CTAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTACG [2700]
 : [2650] -------------------------------------------------- [2700]

S: [2700] ACATTATGAAATCACGGAAGATTATGAGCACTTCGCGGCGTTCTGATTTT [2750]
 : [2700] --------------------------------------><---------- [2750]

S: [2750] ACTTGCCATGCAACCACGCACCTATATAATAGAGAACACCCCTCCATGTG [2800]
 : [2750] -------------------------------------------------- [2800]

S: [2800] CCTAGAGTAAGTTAGATTTTCCAAGTACCTTAGTACCTCATGGCTGTTAC [2850]
 : [2800] -------------------------------------------------- [2850]

S: [2850] GACATTATGAAATCACGGAAGATTATGAGCACTTCGCGGCTGCATGCCC [2898]
 : [2850] --------------------------------------->          [2898]
===================================================================


Model: Suffix Array
Number of matches: 13
Query: TAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGC
===================================================================
S: [0000] GTCGCACCAATAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGC [0050]
 : [0000]           <--------------------------------------- [0050]

S: [0050] GGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGG [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] GCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCG [0150]
 : [0100] -------------------------------------------------- [0150]

S: [0150] TACCTGATGCATCACTGATTGATTGGAGAGCTAAACACATTCACGGCATA [0200]
 : [0150] ------------------------------><------------------ [0200]

S: [0200] GCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCA [0250]
 : [0200] -------------------------------------------------- [0250]

S: [0250] GGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACA [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] ATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGA [0350]
 : [0300] -------------------------------------------------- [0350]

S: [0350] GCTAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAA [0400]
 : [0350] -><----------------------------------------------- [0400]

S: [0400] GACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAAC [0450]
 : [0400] -------------------------------------------------- [0450]

S: [0450] AGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGAT [0500]
 : [0450] -------------------------------------------------- [0500]

S: [0500] GCATCACTGATTGATTGGAGAGCTAAACACATTCACGGCATAGCATTTAA [0550]
 : [0500] ----------------------><-------------------------- [0550]

S: [0550] AGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTA [0600]
 : [0550] -------------------------------------------------- [0600]

S: [0600] GAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTAC [0650]
 : [0600] -------------------------------------------------- [0650]

S: [0650] CATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGCTAAACA [0700]
 : [0650] -------------------------------------------><----- [0700]

S: [0700] CATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCC [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] GCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACG [0800]
 : [0750] -------------------------------------------------- [0800]

S: [0800] CAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACT [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] GATTGATTGGAGAGCTAAACACATTCACGGCATAGCATTTAAAGTGACCG [0900]
 : [0850] --------------><---------------------------------- [0900]

S: [0900] CGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCG [0950]
 : [0900] -------------------------------------------------- [0950]

S: [0950] GGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTC [1000]
 : [0950] -------------------------------------------------- [1000]

S: [1000] TGGCGTACCTGATGCATCACTGATTGATTGGAGAGCTAAACACATTCACG [1050]
 : [1000] -----------------------------------><------------- [1050]

S: [1050] GCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGC [1100]
 : [1050] -------------------------------------------------- [1100]

S: [1100] GTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCA [1150]
 : [1100] -------------------------------------------------- [1150]

S: [1150] CTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATT [1200]
 : [1150] -------------------------------------------------- [1200]

S: [1200] GGAGAGCTAAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGC [1250]
 : [1200] ------><------------------------------------------ [1250]

S: [1250] AGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCC [1300]
 : [1250] -------------------------------------------------- [1300]

S: [1300] AGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTAC [1350]
 : [1300] -------------------------------------------------- [1350]

S: [1350] CTGATGCATCACTGATTGATTGGAGAGCTAAACACATTCACGGCATAGCA [1400]
 : [1350] ---------------------------><--------------------- [1400]

S: [1400] TTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGA [1450]
 : [1400] -------------------------------------------------- [1450]

S: [1450] TGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATA [1500]
 : [1450] -------------------------------------------------- [1500]

S: [1500] GCTACCATCCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGCT [1550]
 : [1500] ------------------------------------------------>< [1550]

S: [1550] AAACACATTCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGAC [1600]
 : [1550] -------------------------------------------------- [1600]

S: [1600] TTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGA [1650]
 : [1600] -------------------------------------------------- [1650]

S: [1650] GCACGCAACCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCA [1700]
 : [1650] -------------------------------------------------- [1700]

S: [1700] TCACTGATTGATTGGAGAGCTAAACACATTCACGGCATAGCATTTAAAGT [1750]
 : [1700] -------------------><----------------------------- [1750]

S: [1750] GACCGCGGGCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAA [1800]
 : [1750] -------------------------------------------------- [1800]

S: [1800] GACCGGGCGGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCAT [1850]
 : [1800] -------------------------------------------------- [1850]

S: [1850] CCGTCTGGCGTACCTGATGCATCACTGATTGATTGGAGAGCTAAACACAT [1900]
 : [1850] ----------------------------------------><-------- [1900]

S: [1900] TCACGGCATAGCATTTAAAGTGACCGCGGGCGGCAGTAAGACTTGCCGCG [1950]
 : [1900] -------------------------------------------------- [1950]

S: [1950] CTTGCGTGCAGGATGGTAGAAGACCGGGCGGGCCAGAACAGAGCACGCAA [2000]
 : [1950] -------------------------------------------------- [2000]

S: [2000] CCCCACTACAATAGCTACCATCCGTCTGGCGTACCTGATGCATCACTGAT [2050]
 : [2000] -------------------------------------------------- [2050]

S: [2050] TGATTGGAGAGCTAAACACATTCACGGCATAGCATTTAAAGTGACCGCGG [2100]
 : [2050] -----------><------------------------------------- [2100]

S: [2100] GCGGCAGTAAGACTTGCCGCGCTTGCGTGCAGGATGGTAGAAGACCGGGC [2150]
 : [2100] -------------------------------------------------- [2150]

S: [2150] GGGCCAGAACAGAGCACGCAACCCCACTACAATAGCTACCATCCGTCTGG [2200]
 : [2150] -------------------------------------------------- [2200]

S: [2200] CGTACCTGATGCATCACTGATTGATTGGAGAGCCACAATACTCTTGCTAA [2250]
 : [2200] -------------------------------->                  [2250]

S: [2250] GGA [2252]
 : [2250]     [2252]
===================================================================


Model: Suffix Array
Number of matches: 17
Query: CATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGAC
===================================================================
S: [0000] CGCGCATGCCTGGGCCTGAAGGCAGTCCCCATTTGGCTGACATGGAGCTT [0050]
 : [0000]                                         <--------- [0050]

S: [0050] GTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGGCT [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] AACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCACA [0150]
 : [0100] -------------------------------------------------- [0150]

S: [0150] AACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGACC [0200]
 : [0150] -------------------------------------------------- [0200]

S: [0200] CGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAATGC [0250]
 : [0200] -------------------------------------------------- [0250]

S: [0250] CGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCGGT [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] GCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGACA [0350]
 : [0300] -------------------------------------------------- [0350]

S: [0350] CGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTCGT [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] AGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAAAG [0450]
 : [0400] -------------------------------------------------- [0450]

S: [0450] ACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTTCC [0500]
 : [0450] -------------------------------------------------- [0500]

S: [0500] ATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATGGAGCT [0550]
 : [0500] ----------------------------------------><-------- [0550]

S: [0550] TGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGGC [0600]
 : [0550] -------------------------------------------------- [0600]

S: [0600] TAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCAC [0650]
 : [0600] -------------------------------------------------- [0650]

S: [0650] AAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGAC [0700]
 : [0650] -------------------------------------------------- [0700]

S: [0700] CCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAATG [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] CCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCGG [0800]
 : [0750] -------------------------------------------------- [0800]

S: [0800] TGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGAC [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] ACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTCG [0900]
 : [0850] -------------------------------------------------- [0900]

S: [0900] TAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAAA [0950]
 : [0900] -------------------------------------------------- [0950]

S: [0950] GACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTTC [1000]
 : [0950] -------------------------------------------------- [1000]

S: [1000] CATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATGGAGC [1050]
 : [1000] -----------------------------------------><------- [1050]

S: [1050] TTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAGG [1100]
 : [1050] -------------------------------------------------- [1100]

S: [1100] CTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATCA [1150]
 : [1100] -------------------------------------------------- [1150]

S: [1150] CAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAGA [1200]
 : [1150] -------------------------------------------------- [1200]

S: [1200] CCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAAT [1250]
 : [1200] -------------------------------------------------- [1250]

S: [1250] GCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTCG [1300]
 : [1250] -------------------------------------------------- [1300]

S: [1300] GTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTGA [1350]
 : [1300] -------------------------------------------------- [1350]

S: [1350] CACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTTC [1400]
 : [1350] -------------------------------------------------- [1400]

S: [1400] GTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATAA [1450]
 : [1400] -------------------------------------------------- [1450]

S: [1450] AGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCTT [1500]
 : [1450] -------------------------------------------------- [1500]

S: [1500] CCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATGGAG [1550]
 : [1500] ------------------------------------------><------ [1550]

S: [1550] CTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGAG [1600]
 : [1550] -------------------------------------------------- [1600]

S: [1600] GCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACATC [1650]
 : [1600] -------------------------------------------------- [1650]

S: [1650] ACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAAG [1700]
 : [1650] -------------------------------------------------- [1700]

S: [1700] ACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTAA [1750]
 : [1700] -------------------------------------------------- [1750]

S: [1750] TGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACTC [1800]
 : [1750] -------------------------------------------------- [1800]

S: [1800] GGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTTG [1850]
 : [1800] -------------------------------------------------- [1850]

S: [1850] ACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCTT [1900]
 : [1850] -------------------------------------------------- [1900]

S: [1900] CGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTATA [1950]
 : [1900] -------------------------------------------------- [1950]

S: [1950] AAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCCT [2000]
 : [1950] -------------------------------------------------- [2000]

S: [2000] TCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATGGA [2050]
 : [2000] -------------------------------------------><----- [2050]

S: [2050] GCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTGA [2100]
 : [2050] -------------------------------------------------- [2100]

S: [2100] GGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACAT [2150]
 : [2100] -------------------------------------------------- [2150]

S: [2150] CACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCAA [2200]
 : [2150] -------------------------------------------------- [2200]

S: [2200] GACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCTA [2250]
 : [2200] -------------------------------------------------- [2250]

S: [2250] ATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAACT [2300]
 : [2250] -------------------------------------------------- [2300]

S: [2300] CGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGTT [2350]
 : [2300] -------------------------------------------------- [2350]

S: [2350] GACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGCT [2400]
 : [2350] -------------------------------------------------- [2400]

S: [2400] TCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTAT [2450]
 : [2400] -------------------------------------------------- [2450]

S: [2450] AAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCCC [2500]
 : [2450] -------------------------------------------------- [2500]

S: [2500] TTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATGG [2550]
 : [2500] --------------------------------------------><---- [2550]

S: [2550] AGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCTG [2600]
 : [2550] -------------------------------------------------- [2600]

S: [2600] AGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGACA [2650]
 : [2600] -------------------------------------------------- [2650]

S: [2650] TCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGCA [2700]
 : [2650] -------------------------------------------------- [2700]

S: [2700] AGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACCT [2750]
 : [2700] -------------------------------------------------- [2750]

S: [2750] AATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAAC [2800]
 : [2750] -------------------------------------------------- [2800]

S: [2800] TCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGGT [2850]
 : [2800] -------------------------------------------------- [2850]

S: [2850] TGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGGC [2900]
 : [2850] -------------------------------------------------- [2900]

S: [2900] TTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGTA [2950]
 : [2900] -------------------------------------------------- [2950]

S: [2950] TAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGCC [3000]
 : [2950] -------------------------------------------------- [3000]

S: [3000] CTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCATG [3050]
 : [3000] ---------------------------------------------><--- [3050]

S: [3050] GAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTCT [3100]
 : [3050] -------------------------------------------------- [3100]

S: [3100] GAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGAC [3150]
 : [3100] -------------------------------------------------- [3150]

S: [3150] ATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCGC [3200]
 : [3150] -------------------------------------------------- [3200]

S: [3200] AAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGACC [3250]
 : [3200] -------------------------------------------------- [3250]

S: [3250] TAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACAA [3300]
 : [3250] -------------------------------------------------- [3300]

S: [3300] CTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAGG [3350]
 : [3300] -------------------------------------------------- [3350]

S: [3350] TTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTGG [3400]
 : [3350] -------------------------------------------------- [3400]

S: [3400] CTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAGT [3450]
 : [3400] -------------------------------------------------- [3450]

S: [3450] ATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCGC [3500]
 : [3450] -------------------------------------------------- [3500]

S: [3500] CCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCAT [3550]
 : [3500] ----------------------------------------------><-- [3550]

S: [3550] GGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCTC [3600]
 : [3550] -------------------------------------------------- [3600]

S: [3600] TGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGGA [3650]
 : [3600] -------------------------------------------------- [3650]

S: [3650] CATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGCG [3700]
 : [3650] -------------------------------------------------- [3700]

S: [3700] CAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGAC [3750]
 : [3700] -------------------------------------------------- [3750]

S: [3750] CTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAACA [3800]
 : [3750] -------------------------------------------------- [3800]

S: [3800] ACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCAG [3850]
 : [3800] -------------------------------------------------- [3850]

S: [3850] GTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCTG [3900]
 : [3850] -------------------------------------------------- [3900]

S: [3900] GCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAAG [3950]
 : [3900] -------------------------------------------------- [3950]

S: [3950] TATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTCG [4000]
 : [3950] -------------------------------------------------- [4000]

S: [4000] CCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACCA [4050]
 : [4000] -----------------------------------------------><- [4050]

S: [4050] TGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTCT [4100]
 : [4050] -------------------------------------------------- [4100]

S: [4100] CTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTGG [4150]
 : [4100] -------------------------------------------------- [4150]

S: [4150] ACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAGC [4200]
 : [4150] -------------------------------------------------- [4200]

S: [4200] GCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAGA [4250]
 : [4200] -------------------------------------------------- [4250]

S: [4250] CCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAAC [4300]
 : [4250] -------------------------------------------------- [4300]

S: [4300] AACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCCA [4350]
 : [4300] -------------------------------------------------- [4350]

S: [4350] GGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGCT [4400]
 : [4350] -------------------------------------------------- [4400]

S: [4400] GGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGAA [4450]
 : [4400] -------------------------------------------------- [4450]

S: [4450] GTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCTC [4500]
 : [4450] -------------------------------------------------- [4500]

S: [4500] GCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGACC [4550]
 : [4500] ------------------------------------------------>< [4550]

S: [4550] ATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTTC [4600]
 : [4550] -------------------------------------------------- [4600]

S: [4600] TCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCTG [4650]
 : [4600] -------------------------------------------------- [4650]

S: [4650] GACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAAG [4700]
 : [4650] -------------------------------------------------- [4700]

S: [4700] CGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCAG [4750]
 : [4700] -------------------------------------------------- [4750]

S: [4750] ACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACAA [4800]
 : [4750] -------------------------------------------------- [4800]

S: [4800] CAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCCC [4850]
 : [4800] -------------------------------------------------- [4850]

S: [4850] AGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAGC [4900]
 : [4850] -------------------------------------------------- [4900]

S: [4900] TGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACGA [4950]
 : [4900] -------------------------------------------------- [4950]

S: [4950] AGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGCT [5000]
 : [4950] -------------------------------------------------- [5000]

S: [5000] CGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGAC [5050]
 : [5000] -------------------------------------------------> [5050]

S: [5050] CATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTTT [5100]
 : [5050] <------------------------------------------------- [5100]

S: [5100] CTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTCT [5150]
 : [5100] -------------------------------------------------- [5150]

S: [5150] GGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAAA [5200]
 : [5150] -------------------------------------------------- [5200]

S: [5200] GCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCCA [5250]
 : [5200] -------------------------------------------------- [5250]

S: [5250] GACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGACA [5300]
 : [5250] -------------------------------------------------- [5300]

S: [5300] ACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGCC [5350]
 : [5300] -------------------------------------------------- [5350]

S: [5350] CAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATAG [5400]
 : [5350] -------------------------------------------------- [5400]

S: [5400] CTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCACG [5450]
 : [5400] -------------------------------------------------- [5450]

S: [5450] AAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCGC [5500]
 : [5450] -------------------------------------------------- [5500]

S: [5500] TCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTGA [5550]
 : [5500] -------------------------------------------------- [5550]

S: [5550] CCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGTT [5600]
 : [5550] ><------------------------------------------------ [5600]

S: [5600] TCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTTC [5650]
 : [5600] -------------------------------------------------- [5650]

S: [5650] TGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATAA [5700]
 : [5650] -------------------------------------------------- [5700]

S: [5700] AGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGCC [5750]
 : [5700] -------------------------------------------------- [5750]

S: [5750] AGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGAC [5800]
 : [5750] -------------------------------------------------- [5800]

S: [5800] AACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAGC [5850]
 : [5800] -------------------------------------------------- [5850]

S: [5850] CCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCATA [5900]
 : [5850] -------------------------------------------------- [5900]

S: [5900] GCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCAC [5950]
 : [5900] -------------------------------------------------- [5950]

S: [5950] GAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGCG [6000]
 : [5950] -------------------------------------------------- [6000]

S: [6000] CTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGTG [6050]
 : [6000] -------------------------------------------------- [6050]

S: [6050] ACCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACGT [6100]
 : [6050] -><----------------------------------------------- [6100]

S: [6100] TTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGTT [6150]
 : [6100] -------------------------------------------------- [6150]

S: [6150] CTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAATA [6200]
 : [6150] -------------------------------------------------- [6200]

S: [6200] AAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGGC [6250]
 : [6200] -------------------------------------------------- [6250]

S: [6250] CAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTGA [6300]
 : [6250] -------------------------------------------------- [6300]

S: [6300] CAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAAG [6350]
 : [6300] -------------------------------------------------- [6350]

S: [6350] CCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCAT [6400]
 : [6350] -------------------------------------------------- [6400]

S: [6400] AGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGCA [6450]
 : [6400] -------------------------------------------------- [6450]

S: [6450] CGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGGC [6500]
 : [6450] -------------------------------------------------- [6500]

S: [6500] GCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGGT [6550]
 : [6500] -------------------------------------------------- [6550]

S: [6550] GACCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACACG [6600]
 : [6550] --><---------------------------------------------- [6600]

S: [6600] TTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGGT [6650]
 : [6600] -------------------------------------------------- [6650]

S: [6650] TCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAAT [6700]
 : [6650] -------------------------------------------------- [6700]

S: [6700] AAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTGG [6750]
 : [6700] -------------------------------------------------- [6750]

S: [6750] CCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATTG [6800]
 : [6750] -------------------------------------------------- [6800]

S: [6800] ACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCAA [6850]
 : [6800] -------------------------------------------------- [6850]

S: [6850] GCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATCA [6900]
 : [6850] -------------------------------------------------- [6900]

S: [6900] TAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTGC [6950]
 : [6900] -------------------------------------------------- [6950]

S: [6950] ACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTGG [7000]
 : [6950] -------------------------------------------------- [7000]

S: [7000] CGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAGG [7050]
 : [7000] -------------------------------------------------- [7050]

S: [7050] TGACCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACAC [7100]
 : [7050] ---><--------------------------------------------- [7100]

S: [7100] GTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACGG [7150]
 : [7100] -------------------------------------------------- [7150]

S: [7150] TTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGAA [7200]
 : [7150] -------------------------------------------------- [7200]

S: [7200] TAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTTG [7250]
 : [7200] -------------------------------------------------- [7250]

S: [7250] GCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAATT [7300]
 : [7250] -------------------------------------------------- [7300]

S: [7300] GACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGCA [7350]
 : [7300] -------------------------------------------------- [7350]

S: [7350] AGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGATC [7400]
 : [7350] -------------------------------------------------- [7400]

S: [7400] ATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTTG [7450]
 : [7400] -------------------------------------------------- [7450]

S: [7450] CACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCTG [7500]
 : [7450] -------------------------------------------------- [7500]

S: [7500] GCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGAG [7550]
 : [7500] -------------------------------------------------- [7550]

S: [7550] GTGACCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAACA [7600]
 : [7550] ----><-------------------------------------------- [7600]

S: [7600] CGTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAACG [7650]
 : [7600] -------------------------------------------------- [7650]

S: [7650] GTTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTGA [7700]
 : [7650] -------------------------------------------------- [7700]

S: [7700] ATAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGTT [7750]
 : [7700] -------------------------------------------------- [7750]

S: [7750] GGCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAAT [7800]
 : [7750] -------------------------------------------------- [7800]

S: [7800] TGACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCGC [7850]
 : [7800] -------------------------------------------------- [7850]

S: [7850] AAGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGAT [7900]
 : [7850] -------------------------------------------------- [7900]

S: [7900] CATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCTT [7950]
 : [7900] -------------------------------------------------- [7950]

S: [7950] GCACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCCT [8000]
 : [7950] -------------------------------------------------- [8000]

S: [8000] GGCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGGA [8050]
 : [8000] -------------------------------------------------- [8050]

S: [8050] GGTGACCATGGAGCTTGTCAAGTAGTCCGGGGGCTACCAACGCATCCAAC [8100]
 : [8050] -----><------------------------------------------- [8100]

S: [8100] ACGTTTCTCTGAGGCTAACCCCGTAAACAGTCTAGCGCAAAAAACGTAAC [8150]
 : [8100] -------------------------------------------------- [8150]

S: [8150] GGTTCTGGACATCACAAACACAACACACTTAGAGCCAGTTTCTGGACTTG [8200]
 : [8150] -------------------------------------------------- [8200]

S: [8200] AATAAAGCGCAAGACCCGTATATCCAACTTCTCCCGCTCAGTCGCGGAGT [8250]
 : [8200] -------------------------------------------------- [8250]

S: [8250] TGGCCAGACCTAATGCCGTGTGATGTCTAGCTTAGCATACTAAGTTCCAA [8300]
 : [8250] -------------------------------------------------- [8300]

S: [8300] TTGACAACAACTCGGTGCCCGCCTACTATTCCCCGTCACTTAGCTACCCG [8350]
 : [8300] -------------------------------------------------- [8350]

S: [8350] CAAGCCCAGGTTGACACGTAGGTGGCGGGCGGCGGGACCGGGGACAGCGA [8400]
 : [8350] -------------------------------------------------- [8400]

S: [8400] TCATAGCTGGCTTCGTAGAATGCAGGCTGTCGGACCGTATAGCGTTGCCT [8450]
 : [8400] -------------------------------------------------- [8450]

S: [8450] TGCACGAAGTATAAAGACGGTTCGGAAAAAATTAACGTGGGAGATGATCC [8500]
 : [8450] -------------------------------------------------- [8500]

S: [8500] TGGCGCTCGCCCTTCCATCATTCTGTGAGAATCGAGGTTTATACTATAGG [8550]
 : [8500] -------------------------------------------------- [8550]

S: [8550] AGGTGACTACCGAAATGCTCGTAGGATCAAATGTGTG [8586]
 : [8550] ------>                               [8586]
===================================================================


Model: Suffix Array
Number of matches: 7
Query: ATTAGTACCCCATCGGTTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTGCTACAGTCTTACTCCCCCTATAGCTGGGTGACT
===================================================================
S: [0000] ATGGCAAATCTTGGGCAAGTAGCTGCGCGCATTAGTACCCCATCGGTTTT [0050]
 : [0000]                               <------------------- [0050]

S: [0050] ATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCGGA [0100]
 : [0050] -------------------------------------------------- [0100]

S: [0100] GATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGTAC [0150]
 : [0100] -------------------------------------------------- [0150]

S: [0150] TAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTTCA [0200]
 : [0150] -------------------------------------------------- [0200]

S: [0200] ACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGACCG [0250]
 : [0200] -------------------------------------------------- [0250]

S: [0250] AGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAAAG [0300]
 : [0250] -------------------------------------------------- [0300]

S: [0300] CCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGCGC [0350]
 : [0300] -------------------------------------------------- [0350]

S: [0350] TATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGGCT [0400]
 : [0350] -------------------------------------------------- [0400]

S: [0400] TATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTTGT [0450]
 : [0400] -------------------------------------------------- [0450]

S: [0450] TAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCAGT [0500]
 : [0450] -------------------------------------------------- [0500]

S: [0500] TTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTGTT [0550]
 : [0500] -------------------------------------------------- [0550]

S: [0550] CAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATCCC [0600]
 : [0550] -------------------------------------------------- [0600]

S: [0600] TTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTATT [0650]
 : [0600] -------------------------------------------------- [0650]

S: [0650] TGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACCAA [0700]
 : [0650] -------------------------------------------------- [0700]

S: [0700] GGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACACGG [0750]
 : [0700] -------------------------------------------------- [0750]

S: [0750] ACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGTCA [0800]
 : [0750] -------------------------------------------------- [0800]

S: [0800] TTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAGAT [0850]
 : [0800] -------------------------------------------------- [0850]

S: [0850] CTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGCTG [0900]
 : [0850] -------------------------------------------------- [0900]

S: [0900] GCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAACT [0950]
 : [0900] -------------------------------------------------- [0950]

S: [0950] ATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTGCT [1000]
 : [0950] -------------------------------------------------- [1000]

S: [1000] ACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATCGGTTT [1050]
 : [1000] ------------------------------><------------------ [1050]

S: [1050] TATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCGG [1100]
 : [1050] -------------------------------------------------- [1100]

S: [1100] AGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGTA [1150]
 : [1100] -------------------------------------------------- [1150]

S: [1150] CTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTTC [1200]
 : [1150] -------------------------------------------------- [1200]

S: [1200] AACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGACC [1250]
 : [1200] -------------------------------------------------- [1250]

S: [1250] GAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAAA [1300]
 : [1250] -------------------------------------------------- [1300]

S: [1300] GCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGCG [1350]
 : [1300] -------------------------------------------------- [1350]

S: [1350] CTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGGC [1400]
 : [1350] -------------------------------------------------- [1400]

S: [1400] TTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTTG [1450]
 : [1400] -------------------------------------------------- [1450]

S: [1450] TTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCAG [1500]
 : [1450] -------------------------------------------------- [1500]

S: [1500] TTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTGT [1550]
 : [1500] -------------------------------------------------- [1550]

S: [1550] TCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATCC [1600]
 : [1550] -------------------------------------------------- [1600]

S: [1600] CTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTAT [1650]
 : [1600] -------------------------------------------------- [1650]

S: [1650] TTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACCA [1700]
 : [1650] -------------------------------------------------- [1700]

S: [1700] AGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACACG [1750]
 : [1700] -------------------------------------------------- [1750]

S: [1750] GACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGTC [1800]
 : [1750] -------------------------------------------------- [1800]

S: [1800] ATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAGA [1850]
 : [1800] -------------------------------------------------- [1850]

S: [1850] TCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGCT [1900]
 : [1850] -------------------------------------------------- [1900]

S: [1900] GGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAAC [1950]
 : [1900] -------------------------------------------------- [1950]

S: [1950] TATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTGC [2000]
 : [1950] -------------------------------------------------- [2000]

S: [2000] TACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATCGGTT [2050]
 : [2000] -------------------------------><----------------- [2050]

S: [2050] TTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTCG [2100]
 : [2050] -------------------------------------------------- [2100]

S: [2100] GAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCGT [2150]
 : [2100] -------------------------------------------------- [2150]

S: [2150] ACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACTT [2200]
 : [2150] -------------------------------------------------- [2200]

S: [2200] CAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGAC [2250]
 : [2200] -------------------------------------------------- [2250]

S: [2250] CGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAAA [2300]
 : [2250] -------------------------------------------------- [2300]

S: [2300] AGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGGC [2350]
 : [2300] -------------------------------------------------- [2350]

S: [2350] GCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGGG [2400]
 : [2350] -------------------------------------------------- [2400]

S: [2400] CTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTTT [2450]
 : [2400] -------------------------------------------------- [2450]

S: [2450] GTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGCA [2500]
 : [2450] -------------------------------------------------- [2500]

S: [2500] GTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCTG [2550]
 : [2500] -------------------------------------------------- [2550]

S: [2550] TTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGATC [2600]
 : [2550] -------------------------------------------------- [2600]

S: [2600] CCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGTA [2650]
 : [2600] -------------------------------------------------- [2650]

S: [2650] TTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAACC [2700]
 : [2650] -------------------------------------------------- [2700]

S: [2700] AAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACAC [2750]
 : [2700] -------------------------------------------------- [2750]

S: [2750] GGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGGT [2800]
 : [2750] -------------------------------------------------- [2800]

S: [2800] CATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAAG [2850]
 : [2800] -------------------------------------------------- [2850]

S: [2850] ATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTGC [2900]
 : [2850] -------------------------------------------------- [2900]

S: [2900] TGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGAA [2950]
 : [2900] -------------------------------------------------- [2950]

S: [2950] CTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGTG [3000]
 : [2950] -------------------------------------------------- [3000]

S: [3000] CTACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATCGGT [3050]
 : [3000] --------------------------------><---------------- [3050]

S: [3050] TTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTTC [3100]
 : [3050] -------------------------------------------------- [3100]

S: [3100] GGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACCG [3150]
 : [3100] -------------------------------------------------- [3150]

S: [3150] TACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTACT [3200]
 : [3150] -------------------------------------------------- [3200]

S: [3200] TCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTGA [3250]
 : [3200] -------------------------------------------------- [3250]

S: [3250] CCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTAA [3300]
 : [3250] -------------------------------------------------- [3300]

S: [3300] AAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGGG [3350]
 : [3300] -------------------------------------------------- [3350]

S: [3350] CGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTGG [3400]
 : [3350] -------------------------------------------------- [3400]

S: [3400] GCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGTT [3450]
 : [3400] -------------------------------------------------- [3450]

S: [3450] TGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTGC [3500]
 : [3450] -------------------------------------------------- [3500]

S: [3500] AGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTCT [3550]
 : [3500] -------------------------------------------------- [3550]

S: [3550] GTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGAT [3600]
 : [3550] -------------------------------------------------- [3600]

S: [3600] CCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTGT [3650]
 : [3600] -------------------------------------------------- [3650]

S: [3650] ATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAAC [3700]
 : [3650] -------------------------------------------------- [3700]

S: [3700] CAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAACA [3750]
 : [3700] -------------------------------------------------- [3750]

S: [3750] CGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCGG [3800]
 : [3750] -------------------------------------------------- [3800]

S: [3800] TCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTAA [3850]
 : [3800] -------------------------------------------------- [3850]

S: [3850] GATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCTG [3900]
 : [3850] -------------------------------------------------- [3900]

S: [3900] CTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCGA [3950]
 : [3900] -------------------------------------------------- [3950]

S: [3950] ACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAGT [4000]
 : [3950] -------------------------------------------------- [4000]

S: [4000] GCTACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATCGG [4050]
 : [4000] ---------------------------------><--------------- [4050]

S: [4050] TTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCTT [4100]
 : [4050] -------------------------------------------------- [4100]

S: [4100] CGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTACC [4150]
 : [4100] -------------------------------------------------- [4150]

S: [4150] GTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTAC [4200]
 : [4150] -------------------------------------------------- [4200]

S: [4200] TTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCTG [4250]
 : [4200] -------------------------------------------------- [4250]

S: [4250] ACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGTA [4300]
 : [4250] -------------------------------------------------- [4300]

S: [4300] AAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACGG [4350]
 : [4300] -------------------------------------------------- [4350]

S: [4350] GCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGTG [4400]
 : [4350] -------------------------------------------------- [4400]

S: [4400] GGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAGT [4450]
 : [4400] -------------------------------------------------- [4450]

S: [4450] TTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCTG [4500]
 : [4450] -------------------------------------------------- [4500]

S: [4500] CAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATTC [4550]
 : [4500] -------------------------------------------------- [4550]

S: [4550] TGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTGA [4600]
 : [4550] -------------------------------------------------- [4600]

S: [4600] TCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTTG [4650]
 : [4600] -------------------------------------------------- [4650]

S: [4650] TATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGAA [4700]
 : [4650] -------------------------------------------------- [4700]

S: [4700] CCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAAC [4750]
 : [4700] -------------------------------------------------- [4750]

S: [4750] ACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCCG [4800]
 : [4750] -------------------------------------------------- [4800]

S: [4800] GTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGTA [4850]
 : [4800] -------------------------------------------------- [4850]

S: [4850] AGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTCT [4900]
 : [4850] -------------------------------------------------- [4900]

S: [4900] GCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCCG [4950]
 : [4900] -------------------------------------------------- [4950]

S: [4950] AACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAAG [5000]
 : [4950] -------------------------------------------------- [5000]

S: [5000] TGCTACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATCG [5050]
 : [5000] ----------------------------------><-------------- [5050]

S: [5050] GTTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGCT [5100]
 : [5050] -------------------------------------------------- [5100]

S: [5100] TCGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTAC [5150]
 : [5100] -------------------------------------------------- [5150]

S: [5150] CGTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCTA [5200]
 : [5150] -------------------------------------------------- [5200]

S: [5200] CTTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTCT [5250]
 : [5200] -------------------------------------------------- [5250]

S: [5250] GACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCGT [5300]
 : [5250] -------------------------------------------------- [5300]

S: [5300] AAAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGACG [5350]
 : [5300] -------------------------------------------------- [5350]

S: [5350] GGCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGGT [5400]
 : [5350] -------------------------------------------------- [5400]

S: [5400] GGGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGAG [5450]
 : [5400] -------------------------------------------------- [5450]

S: [5450] TTTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACCT [5500]
 : [5450] -------------------------------------------------- [5500]

S: [5500] GCAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAATT [5550]
 : [5500] -------------------------------------------------- [5550]

S: [5550] CTGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGTG [5600]
 : [5550] -------------------------------------------------- [5600]

S: [5600] ATCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCTT [5650]
 : [5600] -------------------------------------------------- [5650]

S: [5650] GTATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAGA [5700]
 : [5650] -------------------------------------------------- [5700]

S: [5700] ACCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCAA [5750]
 : [5700] -------------------------------------------------- [5750]

S: [5750] CACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGCC [5800]
 : [5750] -------------------------------------------------- [5800]

S: [5800] GGTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGGT [5850]
 : [5800] -------------------------------------------------- [5850]

S: [5850] AAGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCTC [5900]
 : [5850] -------------------------------------------------- [5900]

S: [5900] TGCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGCC [5950]
 : [5900] -------------------------------------------------- [5950]

S: [5950] GAACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGAA [6000]
 : [5950] -------------------------------------------------- [6000]

S: [6000] GTGCTACAGTCTTACTCCCCCTATAGCTGGGTGACTATTAGTACCCCATC [6050]
 : [6000] -----------------------------------><------------- [6050]

S: [6050] GGTTTTATCAGCGACGTCCACAGACTGCTCATGACATCAGGTGAACTAGC [6100]
 : [6050] -------------------------------------------------- [6100]

S: [6100] TTCGGAGATAGTTTTGGAACATTTTGCACACACCATAAGACGTTTATTTA [6150]
 : [6100] -------------------------------------------------- [6150]

S: [6150] CCGTACTAGGTAGGCATAGCCGGCAATCGGCGGCGGAGGAGATCGCCTCT [6200]
 : [6150] -------------------------------------------------- [6200]

S: [6200] ACTTCAACCCCAGTCGCTTAGACCATGGTTTCACTGTTTCCGGCAAGTTC [6250]
 : [6200] -------------------------------------------------- [6250]

S: [6250] TGACCGAGATAGATCATCCATATTTTGTAGTAGTTGTAAGAAGCGTAGCG [6300]
 : [6250] -------------------------------------------------- [6300]

S: [6300] TAAAAGCCGGAGGCGGCTCTGCGCATGCCGGAGTGCACTTAGAAAGTGAC [6350]
 : [6300] -------------------------------------------------- [6350]

S: [6350] GGGCGCTATTGACAAGGCTACGTATCCCAATGTATGTTAAGTCCCCTGGG [6400]
 : [6350] -------------------------------------------------- [6400]

S: [6400] TGGGCTTATCTATTATGAGAGACTAGTTACGGACCACATATTGGTGGCGA [6450]
 : [6400] -------------------------------------------------- [6450]

S: [6450] GTTTGTTAGGAGGAAGTAACCGCCGGGGCGACTATCGGATCGAGTCTACC [6500]
 : [6450] -------------------------------------------------- [6500]

S: [6500] TGCAGTTTGAACGTCTTAAAGTGCACTAGGAGAATAACTAATAGGCAAAT [6550]
 : [6500] -------------------------------------------------- [6550]

S: [6550] TCTGTTCAGGATTAAAATGTCTAGAATTCGTACAACCGATTACTTCCCGT [6600]
 : [6550] -------------------------------------------------- [6600]

S: [6600] GATCCCTTGCTATGCTATTGAAGCAGATGCTGCTAAGTCCATTGGTTGCT [6650]
 : [6600] -------------------------------------------------- [6650]

S: [6650] TGTATTTGACGGAACCCAGCGGGAGGCGAATTAGCTTTCTGTATTGATAG [6700]
 : [6650] -------------------------------------------------- [6700]

S: [6700] AACCAAGGCCCCATCTGCCTCCGACGCTGTTTCCCCGGACAATGATTCCA [6750]
 : [6700] -------------------------------------------------- [6750]

S: [6750] ACACGGACACTTGCGCCACATGCGCATTGGAGCTCCAGCGACACGGACGC [6800]
 : [6750] -------------------------------------------------- [6800]

S: [6800] CGGTCATTGAACCGACTGGCTGGACGTCCCGTCGCCAAACGATCGTATGG [6850]
 : [6800] -------------------------------------------------- [6850]

S: [6850] TAAGATCTCTCCCCTGTCTCTGAGGAGGTTCCTGGGTATAGCATACGCCT [6900]
 : [6850] -------------------------------------------------- [6900]

S: [6900] CTGCTGGCGGCCGTAGATTTCTTTTGACCTCCAAGATACGGTAACAATGC [6950]
 : [6900] -------------------------------------------------- [6950]

S: [6950] CGAACTATGATGTACGACCAGTGGTACTCACATCTTAGCCCCTACGATGA [7000]
 : [6950] -------------------------------------------------- [7000]

S: [7000] AGTGCTACAGTCTTACTCCCCCTATAGCTGGGTGACTGTTCCACTCCTGA [7050]
 : [7000] ------------------------------------>              [7050]

S: [7050] AGTCAGAGGGATCCCAAGTCTGGTGAT [7076]
 : [7050]                             [7076]
===================================================================

演習問題2(オプション)

　記号列 $q$ が正確なタンデムリピート $p^k$ を部分列として含むとする。$q$ に対して、ランダムに1塩基の変異・挿入・削除を加えた配列を $q^{\prime}$ とする。$q^{\prime}$ から $p^k$ を予測するアルゴリズムを考えよ。

　例題として $10$ 個の配列を与えるので、各配列から $p$ と $k$ を決定せよ（$k$ を正確に予測するのは難しいので問わない）

解答

In [9]:

raise NotImplementedError("Not Implemented.")

---------------------------------------------------------------------------
NotImplementedError                       Traceback (most recent call last)
<ipython-input-14-7c6745d17b09> in <module>
----> 1 raise NotImplementedError("Not Implemented.")

NotImplementedError: Not Implemented.

In [ ]:

1.繰返し構造

演習問題1(必須)

解答

Results¶

find tandem repeat DP¶

find tandem repeat using Suffix Array¶

display tandem repeat¶

演習問題2(オプション)

解答

Table of Contents

Published

Last Updated

Category

Tags

Contact

Other contents