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

GAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGC  >  NZ_CP009273/1548172‑1548339
                                                                                     |                                                                                  
gagaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTc                                                                                >  1:821469/1‑90 (MQ=255)
 agaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTAtttg                                                                               >  2:891835/1‑88 (MQ=255)
  gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg                                                                              >  1:1227232/1‑90 (MQ=255)
  gaTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCgg                                                                              >  1:670164/1‑90 (MQ=255)
   aTCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGc                                                                             <  2:837912/90‑1 (MQ=255)
       gAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCTGCTTGc                                                                         >  1:1447394/1‑90 (MQ=255)
            aGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGCTTTCCAGaa                                                                    >  1:897198/1‑90 (MQ=255)
                 gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAAatat                                                               >  1:241378/1‑90 (MQ=255)
                 gagGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGCGGGATTTGGCTTGCCAGAAAatat                                                               >  1:888202/1‑90 (MQ=255)
                        aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAAt                                                        <  2:1227232/90‑1 (MQ=255)
                        aTAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGCGGGATTCGGCTTGCCCGAAATTATCCTCaaa                                                        >  2:818839/1‑89 (MQ=255)
                           gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCCGAAATTTTCGTCAATTGc                                                     >  1:326569/1‑90 (MQ=255)
                           gCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGGATCGGGGGGTTTGGGTTGGCCGAAATTATCTTCAATggg                                                     >  2:1204672/1‑89 (MQ=255)
                             ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt                                                   <  1:1126416/90‑1 (MQ=255)
                             ccGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGACGCGGTATTCGGGTTGGCCGGAATTATCGGCAATTGCCt                                                   >  1:653731/1‑90 (MQ=255)
                             ccGGTAGCCTTCACTACCGGGCGCAGGCGTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCt                                                   <  2:1295952/90‑1 (MQ=255)
                              cGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGCGGGATTTGGGTTGTCCGAAAATATCTTCAATGGCCTg                                                  >  2:731372/1‑90 (MQ=255)
                                        cACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCt                                        <  2:1013689/90‑1 (MQ=255)
                                               gggCGCAGGCTTAGATGGAGGTACGGCGGGAATCGGGGTATTCGGCTTGCCCGGAAATAAAGTCAAGGGGCTGGTGGAAGGGTTCGGCaa                                 >  1:139814/1‑89 (MQ=255)
                                                 gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag                                 >  1:1547677/1‑88 (MQ=255)
                                                 gcgcAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCag                                 <  2:1547677/88‑1 (MQ=255)
                                                               ggAGGTACGGGGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCAc                 <  1:478250/90‑1 (MQ=255)
                                                                gAGGTACGGCGGTAGTCGCGGGATTTGGCTTTCCAGAAAATATCGTCAATGGGCTGTTGGAGGGGTTCGGCAGAGGGTTTTACCGCCCCg                >  2:1253471/1‑90 (MQ=255)
                                                                gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATTGTCAATGGCCCGTTGCAGGGGTTCGGGAGGGGGTTTTACCCCCCCg                >  2:478250/1‑90 (MQ=255)
                                                                gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTTACCGCCACg                >  1:808700/1‑90 (MQ=255)
                                                                gAGGTACGGCGGTAGTCGCGGGATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCGCCAcc                >  1:578410/1‑89 (MQ=255)
                                                                 aGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGc               <  2:241378/90‑1 (MQ=255)
                                                                         cggTAGTCGCGTTATTCGGCTTGCCAGAAATTATCGTCAATGGCGTGTTGCAGGGCTTCGGCAGGGGTTTTCACCGCCACGCCTtgctgc       <  1:495898/90‑1 (MQ=255)
                                                                          ggTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGGTTTCACCCCCACGCCCtgctgct      >  1:576710/1‑90 (MQ=255)
                                                                          ggTAGTCGCGGGATTCGGCTTGCCAGAAAATATCGTCAATGGCCTGTTGGCGGGCTTCGGCAGAGGGTTTCACCGCCCCGCCTtgctgct      >  2:85371/1‑90 (MQ=255)
                                                                             aGTCGCGGGATTCGGCTTGCCAGAAATTATCGGCAATGGCCCGTTGGCGGGCTTCGGCAGAGGTTTTCACCCCCACGCCCTGCTGCTgcg   >  1:105313/1‑90 (MQ=255)
                                                                              gTCGCGGTATTGGGCTTGCCAAAAATTATCGTCAATGGCCTGTTGCGGGGCTTCGGAAGAGTTTTTCACCGCCACGCCTTGCTGCTgcgc  <  2:1552897/90‑1 (MQ=255)
                                                                                     |                                                                                  
GAGATCTGAAAAAGGGAGAGGGAAATAGCCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGTATTCGGCTTGCCAGAAATTATCGTCAATGGCCTGTTGCAGGGCTTCGGCAGAGGTTTTCACCGCCACGCCTTGCTGCTGCGC  >  NZ_CP009273/1548172‑1548339

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: