Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,348,8160AC59.1% 17.0 / 25.1 22T245P (ACC→CCC) gltBglutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (9/0);  new base C (0/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 2.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC  >  NZ_CP009273/3348730‑3348884
                                                                                      |                                                                    
aGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCg                                                                   >  1:39135/1‑90 (MQ=255)
 gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCgg                                                                  >  2:497257/1‑90 (MQ=255)
                cacCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgc                                                   >  1:337796/1‑90 (MQ=255)
                  ccGTCCCGCGCTGGCCGCGGGCCCACCCGTCCCCCTTTCGGGCGCAAACGGGTAAAACCACCCCCATCCCCGGTAACCGCCAATGGgcgc                                                 <  2:455078/90‑1 (MQ=255)
                   cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg                                                >  1:569633/1‑90 (MQ=255)
                        cgcgCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgc                                           >  1:455078/1‑90 (MQ=255)
                         gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg                                          >  1:492460/1‑90 (MQ=255)
                          cgcgggccgcgggcGCAACCGTCCCGCTATCGGGCGCATAACGGTAAAATCAACCCCATCCCGGGTAACCGCCAATGGgcgcgtgcgcgt                                         <  1:544713/90‑1 (MQ=255)
                          cgcgggccgcggccGCACCCGTCCCGTTTTCGGGCGCATAAGGGTAAAATCACCCCCACCCCCGGTAACCGCCAATGGgcgcgtgcgcgt                                         <  2:238123/90‑1 (MQ=255)
                            tgggccgcgggcGCACCCTTCCCGCTGTTTGGGGAATAGGGGTAAAATCAACCCCACCCCGGGTAACCGCCAATGGGCGCGTGCGCGTAc                                       <  1:470237/88‑1 (MQ=255)
                            cgggccgcgggcCCACCCTTCCCGCTTTCGGGCGCAAAAGGGTAAAATCACCCCCCTCCCCGGTAACCGCCAATGGGCGCGTGCGCGTAc                                       <  2:163577/88‑1 (MQ=255)
                                 cgcgggcCCACCCGTCCCGCTTTCGGGCGAATAAGGGTAAAATAAACCCCATCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTtata                                  <  2:243058/90‑1 (MQ=255)
                                   ctggcGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATaaa                                >  2:544713/1‑90 (MQ=255)
                                             cGTTCCGCTATGGGGCGCATAGGGGTAAAATAACCACCTCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                      <  2:309244/90‑1 (MQ=255)
                                             cGTCCCGCTATCGGGCGCATAGGGGAAAAACAAACCCCCCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                      <  2:588988/90‑1 (MQ=255)
                                               ttCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCg                    >  1:29705/1‑90 (MQ=255)
                                                  cGCTTTCGGGCGAAAAACGGTAAAAACACCCCCACCCCCGTTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                 <  2:47396/90‑1 (MQ=255)
                                                  cGCTGTCTGGCGCATAAGGGTAAATTCAACCCCCTCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                 <  2:41643/90‑1 (MQ=255)
                                                  cGCTGTCGGCCGCAAAACGGAAAAATCACCCCCCTCCCGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                 <  1:497257/90‑1 (MQ=255)
                                                        cTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCt           >  1:431296/1‑90 (MQ=255)
                                                            cgcATAACGGAAAAACCACCCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc       <  1:269485/90‑1 (MQ=255)
                                                                 aagggTAAAACCAACACCACCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac  <  2:569633/87‑1 (MQ=255)
                                                                                      |                                                                    
AGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCAC  >  NZ_CP009273/3348730‑3348884

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