Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,775,2230TG60.0% 13.0 / 21.2 25T234P (ACC→CCC) cysEserine O‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/6);  new base G (15/0);  total (19/6)
Fisher's exact test for biased strand distribution p-value = 1.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTCGCGCCC  >  NZ_CP009273/3775136‑3775311
                                                                                       |                                                                                        
tACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGAt                                                                                        >  2:536926/1‑90 (MQ=255)
   cGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATgcg                                                                                     >  1:218643/1‑90 (MQ=255)
    gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATgcgg                                                                                    >  2:195876/1‑90 (MQ=255)
    gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATgcgg                                                                                    >  2:107121/1‑90 (MQ=255)
    gTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTATCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATgcgg                                                                                    >  1:287075/1‑90 (MQ=255)
         aaTGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGGgggggcc                                                                               >  1:728761/1‑89 (MQ=255)
            gCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCg                                                                            <  1:172332/90‑1 (MQ=255)
              tGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGt                                                                          >  1:667466/1‑90 (MQ=255)
                    atatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTgcagc                                                                    >  2:49333/1‑90 (MQ=255)
                    atatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCCCCGGGTgccgc                                                                    >  1:621330/1‑90 (MQ=255)
                     tatCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTgcagca                                                                   <  2:667466/90‑1 (MQ=255)
                      atCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGGGGCCCCCGGTGCAGcac                                                                  >  1:134527/1‑90 (MQ=255)
                         cATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGcaccac                                                               <  1:243878/90‑1 (MQ=255)
                         cATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGGACCGGGTTCAGcagcac                                                               >  2:853411/1‑90 (MQ=255)
                          aTTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGCGGCACCCGGTGCAGCACCACg                                                              >  1:53099/1‑90 (MQ=255)
                               tGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCCGGTGCAGCACCACGGAAcc                                                         >  2:189909/1‑90 (MQ=255)
                                    tATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGGTGGGGCGGCCCCGGGTGCAGCACCACGGAACCCgcgc                                                    >  1:633170/1‑90 (MQ=255)
                                          tGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCCGGTGCAGCACCCCGGAACCCGCCCCCATCt                                              >  1:19729/1‑90 (MQ=255)
                                           gTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCGGTTGCAGCACCACGGAACCCGCGCCAATCtt                                             >  2:463315/1‑90 (MQ=255)
                                           gTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGGATGCGGCGGCACCGGGTGCAGCACCACGGAACCTGCGCCCATCtt                                             >  2:417164/1‑90 (MQ=255)
                                                gTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGGGGGATGCGGCGGCACCGGGTGCAGCACCCCGGGAACCGCGCCCATCTTcgcgc                                        >  1:829651/1‑90 (MQ=255)
                                                        aCAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCa                                <  1:623112/90‑1 (MQ=255)
                                                                     gAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGtt                                                                         <  1:662375/36‑1 (MQ=255)
                                                                     gAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGtt                                                                         >  2:662375/1‑36 (MQ=255)
                                                                         gCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAg               <  2:880496/90‑1 (MQ=255)
                                                                                  ggtgggATGCGGCGGCACCGGGTGCAGCACCACCGAACCCGCGCCCATCTTCGCGCCGCGCCCCACTTCAATATTGCCGAGGGTTTTcgc      >  1:856423/1‑90 (MQ=255)
                                                                                  ggtgggATGCGGCGGCACCCGGTGCAGCACCACGGGACCTGCGCCCATCTTCTCGCCCCGCCCCACTTCAATATTTCCGAGGATTTTtgc      >  1:880496/1‑90 (MQ=255)
                                                                                     ggtATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTcgcgcc   <  1:437415/90‑1 (MQ=255)
                                                                                      gtATGCGGCGGGACCGGGTTCAAAACCACGGAACCTGCGCCCATCTTCGCGCCGCGCCCCACTTCAAAATTGCCCAGGATTTTcgcgccc  >  1:46637/1‑90 (MQ=255)
                                                                                       |                                                                                        
TACCGTTGAAATGCTGGTCCATATCCATTGATGGCTTATCGCTGTCTGGTTTACCGACAATACGAGCCGGAACGCCAGCGGCGGTGGTATGCGGCGGCACCGGTTGCAGCACCACGGAACCTGCGCCAATCTTCGCGCCGCGCCCAACTTCAATATTGCCGAGGATTTTCGCGCCC  >  NZ_CP009273/3775136‑3775311

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