Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,429,9840GA50.0% 7.5 / 12.2 18intergenic (‑108/+33)uspF/ompNuniversal stress protein UspF/porin OmpN
Reads supporting (aligned to +/- strand):  ref base G (9/0);  new base A (0/9);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 4.11e-05
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.

TAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGG  >  NZ_CP009273/1429902‑1430072
                                                                                  |                                                                                        
aaaaGTGAAACGCCAGATGACTTTTTTTGTGAAATGACGGAGTTCACTTTTTTATTTTAGAAAAAAGGAGAAAAAAAAAGAAAAAAAAGc                                                                                   <  1:151513/89‑1 (MQ=255)
      gtAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCa                                                                             >  2:503031/1‑90 (MQ=255)
       tAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCAt                                                                            >  1:60461/1‑90 (MQ=255)
            caaGATGACTTTTTTTGTGAAAGGACGGATTTCATTTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAg                                                                       <  1:918346/88‑1 (MQ=255)
              aGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCg                                                                     >  1:10974/1‑90 (MQ=255)
                    cacTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                               >  1:870517/1‑90 (MQ=255)
                      ttttttGTGAAAGGAGGGATTTCACTTTTTTAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCaggag                                                             <  1:88783/89‑1 (MQ=255)
                      ttttttGTAAAATGAGGGATTTCACTTTTTAAATTTAGAAAAAAGGAGAAAAAAAAACAAAAAAAACCCCGCCATAACAGCGGGCaggag                                                             <  1:598057/89‑1 (MQ=255)
                        ttttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAggagga                                                           >  1:110704/1‑90 (MQ=255)
                        ttttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAggagga                                                           >  1:429048/1‑90 (MQ=255)
                                          ttCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCa                                         >  1:169000/1‑90 (MQ=255)
                                             acaTTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAc                                      >  2:858483/1‑90 (MQ=255)
                                                  tttAATTTAGAAAAAGGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                 <  2:110704/90‑1 (MQ=255)
                                                  tttAATTAAGAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                 <  2:405747/90‑1 (MQ=255)
                                                                aagaagAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc                   <  2:745284/90‑1 (MQ=255)
                                                                aaGGAGGAAAAAAAAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc                   <  2:563500/90‑1 (MQ=255)
                                                                 aggagaAAAAATAACAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                  <  2:904364/90‑1 (MQ=255)
                                                                                 aGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTgg  >  2:177748/1‑90 (MQ=255)
                                                                                  |                                                                                        
TAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGG  >  NZ_CP009273/1429902‑1430072

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