Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,734,2690TG55.2% 23.2 / 37.2 29V402G (GTG→GGG) avtAvaline‑‑pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base T (3/10);  new base G (16/0);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 1.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-05
Rejected as consensus: Frequency below/above cutoff threshold.

AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT  >  NZ_CP009273/3734186‑3734358
                                                                                   |                                                                                         
aGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGAt                                                                                     >  2:158004/1‑90 (MQ=255)
      ggATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGc                                                                               >  2:42633/1‑90 (MQ=255)
           aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag                                                                          <  1:17189/90‑1 (MQ=255)
           aaCCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag                                                                          <  1:197001/90‑1 (MQ=255)
                tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGGTTTCGGCGGGAGAGGTa                                                                     >  1:506239/1‑90 (MQ=255)
                tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGTGGGCGAGATa                                                                     >  1:445032/1‑90 (MQ=255)
                tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATa                                                                     >  2:493756/1‑90 (MQ=255)
                tGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATa                                                                     >  2:439080/1‑90 (MQ=255)
                 ggCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGGGGAAGAGATaa                                                                    >  1:230220/1‑89 (MQ=255)
                      cATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag                                                                          <  1:134932/79‑1 (MQ=255)
                      cATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAag                                                                          >  2:134932/1‑79 (MQ=255)
                                    tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGTCTGGGCTGAAAg                                                 <  2:173341/90‑1 (MQ=255)
                                    tATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAg                                                 >  2:115223/1‑90 (MQ=255)
                                         gcATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACt                                            >  2:173151/1‑90 (MQ=255)
                                          cATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTa                                           <  2:349260/90‑1 (MQ=255)
                                           aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTTTGGCGGGAGAGATAGAAAAAGCCTGGGGTGGAAGGCACTaa                                          >  1:115305/1‑90 (MQ=255)
                                           aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGAAAGAAAAAGGCCGGGGTGGAAGGCaacaa                                          >  1:244270/1‑86 (MQ=255)
                                           aTGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGGTGAAAGGCACTaa                                          >  1:178655/1‑90 (MQ=255)
                                               aCTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAAt                                      <  2:245375/90‑1 (MQ=255)
                                                  aCGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGGAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGc                                   >  2:191484/1‑90 (MQ=255)
                                                    gTACCAGAGCCGGAGAAAATTGAGGCGGGGGGGGAGATTCTGGCGGAAGAGATAGAAAGAGCCTTGGCCGAAAGTCACTAAAAAATTCCt                                 >  2:469875/1‑90 (MQ=255)
                                                       ccAGAGCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGGAGAGATAGAAAAAGCCCGGGCTGAAAGGCACTAAAAAATTCCTGAt                              >  1:317397/1‑90 (MQ=255)
                                                           agCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATgcgc                          <  2:244270/90‑1 (MQ=255)
                                                            gCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCGGAAAGGCACTAAAAATTGCCTGATGCGCt                         >  1:298783/1‑90 (MQ=255)
                                                            gCCGGAGAAAATTGAGGCGGGGGGGAAGATTCTGGCGGAAGAGATAGAAAAAGCCTGGGCTGAAAGGCACTAAAAAATGCCTGATGCGCt                         >  1:172047/1‑90 (MQ=255)
                                                                          aGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGcc           <  2:179325/90‑1 (MQ=255)
                                                                              gggggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTTGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTAc                       >  2:6364/1‑74 (MQ=255)
                                                                              gggggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTTGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTAc                       <  1:6364/74‑1 (MQ=255)
                                                                              gggggTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa       <  2:481563/90‑1 (MQ=255)
                                                                                   tGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAtt  <  2:230220/90‑1 (MQ=255)
                                                                                   tGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGAtt  <  2:172298/90‑1 (MQ=255)
                                                                                   |                                                                                         
AGGGCTGGATAAACCGTGGCCGCATACGCATCAATGTATGCGCATGAACTACGTACCAGAGCCGGAGAAAATTGAGGCGGGGGTGAAGATTCTGGCGGAAGAGATAGAAAGAGCCTGGGCTGAAAGTCACTAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATGATT  >  NZ_CP009273/3734186‑3734358

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