Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,719,058 G→A P182L (CCG→CTG)  anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,719,0580GA100.0% 229.9 / NA 70P182L (CCG→CTG) anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/31);  total (39/31)

TTGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCAAT  >  NC_000913/1718928‑1719191
                                                                                                                                  |                                                                                                                                     
ttGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCa                                                                                                                               <  1:366102/139‑1 (MQ=255)
    agcagTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGaa                                                                                                                           >  1:965308/1‑139 (MQ=255)
      cagTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATg                                                                                                                         <  2:965308/139‑1 (MQ=255)
      cagTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATg                                                                                                                         <  2:1886919/139‑1 (MQ=255)
        gTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                       <  1:476540/139‑1 (MQ=255)
        gTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                       <  1:718066/139‑1 (MQ=255)
                 tAACTTTACCCGGCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAg                                                                                                              >  2:291838/1‑139 (MQ=255)
                         cccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCTTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGa                                                                                                      >  2:1728617/1‑139 (MQ=255)
                         cccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGa                                                                                                      >  2:1098367/1‑139 (MQ=255)
                          ccgcccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                     <  1:1098367/139‑1 (MQ=255)
                             cccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgcc                                                                                                  >  2:195665/1‑139 (MQ=255)
                             cccgTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgcc                                                                                                  >  2:149835/1‑139 (MQ=255)
                                        tCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGa                                                                                       >  2:532401/1‑139 (MQ=255)
                                          ggCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCaa                                                                                             <  1:2013221/131‑1 (MQ=255)
                                          ggCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCaa                                                                                             >  2:2013221/1‑131 (MQ=255)
                                                  tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGa                                                                                                                 <  2:1000191/103‑1 (MQ=255)
                                                  tATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGa                                                                                                                 >  1:1000191/1‑103 (MQ=255)
                                                          ggTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                     <  1:1956396/139‑1 (MQ=255)
                                                          ggTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCgg                                                                     >  1:116512/1‑139 (MQ=255)
                                                              tACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTgg                                                                 >  1:1889851/1‑139 (MQ=255)
                                                               aCCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTggg                                                                <  2:116512/139‑1 (MQ=255)
                                                                      tGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCca                                                         >  1:845055/1‑139 (MQ=255)
                                                                        aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGaa                                                                                                                           >  2:904579/1‑71 (MQ=255)
                                                                        aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGaa                                                                                                                           <  1:904579/71‑1 (MQ=255)
                                                                        aCGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCa                                                                                    >  2:1217419/1‑110 (MQ=255)
                                                                        aCGCCAGATCCAGGCATCCATCAGCATGGTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGTCAGATTGGCGATGCCCCCAATATTGAGCa                                                                                    <  1:1217419/110‑1 (MQ=255)
                                                                              gATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCa                                                 >  2:17237/1‑139 (MQ=255)
                                                                                    ggCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcg                                                   >  1:388677/1‑131 (MQ=255)
                                                                                    ggCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAgcg                                                   <  2:388677/131‑1 (MQ=255)
                                                                                    ggCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATgg                                           <  1:17237/139‑1 (MQ=255)
                                                                                    ggCATCCATCAGCATGTTACCAGAACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                     <  1:431612/81‑1 (MQ=255)
                                                                                    ggCATCCATCAGCATGTTACCAGAACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                     >  2:431612/1‑81 (MQ=255)
                                                                                        tccatcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACACATTGGCGAt                                                                                                     <  2:1378136/77‑1 (MQ=255)
                                                                                        tccatcAGCAGGTAACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGAt                                                                                                     >  1:1378136/1‑77 (MQ=255)
                                                                                          catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGa                                          >  2:751919/1‑134 (MQ=255)
                                                                                          catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGa                                          <  1:751919/134‑1 (MQ=255)
                                                                                          catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                                     >  1:2029650/1‑139 (MQ=255)
                                                                                          catcAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCa                                     >  2:1538239/1‑139 (MQ=255)
                                                                                                    ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACg                                       <  1:1715493/127‑1 (MQ=255)
                                                                                                    ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACg                                       >  2:1715493/1‑127 (MQ=255)
                                                                                                    ttACCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCg                           >  1:19854/1‑139 (MQ=255)
                                                                                                      aCCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccgc                                                                                                >  1:1737743/1‑68 (MQ=255)
                                                                                                      aCCAGGACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATgccgc                                                                                                <  2:1737743/68‑1 (MQ=255)
                                                                                                          ggACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGc                                                           >  2:1864448/1‑101 (MQ=255)
                                                                                                          ggACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGc                                                           <  1:1864448/101‑1 (MQ=255)
                                                                                                           gACCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGcc                    >  2:1811217/1‑139 (MQ=255)
                                                                                                            aCCGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCt                   >  1:1994546/1‑139 (MQ=255)
                                                                                                              cGGTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTg                 >  1:1425227/1‑139 (MQ=255)
                                                                                                               ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                       <  2:252149/36‑1 (MQ=255)
                                                                                                               ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAg                                                                                                                       >  1:252149/1‑36 (MQ=255)
                                                                                                               ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTgcg                      <  2:1830636/133‑1 (MQ=255)
                                                                                                               ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTgcg                      >  1:1830636/1‑133 (MQ=255)
                                                                                                               ggTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGt                <  1:195665/139‑1 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                              >  1:1979952/1‑60 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCa                                                                                              <  2:1979952/60‑1 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGa                                                             <  1:1670831/93‑1 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGa                                                             >  2:1670831/1‑93 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAc                                <  2:1813230/122‑1 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTAc                                >  1:1813230/1‑122 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtc               <  2:845055/139‑1 (MQ=255)
                                                                                                                gTATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtc               >  2:1832956/1‑139 (MQ=255)
                                                                                                                 tATCGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcc              >  1:578902/1‑139 (MQ=255)
                                                                                                                    cGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcctcc           <  2:1994546/139‑1 (MQ=255)
                                                                                                                    cGTAGCCCCCAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGtcctcc           <  1:1832956/139‑1 (MQ=255)
                                                                                                                         cccccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAgg      <  1:149835/139‑1 (MQ=255)
                                                                                                                           cccAACCAGCTGCCCAGGAATTAGCAGTGACAGATTGGCGATgccgc                                                                                                >  1:1063313/1‑47 (MQ=255)
                                                                                                                           cccAACCAGCTGCCCAGGAATTAGCAGTGACAGATTGGCGATgccgc                                                                                                <  2:1063313/47‑1 (MQ=255)
                                                                                                                           cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCa    >  1:286557/1‑139 (MQ=255)
                                                                                                                           cccAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCa    >  1:1069761/1‑139 (MQ=255)
                                                                                                                             cAACCAGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGAGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCAat  >  2:1302727/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                     
TTGCAGCAGTGGGAGAATAACTTTACCCGCCCGTGCCCACTCGGCATCTTTATCGTAAGGTTTACCGGCCTGACGCCAGATCCAGGCATCCATCAGCATGTTACCAGGACCGGTATCGTAGCCCCCAACCGGCTGCCCAGGAATGAGCAGTGACAGATTGGCGATGCCGCCAATATTGAGCACCATTCGTCGCTCGGTTGGGTGAGCCAGCAGCGCATGATGGAACGCAGGTACCAGCGGTGCGCCTTGTCCTCCCAAGGCAAT  >  NC_000913/1718928‑1719191

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