Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,371,3140GA44.6% 40.1 / 32.3 56F475F (TTC→TTTmenD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase
Reads supporting (aligned to +/- strand):  ref base G (25/6);  new base A (0/25);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.57e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC  >  NZ_CP009273/2371227‑2371397
                                                                                       |                                                                                   
gAACATCTCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGaa                                                                                   >  2:963644/1‑90 (MQ=255)
gAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGaa                                                                                   >  1:339109/1‑90 (MQ=255)
     tCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAAttt                                                                              >  2:171905/1‑90 (MQ=255)
         ggcggcgTGCTAAAAAGGGACTTTTGGCGGCATCAGATAGAACCGCTCACGCTCGTTTTGCGGGGTTGGCAACAGCGAGAAATTTTGccc                                                                          <  2:1618058/90‑1 (MQ=255)
         ggcggcgTGCAAAAAATGGAGTTTTTGGGGAATAAGAGGAAAACGCGCACGCTGGCTTTGGGGTGTTGGGAACAGAGAGAAAATTTGccc                                                                          <  1:1459632/90‑1 (MQ=255)
           cggcgTGCAAAAAATGGACGTTTTGCGGCATAAGATAGAAACCCCCCCGCTCGTTTTGGGGCGTGGGCAACAGCGAGAAAATTTGCccgc                                                                        <  2:928363/90‑1 (MQ=255)
            ggcgTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCccgcc                                                                       <  1:3315459/90‑1 (MQ=255)
              cgTGCTAAAAAGGGACGTTTGGCGGAATAAAAAAGAACCCCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTCCCCGCCgt                                                                     <  2:2034316/90‑1 (MQ=255)
               gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt                                                                    >  2:1522464/1‑90 (MQ=255)
               gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt                                                                    >  2:3315873/1‑90 (MQ=255)
                      aaaTGGTCTTTTTGCGGCATAAGAAAAAACCCCTCGCGCTCGTTTTGGGGCGTTGGCAACGGGGAAAAAATTTCCCCGCCGTTGTTGTTc                                                             <  2:1451211/90‑1 (MQ=255)
                      aaaTGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTc                                                             <  1:1257223/90‑1 (MQ=255)
                      aaaTGGACGTTTTGCGGCATCAGAAAGAACCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTc                                                             <  1:2256374/90‑1 (MQ=255)
                         tGGACTTTTTGGGGCATCAGATAAAACCGCCCCCCCTCGTTTGGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacc                                                          <  1:1522208/90‑1 (MQ=255)
                         tGGACGTTTTGGGGCATCATAAAAAAACGCTCACGCTCGCCTTGCGGAGTTGGCAACAGCGTGAAAATTTACCCGCCGTTGCCGTTgacc                                                          >  2:1616087/1‑90 (MQ=255)
                         tGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacc                                                          >  1:3402443/1‑90 (MQ=255)
                         tGGACGTTTGGCGACATAAGAAAAAACCCCCCCCGCTCGTTTTGCGGCGTTGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacc                                                          <  2:2838274/90‑1 (MQ=255)
                         gggACGTTTGGGGGCATCAGATAAAACCCCCCACGCTCGTTTTGGGGCGTGGGCAACAGAGAAAAATTTTCCCCGCCGTTGTTGTTcacc                                                          <  2:2384435/89‑1 (MQ=255)
                          gggggTTTTGCGGAATAAGATAAAACCCCCCCCGCTCGTTTTGGGGCGTTGGAAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcacca                                                         <  2:3281235/86‑1 (MQ=255)
                           gccGTTTGGCGGCATCAGATAAAACCGCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  2:1890583/88‑1 (MQ=255)
                           gatgTTTTGCGGCATCAAAAAGAACCGCTCCCCCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  2:496422/87‑1 (MQ=255)
                           gACTTTTGGCGGCACCAGATAAAAACCCTCCCGCTCGCTTGGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac                                                        <  1:2841724/90‑1 (MQ=255)
                           gACGTTTTGCGGCATCAGAAAAAAACGCTCCCGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  1:2966977/90‑1 (MQ=255)
                           gACGTTTGGCGGCATCAGATAGAACCGCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac                                                        <  2:356107/90‑1 (MQ=255)
                            aCGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACcaca                                                       >  2:3001615/1‑90 (MQ=255)
                            aCGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACcaca                                                       >  1:2235392/1‑90 (MQ=255)
                                tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa                                                   >  2:497167/1‑90 (MQ=255)
                                ttggCGGCATAAGAAAAAACCCCCCCCGCTCGTTTGGGGGGGTTGGCAACAGGGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTa                                                   <  2:3402443/87‑1 (MQ=255)
                                        aTCAAATAAAAACCCTCCCCCTCGTTTTGCGGAGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgc                                           <  1:781198/90‑1 (MQ=255)
                                         tCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgcg                                          >  1:1547110/1‑90 (MQ=255)
                                                  aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa                                 >  1:2417546/1‑90 (MQ=255)
                                                    cgctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAAcc                               >  1:3241490/1‑90 (MQ=255)
                                                     gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacca                                                         <  1:2848577/63‑1 (MQ=255)
                                                     gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacca                                                         >  2:2848577/1‑63 (MQ=255)
                                                     gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt                              >  1:1707018/1‑90 (MQ=255)
                                                     gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt                              >  1:2781257/1‑90 (MQ=255)
                                                      ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg                             >  1:1445136/1‑90 (MQ=255)
                                                      ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg                             >  1:911523/1‑90 (MQ=255)
                                                       tcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGa                            >  2:85193/1‑90 (MQ=255)
                                                        cccgctcgttTTGGGGCGTGGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc                           <  1:3076586/88‑1 (MQ=255)
                                                        cccgctcgttTGGCGGCGTTGAAAACAGAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc                           <  2:1707018/88‑1 (MQ=255)
                                                        cacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc                           >  2:2777608/1‑90 (MQ=255)
                                                         ccgctcgttTTGCGGGGTGGGAAACAGGAAAAAAATTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg                          <  2:598749/89‑1 (MQ=255)
                                                          cgctcgctTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGc                         >  1:1700322/1‑90 (MQ=255)
                                                                  ttGCGGCTTGGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGc                 <  1:2599828/90‑1 (MQ=255)
                                                                    gcggcggcgGCAACAGAGAAAATCTTGCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa               <  2:212462/90‑1 (MQ=255)
                                                                    gcggcgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa               >  2:2525019/1‑90 (MQ=255)
                                                                      ggcgTTGCCAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgc             <  2:2514168/90‑1 (MQ=255)
                                                                        cgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgcgc           >  1:3345877/1‑90 (MQ=255)
                                                                          ttGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGt         >  2:1337218/1‑90 (MQ=255)
                                                                          ttGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGt         >  2:2711263/1‑90 (MQ=255)
                                                                             gaaaCGGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  2:784734/88‑1 (MQ=255)
                                                                             gCAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  2:48544/90‑1 (MQ=255)
                                                                             gCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  1:2711263/90‑1 (MQ=255)
                                                                             gCAACACAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  2:2933896/90‑1 (MQ=255)
                                                                                 aaGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATc  <  2:339109/89‑1 (MQ=255)
                                                                                       |                                                                                   
GAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC  >  NZ_CP009273/2371227‑2371397

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: