Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,371,3140GA56.0% 4.4 / 26.8 25F475F (TTC→TTTmenD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base A (0/14);  total (9/16)
Fisher's exact test for biased strand distribution p-value = 2.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC  >  NZ_CP009273/2371249‑2371397
                                                                 |                                                                                   
aaaGGGACGTTTTGCGGCATCAGATAGAACCCCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTc                                                             <  2:1487084/90‑1 (MQ=255)
    ggACTTTTGGCGGCATAAAAAAGAACCGCCCCCGCTCCTTTTGCGGGGTTGGCAACAGCGAAAAAATTTCCCCGCCGTTGTTGTTcacca                                                         <  1:541184/90‑1 (MQ=255)
    ggACGTTTTGCGGCATCAGATAGAAACGCTCCCGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacca                                                         <  2:1107715/90‑1 (MQ=255)
    ggACGTTTTGCGGCATCAGAAAAAACCCCCCCCCCTCGTTTTGCGGCGTTGGCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcacca                                                         <  1:1328272/90‑1 (MQ=255)
     ggcGTTTTGCGGCATCAAATAAAACCGCCCCCGCTCGCTTTGCGGGGTTGGAAAAAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac                                                        <  1:705123/88‑1 (MQ=255)
     gatttttcGCGGCATAAAATAAAACCGCCCCCGCTCGTTTGGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac                                                        <  1:1369936/86‑1 (MQ=255)
     gACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        >  2:946498/1‑90 (MQ=255)
     gACGTTTTGCGGCATCAGATAAAACCCCCCCCCCTCGTTTGGGGGCGTTGGCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  1:93013/90‑1 (MQ=255)
     gACGTTTTGCGGCATCAGATAAAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  2:1050204/90‑1 (MQ=255)
       ctttttGCGGCATCAAATAAAACCGCTCCCGCTCGCTTTGCGGGGTTGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACaa                                                      <  2:1168887/88‑1 (MQ=255)
          tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa                                                   >  1:945914/1‑90 (MQ=255)
          tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa                                                   >  1:1217561/1‑90 (MQ=255)
                            aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa                                 >  1:429511/1‑90 (MQ=255)
                            aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa                                 >  2:786757/1‑90 (MQ=255)
                             acgctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAAc                                >  1:1483991/1‑90 (MQ=255)
                                ctcacgctcGCTTTGCGGCGTTGGCAACAGCgaga                                                                                    >  1:679362/1‑35 (MQ=255)
                                ctcacgctcGCTTTGCGGCGTTGGCAACAGCgaga                                                                                    <  2:679362/35‑1 (MQ=255)
                                ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        >  1:459613/1‑63 (MQ=255)
                                ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac                                                        <  2:459613/63‑1 (MQ=255)
                                  cccgctcgttTTGGGGCGTGGGCAACGGCGAAAAAATTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc                           <  2:642212/88‑1 (MQ=255)
                                   ccgctcgttTGGCGGCGTGGCAAACCGAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg                          <  1:473677/89‑1 (MQ=255)
                                            ttGGGGTGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGc                 <  2:944842/90‑1 (MQ=255)
                                              gcggcgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa               >  1:1584442/1‑90 (MQ=255)
                                              gcggcgTGGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa               <  2:945914/90‑1 (MQ=255)
                                                       gCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  2:968580/90‑1 (MQ=255)
                                                       ccAACAGAGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga      <  2:803283/89‑1 (MQ=255)
                                                           cAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATc  >  2:478086/1‑90 (MQ=255)
                                                                 |                                                                                   
AAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC  >  NZ_CP009273/2371249‑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: