Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,332,7910AC43.9% 100.8 / 13.6 57V161G (GTG→GGG) basSsensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (30/2);  new base C (0/25);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 3.53e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAT  >  NC_000913/4332645‑4332930
                                                                                                                                                  |                                                                                                                                           
gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAcccc                                                                                                                                         >  1:792976/1‑151 (MQ=255)
gTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAccc                                                                                                                                          <  2:1064977/150‑1 (MQ=255)
  tGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCAccccc                                                                                                                                        <  2:79753/150‑1 (MQ=255)
    aTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                     >  1:280948/1‑151 (MQ=255)
     tAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                     >  1:1106719/1‑150 (MQ=255)
     tAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGcc                                                                                                                                     >  2:1069477/1‑150 (MQ=255)
          aCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGt                                                                                                                                  >  1:612352/1‑148 (MQ=255)
             ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                            >  1:780499/1‑151 (MQ=255)
             ggAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGtt                                                                                                                            >  1:506912/1‑151 (MQ=255)
             gaaaaaaaaTGACTGTCGGCCACGCGCAAGTGGCAGCAGCGGGGAGCCCCTCTCCATCATCTAACCACGCCGTGCCACTAAGGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGAAAACCCCCCCCCGCCAGTgg                                                                                                                                <  2:868874/142‑1 (MQ=255)
               aagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCg                                                                                                                          >  1:499792/1‑151 (MQ=255)
                 aaaaaTGACTGCCGGGCCCGCGCCGGTGGCGCGAGCGGGGAGACCCTCCCCATCATCTGACCACCCCGTCCCACAACCGGTCCAACAAAATGGTGTGTCTTTTCCGCCAGAGGTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCa                                                                                                                        <  2:888198/150‑1 (MQ=255)
                    aaTGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGtt                                                                                                                     >  1:1181359/1‑151 (MQ=255)
                                    gcgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGc                                                                                                     >  2:953490/1‑151 (MQ=255)
                                      gcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCgg                                                                                                   >  2:513006/1‑151 (MQ=255)
                                       ccAGTTGCAGCAGCGGGGGGCCCTTCTCCACCTTCTGATCAGCCCGTGCCACTAAGGGTTCTACATCAATGTGATGCTTTTCCGCCAGCAGTCCAAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                        <  1:1013205/145‑1 (MQ=255)
                                        cAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                  >  2:1024295/1‑150 (MQ=255)
                                         aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTa                                                                                                 >  1:804318/1‑150 (MQ=255)
                                            tGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                              >  1:379010/1‑150 (MQ=255)
                                                gcagcGGGGAGCCCCTCTCCTTCTTCTGATCAAGCCGTGCCACTAACGGTGCTACACAAATGTGAGGCGTTTCCGCCAGCATTTCCAGATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaacaac                                                                                           <  1:167102/149‑1 (MQ=255)
                                                  ggcggggAGCCCCTCCCCATCTTCTAATCAAGCCGGCCCAAAACGGGTCTTACATCAATGTGAGGGGTTTCCCCCAGCTGTCCCAGAGGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                             <  2:650630/141‑1 (MQ=255)
                                                  agcTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                       >  2:1082295/1‑151 (MQ=255)
                                                   gcggggAGCCGCTCCCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTTCCAGAGGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                             <  1:475153/141‑1 (MQ=255)
                                                      gggggcccccccccATCATCTGACCACCCCGTCCCACTACCGGTCCTACAACATTGTGATGTTTTTTCCCCAGCTGTTCCAGATCCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAAcc                                                                                          <  2:299616/132‑1 (MQ=255)
                                                           aCGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                               >  1:583952/1‑150 (MQ=255)
                                                                ctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTg                                                                         >  1:609184/1‑151 (MQ=255)
                                                                   catcatCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                      >  2:109951/1‑151 (MQ=255)
                                                                               aggCCGGGCCACAAGCGGTGCTACACAATTGTGTTGTTTTTCCGCCAGTAGTTCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAg                                                                 <  2:7854/142‑1 (MQ=255)
                                                                               acgCCGTCCCCCAAACGGTCCTACACAATTGGGTTGCTTTTTCCCAAGCATTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                      <  1:975185/137‑1 (MQ=255)
                                                                                 gCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                        >  1:561967/1‑151 (MQ=255)
                                                                                  cctgcccACTACCGGTGCTACACAAATGGGTGGCTTTTCCGCCGGCAGTCCCAAATACAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                   <  2:417889/134‑1 (MQ=255)
                                                                                  ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                       >  1:815531/1‑151 (MQ=255)
                                                                                  ccGTCCCACAAACGGTGCTACATCAATGTGATGCTTTTTCGCCGGCAGTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                       <  1:205608/151‑1 (MQ=255)
                                                                                       ccACTAAGGGTCCTACATAAAGGTGATGCGTTTTCCCCAGCGTTTCAAGAGGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCa                                                      <  2:801541/147‑1 (MQ=255)
                                                                                         aCTAACGGTGCTACATAAAGGTGATGCTTTTCCGCCAGAATTCCCAGATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCg                                                 <  2:479703/150‑1 (MQ=255)
                                                                                           tAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                              >  2:781446/1‑151 (MQ=255)
                                                                                           aaagggTCCTACATCATTGGGATGCTTTTCCCCCAGCGGTCCAAGACGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                <  2:665482/145‑1 (MQ=255)
                                                                                                gTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                         >  1:458967/1‑151 (MQ=255)
                                                                                                 tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                        >  1:684088/1‑151 (MQ=255)
                                                                                                 tGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                        >  2:67603/1‑151 (MQ=255)
                                                                                                  gCTACAACAAGGTGTTGGTTTGGCCCCAGCTGTTCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGcc                                        <  2:846987/150‑1 (MQ=255)
                                                                                                     aCATAAATGTGATGCTTTTCCGCCAGCAGTCCCAGATCCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaa                                     <  2:583935/150‑1 (MQ=255)
                                                                                                       aTCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                  >  2:221094/1‑151 (MQ=255)
                                                                                                          agtgtgGTGCTTTTTCCCCACCGGTCCAAGAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAac                                   <  2:887021/145‑1 (MQ=255)
                                                                                                          aaTGTGATGCTTTTCCGCCAGCAGTCCCAGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                  <  2:765981/148‑1 (MQ=255)
                                                                                                          aaTGTGATGCTTTTCCGCCAGCAGTCCCAGATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacca                                 <  2:1088046/149‑1 (MQ=255)
                                                                                                              tgATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCtt                           >  1:1006725/1‑151 (MQ=255)
                                                                                                                    gTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt                     >  1:1103007/1‑151 (MQ=255)
                                                                                                                         ccccAGCAGTCCCAGATGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGa                  <  1:94134/147‑1 (MQ=255)
                                                                                                                            cagctttccaaGATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg             <  1:274600/141‑1 (MQ=255)
                                                                                                                                agTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc         >  1:990260/1‑151 (MQ=255)
                                                                                                                                  ttCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg       >  1:347344/1‑151 (MQ=255)
                                                                                                                                    caaGATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg       <  1:849116/147‑1 (MQ=255)
                                                                                                                                    caaGATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg       <  2:1097883/147‑1 (MQ=255)
                                                                                                                                       gATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAt  >  1:1180201/1‑151 (MQ=255)
                                                                                                                                        agacAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa    <  2:780480/145‑1 (MQ=255)
                                                                                                                                        aTGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAt  <  1:40744/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAAT  >  NC_000913/4332645‑4332930

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