Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927371,5610TG57.1% 4.9 / 13.3 14N136T (AAC→ACC) thiBthiamine ABC transporter substrate binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (8/0);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 3.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGA  >  NZ_CP009273/71475‑71644
                                                                                      |                                                                                   
gCCCCGGTGTACTGGGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGttttt                                                                                 <  2:700295/90‑1 (MQ=255)
 ccccGGTGTACTGGTGCGCGGTTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGttttt                                                                                 <  2:98938/90‑1 (MQ=255)
    cGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGgttttttt                                                                              >  1:537559/1‑87 (MQ=255)
          aCTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGtttttt                                                                        >  1:161593/1‑90 (MQ=255)
                gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTttttttttttt                                                                  >  2:101698/1‑87 (MQ=255)
                gcgcGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTGTttttttt                                                                  >  2:138942/1‑87 (MQ=255)
                          tGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGAGGTTTTCAGTTTGTTCTTGTCATAAACGAAg                                                        >  1:587846/1‑90 (MQ=255)
                            aTAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGc                                                      <  1:29428/90‑1 (MQ=255)
                                  cACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTAGTTTTTTTTTTTTATAAAAGAAGGCGGAgga                                                >  1:193547/1‑88 (MQ=255)
                                     ccGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGcc                                             <  2:724004/90‑1 (MQ=255)
                                             tttGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCaa                                     <  2:560294/90‑1 (MQ=255)
                                               tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAc                                   <  1:20949/90‑1 (MQ=255)
                                               tGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTGGTTCCTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAc                                   <  1:482574/90‑1 (MQ=255)
                                                     ctctCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTTGTTTTTTTTTTTTTTAAAAAAAAGGGAAAAAGGCATAATAAAAAGGGACg                             >  2:224664/1‑90 (MQ=255)
                                                                   tttCAGGCTTTGTGGCGGGGTTTTTTTTTTTTTTTTTTTTTTAAAAAAGGCGGAGAAACCATAATTAAAAAGTACGCAAAAGTCgttatt               >  2:614649/1‑90 (MQ=255)
                                                                   tttCAGGCTTTGTGGCGGGGTTTTTAGTTTGTTTTTTTCATAAAAGAAAGCGAAGAAGCCATAATCAAAAGGGACGAAAGAGTCattatt               >  1:84053/1‑90 (MQ=255)
                                                                                ggCGGTTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGa  <  2:235073/90‑1 (MQ=255)
                                                                                      |                                                                                   
GCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGA  >  NZ_CP009273/71475‑71644

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