Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273799,3250GT66.7% ‑3.7 / 18.2 18V123L (GTG→TTG) ybhJhydratase
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (12/0);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 7.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA  >  NZ_CP009273/799237‑799412
                                                                                        |                                                                                       
tAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgt                                                                                        <  1:796055/90‑1 (MQ=255)
tAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTtgt                                                                                        <  1:1487995/90‑1 (MQ=255)
      actcTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTGGCTcc                                                                                  >  2:567723/1‑90 (MQ=255)
                 tCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCATAttttgg                                                                       >  1:1042770/1‑86 (MQ=255)
                      ggcACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCTGTcatc                                                                  <  1:1466262/90‑1 (MQ=255)
                        cACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatcca                                                                <  2:995069/90‑1 (MQ=255)
                          cTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTcatccatc                                                              <  1:1546882/90‑1 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTGTTTTTTTGTCTCCCCCTATTTGGGGCATTCACCAAAATATggg                                                    >  2:1528023/1‑88 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTGTTTTTTTGCCCCCGCATATTTGGGTTCTTCCTCAatatatttg                                                    >  1:870288/1‑87 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTTTGCCTCCGCATATTGCGGGCATCCATCAATATATGCg                                                    >  1:122034/1‑90 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTTTCCCCCCCTATTTCGGGCCTCCCTCAatataatcg                                                    >  2:1515532/1‑86 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCATATTTCGGGCATCCCTCAATATATTCg                                                    >  1:1149567/1‑90 (MQ=255)
                                    tgatgaCCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCGCATATTGCGGTCATCCATCAATATATGCg                                                    >  1:196711/1‑90 (MQ=255)
                                             tGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGGGTTTTTTTGTCCCCCCCTATTTGGGGCCTCCCCCAATATATGGGTGGgatggt                                           >  2:722185/1‑88 (MQ=255)
                                                        tATCGGCGGCCCAGCGTTATGGCGGGATTTTTTTGCCTCCCCCTATTTCGGTCATCCATCAATATATTCGGGAGATGATGGGAggggggg                                >  1:1549272/1‑88 (MQ=255)
                                                               ggcCCAGCGTTATGGCGGGATTTTTTTGCCCCCCCATATTGCGGGCCTCCATCAATATATGGGGGAGGTGATGGCAGGGGGCGGGaaaaa                         >  1:1260611/1‑89 (MQ=255)
                                                                               gggaTTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATTCGTGAGGTGATGGCAGGGGGCGGCCAAATGATCCTCGGGAAAGAc         >  2:518328/4‑90 (MQ=255)
                                                                               ggTATTTTTTTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                  >  1:1452178/1‑49 (MQ=255)
                                                                               ggTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTg                                                  <  2:1452178/49‑1 (MQ=255)
                                                                                      ttGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCcaca  <  1:1528023/90‑1 (MQ=255)
                                                                                       tgtgCCTCCACATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCca    <  1:890191/87‑1 (MQ=255)
                                                                                       tgtgCCTCCACATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCca    >  2:890191/1‑87 (MQ=255)
                                                                                        |                                                                                       
TAACTCACTCTGCGCCGTCGGCGGCACTATTAACGGTGATGACCATGTTTTTGGTTTATCGGCGGCCCAGCGTTATGGCGGTATTTTTGTGCCTCCGCATATTGCGGTCATCCATCAATATATGCGTGAGATGATGGCAGGCGGCGGCAAAATGATCCTCGGGTCAGACAGCCACA  >  NZ_CP009273/799237‑799412

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