Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9550TG65.4% 5.1 / 21.1 26T495P (ACC→CCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (2/7);  new base G (17/0);  total (19/7)
Fisher's exact test for biased strand distribution p-value = 5.47e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822873‑823037
                                                                                  |                                                                                  
tCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggg                                                                             >  1:467566/1‑90 (MQ=255)
     cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAAc                                                                        <  1:577732/90‑1 (MQ=255)
             gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTGGACGCCGGGAg                                                                >  2:944018/1‑90 (MQ=255)
             gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg                                                                >  1:1105524/1‑90 (MQ=255)
             gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg                                                                >  1:995971/1‑90 (MQ=255)
             gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAg                                                                >  1:13543/1‑90 (MQ=255)
              ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCCGAAgc                                                               >  1:1110047/1‑89 (MQ=255)
                    gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGAAGGCGGGtc                                                         >  2:814400/1‑90 (MQ=255)
                    gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCgggtt                                                         >  1:1224806/1‑87 (MQ=255)
                     aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAGTCGGtttg                                                        >  2:721950/1‑88 (MQ=255)
                     aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCCGAAGTCGGGtcg                                                        >  2:69016/1‑90 (MQ=255)
                         tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGTCGGGtcgtcgg                                                    >  2:162776/1‑89 (MQ=255)
                                      tttCTCTACCATGCTGTTGACGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTTGTCGGGAAACAGAATggc                                       >  2:809277/1‑88 (MQ=255)
                                       ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  1:431545/90‑1 (MQ=255)
                                       ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  1:956956/90‑1 (MQ=255)
                                       ttCTCTACCATGCTGTTGATGGGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  2:1092122/90‑1 (MQ=255)
                                        tctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCg                                     <  1:78309/90‑1 (MQ=255)
                                                aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGGGAAAAAAAATGTCCGGGTCATGc                             >  1:1083073/1‑90 (MQ=255)
                                                   cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATc                          >  1:444899/1‑90 (MQ=255)
                                                     gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGACGGTTCGTTGAAAAAAAAAAATTCCGGGTTATTGATCAg                        >  1:1307449/1‑90 (MQ=255)
                                                     gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAAAAAGTCCGGGTCCTGCATCAg                        >  1:773809/1‑90 (MQ=255)
                                                        gATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCGGGTCGTCGAGAAACAAAATGTCCGGGTTATGCATCAGCgg                     >  1:1218416/1‑89 (MQ=255)
                                                             gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGGCCGGGTCATGCATCAGCGAACAgg                >  1:284908/1‑90 (MQ=255)
                                                             gCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCCTCAGCGAACAgg                >  2:1124204/1‑90 (MQ=255)
                                                               aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgggc              >  2:1057938/1‑88 (MQ=255)
                                                                      aTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca       <  2:1048863/90‑1 (MQ=255)
                                                                           cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGAGAACAGGCCAGCGCCAGCCGc  <  1:809277/90‑1 (MQ=255)
                                                                                  |                                                                                  
TCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822873‑823037

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