Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA56.4% ‑1.8 / 49.9 39intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (1/16);  new base A (22/0);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 4.51e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-08
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380260‑2380429
                                                                                 |                                                                                        
aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                  <  2:466565/90‑1 (MQ=255)
aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                  <  1:1138791/90‑1 (MQ=255)
aTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGc                                                                                  <  2:170188/90‑1 (MQ=255)
                aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCTGAAAAAAAAACAAGCACaaattttttttctt                                                                  >  2:1651353/1‑86 (MQ=255)
                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGACAAACACAATTTTTTTTTTCTCAt                                                               >  2:1584864/1‑90 (MQ=255)
                   gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTcct                                                               >  1:670656/1‑88 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                              <  1:1060488/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                              <  2:37078/90‑1 (MQ=255)
                    tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                              <  2:1507131/90‑1 (MQ=255)
                      ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGGAAAAAAGAAATGCACATTTTTTTTTTTCTCATCtt                                                            >  2:750746/1‑89 (MQ=255)
                      ttATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACCCTTTTTTTTTTCTCctctt                                                            >  1:778196/1‑89 (MQ=255)
                            cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTaa                                                      <  1:1584864/90‑1 (MQ=255)
                            cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACACAATTTTTTTTTTCTCATCtttttttt                                                      >  2:352976/1‑88 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAg                                                     <  1:1650453/90‑1 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGCAAGAACCTTTTTTTTTTCCTCTTCTCttttttg                                                     >  1:1313346/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAAAACAATTTTTTTTTCCTCCTCTGtttttaa                                                     >  2:246077/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAAAAAAAAAGAAAAAAACATTTTTTTTTTTCTCCTCTTTTTTTAg                                                     >  2:689045/1‑90 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAACAAATTTTTTTTTTTCTCTTCtttttttaa                                                     >  2:593597/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCTTCtttttttta                                                     >  2:1060488/1‑87 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCCTCTGTTTTTAg                                                     >  1:1191226/1‑90 (MQ=255)
                             tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGAACAATTTTTTTTTTCTCATCtttttttaa                                                     >  2:1542684/1‑87 (MQ=255)
                              ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAAACAAGCACAATTTTTTTTTCCTCCTCtttttttaaa                                                    >  1:662329/1‑86 (MQ=255)
                              ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGAAAAGAACAATTTTTTTTTTTTTATCTGttttttga                                                    >  2:1168071/1‑86 (MQ=255)
                                   ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                                               <  1:750746/90‑1 (MQ=255)
                                   ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCTCATTTTATTTTCCCTCATCTGTTTTAAGATGCGa                                               <  1:569437/90‑1 (MQ=255)
                                     acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCTACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGaaa                                             <  1:273856/90‑1 (MQ=255)
                                        acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGTTTTTAAATGCGaaaaaa                                          >  1:1448437/1‑87 (MQ=255)
                                        acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTTTTTTTTTAAATGCGaaaaaa                                          >  1:799542/1‑87 (MQ=255)
                                        acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAACACATTTTTTTTTTTCTCCTCCTTTTTTAAATGGGaaaaaa                                          >  2:622352/1‑87 (MQ=255)
                                          aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa                                        <  1:540878/90‑1 (MQ=255)
                                          aGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGa                                        <  1:1020455/90‑1 (MQ=255)
                                             aaaTAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAAt                                     <  1:743566/90‑1 (MQ=255)
                                             aaaTAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAAAGAAATGAACAGTTTATTTTTCCTCCTCTTTTTTTAAATGCGAAACATGaaaa                                     >  1:366562/1‑89 (MQ=255)
                                             aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGTACAATTTTTTTTTTCTCCTCTTTTTTTAAATGCGGAAAATaaaaa                                     >  1:1718502/1‑89 (MQ=255)
                                             aaaTAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGCACAATTTTTTTTTCCTCCTCTTTTTTTAAAAGGGGAAAATGaaaa                                     >  1:1102002/1‑89 (MQ=255)
                                                  ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCCTTTGTTTTTAGATGGGAAAAAAGAAAATTAAc                                >  1:40163/1‑90 (MQ=255)
                                                   taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg                               <  2:1448437/90‑1 (MQ=255)
                                                                       tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt           <  1:1370657/90‑1 (MQ=255)
                                                                                aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa  >  1:30264/1‑90 (MQ=255)
                                                                                 |                                                                                        
ATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380260‑2380429

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