Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,380,3410GA52.5% 12.8 / 54.1 40intergenic (+31/‑72)yfbL/yfbMM28 family metallopeptidase/YfbM family protein
Reads supporting (aligned to +/- strand):  ref base G (2/17);  new base A (21/0);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.93e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-07
Rejected as consensus: Frequency below/above cutoff threshold.

CAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380254‑2380429
                                                                                        |                                                                                          
cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCg                                                                                          <  1:393433/90‑1 (MQ=255)
                   tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTc                                                                     <  1:40144/90‑1 (MQ=255)
                     gATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                   <  2:302487/90‑1 (MQ=255)
                      aTGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCACAAAAAACGGCAAGCACCATTTTTTTTTcct                                                                   >  2:1258824/1‑89 (MQ=255)
                        ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGACCGGAAAAAGGACATGGGCAATTTTTTTTTCCTct                                                                >  2:1067721/1‑89 (MQ=255)
                        ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACGACAAGCACAATTTTTTTTTccccc                                                                 >  2:1439630/1‑88 (MQ=255)
                        ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAAAAAAAAAGAAAAGCACATTTTTTTTTTCCTct                                                                >  2:1351414/1‑89 (MQ=255)
                          tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                              <  1:302493/90‑1 (MQ=255)
                          tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                              <  1:714435/90‑1 (MQ=255)
                              ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAAGAAAAGCACAATTTTTTTTTTCTCCTTTGtt                                                          >  1:999890/1‑90 (MQ=255)
                              ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGGACATTTTTTTTTTCCTCTTCtttt                                                          >  2:79543/1‑90 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAAAAAAAAAAAAAAGGACAATTTATTTTTTCTCTTCtttttttaa                                                     >  1:1207315/1‑87 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCATCTGtttttaa                                                     >  2:639148/1‑87 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGGACATTTTTTTTTTCCTCATCtttttttaa                                                     >  2:292079/1‑87 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCCTCTTTTTTTAg                                                     >  2:793245/1‑90 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATGCACATTTTTTTTTTTCCCATCtttttttta                                                     >  2:917832/1‑87 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGAAAAAAGGACAAGCACAATTTTTTTTTTCCCATCTTTTTTTAg                                                     >  2:576545/1‑90 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACATGCACAATTTATTTTTTCTCATCTGTTTTTAg                                                     >  2:525254/1‑90 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGGACATTTTTTTTTTCCTCCtttttttttat                                                     >  2:268197/1‑87 (MQ=255)
                                   tttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGAAAAGGAAATTTTTTTTTTTCTCCTCtttttttaa                                                     >  1:498431/1‑87 (MQ=255)
                                    ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAGCACAATTTTTTTTTCCTCCTCTTTTTTTAGa                                                    >  2:450818/1‑90 (MQ=255)
                                        tataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTTTTTTTCCTCCTCTTTTTTTAGATGCg                                                >  2:126103/1‑90 (MQ=255)
                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCGGAAAAAAAGACAAGCACATTTTTTTTTTTCTCCTCTGTTTTTAGATTCGgaaaaa                                          >  2:639077/1‑87 (MQ=255)
                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt                                          <  1:126103/90‑1 (MQ=255)
                                              acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACAt                                          <  1:918360/90‑1 (MQ=255)
                                                 gTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGACAAGGACATTTTTTTTTTCCTCCTCTGTTTTTAGATGCGAAAAAAGaa                                       >  2:186728/1‑90 (MQ=255)
                                                   aaaTAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGAACATTTTTTTTTTTCTCCTCTTTTTTTAGATGCGAAAAAAGaaaa                                     >  1:363678/1‑89 (MQ=255)
                                                        ataaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                <  2:919578/90‑1 (MQ=255)
                                                        ataaAAATTTCTCGAGGCGGGGCCAGAAAAAACGACATGCACATTTTTTTTTTCCTCATCTTTTTTTAGATGCGAAAAAAAAAAATCaaa                                >  2:440046/1‑89 (MQ=255)
                                                        ataaAAATTTCTCGAGGCGGGACCAGAAAAAAAGACAAGCACCTTTTTTTTTTCCTCTTCTGTTTTAAGAAGAGAAAAATGAAATTCaaa                                >  2:668367/1‑89 (MQ=255)
                                                         taaAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACg                               <  1:1333998/90‑1 (MQ=255)
                                                           aaaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACgga                             <  1:1183552/90‑1 (MQ=255)
                                                            aaaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggag                            <  1:1212665/90‑1 (MQ=255)
                                                             aaTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACggagg                           <  1:435377/90‑1 (MQ=255)
                                                                            ggTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGa             <  1:19988/90‑1 (MQ=255)
                                                                             gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa            <  2:498431/90‑1 (MQ=255)
                                                                             gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGaa            <  1:79543/90‑1 (MQ=255)
                                                                              tCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAAt           <  2:971274/90‑1 (MQ=255)
                                                                                     aaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGc    <  2:1391883/90‑1 (MQ=255)
                                                                                       aGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa  >  1:1218844/1‑90 (MQ=255)
                                                                                        |                                                                                          
CAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTA  >  NZ_CP009273/2380254‑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: