Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,234,3440TG60.0% 2.9 / 20.2 20S348A (TCT→GCT) sstTserine/threonine transporter SstT
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base G (12/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.94e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.14e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGC  >  NZ_CP009273/3234264‑3234433
                                                                                |                                                                                          
cGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgc                                                                                  >  2:61826/1‑90 (MQ=255)
  gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgctg                                                                                >  1:163784/1‑90 (MQ=255)
  gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgctg                                                                                >  1:93370/1‑90 (MQ=255)
  gTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTctgctgctg                                                                                >  1:63820/1‑90 (MQ=255)
    cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGAt                                                                              <  2:9818/90‑1 (MQ=255)
    cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGAt                                                                              >  1:330778/1‑90 (MQ=255)
    cGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCCctgctgctgtg                                                                              >  1:59049/1‑88 (MQ=255)
         tGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                         <  2:208675/90‑1 (MQ=255)
         tGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCAc                                                                         <  2:330778/90‑1 (MQ=255)
            ccACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTgg                                                                      <  1:143502/90‑1 (MQ=255)
              aCGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGcc                                                                    <  2:187105/90‑1 (MQ=255)
                    ctgctgTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCACTGGCCTGTAat                                                              >  1:21151/1‑90 (MQ=255)
                                gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCCGGCCTGTAATATGTTTggggtt                                                  >  2:119923/1‑86 (MQ=255)
                                gtggtgGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCCGGCCCGTAATATGTTTGGGAtt                                                  >  2:157760/1‑90 (MQ=255)
                                      gCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAAc                                            <  2:93370/90‑1 (MQ=255)
                                            ctGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTTGCGGTTGTTCCCCCTGGCCTGAAAAATTTTGGGATTTTCGAACGaatt                                       >  1:94909/1‑87 (MQ=255)
                                              gtgtgCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGc                                    <  2:251979/90‑1 (MQ=255)
                                                   ccTGTGGCGCATCCGGCGTGGCGGGGGGGGCTCTTCTGCTGATCCCCCTGGCCCGTAATATGTTTGGGGTTTCGGACGGTATCGCCACGc                               >  2:245984/1‑90 (MQ=255)
                                                             aTCCGGCGTGGCGGGGGGGGCTCTGCTGCTGATCCCCCCGGCCCGTAATATGTTTGGGGTTTTGAACGATATCGCCATGCAGGGGGGTGc                     >  2:344060/1‑90 (MQ=255)
                                                                          ggggggTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTatc        <  1:327034/90‑1 (MQ=255)
                                                                                tctcTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGCTGCCGTCGGCTTTATCATCggt  >  1:191467/1‑89 (MQ=255)
                                                                                |                                                                                          
CGGTCGATCTGCCCACGGCGCTGCTGTTGAGCGTGGTGGCTTCTCTGTGTGCCTGTGGCGCATCCGGCGTGGCGGGGGGGTCTCTGCTGCTGATCCCACTGGCCTGTAATATGTTCGGTATTTCGAACGATATCGCCATGCAGGTGGTTGCCGTCGGCTTTATCATCGGC  >  NZ_CP009273/3234264‑3234433

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