Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,367,4000AC60.0% 7.5 / 12.5 15G361G (GGT→GGGmenEo‑succinylbenzoate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (5/1);  new base C (0/9);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 2.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367313‑2367485
                                                                                       |                                                                                     
aaaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAcc                                                                                     >  2:391090/1‑90 (MQ=255)
aaaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAcc                                                                                     >  2:509325/1‑90 (MQ=255)
               aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTgaa                                                                      <  1:606575/90‑1 (MQ=255)
                                  tGCAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                   <  1:10379/90‑1 (MQ=255)
                                   gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:135149/78‑1 (MQ=255)
                                   gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGAt                                                              >  2:135149/1‑78 (MQ=255)
                                              ggATGTGCAGCAATTACGCGCTCGACTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       <  1:186171/90‑1 (MQ=255)
                                                  gtgCAGCAATTACGCGCTCGCCTTTTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                   <  1:509325/90‑1 (MQ=255)
                                                   tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                  <  2:128177/90‑1 (MQ=255)
                                                          aTTACGCGCTCGATTTCTCCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                           <  1:75267/90‑1 (MQ=255)
                                                                 gcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  1:555129/61‑1 (MQ=255)
                                                                 gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 >  2:555129/1‑61 (MQ=255)
                                                                                 cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt    <  1:555682/90‑1 (MQ=255)
                                                                                   gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  >  1:117298/1‑90 (MQ=255)
                                                                                   aaaTCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  1:310412/89‑1 (MQ=255)
                                                                                       |                                                                                     
AAACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367313‑2367485

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