Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 12,384 T→C G74G (GGT→GGC dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,3840TC97.1% 106.7 / NA 35G74G (GGT→GGCdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/11);  minor base A (1/0);  total (24/11)
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.85e-01

ACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAA  >  NC_000913/12137‑12626
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aCCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCcgc                                                                                                                                                                                                                                                   >  1:130462/1‑251 (MQ=255)
      tatataGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTc                                                                                                                                                                                                                                             >  1:317821/1‑251 (MQ=255)
        tataGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTcc                                                                                                                                                                                                                                           >  1:403294/1‑250 (MQ=255)
            gTGGAGATGTTTAGATGGGTAAAATAATTGGTATCGAACTGGGTACTAACAAATATTGTGTAGCGATTATGGATGGCACCACTCATCGTGTGCTGGAGAACGCCGAAGGCGATCGCAACAGGCCTTCTATCATTGACTATAACCAGGATGGTGAAACTCTGGTTGGTCAGACGGCTAAACGTCAGGCAGTGAAGAAACCGAAAAAAACCCTGTTTGCGATTAAACGACTGATTGGCGGCCGCTTCCAGGAc                                                                                                                                                                                                                                      >  2:211800/1‑251 (MQ=255)
                          aTGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAAGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGACCCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTACAGGACGAAGAAGTACAGCg                                                                                                                                                                                                                        >  1:85784/1‑251 (MQ=255)
                              gTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGAt                                                                                                                                                                                                                    <  1:160091/251‑1 (MQ=255)
                                   ataatTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCGGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTc                                                                                                                                                                                                              <  1:128115/251‑1 (MQ=255)
                                      agtGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCcat                                                                                                                                                                                                           <  2:94134/249‑1 (MQ=255)
                                            aTCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGcc                                                                                                                                                                                                     >  2:202140/1‑251 (MQ=255)
                                                                  ctTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                                                                                                               >  2:331539/1‑251 (MQ=255)
                                                                                        caccacTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGa                                                                                                                                                           >  1:119099/1‑249 (MQ=255)
                                                                                                     gTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATAATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCCTGGGTCGACGTTAAAGGCCAGa                                                                                                                                             >  1:272821/1‑250 (MQ=255)
                                                                                                         tGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGCCGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATg                                                                                                                                        >  2:266998/1‑251 (MQ=255)
                                                                                                              aaCGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCATCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCAGCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAcc                                                                                                                                   <  2:272821/251‑1 (MQ=255)
                                                                                                                            ccGCACCACGCTTTCAATCATTGCCTATACCCAGGAGGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACCCCCTGTTTGGGATAAACCGCCTGATTGGCCGCCGCTCCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTctg                                                                                                                     <  1:76926/250‑1 (MQ=255)
                                                                                                                              gCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACAGCCGCAGATCTctgct                                                                                                                   >  1:168070/1‑251 (MQ=255)
                                                                                                                                                   ccTATACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGACCCCGCAAACCCCCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCGGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATgaag                                                                                              <  1:202140/251‑1 (MQ=255)
                                                                                                                                                     tataCCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGGGCTGAAAAAAATgaaga                                                                                             >  2:30164/1‑250 (MQ=255)
                                                                                                                                                         cccAGGATGTGGAAACTCTGGTGGGTCGGCCGGCTAACCGTCAGGCATTGACGAACCCGCAAAACCCCCTGTTTGCGATAAAACGCCTGATTGGCCGCCGTTCCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGTTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAAcc                                                                                        <  1:278785/251‑1 (MQ=255)
                                                                                                                                                              gATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTg                                                                                    >  1:153041/1‑250 (MQ=255)
                                                                                                                                                                 ggTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGCAGTGCTGAAAAAAATGAAGAAAACCGCTGa                                                                                   >  2:52885/1‑248 (MQ=255)
                                                                                                                                                                  ctcccagtcggggtggccggCCGGCAAAACGGCAGGCAGGGACAACCCCGAAAACCACCCGGTTGGCGATTAAACGCCTGATTGGCCGCCGCTTCAAGGCCGAAAAAGTACAGCGTGATGTTTCCATCATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGAATGGGTCGAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAg                                                                                 <  2:477945/230‑1 (MQ=255)
                                                                                                                                                                     aaaCTCTGGTTGGTCATCCGGCTAAACGTCACGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGAGATGTTTCCATCATGCCGTTCAGAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATTGCCACGCCGCTGATCTCTGCTGAAGAGCTGAAAAAAATGAAGAAACACGCTTAAGATTa                                                                             >  1:376084/1‑250 (MQ=255)
                                                                                                                                                                        ctctGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTg                                                                         >  1:112804/1‑251 (MQ=255)
                                                                                                                                                                        ctctGGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGACGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTg                                                                         >  1:268547/1‑251 (MQ=255)
                                                                                                                                                                            ggTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAg                                                                                 >  1:477945/2‑239 (MQ=255)
                                                                                                                                                                               gggccAGCCGGCTAAACTTCAGGAAGTGACGAACCCGCAAAACCCCCTGTTTCCGATTAAACGCCTGATGGGCCGCCGCTTCCAGGACAAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCAGGGTTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                  <  1:81940/247‑1 (MQ=255)
                                                                                                                                                                                ggTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAcc                                                                 >  1:431102/1‑251 (MQ=255)
                                                                                                                                                                                           tAAACGTCAGGCAGTGACGAACCCGCAAAACACCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGa                                                        >  2:128294/1‑249 (MQ=255)
                                                                                                                                                                                                              aaCCCGCAAAACACACTGTTTGCGATTAAACGCCTGATTGGCCGCCGATTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCt                                                                                                           >  2:442692/1‑179 (MQ=255)
                                                                                                                                                                                                                              gtgtggtgtatAAACGCCTGAGGGGCCGCGCTTCGCAAGACGAAGAGGTACAGGGTGTGTTTTCAAGCATGCCGTTAAAAATTTTTGATGGTGATAACGGCGCCGCAGGGGGCGAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGTTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCt                   <  2:479909/241‑1 (MQ=255)
                                                                                                                                                                                                                                    cGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGCAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGATCAGCGt             >  2:130149/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCgg             >  1:119933/1‑250 (MQ=255)
                                                                                                                                                                                                                                         aaaCCCCTTATTGGCCCCCGTTTCCAGGACGAAAAATTACGGCGTGTTGTTTCCATCGTGCCTTTCAAAATTTGTGCTGCTGATAACGGCGACGCATGGGTCGAATTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAgg         <  1:265860/250‑1 (MQ=255)
                                                                                                                                                                                                                                                tGATTGGCCGCCGCTTCCAGGACGAAGAAGCACAGAGTGATGTTTCCATCATGCCGTTCAAAATTATTGATGCTGATAGCGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGCAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaa  >  2:439904/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAA  >  NC_000913/12137‑12626

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