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

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822886‑823044
                                                                     |                                                                                         
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAg                                                                       >  2:450366/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCgggaa                                                                       >  1:77407/1‑87 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg                                                                       >  1:762082/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCCGAAg                                                                       >  1:377527/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGt                                                                      >  2:1095496/1‑90 (MQ=255)
 ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAgg                                                                      >  1:525936/1‑89 (MQ=255)
       gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                <  2:1227113/90‑1 (MQ=255)
       gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTtc                                                                >  2:355753/1‑90 (MQ=255)
            tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGtcgtcgg                                                           >  2:1147859/1‑89 (MQ=255)
               ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaa                                                        <  1:1229219/90‑1 (MQ=255)
               ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGaa                                                        >  1:406568/1‑90 (MQ=255)
                         tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTTGTCGAGAAACaaaaaatc                                              >  2:836815/1‑86 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  1:831310/90‑1 (MQ=255)
                            ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg                                           <  2:137253/90‑1 (MQ=255)
                                  cATGCTGTTGATGTGCGGCCAAAATTCACGGCGGGGGGGGGGGGGCACGCCGGGAGTCGGGTTGTCGGGAAAAAGAAAGGCCGGGTCATg                                     >  1:1091452/1‑90 (MQ=255)
                                  cATGCTGTTGATGTGCAGCCAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  1:1095496/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg                                                                                       >  2:245632/1‑39 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg                                                                                       <  1:245632/39‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGGGGCGGGTCGTCGGGAAACAGAAAGTCCGGGTCATGc                                    >  1:972731/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGc                                    >  2:1360982/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGTTTGTCGGGAAAAAAAAAGTCCGGGTCCTGc                                    >  1:147849/1‑90 (MQ=255)
                                     gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCAt                                  >  2:618580/1‑90 (MQ=255)
                                      cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATc                                 >  1:549406/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  1:280952/76‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             >  2:280952/1‑76 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                               <  2:1390988/90‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                               <  2:472778/90‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAAGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAg                               >  2:1291418/1‑90 (MQ=255)
                                          tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGTCCACGCCCGAAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTTATGCATCAACg                             >  1:118354/1‑90 (MQ=255)
                                          tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAACg                             >  2:1140148/1‑90 (MQ=255)
                                          tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTTCTGCATCAGCg                             >  2:1117332/1‑90 (MQ=255)
                                          tGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCg                             >  2:60335/1‑90 (MQ=255)
                                           gATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGGAGGCGGGTTGTCGGGGAAAAAAAAGTGCGGGTCATGCATCAGCGa                            >  1:1198033/1‑90 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:1079133/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  2:1198033/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:76206/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:549755/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                          <  1:1147859/90‑1 (MQ=255)
                                               tgCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGGACAg                        >  2:606726/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATTCATCAGCGAACAgg                       >  1:466654/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAgg                       >  2:601354/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAAGCCGGGTCATGCATCAGCGAACAAgcc                     >  1:1220696/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAAAAgggc                     >  2:727782/1‑88 (MQ=255)
                                                     cAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGCGGACAGgccagg                  >  1:1056985/1‑89 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  2:814233/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc         <  2:95134/90‑1 (MQ=255)
                                                                     tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa  <  2:466654/90‑1 (MQ=255)
                                                                     |                                                                                         
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822886‑823044

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