Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,781,497 C→T 6.7% S284S (AGC→AGT sotB → sugar efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,781,4970CT 6.7% 192.9 / 3.1 60S284S (AGC→AGTsotBsugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base T (2/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGTTCACGGCCTACAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTGAGCATGTTCTGGGGCGTGGCAA  >  NC_002947/2781364‑2781624
                                                                                                                                     |                                                                                                                               
aGTTCACGGCCTACAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGa                                                                                                                               >  1:203552/1‑136 (MQ=255)
    cACGGCCTACAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGc                                                                                                                           >  1:161852/1‑136 (MQ=255)
     aCGGCCTACAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTGTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCt                                                                                                                          >  1:76201/1‑136 (MQ=255)
             cAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCAt                                                                                                                  >  1:295677/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGTGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:534275/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGTGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:536742/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  2:137160/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  2:254515/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:480751/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:270831/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  2:280450/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  2:401841/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:210012/1‑136 (MQ=255)
                      cGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTGTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTg                                                                                                         >  1:394268/1‑136 (MQ=255)
                           cGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAgg                                                                                                    >  1:537319/1‑136 (MQ=255)
                           cGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAgg                                                                                                    >  2:337147/1‑136 (MQ=255)
                           cGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAgg                                                                                                    >  1:425516/1‑136 (MQ=255)
                           cGTTAGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCACGCATGGCTTCCTGGTAgg                                                                                                    >  1:243269/1‑136 (MQ=255)
                               cGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCg                                                                                                <  2:425516/136‑1 (MQ=255)
                                gCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGa                                                                                               >  1:532531/1‑136 (MQ=255)
                                gCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCAGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGa                                                                                               >  1:359886/1‑136 (MQ=255)
                                     gCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCggg                                                                                          >  2:282224/1‑136 (MQ=255)
                                      cACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCgggg                                                                                         <  1:401416/136‑1 (MQ=255)
                                       aCGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGt                                                                                        <  1:280450/136‑1 (MQ=255)
                                        cGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTa                                                                                       >  1:241558/1‑136 (MQ=255)
                                              gcAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTcgcc                                                                                 >  1:267840/1‑136 (MQ=255)
                                              gcAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTAAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTcccc                                                                                 >  1:547954/1‑136 (MQ=255)
                                                  aTCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGt                                                                             <  2:295677/136‑1 (MQ=255)
                                                     ggcggcGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGcc                                                                          >  2:269512/1‑136 (MQ=255)
                                                          cgACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTg                                                                     >  1:492329/1‑136 (MQ=255)
                                                                           tgctgcTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTc                                                    <  1:337147/136‑1 (MQ=255)
                                                                            gctggggttgttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCg                                                   <  1:222985/130‑1 (MQ=255)
                                                                             gtgctgTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCgg                                                  <  1:71217/135‑1 (MQ=255)
                                                                             ctgctgTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCgg                                                  >  1:506482/1‑136 (MQ=255)
                                                                               gctgTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTTCCGCTGTCGGGc                                                >  1:549848/1‑136 (MQ=255)
                                                                                  gttgttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGCTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAAc                                             >  1:225768/1‑136 (MQ=255)
                                                                                   ttgttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACt                                            >  2:357685/1‑136 (MQ=255)
                                                                                    ttttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACtt                                           <  1:156446/134‑1 (MQ=255)
                                                                                    tgttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACtt                                           <  2:210012/136‑1 (MQ=255)
                                                                                    tgttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACtt                                           <  2:270831/136‑1 (MQ=255)
                                                                                     gttCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTc                                          >  2:505998/1‑136 (MQ=255)
                                                                                       tCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTa                                        <  1:505998/136‑1 (MQ=255)
                                                                                       tCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTa                                        <  2:53503/136‑1 (MQ=255)
                                                                                         ggcggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGTGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTAtg                                      <  2:534275/136‑1 (MQ=255)
                                                                                         ggcggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGTGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTAtg                                      <  2:536742/136‑1 (MQ=255)
                                                                                         ggcggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTAtg                                      <  2:549848/136‑1 (MQ=255)
                                                                                           cggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCgg                                                  >  2:409609/1‑122 (MQ=255)
                                                                                           cggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCgg                                                  <  1:409609/122‑1 (MQ=255)
                                                                                            ggTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATgtgt                                   <  1:232471/136‑1 (MQ=255)
                                                                                                  ggggTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAg                             <  2:359886/136‑1 (MQ=255)
                                                                                                  ggggTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAg                             <  2:532531/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGTTGTTCAGCCGTTTCAGCGAACGCTTCCCGCATGGCTTGCTGGTAGGTTCGATCGGGGTTCTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTg                        <  1:226794/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTg                        <  2:1173/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTg                        <  2:537319/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTg                        <  1:44487/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGCGTTTGCAGTGCTg                        <  2:480751/136‑1 (MQ=255)
                                                                                                       ttttGGTTCATTGGTGTTCAGACGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTTGGTTCGTTCGGGGTACTCGCCGCGTGCCTGTTTCTGATGTTGCCGCTGTCTGGCAACTTCTATGTGTTTGCAGTGCTg                        <  2:24078/136‑1 (MQ=255)
                                                                                                                          aGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTGAGCATGTTCTGGGGCGTgg     >  1:236583/1‑136 (MQ=255)
                                                                                                                          aGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGGTCGATCGGTGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGTGCAACTTCTATGTGTTTGCAGTGCTGAGCATGTTCTGGGGCGTgg     <  2:206409/136‑1 (MQ=255)
                                                                                                                             cGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTGAGCATGTTCTGGGGCGTGGCaa  >  2:44741/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                               
AGTTCACGGCCTACAGCTATATCGAGCCGTTCGCCCTGCACGTAGCGCAGATCGGCGGCGACCGCACCACGCTGTTGCTGCTGTTGTTCGGCGGTGCCGGGGTTTTTGGTTCATTGTTGTTCAGCCGTTACAGCGAACGCTTCCCGCATGGCTTCCTGGTAGGTTCGATCGGGGTACTCGCCGCGTGCCTGTTGCTGATGTTGCCGCTGTCGGGCAACTTCTATGTGTTTGCAGTGCTGAGCATGTTCTGGGGCGTGGCAA  >  NC_002947/2781364‑2781624

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