Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,348,8100AC52.9% 17.4 / 21.6 17T243P (ACC→CCC) gltBglutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base C (0/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCC  >  NZ_CP009273/3348731‑3348890
                                                                               |                                                                                
gcgcTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGTGCATAACGGTGAAATCAACACCATCACCgg                                                                        >  2:172408/1‑90 (MQ=255)
                 ccGTCCCGCGCTGCCCGCGGGCCCACCCTTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGgcgc                                                       <  2:513698/90‑1 (MQ=255)
                  cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac                                                                                 <  2:473374/63‑1 (MQ=255)
                  cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac                                                                                 >  1:473374/1‑63 (MQ=255)
                  cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg                                                      >  2:422906/1‑90 (MQ=255)
                  cGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcg                                                      >  2:337319/1‑90 (MQ=255)
                         cgcgggccgcgggcGCACCCGTTCCGCTATGTGGCTAATAACGGGAAAATCAACCCCACCCCCGGAAACCGCCAATGGgcgcgtgcgcgt                                               <  2:376959/90‑1 (MQ=255)
                                          cccGTCCCGCTATCTGGCGCATAACGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac                              <  2:127472/89‑1 (MQ=255)
                                          cccGTCCCGCTATCGGGCGAATAGGGGAAAAAACACCCCCACCCCCGTTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGac                              <  1:435190/89‑1 (MQ=255)
                                            cttcccGCTATCTGGCGCATAGCGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                            <  2:252491/86‑1 (MQ=255)
                                            cGTCCCTCTATCTGCCGCAAAAGGGTAAAATAAACCCCCCCCCCGGTACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                            <  1:436017/90‑1 (MQ=255)
                                                 cGCTATCTGGCGCATAAGGGTAAAATCAACCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                       <  2:369098/90‑1 (MQ=255)
                                                 cGCTATCGGGCGCATAAGGGTAAAACCACCCCCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                       <  2:399561/90‑1 (MQ=255)
                                                    tATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATc                    >  1:136614/1‑90 (MQ=255)
                                                     aTCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATcc                   >  1:145130/1‑90 (MQ=255)
                                                           cgcATAACGGTAAAACCACCCCCACCCCCGGAAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc             <  2:145130/90‑1 (MQ=255)
                                                                aaCGGTAAAATAAACACCATCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCac        <  2:115202/90‑1 (MQ=255)
                                                                    gTGAAATCAAAACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacgacg    >  1:485003/1‑90 (MQ=255)
                                                                      gAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGAcgcc  >  1:132830/1‑90 (MQ=255)
                                                                               |                                                                                
GCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCC  >  NZ_CP009273/3348731‑3348890

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