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

CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822871‑823037
                                                                                    |                                                                                  
cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggg                                                                               >  2:39604/1‑90 (MQ=255)
    aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGgggggggggggg                                                                           >  2:104960/1‑88 (MQ=255)
    aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAggggggc                                                                           >  1:813372/1‑88 (MQ=255)
         tAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGc                                                                      <  2:324606/90‑1 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                >  2:358795/1‑90 (MQ=255)
               gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAg                                                                >  1:160130/1‑90 (MQ=255)
                      gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGtc                                                         >  2:209108/1‑90 (MQ=255)
                           tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGttgtcgg                                                    >  1:246220/1‑89 (MQ=255)
                                gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGGCGGTTTGTCGAGaaaa                                               >  1:512839/1‑89 (MQ=255)
                                     gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAt                                          >  1:310949/1‑90 (MQ=255)
                                      ccTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGAAGCCCGGAGGCGGGTCGTCGAGAAACAGaaag                                         >  2:698826/1‑88 (MQ=255)
                                        tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGCGAAACAGAAAGTc                                       >  1:661896/1‑90 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAGCGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  2:1089292/90‑1 (MQ=255)
                                         ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  1:104960/90‑1 (MQ=255)
                                         acagctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                      <  2:635770/86‑1 (MQ=255)
                                           ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg                                    <  1:594950/90‑1 (MQ=255)
                                             ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                  <  2:677528/90‑1 (MQ=255)
                                                ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                               <  1:135592/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                             <  1:21970/90‑1 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTGATGc                             >  2:554287/1‑90 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCGGAAGTCGGGTCGTCGGGAAAAAGAATGTCCGGGGCATgg                             >  2:1062444/1‑89 (MQ=255)
                                                  aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGc                             >  2:280796/1‑90 (MQ=255)
                                                     cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATc                          <  1:280796/90‑1 (MQ=255)
                                                       gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                        <  1:602910/90‑1 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGGCCGGGTCATGCATCAGCGGACaag                >  2:934680/1‑88 (MQ=255)
                                                               gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgg                >  1:132913/1‑90 (MQ=255)
                                                                 aGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc              <  1:1114260/90‑1 (MQ=255)
                                                                            aCGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGGCGAGAAAAAAAAAGGCCGGGTCGTGCAACCGCGGACAGGCCAGCGCccgccg   >  2:677713/1‑90 (MQ=255)
                                                                             cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc  <  1:554287/90‑1 (MQ=255)
                                                                                    |                                                                                  
CTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822871‑823037

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