Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG56.2% 26.5 / 64.8 48D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (3/18);  new base G (27/0);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 1.82e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-07
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT  >  NZ_CP009273/822886‑823052
                                                                             |                                                                                         
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)
                                   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                                                     >  2:280952/1‑76 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                     <  1:280952/76‑1 (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:95134/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  2:814233/90‑1 (MQ=255)
                                                                     tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa          <  2:466654/90‑1 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa         <  1:131978/90‑1 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa         <  2:1091452/90‑1 (MQ=255)
                                                                            gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa   >  2:121935/1‑90 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt  >  2:459287/1‑90 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt  >  2:216211/1‑90 (MQ=255)
                                                                             |                                                                                         
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT  >  NZ_CP009273/822886‑823052

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