Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,548,2570TG53.3% 10.5 / 18.7 16Y557S (TAC→TCC) maeAmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (8/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 6.99e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCA  >  NZ_CP009273/1548174‑1548341
                                                                                   |                                                                                       
gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg                                                                                  >  2:326584/1‑90 (MQ=255)
gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTgg                                                                                  >  1:114804/1‑90 (MQ=255)
gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTCgg                                                                                  >  2:370708/1‑90 (MQ=255)
            ggAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAAt                                                                     <  2:573332/90‑1 (MQ=255)
               gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGCGGGTGCGGCGTTTGGCTTGCCGAAAAATTa                                                                   >  1:329658/1‑90 (MQ=255)
                                              ggCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCaga                                  <  2:329658/90‑1 (MQ=255)
                                               gcgcAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTCGGTTTGCCGAAAATAACTGCAAAGGGCTGGTTGACGGGTTCGGCaggg                                 >  1:577455/1‑88 (MQ=255)
                                                        tAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAATTATCGTCAATGGCCTGTTTCAGGGGTTTGGGGGAGGGTTTTAcc                        >  1:98460/1‑90 (MQ=255)
                                                              gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag                                   >  1:91146/1‑73 (MQ=255)
                                                              gAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag                                   <  2:91146/73‑1 (MQ=255)
                                                              gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCCGTTGCAGGGGCTTCGGCAGAGGGTTTCACCGcccc                   >  1:350143/1‑88 (MQ=255)
                                                               aGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGc                 <  2:239588/90‑1 (MQ=255)
                                                                ggTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa                                                <  1:51925/59‑1 (MQ=255)
                                                                ggTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCa                                                >  2:51925/1‑59 (MQ=255)
                                                                    cggcggTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCACGCCTtgc            >  2:89169/1‑90 (MQ=255)
                                                                              cgcgGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACGCCTTGCTGCTGCGCCa  >  1:529874/1‑90 (MQ=255)
                                                                                   |                                                                                       
GATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGCCA  >  NZ_CP009273/1548174‑1548341

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: