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

AAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAT  >  NZ_CP009273/2380258‑2380423
                                                                                   |                                                                                   
aaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAt                                                                               <  2:310996/90‑1 (MQ=255)
 aaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATg                                                                              <  1:1045060/90‑1 (MQ=255)
     gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGcaca                                                                          <  1:966710/90‑1 (MQ=255)
            tAGTGCATGGTGTTATAAATTTATTATAAAATAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATCCACattttatt                                                                   <  1:605031/90‑1 (MQ=255)
                  aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACaattttttttttct                                                             >  1:183605/1‑86 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                         <  2:806919/90‑1 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                         <  1:522876/90‑1 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                         <  2:141720/90‑1 (MQ=255)
                      tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                         <  2:114973/90‑1 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGAAAATTTTATTTTTCCTTTTCTGTTTTTAg                                               >  1:512627/1‑90 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGAACATTTTATTTTTTCTCCTCTGtttttaa                                               >  1:1031404/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCtttttttaa                                               >  2:1084658/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGGACATTTTTTTTTTCCTCATCTTTTTTTAg                                               >  2:700585/1‑90 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCATCTGTTTTTAg                                               >  2:835294/1‑90 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAACAAAAGAAAATTTTATTTTTTCTCAGCTGGTTTTAg                                               >  1:979837/1‑90 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCtttttttaa                                               >  2:432432/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTATTTTTCCTCATCTGtttttaa                                               >  2:446295/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACATTTTATTTTTCCTCCtttttttttaa                                               >  2:483802/1‑87 (MQ=255)
                               tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTTCTCCTCTGtttttaa                                               >  1:314224/1‑87 (MQ=255)
                                 tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAGAt                                             >  2:106107/1‑90 (MQ=255)
                                        caacaGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAAACAAGATCAATTTTTTTTTTCTCCTCTGTTTTTAAAATGGGaaa                                       >  1:633039/1‑90 (MQ=255)
                                               aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTCCCCCTCTTTTTTTAAAAGCGGAAAAAGaaaa                               >  2:240837/1‑89 (MQ=255)
                                                  taataaAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCATCTTTTTTTAAATGCGAAAAAAGaaaatta                            >  1:69546/1‑86 (MQ=255)
                                                    ataaAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTTTTTTTCCCCCTCTGTTTTTAGATGCGAAAAATGAAAATCAAc                          >  2:281878/1‑90 (MQ=255)
                                                    ataaAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACATTTTTTTTTTCCTCATCTTTTTTTTGATGGGAAAAAAGAAAATTAAc                          >  1:132677/1‑90 (MQ=255)
                                                                            aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt  <  1:206830/90‑1 (MQ=255)
                                                                                   |                                                                                   
AAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAT  >  NZ_CP009273/2380258‑2380423

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