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

TCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGA  >  NZ_CP009273/2371232‑2371393
                                                                                  |                                                                               
tCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAAttt                                                                          >  2:758871/1‑90 (MQ=255)
       ggcgTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCccgcc                                                                   <  1:563301/90‑1 (MQ=255)
          gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt                                                                >  1:1711978/1‑90 (MQ=255)
            cctaaaaaGGGACGTTTTGCGGCATAAAAAAAAACCCCCCCCGCTCGTTTTGGGGCGTTGGCAACAGGGAAAAAATTTGCCCGCCgttgt                                                              <  2:1479097/86‑1 (MQ=255)
                    ggggcGTTTTGCGGCATAAGATAGAACCGCCCCCGCTCGTTTTGCGGCGTGGGAAACAGCAAAAAATTTTGCCCGCCGTTGTTGTTcacc                                                      <  2:1812564/86‑1 (MQ=255)
                     ggccGTTTGGCGGCATAAGAAAGAACCCCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacca                                                     <  2:37120/87‑1 (MQ=255)
                     ggccGTTTGGCGGAATCAGAAAAAAACCCTCACGCTCGTTTTGCGGCGTTGGAAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTcacca                                                     <  2:293407/87‑1 (MQ=255)
                     ggAGGTTTTGCGGCATCAGATAGAACCGCTCCCCCTCGCTTTGCGGGGTTGGCAAAAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacca                                                     <  1:1349373/90‑1 (MQ=255)
                     ggACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacca                                                     <  1:303280/90‑1 (MQ=255)
                      gACTTTTTGCGGCATAAGATAGAACCCCCCCCCCTCGTTTTGCGGCGTGGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                    <  2:1711978/90‑1 (MQ=255)
                      gACGTTTTGCGGCATCAGATAGAAACGCTCCCGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTCCCCGCCGTTGTTGTTcaccac                                                    <  1:2252080/90‑1 (MQ=255)
                      gACGTTTTGCGGCATCAGATAAAACCGCTCCCGCTCGCTTTGCGGCGTTGGAAAAAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                    <  1:1293333/90‑1 (MQ=255)
                      gACGTTTTGCGGCATAAGATAAAACCCCCCCCCCTCGTTTTGCGGCGTTGGCAACAGCGAAAAATTTCGCCCGCCGTTGTTGTTcaccac                                                    <  2:1817135/90‑1 (MQ=255)
                      gACGTTTTGCGGAATAAGATAAAACCGCTCCCCCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                    <  2:643924/90‑1 (MQ=255)
                                   aTCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgc                                       >  2:353170/1‑90 (MQ=255)
                                       gATAAAACCGCTCCCGCCCCCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgcggcg                                   <  1:1802395/90‑1 (MQ=255)
                                           gAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt                                                                <  1:512/57‑1 (MQ=255)
                                           gAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt                                                                >  2:512/1‑57 (MQ=255)
                                                gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt                          >  1:392195/1‑90 (MQ=255)
                                                gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt                          >  1:1624862/1‑90 (MQ=255)
                                                 ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCTGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg                         >  2:202162/1‑90 (MQ=255)
                                                 ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac                                                    >  2:1307683/1‑63 (MQ=255)
                                                 ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac                                                    <  1:1307683/63‑1 (MQ=255)
                                                 ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg                         >  2:1891072/1‑90 (MQ=255)
                                                   cccgctcgttTTGCGGGGTTGCCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc                       <  2:40124/88‑1 (MQ=255)
                                                    ccgctcgttTGGGGGTGTGGGCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg                      <  2:227422/89‑1 (MQ=255)
                                                    ccgctcgttTGGGGGCGTGGCAAACAGCGAAAAATTTCCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg                      <  2:2263806/89‑1 (MQ=255)
                                                     cgctcgctTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGc                     >  1:659594/1‑90 (MQ=255)
                                                           ctTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAAc               >  1:227422/1‑90 (MQ=255)
                                                               ggggcgTGGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa           <  1:1501119/88‑1 (MQ=255)
                                                               gcggcgTTGGCAACAGCGAGAGAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa           >  1:310976/1‑90 (MQ=255)
                                                                 ggcgTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgc         <  1:322372/90‑1 (MQ=255)
                                                                   cgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgcgc       >  2:320352/1‑90 (MQ=255)
                                                                     ttGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGt     >  1:1147262/1‑90 (MQ=255)
                                                                      ggccAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGtt    <  2:2133427/87‑1 (MQ=255)
                                                                        gCAACAGCGAAAAAATTCCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga  <  2:1549095/90‑1 (MQ=255)
                                                                                  |                                                                               
TCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGA  >  NZ_CP009273/2371232‑2371393

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: