Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,558,699 Δ6,790 bp intZ[eutA] intZ, yffL, yffM, yffN, yffO, yffP, yffQ, yffR, yffS, [eutA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2558699 2565488 6790 34 [0] [4] 38 intZ–[eutA] intZ,yffL,yffM,yffN,yffO,yffP,yffQ,yffR,yffS,[eutA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25586980 (0.000)34 (1.400) 30/476 0.1 100% intergenic (‑19/‑160) eutB/intZ ethanolamine ammonia‑lyase subunit alpha/putative phage integrase IntZ
?NC_000913 2565489 = 0 (0.000)coding (1396/1404 nt) eutA ethanolamine ammonia‑lyase reactivase EutA

AGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGTACCGCGGC  >  NC_000913/2565239‑2565738
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
ctgaatcagccgcgtcacgcagtcatcttcataggcaataaccggattattgcggatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagg                                                                                                                                                                                                                                                           >  1:106208‑M2/251‑251 (MQ=255)
        gccgcgtcacgcagtcatcttcataggcaatcaccggattattgcggatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCg                                                                                                                                                                                                                                                    >  1:183821‑M2/243‑250 (MQ=255)
                ccgcagtcatcttcataggcaatccccggattatggcggatgcccgctacggtcattcccgacaacccctgctttcccgccccgcgcccctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtc                                                                                                                                                                                                                                           <  1:205950‑M2/17‑1 (MQ=255)
                    agtcatcttcataggcaataaccggattaatgcggttgtccgctacggtcatttccgacaacacctgctttaccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcccatcccccgaacgcagttcgctggctttccccagcccctcttttacctcctaaaactgacataccttgccgaaccatgtggtctttagtttcctaagtcgttcactcgggaagGAAATGCGAGTGATTtcaccg                                                                                                                                                                                                                                       >  2:244187‑M2/231‑251 (MQ=255)
                            taataggcaatcacgggattattgcggatgtccgctacggttattccccacaacccctgcttggccgccacgctctcctgtgagtttgccgctgcggcgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgtgccctcaggaagGAAATGCGAGTGATTtcaccgtcaccg                                                                                                                                                                                                                                 <  2:23845‑M2/27‑1 (MQ=255)
                               taggcaatcaccggattattgcggatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctatgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGcac                                                                                                                                                                                                                             >  2:119040‑M2/220‑250 (MQ=255)
                                            gattattgcggatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCa                                                                                                                                                                                                                >  2:162417‑M2/207‑250 (MQ=255)
                                            gattattgcggatgtccgctaaggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCa                                                                                                                                                                                                                >  1:80667‑M2/207‑250 (MQ=255)
                                                    cggatgcccgctacggtcatttccgacaccccctgctttgccgccacgcgctcctgtgagcgtgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGc                                                                                                                                                                                                       <  1:243277‑M2/53‑1 (MQ=255)
                                                      gatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGt                                                                                                                                                                                                     >  1:242867‑M2/197‑251 (MQ=255)
                                                      gatgtccgctacggtcatttccgacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGt                                                                                                                                                                                                     >  1:242873‑M2/197‑251 (MQ=255)
                                                                            gacaacacctgctttgccgccacgcgctcctgtgagcttgccgctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTccccc                                                                                                                                                                               >  2:196040‑M2/175‑251 (MQ=255)
                                                                                                                      gctgcgacgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAgttgt                                                                                                                                     >  1:167335‑M2/133‑251 (MQ=255)
                                                                                                                            ccgcccgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGc                                                                                                                               <  1:86151‑M2/125‑1 (MQ=255)
                                                                                                                                 cgcccgcacatcccccgaacgcagttcgttggctttagccagcacctcgtttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCCCCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                                                                                                                          <  2:156595‑M2/130‑1 (MQ=255)
                                                                                                                                 cgccagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTgg                                                                                                                          <  1:20601‑M2/130‑1 (MQ=255)
                                                                                                                                    cagcacatcccccgaacgcagttcgttggctttagccagcacctcttttacgtccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCCCAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCtgtggccg                                                                                                                       >  2:149696‑M2/119‑248 (MQ=255)
                                                                                                                                                          tcgttggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCttt                                                                                                  <  2:187414‑M2/154‑1 (MQ=255)
                                                                                                                                                               ggctttagccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCaa                                                                                             >  2:59945‑M2/92‑250 (MQ=255)
                                                                                                                                                                     agccagcacctcttttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTg                                                                                      >  1:167928‑M2/86‑251 (MQ=255)
                                                                                                                                                                                  tttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGccaccac                                                                          <  1:196040‑M2/178‑1 (MQ=255)
                                                                                                                                                                                   ttacatccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGccaccacc                                                                         >  2:91069‑M2/72‑250 (MQ=255)
                                                                                                                                                                                       atccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACcagcag                                                                    <  2:108414‑M2/184‑1 (MQ=255)
                                                                                                                                                                                        tccttaaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAgg                                                                   >  1:193227‑M2/67‑251 (MQ=255)
                                                                                                                                                                                            taaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGa                                                                <  1:184638‑M2/188‑1 (MQ=255)
                                                                                                                                                                                            taaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGAt                                                               >  1:89076‑M2/63‑251 (MQ=255)
                                                                                                                                                                                            taaactgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGAt                                                               >  1:3147‑M2/63‑251 (MQ=255)
                                                                                                                                                                                             aaactgatatacattgccgaacaatgtgggctttattttcctaagtcgttccctcaggaagTAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCTGACAATGACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAAGAACATGCCCAGAGCTTCGCCAAAGTCCAGCCCGGCCCCCGACAGCAGGGGAt                                                               <  2:96127‑M2/189‑1 (MQ=255)
                                                                                                                                                                                                ctgatatacattgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCgg                                                           >  2:245190‑M2/59‑251 (MQ=255)
                                                                                                                                                                                                                    atgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGcgacga                                       <  2:136342‑M2/213‑1 (MQ=255)
                                                                                                                                                                                                                        ggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAAt                                    <  1:104044‑M2/216‑1 (MQ=255)
                                                                                                                                                                                                                        ggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATc                                   <  2:204091‑M2/217‑1 (MQ=255)
                                                                                                                                                                                                                                 tttcataagtcgttcccgcaggaagGAAATGCGGATGATTTTACCGTCACCGGAACAACAGATCCGCCAAAAAGAGGCGTACCCATGTCGATATAGTCCCCCCCGCGGACAATCAATTCGTCAATGACTGCCACCGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCAAGAGCTTTACCAAAGTCCTGCCCGGCCCCCACCAGCAGGGGATGCGGATTCGGAAAACGAGCGACGAAATCGACCAAcgc                          >  2:98233‑M2/26‑251 (MQ=255)
                                                                                                                                                                                                                                  ttcataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAAcgcg                         >  1:135233‑M2/25‑251 (MQ=255)
                                                                                                                                                                                                                                                   cAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGt          <  2:123144/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTCCCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGTAcc       >  1:243874/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    aggaaggaAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGTAcc       >  1:102190/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         ggaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGTACCgcggc  >  1:228819/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
AGACGAAAACGGCAACCAGATCGGCAAAGGCAAGAAAGGTGACAAGCTGATAACTTTGCTTCCCAACTGGATTGATAAGCTGGGAGACGAATAAACCAGCCTTCAACCCCATCTCATCAATCAACGCCCGGCCCCGCTGCCGGGTTTTTGCTATGCACCACAATTACCCCAACCGGATACACAGCCGGATACAATTCCACCAGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAGAGGCGTACCAATGTCGATATAGTCCCCCGCGCGGACAATCACTTCGTCAATGACTGCCAACGGGAGTTGTTGTAGCTGTGGGCGCAACAACATGCCCAGAGCTTTACCAAAGTCCTGCCCGGCCACCACCAGCAGGGGATGCGGATTCGGAAAACGCGCGACGAAATCGACCAACGCGTTGATGACCGTCAGTACCGCGGC  >  NC_000913/2565239‑2565738

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: