Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,431,7690AG23.8% 70.0 / 24.4 42K340R (AAA→AGA) lysC2aspartate kinase
Reads supporting (aligned to +/- strand):  ref base A (10/22);  new base G (5/5);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGATTTCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGATTTTAA  >  CP000730/1431624‑1431910
                                                                                                                                                 |                                                                                                                                             
atGATTTCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCaaaa                                                                                                                                             >  1:376193/1‑148 (MQ=255)
atGATTTCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCaaaa                                                                                                                                             >  2:244301/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCANAATTTCa                                                                                                                                        >  2:287509/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAGAATTTCa                                                                                                                                        >  1:351832/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAGAATTTCa                                                                                                                                        >  1:233642/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAGAATTTCa                                                                                                                                        >  1:237064/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAGAATTTCa                                                                                                                                        >  1:290918/1‑148 (MQ=255)
     ttCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAGAATTTCa                                                                                                                                        >  1:9896/1‑148 (MQ=255)
                          tGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGt                                                                                                                   >  2:90314/1‑148 (MQ=255)
                          tGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGAGt                                                                                                                   >  2:291815/1‑148 (MQ=255)
                           gggCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTg                                                                                                                  <  1:31465/148‑1 (MQ=255)
                           gggCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTg                                                                                                                  <  2:188474/148‑1 (MQ=255)
                                aCAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATg                                                                                                             <  1:324975/148‑1 (MQ=255)
                                aCAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATg                                                                                                             >  1:175257/1‑148 (MQ=255)
                                aCAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATg                                                                                                             >  1:359341/1‑148 (MQ=255)
                                           tCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGc                                                                                                  >  2:60431/1‑148 (MQ=255)
                                            cACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCa                                                                                                 <  2:234583/148‑1 (MQ=255)
                                                           tGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTg                                                                                  >  1:60354/1‑148 (MQ=255)
                                                               tttCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACAt                                                                              >  1:25230/1‑148 (MQ=255)
                                                                  catcaAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTaa                                                                           <  1:304061/148‑1 (MQ=255)
                                                                  catcaAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTaa                                                                           <  2:283425/148‑1 (MQ=255)
                                                                       aaTTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGaa                                                                      >  2:148163/1‑148 (MQ=255)
                                                                          ttCTATGATTCTTTAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGTCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAaat                                                                   <  2:227763/148‑1 (MQ=255)
                                                                          ttCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAaat                                                                   <  2:366905/148‑1 (MQ=255)
                                                                               tGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCGAAAGCATTTTTGACATTAATTGAAAATAatat                                                              <  1:9120/148‑1 (MQ=255)
                                                                               tGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAatat                                                              <  1:174937/148‑1 (MQ=255)
                                                                               tGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAatat                                                              <  1:213117/148‑1 (MQ=255)
                                                                                                 aGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAaacaac                                            >  1:358750/1‑148 (MQ=255)
                                                                                                                   tAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTAt                          <  2:323372/148‑1 (MQ=255)
                                                                                                                   tAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTAt                          <  1:244301/148‑1 (MQ=255)
                                                                                                                   tAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTAt                          <  2:257272/148‑1 (MQ=255)
                                                                                                                     gCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAGATAAGTAttt                        <  2:199797/148‑1 (MQ=255)
                                                                                                                     gCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTAttt                        <  2:89796/148‑1 (MQ=255)
                                                                                                                                    gAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAGCATCAACATCTGAAATAAGTATTTCATACGTCATtgatg         <  2:96496/148‑1 (MQ=255)
                                                                                                                                     accATTACGTCAAACTTTCATTAATTGGTTCACGCCTGCGTGATATGTCAGGTTTGTCATCAAAAGCATTTTTGACATTAATTGAAAATAATATTCCTTTCTTCCTAACAACATCATCTTACATAAGTGGTTAATACGTCATtgacgc        <  2:262664/146‑4 (MQ=255)
                                                                                                                                     aGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATtgatga        <  2:279925/148‑1 (MQ=255)
                                                                                                                                     aGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATtgatga        <  1:60431/148‑1 (MQ=255)
                                                                                                                                      gCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCGAACAACAACATCTGAAATAAGTATTTCATACGTCATTgatgat       >  1:300866/1‑148 (MQ=255)
                                                                                                                                         ttATGTCAGAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATACTATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  2:237064/148‑1 (MQ=255)
                                                                                                                                         ttATGTCAGAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  2:351832/148‑1 (MQ=255)
                                                                                                                                         ttATGTCAGAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  2:290918/148‑1 (MQ=255)
                                                                                                                                         ttATGTCAGAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  2:233642/148‑1 (MQ=255)
                                                                                                                                         ttATGTCAGAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  2:9896/148‑1 (MQ=255)
                                                                                                                                         ttATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGAtttt    <  1:287509/148‑1 (MQ=255)
                                                                                                                                           aTGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGATTTTaa  <  1:148163/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ATGATTTCACAAATCGTCAACTTGGATGGGCTACAACTATCCTTCACGATTAAAGATAGTGATTTTCATCAAATTTCTATGATTCTTGAAACATTAAAGAATCAATATGAAGCATTAGCTTATAAAATCAATGAGCATTATGTCAAAATTTCATTAATTGGCTCAGGCATGCGTGATATGTCAGGTGTGGCATCAAAAGCATTTTTGACATTAATTGAAAATAATATACCTTTCTACCAAACAACAACATCTGAAATAAGTATTTCATACGTCATTGATGATTTTAA  >  CP000730/1431624‑1431910

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