Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,054,3090AG57.6% 6.8 / 43.5 33E168G (GAA→GGA) yihLGntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base G (11/8);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCA  >  NZ_CP009273/4054182‑4054447
                                                                                                                               |                                                                                                                                          
tATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTAc                                                                                                                                  >  2:761937/1‑138 (MQ=255)
  tCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTAc                                                                                                                                  <  2:67488/136‑1 (MQ=255)
         aCCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAgaga                                                                                                                        <  1:494218/139‑1 (MQ=255)
                        cTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAAt                                                                                                         >  2:354972/1‑139 (MQ=255)
                                          ttCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTg                                                                                                                                      <  1:259492/92‑1 (MQ=255)
                                          ttCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTg                                                                                                                                      >  2:259492/1‑92 (MQ=255)
                                           tCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTg                                                                                                                                      >  2:305579/1‑91 (MQ=255)
                                           tCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTg                                                                                                                                      <  1:305579/91‑1 (MQ=255)
                                                 aCGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGc                                                                                 <  1:107309/138‑1 (MQ=255)
                                                           tatTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGatat                                                                       <  2:591099/138‑1 (MQ=255)
                                                                 tCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAg                                                                >  2:493091/1‑139 (MQ=255)
                                                                 tCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAg                                                                >  2:566957/1‑139 (MQ=255)
                                                                     gAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATc                                                            >  2:763616/1‑139 (MQ=255)
                                                                      aaGTTGTTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCt                                                           >  1:154651/1‑139 (MQ=255)
                                                                            cTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAgaga                                                                                                                        <  2:277838/72‑1 (MQ=255)
                                                                            cTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAgaga                                                                                                                        >  1:277838/1‑72 (MQ=255)
                                                                            cTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCt                                                                                                                                    >  2:666460/1‑60 (MQ=255)
                                                                            cTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCt                                                                                                                                    <  1:666460/60‑1 (MQ=255)
                                                                                      tATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCaaa                                                                                              <  2:667847/88‑1 (MQ=255)
                                                                                      tATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCaaa                                                                                              >  1:667847/1‑88 (MQ=255)
                                                                                       aTTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTgg                                                        <  2:292784/125‑1 (MQ=255)
                                                                                       aTTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTgg                                                        >  1:292784/1‑125 (MQ=255)
                                                                                           aacaacTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAgaga                                                                                                                        <  1:547884/57‑1 (MQ=255)
                                                                                           aacaacTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAgaga                                                                                                                        >  2:547884/1‑57 (MQ=255)
                                                                                               acTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCg                                  <  1:761937/139‑1 (MQ=255)
                                                                                                cTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGa                                 >  1:760591/1‑139 (MQ=255)
                                                                                                cTTGAAAACCACTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGa                                 >  1:15763/1‑139 (MQ=255)
                                                                                                cTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGa                                 >  2:769077/1‑139 (MQ=255)
                                                                                                           cTCATTTTCTGCAGTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGa                       >  1:584600/1‑138 (MQ=255)
                                                                                                                   cTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATcg              >  1:430048/1‑139 (MQ=255)
                                                                                                                        gTCTGGGGAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGcc          >  1:17670/1‑138 (MQ=255)
                                                                                                                             ggAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAg    <  2:430048/139‑1 (MQ=255)
                                                                                                                             ggAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAg    <  1:678947/139‑1 (MQ=255)
                                                                                                                               aaaaGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCa  <  2:179325/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                          
TATCGCATCACCGGCTGGGGGGCGCTGGAAGGACATAAAGTTTTCTATCACGAAACATATATTAATCCTGAAGTTGCTCCGGGCTTTATTGAACAACTTGAAAACCACTCATTTTCTGCAGTCTGGGAAAAGTGCTACCAAAAAGAGACGGTAGTAAAAAAATTGATTTTCAAACCCGTCAGAATGCCGGGCGATATCAGCAAGTATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCA  >  NZ_CP009273/4054182‑4054447

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