Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,329,9480TG54.5% ‑4.1 / 24.4 22intergenic (+24/‑140)BW25113_RS25560/acnAhypothetical protein/aconitate hydratase AcnA
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base G (12/0);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 1.41e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTA  >  NZ_CP009273/1329865‑1330027
                                                                                   |                                                                               
ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAt                                                                           >  1:12833/1‑90 (MQ=255)
      ttttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctt                                                                     <  1:238534/90‑1 (MQ=255)
                              gCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTa                                             <  2:515273/90‑1 (MQ=255)
                                 gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTCTCCTCttttttta                                          >  1:665891/1‑86 (MQ=255)
                                 gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTCTCCTCTTTTATCa                                          >  1:452453/1‑90 (MQ=255)
                                                aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctt                                                                     <  2:451450/48‑1 (MQ=255)
                                                aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAtttgctt                                                                     >  1:451450/1‑48 (MQ=255)
                                                aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTTGGGTGTTTTc                           >  2:126445/1‑90 (MQ=255)
                                                aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCCCAAGGTTTTTCCTCTTTTTTCAATTTGGGTTGTTATc                           >  2:572879/1‑90 (MQ=255)
                                                     tCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTCTCCTCTTTTTTCAATTTGGGTTGTTTTCAAATc                      >  2:242891/1‑90 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGtt                   <  1:126445/90‑1 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGtt                   <  1:572879/90‑1 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCCCAAGGTTTCTCCTCTTTTTTTCATTTGGGGTGTTATCAAATCCtt                   >  2:362357/1‑90 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTTAATTTTGGGTGTTATTAAATCGtt                   >  2:546286/1‑90 (MQ=255)
                                                        cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTTTCCTCTTTTTTTAATTTGGGTTGTTATCAAATCGtt                   >  2:358384/1‑90 (MQ=255)
                                                                tGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCACAAGGTTTTTCCTCTTTTTTCAATTTGGGGTGGTTTCAAATCGTTACGCGATg           >  2:478193/1‑90 (MQ=255)
                                                                    cccGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTCTCCTCTTTTTTCAATTTGGGGTGTTTTCAAATCGTTACGCGGTGtttt       >  2:245339/1‑89 (MQ=255)
                                                                      cGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg     >  1:576890/1‑90 (MQ=255)
                                                                      cGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg     >  1:137148/1‑90 (MQ=255)
                                                                        cTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGtt   <  1:362357/90‑1 (MQ=255)
                                                                        cTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGtt   <  1:358384/90‑1 (MQ=255)
                                                                         tGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa  <  1:245339/90‑1 (MQ=255)
                                                                                   |                                                                               
GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTA  >  NZ_CP009273/1329865‑1330027

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