Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA57.1% ‑4.9 / 22.2 21intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (0/9);  new base A (12/0);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 3.40e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA  >  NZ_CP009273/2380260‑2380419
                                                                                 |                                                                               
aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                         <  1:280812/90‑1 (MQ=255)
            gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACattttttttt                                                            >  1:37148/1‑90 (MQ=255)
            gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAGGACAAGCACattttatttt                                                            >  1:313083/1‑90 (MQ=255)
                aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGGCAAACACAATTTATTTTTCCt                                                        >  1:848259/1‑90 (MQ=255)
                  ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAGAAGCACAATTTTTTTTTcccca                                                       >  1:564174/1‑88 (MQ=255)
                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGGACAATTTTTTTTTCCTCAt                                                     >  1:505132/1‑90 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAATTAATAAAAATTTCTCGAGGCGGTTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                    <  1:760240/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                    <  1:319467/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                    <  1:851258/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                    <  1:4056/90‑1 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACatt                                                                   <  2:89043/66‑1 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACatt                                                                   >  1:89043/1‑66 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGAAAAAGACAATTTTTTTTTCCTCCTCtttttttaa                                           >  2:797635/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCATCTGtttttta                                           >  2:760240/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAAAAAAGAACAATTTTTTTTTTCTCATCttttttttt                                           >  2:233268/1‑87 (MQ=255)
                                      caacaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGGGaaaa                                  >  1:710258/1‑89 (MQ=255)
                                                  ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                         <  2:335051/55‑1 (MQ=255)
                                                  ataaAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTATTTTccc                                                         >  1:335051/1‑55 (MQ=255)
                                                     aaaaTTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTTTTTTTAAATGGGAAAAAAAAAAATCaaaaga                   >  2:354854/1‑86 (MQ=255)
                                                          ttCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaa              <  2:913955/90‑1 (MQ=255)
                                                                      gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa  <  1:797635/90‑1 (MQ=255)
                                                                                 |                                                                               
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAA  >  NZ_CP009273/2380260‑2380419

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