Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,348,8100AC52.9% 14.4 / 18.4 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 = 7.13e-03
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCG  >  NZ_CP009273/3348721‑3348896
                                                                                         |                                                                                      
tGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAaca                                                                                        >  2:346536/1‑90 (MQ=255)
 gTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac                                                                                       >  2:76369/1‑90 (MQ=255)
 gTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAacac                                                                                       >  1:250643/1‑90 (MQ=255)
                       aaCACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATgg                                                                 >  1:135913/1‑90 (MQ=255)
                          acCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcg                                                              >  1:42575/1‑90 (MQ=255)
                           ccGTCCCGCGCGGCCCGCTGGCGCACCCGTCCCGCTATCGGGCGCATAACGGTAAAACCACCCCCACCCCCGGTAACCGCCAATGGgcgc                                                             <  1:76369/90‑1 (MQ=255)
                             gTACGGCGCCGGCCGCCGGCGCAACCGTTGCGTTATGTGGCGCGTAACGGTGAGATCAACACCATCACCGGTAGCCGCCAATGGgcgcgt                                                           >  2:187302/1‑90 (MQ=255)
                                 cgcgcgGGCCGCGGGCCCAACCGTCCCGCTATCGGGCGCATAAGGGAGAAACAAACACCACCCCCGGTAACCGCCAATGGgcgcgtgcgc                                                       <  1:346536/90‑1 (MQ=255)
                                 ccggctggccgctggccCACCCTTTCCTTTACCTGGCGAAAACCGGCAAAAACACCCCCACCCCCGATAACCGCCAATGGgcgcgtgcgc                                                       <  2:215143/87‑1 (MQ=255)
                                  gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGgcgcgtgcgcg                                                      >  1:215143/1‑90 (MQ=255)
                                                 ccacccGTCCCGCTATCGGGCGCAAAAGGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCa                                       <  2:135913/86‑1 (MQ=255)
                                                      cGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                                  >  1:240866/1‑90 (MQ=255)
                                                      cGTCCCGTTTTCTGGCGCAAAAGGGAAAAATCAACCCCACCACGGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGacac                                  <  1:212088/90‑1 (MQ=255)
                                                           cGCTTTCTGGCGCATAACGGTGAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                             <  1:24359/90‑1 (MQ=255)
                                                           cGCTTTCGGGCGCATAACGGTAAAATCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                             <  1:246973/90‑1 (MQ=255)
                                                           cGCTATCGGGCGCATACGGGTAAAACCACCCCCACCCCCGGAACCCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCtt                             <  2:344977/90‑1 (MQ=255)
                                                                tcgggCGCAAAAGGGTAAAACCAACCCCACCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATccc                        <  2:136476/87‑1 (MQ=255)
                                                                     cgcATAACGGAAAAACCAACCCCCCCCCCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGAcc                   <  2:42575/90‑1 (MQ=255)
                                                                              gTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCacgacg          >  1:206459/1‑90 (MQ=255)
                                                                                      aaCACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCg  >  2:206507/1‑90 (MQ=255)
                                                                                         |                                                                                      
TGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCG  >  NZ_CP009273/3348721‑3348896

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