Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,044,6190GT20.6% 62.2 / 11.8 34intergenic (‑494/‑88)yihA/yihIcell division GTP‑binding protein/activator of Der GTPase
Reads supporting (aligned to +/- strand):  ref base G (15/12);  new base T (5/2);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAA  >  CP009273/4044475‑4044761
                                                                                                                                                |                                                                                                                                              
ggAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCCGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                                                            <  2:137860/149‑1 (MQ=255)
      cAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                                                                      <  2:431264/149‑1 (MQ=255)
             ccGGAGACGTAATTAAACGGAAATGGGATCAACCCGGATTGTTCCGGCTAAAGGAAAACCAGGGTGTTTTGGCGGCCTGCAAGGATTGTGAGACCCGTTAAGGGTTATGAGGCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                                                                     <  1:570547/143‑1 (MQ=255)
                     gcaaTTAAACTGAAATGGAATCAACACTGATGGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGCTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                       <  2:402353/147‑1 (MQ=255)
                     gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGCTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGAAATCTGTCGACTTTTTTCTttgct                                                                                                                       >  2:180381/1‑149 (MQ=255)
                     gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGAAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGGGACAGAGGAATCTGTCGCCTTTTGTc                                                                                                                             >  1:543738/1‑143 (MQ=255)
                      tAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                                                                      <  2:604788/149‑1 (MQ=255)
                      tAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                                                                      >  2:679378/1‑149 (MQ=255)
                       aaTTAAACGGAAATGGAATCAACACGGATTGTTCGGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                                                           <  1:51812/143‑1 (MQ=255)
                        aTTAAAGCGAAAGGTGATCAACACGGAGTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTCGCGGCCTGCAAGTATTGTAAGACCCGTTAAGGGTTATTAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                                                                          <  1:360653/143‑1 (MQ=255)
                        aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttgc                                                                                                                    >  2:82790/1‑149 (MQ=255)
                         tcaaaCGGAAATGGAATCAACACGGATTGTTCAGGATAAAGGAAACGCAGGGTCTGCTGGCCGCGCGCCAGGGTGCTAAGGCCCGTTAAGGGTGATGAGTCAGGAAAAAAGGCGACAGAGCAATCTGTCGCCTTTTTTCTTtgcttgct                                                                                                                   <  2:378330/147‑1 (MQ=255)
                                     ggAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGGTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                             >  1:269305/1‑143 (MQ=255)
                                                       tCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCCGTTAGATTCCGCCGTCAATCTATACt                                                                                     >  2:7912/1‑149 (MQ=255)
                                                        ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAATTCt                                                                                          >  1:391334/1‑143 (MQ=255)
                                                         cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTCCTTGCTTTCTGTTAGATTCCGCCGCAAATCCATACtca                                                                                   >  2:522333/1‑147 (MQ=255)
                                                                  gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCCAATCTATACTGAAtc                                                                                >  1:456054/1‑142 (MQ=255)
                                                                            gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTATATTCCGCCGCAAAACTATACTGAATAAAACGGCTAAAAGaag                                                                >  2:8950/1‑147 (MQ=255)
                                                                                         cTGCAAGGATTGTAAGGCCCGTTAAGGGTTATGAGTCAGGAAAAACGGCGACAGAGTAATCTGTCGCCTTTTTTCTTGGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAtt                                                         <  1:259709/143‑1 (MQ=255)
                                                                                           gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAACTATACTGAATAAAACGGCAAAAAGACGAACCATTag                                                       >  1:335632/1‑142 (MQ=255)
                                                                                             aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGACTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATGCTGAATGAAACGGCTAAAAGACGAAACACTATGa                                                     >  1:434893/1‑143 (MQ=255)
                                                                                                gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTCTACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatcat                                            >  2:137654/1‑149 (MQ=255)
                                                                                                           ccGTTAAGGGTTATGGGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                 <  2:188572/149‑1 (MQ=255)
                                                                                                            cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTAACCCGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATATATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaa                                      <  1:250939/143‑1 (MQ=255)
                                                                                                              ttAAGGGTTATGAGTCAGGAAAGAGGGCGACTGAGTAATCTGGCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAGTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAAc                                    <  1:240111/143‑1 (MQ=255)
                                                                                                                aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgcag                            <  2:261823/149‑1 (MQ=255)
                                                                                                                       aTGAGTCAGGAAAAAAGGCGACAGAGTAAGCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGt                     <  2:209524/149‑1 (MQ=255)
                                                                                                                          aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTcct                  >  2:425696/1‑147 (MQ=255)
                                                                                                                               agAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa             >  2:665004/2‑149 (MQ=255)
                                                                                                                                gAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTAACGCAGCAAAGGTCAt                  >  1:443764/1‑143 (MQ=255)
                                                                                                                                     aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGcgcg       >  2:548398/1‑149 (MQ=255)
                                                                                                                                        gCGACAGATTACTCTGTCGCCTTTTTTTTTTTTTTTTTTTCTATTAGATTCCGCCGCAAATCTATATTGAATAATACGGGTAAAAGACGAAGCATTATGAAACCATCATCTTCAAACTCGTCCCGCACATGTCATGCtaaatc          >  1:475228/1‑141 (MQ=255)
                                                                                                                                          gACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCCCCCAGCAACGGCCATGCAAAGgcgc        >  1:284511/1‑143 (MQ=255)
                                                                                                                                          gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAACGGTCATGCAAAAgcgc        >  1:603818/1‑143 (MQ=255)
                                                                                                                                                gTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaa  <  1:74474/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAA  >  CP009273/4044475‑4044761

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