Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG48.1% 53.7 / 33.8 54D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/27);  new base G (26/0);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 2.88e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822881‑823045
                                                                                  |                                                                                  
aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                             <  1:1408504/90‑1 (MQ=255)
     gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                        >  1:1035017/1‑90 (MQ=255)
     gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                        >  2:2070491/1‑90 (MQ=255)
     gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                        >  1:489611/1‑90 (MQ=255)
     gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCGACGCCggtgc                                                                        >  1:2755726/1‑87 (MQ=255)
     gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAg                                                                        >  2:663869/1‑90 (MQ=255)
       gtggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAggc                                                                      >  1:652443/1‑88 (MQ=255)
        tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCCGAAGGCg                                                                     >  1:204438/1‑90 (MQ=255)
               cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAGGCGGGtcggc                                                              >  2:1532589/1‑88 (MQ=255)
                aTCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtcgtcg                                                             <  1:1454599/90‑1 (MQ=255)
                    ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaa                                                         <  1:326852/90‑1 (MQ=255)
                    ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGtcggcgggga                                                         >  1:253878/1‑86 (MQ=255)
                      gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGTCGGGTCTGCGAGaaaa                                                       >  1:2871076/1‑89 (MQ=255)
                      gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAGAAAc                                                       >  2:326852/1‑90 (MQ=255)
                        gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAg                                                     >  2:434743/1‑90 (MQ=255)
                              tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGTTTGTCGAGAAAAAAAATGTc                                               >  1:836709/1‑90 (MQ=255)
                              tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGCCGGGTCGTTGGGAAAAAGAATGTc                                               >  1:2539865/1‑90 (MQ=255)
                               ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  1:1532589/90‑1 (MQ=255)
                               ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  2:2093453/90‑1 (MQ=255)
                               ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  2:2103915/90‑1 (MQ=255)
                                   ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                          <  1:2637725/90‑1 (MQ=255)
                                     aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCa                                        >  1:897351/1‑90 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     <  2:489611/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     <  1:887104/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGTCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCCTGc                                     >  1:59957/1‑90 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                     >  2:446716/1‑90 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGGCGGGTCGTCGGGGAAAAGAAAGTCCCGGTCATGc                                     >  1:2103915/1‑90 (MQ=255)
                                          gCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCAt                                   <  1:1958613/90‑1 (MQ=255)
                                          gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCAt                                   >  2:463716/1‑90 (MQ=255)
                                          gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCAt                                   >  2:1655419/1‑90 (MQ=255)
                                             gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                <  1:2547054/90‑1 (MQ=255)
                                             gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGGGGGTCGGCGGGAAAAAAAAAGTCCGGGTCATGCATCAg                                >  1:2154453/1‑90 (MQ=255)
                                                gATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGGTCGGCGGGAAAAAGAAAGTGCGGGTGATGCATCAGCGa                             >  1:1713649/1‑90 (MQ=255)
                                                  tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  1:2070491/90‑1 (MQ=255)
                                                  tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  1:1005818/90‑1 (MQ=255)
                                                  tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  2:2016837/90‑1 (MQ=255)
                                                     gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTGGGGTCGTCGGGAAAAAAAAAGGCCGGGTCATGCATCAGCGCAAAgg                        >  1:2437964/1‑90 (MQ=255)
                                                     gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAATCCGGGTCATGCATCAGCGAACAgg                        >  2:1548839/1‑90 (MQ=255)
                                                       aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGGCGGGAAAAAGAAAGTCCGGGTTATGCATCAGCGGACAgggc                      >  2:2795500/1‑88 (MQ=255)
                                                            aaaTCCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgc                 <  2:1013213/90‑1 (MQ=255)
                                                                  aCGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAACAGAAAGTCCGGGTCATGCATCAGCGGACAGGCCAGCGCCAGccc           >  1:2408993/1‑89 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAGCGCCGGATGTCGGTTCGTCGAAAACCGGAATGTACGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  1:2064939/90‑1 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:1179490/90‑1 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  1:907942/90‑1 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:2019343/90‑1 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  1:462653/90‑1 (MQ=255)
                                                                   cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc          <  2:936592/90‑1 (MQ=255)
                                                                        ggTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGttt     <  2:514134/90‑1 (MQ=255)
                                                                          tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa   <  1:1548839/90‑1 (MQ=255)
                                                                           gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  2:1022544/90‑1 (MQ=255)
                                                                           gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  2:2365735/90‑1 (MQ=255)
                                                                           gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  2:2408993/90‑1 (MQ=255)
                                                                               gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                       <  1:2460571/49‑1 (MQ=255)
                                                                               gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                       >  2:2460571/1‑49 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  1:1251686/58‑1 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  2:2627623/58‑1 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           >  2:1251686/1‑58 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           >  2:1205489/1‑58 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  1:1205489/58‑1 (MQ=255)
                                                                                  tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           >  1:2627623/1‑58 (MQ=255)
                                                                                  |                                                                                  
AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822881‑823045

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