Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,165,556 T→C T140A (ACG→GCG)  parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,165,5560TC100.0% 191.2 / NA 55T140A (ACG→GCG) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/31);  total (24/31)

CAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATAACGGT  >  NC_000913/3165419‑3165690
                                                                                                                                         |                                                                                                                                      
cAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGccc                                                                                                                                      >  2:671182/1‑140 (MQ=255)
     tATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGc                                                                                                                                 <  1:693658/140‑1 (MQ=255)
      aTGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGcccc                                                                                                                                     <  1:215638/135‑1 (MQ=255)
      aTGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGcccc                                                                                                                                     >  2:215638/1‑135 (MQ=255)
       tGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGccc                                                                                                                               >  2:948859/1‑140 (MQ=255)
         cGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCa                                                                                                                             <  1:76808/140‑1 (MQ=255)
            tGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgct                                                                                                                          >  2:196955/1‑140 (MQ=255)
            tGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgct                                                                                                                          <  2:834988/140‑1 (MQ=255)
             ggAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGcccc                                                                                                                                     <  1:883401/128‑1 (MQ=255)
             ggAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGcccc                                                                                                                                     >  2:883401/1‑128 (MQ=255)
              gAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgctcg                                                                                                                        >  2:376987/1‑140 (MQ=255)
              gAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgctcg                                                                                                                        <  1:376987/140‑1 (MQ=255)
                atatCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgctcgct                                                                                                                      <  2:981563/140‑1 (MQ=255)
                atatCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgctcgct                                                                                                                      <  1:847703/140‑1 (MQ=255)
                 tatCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAgctcgctc                                                                                                                     >  1:291639/1‑140 (MQ=255)
                    cGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAAt                                                                                                                  <  1:975765/140‑1 (MQ=255)
                    cGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAAt                                                                                                                  >  2:161661/1‑140 (MQ=255)
                     ggTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATa                                                                                                                 <  2:955084/140‑1 (MQ=255)
                      gTCGCCATGCCGACGGCAGTACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATag                                                                                                                <  1:555538/140‑1 (MQ=255)
                      gTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATag                                                                                                                <  2:57296/140‑1 (MQ=255)
                        cGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATagca                                                                                                              >  1:807485/1‑140 (MQ=255)
                                   cGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAatat                                                                                                   >  2:135601/1‑140 (MQ=255)
                                         tACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGAc                                                                                             <  2:748249/140‑1 (MQ=255)
                                         tACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGAc                                                                                             <  1:988207/140‑1 (MQ=255)
                                            cGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCg                                                                                               >  1:186838/1‑135 (MQ=255)
                                            cGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCg                                                                                               <  2:186838/135‑1 (MQ=255)
                                              gtggtgCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCg                                                                                        >  2:277133/1‑140 (MQ=255)
                                              gtggtgCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCg                                                                                        >  2:1018495/1‑140 (MQ=255)
                                               tggtgCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCgg                                                                                       >  2:847406/1‑140 (MQ=255)
                                                 gtgCCGTTAAGCAAAATGCTTCGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGa                                                                                     >  2:896531/1‑140 (MQ=255)
                                                  tgCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAt                                                                                    <  1:1018495/140‑1 (MQ=255)
                                                   gCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAtt                                                                                   <  1:135601/140‑1 (MQ=255)
                                                    ccGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAcc                                                                                         <  1:220421/133‑1 (MQ=255)
                                                    ccGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAAcc                                                                                         >  2:220421/1‑133 (MQ=255)
                                                    ccGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAtt                                                                                   <  2:127666/139‑1 (MQ=255)
                                                    ccGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGAtt                                                                                   >  1:127666/1‑139 (MQ=255)
                                                     cGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCg                                                                                 >  1:655454/1‑140 (MQ=255)
                                                     cGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCg                                                                                 <  1:196955/140‑1 (MQ=255)
                                                        tAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGgtg                                                                              <  1:101102/140‑1 (MQ=255)
                                                         aaGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGAc                                                                                             >  2:264475/1‑124 (MQ=255)
                                                         aaGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGAc                                                                                             <  1:264475/124‑1 (MQ=255)
                                                             aaaaTGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACg                                                                         <  1:277133/140‑1 (MQ=255)
                                                                       gCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCcgcg                                                               <  1:1108699/140‑1 (MQ=255)
                                                                           aCGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGa                                                          >  1:109279/1‑140 (MQ=255)
                                                                                     aGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATc                                                 >  2:251584/1‑140 (MQ=255)
                                                                                            tCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgc                                          >  2:467329/1‑140 (MQ=255)
                                                                                             cGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgcg                                         <  1:431606/140‑1 (MQ=255)
                                                                                             cGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgcg                                         >  2:392783/1‑140 (MQ=255)
                                                                                                cTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgcgcgc                                      <  2:576268/140‑1 (MQ=255)
                                                                                                cTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgcgcgc                                      <  1:671182/140‑1 (MQ=255)
                                                                                                cTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGgcgcgc                                      <  1:948859/140‑1 (MQ=255)
                                                                                                             ccGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCt                          <  2:830319/139‑1 (MQ=255)
                                                                                                             ccGTCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCt                          >  1:830319/1‑139 (MQ=255)
                                                                                                                tCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGAcc                      >  2:784364/1‑140 (MQ=255)
                                                                                                                tCGAAGTTTGGCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGAcc                      >  2:149822/1‑140 (MQ=255)
                                                                                                                          gCACCCAGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGc            <  2:807485/140‑1 (MQ=255)
                                                                                                                                aGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATaa      <  1:251584/140‑1 (MQ=255)
                                                                                                                                aGTCAGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATaa      <  1:784364/140‑1 (MQ=255)
                                                                                                                                    aGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATAACGGt  >  2:701678/1‑140 (MQ=255)
                                                                                                                                    aGCCGCCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATAACGGt  <  2:928335/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CAGGTTATGCGGTGGAATATCGGTCGCCATGCCGACGGCAATACCGGTGGTGCCGTTAAGCAAAATGTTTGGCAGACGGGCAGGTAGCATTTTCGGCTCCTGCAAAGTGCCGTCGAAGTTTGGCACCCAGTCAGCCGTCCCCTGCCCCAGCTCGCTCAATAGCAGCTCGGAATATTTCGACAACCGGGATTCGGTGTAACGCATTGCCGCGAACGATTTCGGATCGTCCGGCGCGCCCCAGTTCCCCTGACCATCAACCAGCGGATAACGGT  >  NC_000913/3165419‑3165690

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