Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,577,360 C→T 8.2% R78H (CGC→CAC)  gntK ← gluconate kinase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,577,3600CT 8.2% 154.3 / 4.6 49R78H (CGC→CAC) gntKgluconate kinase 2
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAG  >  NC_000913/3577222‑3577499
                                                                                                                                          |                                                                                                                                           
cgcgcTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                                                                                            <  1:817421/140‑1 (MQ=255)
cgcgcTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTgcg                                                                                                                                            <  1:1524108/140‑1 (MQ=255)
   gcTTTCAGGCGTCTTTAAATCACATCAAAATCGAATTTCTAATATATTAAATATAGATTAGGATTACCTTCACGCAGCAAGTCGAGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTg                                                                                                                                         <  2:180017/140‑1 (MQ=255)
    cTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGTAAGAGAGATTCGGATTAACTTCACGCAGCAAGTCGCGATACTGTTTTTTCAAAGCAGAACAGACGATCAACGAAACTTTATTAGGGCGCTGc                                                                                                                                        >  2:1913511/1‑140 (MQ=255)
    cTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGc                                                                                                                                        >  2:1834647/1‑140 (MQ=255)
      ttCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCAt                                                                                                                                      <  1:1556083/140‑1 (MQ=255)
         aGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGc                                                                                                                                   <  2:1667620/140‑1 (MQ=255)
               cTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAAcgc                                                                                                                             >  2:1348218/1‑140 (MQ=255)
               cTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAAcgc                                                                                                                             >  2:935262/1‑140 (MQ=255)
                   cAATCTCATCAACATCGCCTTTCAATTAGATTAAATTTTGATTGTTTTTACCTTCTCGCTTCAAGGCGCGATTGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACgcggcg                                                                                                                         <  2:617720/140‑1 (MQ=255)
                         caTCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTc                                                                                                                   >  2:1719645/1‑140 (MQ=255)
                          aTCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCa                                                                                                                  >  2:1806246/1‑140 (MQ=255)
                               aaTCTCCTTTCAAATATATGAAATAGAGATTCGTATTAACTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGcc                                                                                                             <  2:208677/140‑1 (MQ=255)
                                              aGATGAAAGAGAGATTCGGATTACCTTCACGCAGCTAGTCGCGTTAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTg                                                                                              <  2:340102/140‑1 (MQ=255)
                                                  gAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGAt                                                                                          <  1:1408753/140‑1 (MQ=255)
                                                              cGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCa                                                                              >  1:101037/1‑140 (MQ=255)
                                                                         cACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCTATCGTCGTCATTCAGTGGTTCGCCa                                                                   >  2:1124569/1‑140 (MQ=255)
                                                                          aCGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAg                                                                  >  1:149679/1‑140 (MQ=255)
                                                                          aCGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAg                                                                  >  1:235486/1‑140 (MQ=255)
                                                                               tcaAGCCGCGATATTGTTTTTTCAATGCATATCAGACGATCAGCGTCACTTTATTAGTGCGCTGCATAGCTTACGCGGCGTCGTTCAGCGCCTGCAACCACGTTTTGAGATCGTCGTCTTTCTGTGGTTCGCCAGACGcc                                                             <  2:1681533/139‑1 (MQ=255)
                                                                               gcacgacgctATAGTGTTTTTTCAATGCAGAACAGACGATCATCGAAACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTAGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGcc                                                             <  2:1681336/134‑1 (MQ=255)
                                                                                caAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACTCGTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGGTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCa                                                            >  1:1829826/1‑140 (MQ=255)
                                                                                caAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCa                                                            >  1:1829525/1‑140 (MQ=255)
                                                                                caAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCa                                                            >  1:680330/1‑140 (MQ=255)
                                                                                caAGTCGCGATAGTGTTTTTTAAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCa                                                            >  1:1743816/1‑140 (MQ=255)
                                                                                      gcgtttgtgtTTTTTCAATGCAGAACAGACGATCAGCGACTCTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAtttttt                                                      <  2:1667046/134‑1 (MQ=255)
                                                                                      gcgATTGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAtttttt                                                      <  2:1141845/140‑1 (MQ=255)
                                                                                      gcgATTGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAtttttt                                                      <  2:1328452/140‑1 (MQ=255)
                                                                                      gcgATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTTTTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAtttttt                                                      <  2:276367/140‑1 (MQ=255)
                                                                                      gcgATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCAtttttt                                                      <  1:978824/140‑1 (MQ=255)
                                                                                        gATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGTGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCg                                                    >  2:1779509/1‑140 (MQ=255)
                                                                                        gATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGTGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCg                                                    >  2:1778832/1‑140 (MQ=255)
                                                                                               ttttttCAATGCAGAACAGAAGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGGCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTgc                                             <  2:1401831/140‑1 (MQ=255)
                                                                                                  tttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGcc                                          <  1:372628/140‑1 (MQ=255)
                                                                                                  tttCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGcc                                          >  2:314749/1‑140 (MQ=255)
                                                                                                                 gACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAAt                           <  1:1385737/140‑1 (MQ=255)
                                                                                                                       cAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCAt                     >  1:90276/1‑140 (MQ=255)
                                                                                                                          cGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCCTTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCaa                  >  2:1336062/1‑140 (MQ=255)
                                                                                                                            acacTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGa                >  2:382889/1‑140 (MQ=255)
                                                                                                                            acacTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGCTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGa                >  2:876893/1‑140 (MQ=255)
                                                                                                                              agtttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCTTTCAGTGGTTCGCCAGACGGCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGaat              >  2:1827883/3‑139 (MQ=255)
                                                                                                                                tttATTAGTGTGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:1779509/140‑1 (MQ=255)
                                                                                                                                tttATTAGTGTGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:1778832/140‑1 (MQ=255)
                                                                                                                                tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:1764137/140‑1 (MQ=255)
                                                                                                                                tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:266850/140‑1 (MQ=255)
                                                                                                                                tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:759052/140‑1 (MQ=255)
                                                                                                                                tttATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcg            <  1:759506/140‑1 (MQ=255)
                                                                                                                                 ttattaGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAAcgc           >  1:877445/1‑140 (MQ=255)
                                                                                                                                  tattaGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGTAAATCGCCATCAAGAAAcgcg          >  2:1111518/1‑140 (MQ=255)
                                                                                                                                        tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGa    >  1:899178/1‑140 (MQ=255)
                                                                                                                                        tGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGa    >  2:1528087/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  2:139850/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:878665/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:1556367/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:69035/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:1510158/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  2:427689/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:223496/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  2:751302/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  2:775536/1‑140 (MQ=255)
                                                                                                                                          cgcTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAg  >  1:223495/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                           
CGCGCTTTCAGGCGGCTTTCAATCACATCAAAATCGCCTTTCAAATAGATGAAAGAGAGATTCGGATTACCTTCACGCAGCAAGTCGCGATAGTGTTTTTTCAATGCAGAACAGACGATCAGCGACACTTTATTAGTGCGCTGCATAGCAAACGCGGCGTCGTTCAGCGCCTGCAACCACGGTTTGCGATCGTCGTCATTCAGTGGTTCGCCAGACGCCATTTTTTCGATATTGCGCCGTGGATGGAGGAAATCGCCATCAAGAAACGCGGCATGAAG  >  NC_000913/3577222‑3577499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑