Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,350,8530AC45.3% 81.9 / 13.9 54S725A (TCA→GCA) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (26/3);  new base C (0/24);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350705‑3351003
                                                                                                                                                       |                                                                                                                                                      
tttGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GAcc                                                                                                                                                      >  1:189465/1‑151 (MQ=255)
       cccccGTTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACTTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCccgccg                                                                                                                                                <  1:916443/150‑1 (MQ=255)
         cacGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCAc                                                                                                                                             <  2:868613/151‑1 (MQ=255)
              tttCAGCCTTTCTACTTCGTGCCGCCACTCTTCTTTCACCCCTCCAACCAATCCACCGACTTTATCTCCCCAGGCCGGAA‑GGTCAGGAGACTGATTTTGCTGCCCCAGCTTTTGCAGATGGCTTAC‑CCCCACT‑GCCCCCGCCGCACCt                                                                                                                                           <  1:42609/148‑1 (MQ=255)
                    cACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAAtttt                                                                                                                                   >  1:169503/1‑150 (MQ=255)
                     aCTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTAt                                                                                                                                 >  1:50445/1‑151 (MQ=255)
                      cTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATg                                                                                                                                >  2:951689/1‑151 (MQ=255)
                       ttCTCCGTCGTGCCGCCCCTCTTCTTTCATCTCTACAACCCATCCACCGACTTTATCTCCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTAt                                                                                                                                 <  2:422777/149‑1 (MQ=255)
                            cgtcgtGCCGCCCCTCTTTTTTCATCTCTCCAACCCATCCACCGACTTTTTCTCCCCAGCCCGGAA‑GGTCAGGAGACTAATTTTGCTGACCCAGCTGTTGCAGATGGCTTAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGtcctt                                                                                                                           <  2:127801/150‑1 (MQ=255)
                               cgtGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                       >  2:620364/1‑151 (MQ=255)
                                gtGACCCCACTCTTCTTTCATCTCTTCAATCCATCCCCCGACTTTATCTCCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAC‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTc                                                                                                                       <  2:76622/150‑1 (MQ=255)
                                    cccccctctTCTTTCATCCCTTAAACCCATCCCCCGACTTTTTCTCCCCAGGCCGGAA‑GGTCGGGAACCTGAATTTGCTCACCCAGTTTTTGCAGAGGGCTTAC‑CCCCCCT‑CCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                   <  1:90701/145‑1 (MQ=255)
                                    cGCCCCTCTTCTTTCATCTCTCCAACCCATCCCCCGCCGTTATCTCCCCAGCCCGGAA‑GGTCAGGAGACTGAATTTCCTGCCCCAGTTGTTGAAGATGGCTTAC‑CCCCACG‑GCCCCCGCCGCACCTTTAATTTTATGtccttcct                                                                                                                        <  2:632314/145‑1 (MQ=255)
                                      cccctctTTTTTCATCTCTTCAACCCATTCCCCGACGTTATCTCCCCAGCCCGAAA‑GGTCAGGAGACTGATTTTGCTGCCCCAGTTGTTGCAGATGGCGTAC‑CCCCCCT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTcaaca                                                                                                                   <  2:184148/145‑1 (MQ=255)
                                       cACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACTCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGc                                                                                                               >  1:933573/1‑151 (MQ=255)
                                        aCTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGcc                                                                                                              >  2:1028597/1‑151 (MQ=255)
                                            ttcttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttt                                                                                                          >  1:343100/1‑151 (MQ=255)
                                              cttTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                        >  1:589273/1‑151 (MQ=255)
                                               tttCTTCTCTAAAACCCATTCCCCGACTTTATCTCCCCAGGCCGAAA‑GGTCAGGAGCCTGATTTTCCTGCCCCATTTTTTGCAGGTGCCATAC‑CCCCCTT‑CCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                         <  1:34310/149‑1 (MQ=255)
                                               ttccacctctCCAATCCATCCCCCGTCGTTTTCTCCCCAGGCCGGAG‑GGTCGGGAGACTAAATTTGCTCCCCCAGCTTTTGCAGATGGCGTAC‑CCCCCCT‑GGCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCtttttt                                                                                                        <  1:309471/144‑1 (MQ=255)
                                                ttCATCTCTCCAACCCATCCCCCGACGTTATCTTCCCAGGCCGGAG‑GGTCAGGAGACTGATTTTGCTGCCCCATTTGTTGCAGATGGCGTAA‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCttttt                                                                                                         <  1:910047/148‑1 (MQ=255)
                                                   aTCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCt                                                                                                   >  2:435/1‑151 (MQ=255)
                                                     ctctTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                                 >  1:254238/1‑151 (MQ=255)
                                                        tccAACCCATTCCCCGACTTTATCTCCCCAGGCGGGAA‑GGCCGGGAGACTAATTTTGCTGCCCCATTTTTTGCAGATGGCTTAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTgc                                                                                                 <  2:192250/146‑1 (MQ=255)
                                                             cccttcccccGTCGTTATCTCCCCAGCCCGGAG‑GGCCAGGAGACGGATTTTGCGGCCCCAGTTGTTGAAGATGGCGAAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCa                                                                                          <  1:882442/142‑1 (MQ=255)
                                                              ccATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTc                                                                                           >  2:379301/1‑148 (MQ=255)
                                                                aTTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAtt                                                                                      >  1:446783/1‑151 (MQ=255)
                                                                  tCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  1:782046/1‑151 (MQ=255)
                                                                   ccccGACTTTATCTCCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGCCCCAGCTGTTGCAGATGGCGTAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    <  2:203054/148‑1 (MQ=255)
                                                                   ccccGACGTTATCTCCCCAGGCCGGAA‑GGTCGGGAGACTGATTTTGCTGCCCCAGCTGTTGCAGAGGGCGAAA‑CCCCACG‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    <  2:589897/148‑1 (MQ=255)
                                                                   cACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  1:570026/1‑150 (MQ=255)
                                                                   cACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCg                                                                                    >  2:319028/1‑150 (MQ=255)
                                                                    cccGACTTTCTCTCCCCA‑GCCGGAAGGGTCGGGAGACGAATTTTGCGCCCCCGGTTGTTGCAGA‑GGCCTAACCCCCCCT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGAt                                                                                       <  1:982174/145‑1 (MQ=255)
                                                                     ccGACGTTATCTCCCCAGGCGGGAA‑GGTCAGGAGACTGAATTTGCTGCCCCAGTTTTTGCAGATGGCGTAC‑CCCCCCT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACt                                                                                 <  2:199181/151‑1 (MQ=255)
                                                                          gTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGc                                                                            >  2:650930/1‑151 (MQ=255)
                                                                                 tCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTg                                                                      >  2:1013841/1‑150 (MQ=255)
                                                                                  cccAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGAc                                                                    >  2:500394/1‑151 (MQ=255)
                                                                                                     gagaCTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtctt                                                  >  1:178557/1‑151 (MQ=255)
                                                                                                       gaCTGATTTCGCGGCCCCAGCTGTTCCAGAGGGCTTAA‑CCCCCCT‑GCCCCCGCCGCTCCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCAtcttc                                                 <  2:859447/150‑1 (MQ=255)
                                                                                                            aaTTTGCTGCCCCAGCTTTTGCAGATGGCGTAA‑CCCCACT‑CCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAac                                           <  2:237677/151‑1 (MQ=255)
                                                                                                              tttGCTGACCCAGCTGTTGCAGATGGCGTAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAaca                                          <  2:159911/150‑1 (MQ=255)
                                                                                                                 gCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGc                                      >  1:843364/1‑151 (MQ=255)
                                                                                                                   tGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCt                                     >  1:88565/1‑150 (MQ=255)
                                                                                                                         aGCTGTTCCGGGGGGCGAAC‑CCCCACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGt                              <  1:577446/151‑1 (MQ=255)
                                                                                                                              ttGCAGGTGGC‑AAA‑CCCCCCTGGCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCgg                           <  1:929739/149‑1 (MQ=255)
                                                                                                                              ttGCAGATGGCATAC‑CCCCCCT‑CCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGt                          <  1:411447/150‑1 (MQ=255)
                                                                                                                                  aGATGGCTAAC‑CCCCCCT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGa                        <  1:1028576/148‑1 (MQ=255)
                                                                                                                                        cGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCg               >  1:66581/1‑151 (MQ=255)
                                                                                                                                         gTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCgg              >  2:212595/1‑151 (MQ=255)
                                                                                                                                         cttacccccACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCtt                 <  2:696320/144‑1 (MQ=255)
                                                                                                                                              ccccccT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCt          <  1:275122/151‑1 (MQ=255)
                                                                                                                                                ccACT‑GCCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTAc        <  2:797817/151‑1 (MQ=255)
                                                                                                                                                    t‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGt    >  1:569984/1‑151 (MQ=255)
                                                                                                                                                    t‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGt    >  1:167332/1‑151 (MQ=255)
                                                                                                                                                       aCCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTc  >  2:514347/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TTTGCCACCCACGCTTTCAGCACTTCTACGTCGTGACGCCACTCTTCTTTCATCTCTTCAATCCATTCACCGACGTTATCTTCCCAGGCCGGAA‑GGTCAGGAGACTGAATTTGCTGACCCAGCTGTTGCAGATGGCGTAA‑CCCCACT‑GACCCCGCCGCACCTTTAATTTTATGTCCTTCCTCAACAATGCCTTTTTTATCCTGCGCCGTCAGATTCGACTCCAGCACGCTGACATAGCCCGGCATCATCTTCTCAAACACCGCTAACCCGTCGGTGATCAGCTTCGGTCCTACAAGTTC  >  NC_000913/3350705‑3351003

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