Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9550TG60.0% 21.2 / 28.8 30T495P (ACC→CCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/12);  new base G (18/0);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.16e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-07
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCC  >  NZ_CP009273/822868‑823035
                                                                                             |                                                                                
ccGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGa                                                                                <  1:1125920/90‑1 (MQ=255)
   cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑Ggggggggg                                                                             >  2:205831/1‑90 (MQ=255)
       aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑Ggggggggggggg                                                                         >  1:1568198/1‑88 (MQ=255)
        tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCa                                                                        <  1:1048690/90‑1 (MQ=255)
          cATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAAc                                                                      <  1:937580/90‑1 (MQ=255)
                 tgggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAcccccgga                                                               >  1:470546/1‑86 (MQ=255)
                  gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGGACGCCGGGAg                                                              >  2:937580/1‑90 (MQ=255)
                  gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGGACGCCGGGAg                                                              >  2:776819/1‑90 (MQ=255)
                      ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGGAGGCgg                                                          >  1:610638/1‑90 (MQ=255)
                         gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                       <  1:742796/90‑1 (MQ=255)
                         gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                       <  1:79760/90‑1 (MQ=255)
                            cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGAAGGCGGgtcgtc                                                    >  1:891919/1‑90 (MQ=255)
                                   gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGaaaa                                             >  1:1143546/1‑89 (MQ=255)
                                            ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                    <  2:35153/90‑1 (MQ=255)
                                                 tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                     <  1:476915/90‑1 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGTGGGGGGGGACGCCGGAGGGCGGGTTGGCGGGAAAAAAAAAGTCCCGGTCAttc                           >  2:1276375/1‑88 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGAACGCCCGAAGGCGGGTTGTCCAGAAAAAGAAAGTGCGGGTGATgg                           >  2:235850/1‑89 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGTGGGGTCGTCGAGAAAAAGAATGTGCGGGGCATGc                           >  2:1027209/1‑90 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGGCGGGTTGTCGAGGAAAAAAAAATCCGGGTTATGc                           >  1:148693/1‑90 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGGGACCCCGGGAAGCGGGTCGGCGAGGAAAAGAAAGTTCGGGTCATGc                           >  2:60360/1‑90 (MQ=255)
                                                     aTGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGc                           >  1:32911/1‑90 (MQ=255)
                                                          gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                      <  1:60360/90‑1 (MQ=255)
                                                          gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAATGTCCGGTTCATGCATccg                      >  2:1490583/1‑88 (MQ=255)
                                                          gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTCATGCCTCAg                      >  1:683374/1‑90 (MQ=255)
                                                          gTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCCTGCATCAg                      >  1:874603/1‑90 (MQ=255)
                                                               tgtgCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                 <  1:1597129/90‑1 (MQ=255)
                                                               tgtgCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                 <  2:514497/90‑1 (MQ=255)
                                                                    aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGGCACGCCCGAAGGCGGGTCGTCGGGAAAAAAAAAGGCCGGGGCATGCATCAGCGAACAGgcc            >  2:436753/1‑90 (MQ=255)
                                                                    aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGGGGGGGGCAAAGCCGGAAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAAgcc            >  1:166591/1‑90 (MQ=255)
                                                                    aGCCAAAATTCACGGC‑‑‑‑‑‑GGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAAAAATTCCGGGTTATGCATCAGCGAAAAGgcc            >  1:402090/1‑90 (MQ=255)
                                                                              cACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc  <  1:1276375/90‑1 (MQ=255)
                                                                              cACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc  <  1:235850/90‑1 (MQ=255)
                                                                                             |                                                                                
CCGCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGC‑‑‑‑‑‑GGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCC  >  NZ_CP009273/822868‑823035

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