Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,630,168 A→G 8.9% A375A (GCA→GCG)  clpX → ATP‑dependent Clp protease ATP‑binding subunit ClpX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,630,1680AG 8.9% 152.8 / 5.6 56A375A (GCA→GCG) clpXATP‑dependent Clp protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (3/2);  total (29/27)
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.69e-01

ACTGGCTACCCTTGACGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAGGGCACCTCG  >  NC_002947/2630021‑2630312
                                                                                                                                                   |                                                                                                                                                
aCTGGCTACCCTTGACGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCAc                                                                                                                                                 >  4:80869/1‑149 (MQ=255)
  tGGCTACCCATGAAGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATCCGCTGACCAAGCAGTATGCCAAGCTGTACGAGATGGAAAGCCTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACg                                                                                                                                                >  2:24288/1‑148 (MQ=255)
         ccTTGACGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGccc                                                                                                                                        >  6:265314/1‑149 (MQ=255)
          cTTGACGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTAAGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCt                                                                                                                                       <  8:25895/149‑1 (MQ=255)
              aCGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCCAGGGCCTTGCGTGCAACGGCCTTGAg                                                                                                                                   <  5:182767/149‑1 (MQ=255)
              aCGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGa                                                                                                                                    <  3:80869/148‑1 (MQ=255)
                 aGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAgcg                                                                                                                                 >  3:281849/1‑148 (MQ=255)
                      gATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGa                                                                                                                            <  6:331617/148‑1 (MQ=255)
                      gATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGAc                                                                                                                           >  1:317362/1‑149 (MQ=255)
                      gATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGAc                                                                                                                           >  2:62357/1‑149 (MQ=255)
                            gCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGgcgc                                                                                                                     >  5:46426/1‑149 (MQ=255)
                              tGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTc                                                                                                                   >  1:119251/1‑149 (MQ=255)
                              tGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTc                                                                                                                   >  5:109455/1‑149 (MQ=255)
                                    gATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGc                                                                                                              <  5:265314/148‑1 (MQ=255)
                                     aTGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTAAGCCAAGCTGTTCGAGATGGACAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCt                                                                                                            <  7:247836/149‑1 (MQ=255)
                                       gCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCCAGGGCCTTGCGTGCAACGGCCtt                                                                                                                                      <  5:92623/121‑2 (MQ=38)
                                       gCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCCAGGGCCTTGCGTGCAACGGCCtt                                                                                                                                      >  6:92623/1‑120 (MQ=38)
                                        cAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGc                                                                                                          <  2:119251/148‑1 (MQ=255)
                                               tCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGa                                                                                                    <  5:88695/147‑1 (MQ=255)
                                               tCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGa                                                                                                    >  6:88695/1‑147 (MQ=255)
                                               tCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATc                                                                                                  >  6:50621/1‑149 (MQ=255)
                                                 aCCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATc                                                                                                  <  4:76396/147‑1 (MQ=255)
                                                    gAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGa                                                                                             <  6:46426/149‑1 (MQ=255)
                                                     aGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGaa                                                                                            <  5:50621/149‑1 (MQ=255)
                                                      gCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAg                                                                                           <  6:109455/149‑1 (MQ=255)
                                                        cGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGc                                                                                         <  5:268642/149‑1 (MQ=255)
                                                          aagaaTGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCgt                                                                                       <  8:135726/149‑1 (MQ=255)
                                                             aaTGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGtt                                                                                    >  8:303133/1‑149 (MQ=255)
                                                               tGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGc                                                                                  <  1:62357/149‑1 (MQ=255)
                                                                  gcTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTc                                                                                >  8:249834/1‑148 (MQ=255)
                                                                     gACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATa                                                                            >  2:22560/1‑149 (MQ=255)
                                                                             aGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCCTGTAt                                                                    <  7:303133/149‑1 (MQ=255)
                                                                               tATGCCAAGTTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTAt                                                                    >  3:19989/1‑147 (MQ=255)
                                                                               tATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACAATGTAt                                                                    >  5:41985/1‑147 (MQ=255)
                                                                                        cTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTc                                                         <  4:281849/149‑1 (MQ=255)
                                                                                              gagaTGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGaaa                                                     >  7:291728/1‑147 (MQ=255)
                                                                                               acaTGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTTCACGCAAGGCCCTGGATCGCAAGACTGGCGCTCGTTGCATGCGTTCGATCATCGAAGGCGTGTTGCTCGATACCATGTATGATATCCCTTCTaaaaa                                                   <  1:24288/146‑1 (MQ=255)
                                                                                                     aaaGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGt                                             >  6:149696/1‑148 (MQ=255)
                                                                                                     aaaGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGt                                             >  2:122725/1‑148 (MQ=255)
                                                                                                      aaGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGTCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCa                                           <  5:238323/149‑1 (MQ=255)
                                                                                                      aaGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCa                                           <  5:331148/149‑1 (MQ=255)
                                                                                                          gTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTc                                                                                >  5:304237/1‑108 (MQ=255)
                                                                                                          gTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTc                                                                                <  6:304237/108‑1 (MQ=255)
                                                                                                          gTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGc                                         >  3:322310/1‑147 (MQ=255)
                                                                                                              aCCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAggtg                                   <  7:249834/149‑1 (MQ=255)
                                                                                                               ccTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAggtgg                                  <  1:22560/149‑1 (MQ=255)
                                                                                                                            aGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGagag                     >  1:292516/1‑149 (MQ=255)
                                                                                                                             gCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGc                    >  8:169480/1‑149 (MQ=255)
                                                                                                                             gCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGc                    <  4:322310/149‑1 (MQ=255)
                                                                                                                                ttGCGCTCCAGGGCCTTGCGTGCAACGGCCTTGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGt                  >  8:40274/2‑148 (MQ=38)
                                                                                                                                    gcTCCAGGGCCTTGCGTGCAACGGCCTTGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCAtt              >  2:159124/1‑148 (MQ=38)
                                                                                                                                      tCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAg           >  6:187703/1‑149 (MQ=255)
                                                                                                                                       cAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAg           <  1:122725/148‑1 (MQ=255)
                                                                                                                                        aaGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAggg         >  2:235562/1‑149 (MQ=255)
                                                                                                                                            ccGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAGGGCAcc     <  1:235562/149‑1 (MQ=255)
                                                                                                                                                tGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTTTCGACGAGAGCGTCATTGAGGGCACCTCg  >  8:51388/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACTGGCTACCCTTGACGAGCTGGATGAGGCTGCGTTGATGCAGATCCTCACCGAGCCGAAGAATGCGCTGACCAAGCAGTATGCCAAGCTGTTCGAGATGGAAAGCGTCGACCTCGAGTTCCGCAGCGATGCGCTCAAGGCCGTTGCACGCAAGGCCCTGGAGCGCAAGACTGGCGCTCGTGGCCTGCGTTCGATCCTCGAAGGCGTGTTGCTCGATACCATGTATGAGATCCCTTCGAAAAAGGATGTCAGCAAGGTGGTTATCGACGAGAGCGTCATTGAGGGCACCTCG  >  NC_002947/2630021‑2630312

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