Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG47.4% 9.5 / 29.0 38D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/19);  new base G (18/0);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTA  >  NZ_CP009273/822881‑823050
                                                                                  |                                                                                       
aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc                                                                                  <  1:1250790/90‑1 (MQ=255)
         ggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGGAGTggg                                                                         >  2:677993/1‑90 (MQ=255)
               cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACCCCGGAAGTCGGgtcgtc                                                                   >  2:1151798/1‑90 (MQ=255)
                  cACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGACGGAAGTCGGTTCGTCgag                                                                <  2:54477/90‑1 (MQ=255)
                  cACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGACGGAAGTCGGTTCGTCgag                                                                <  2:210423/90‑1 (MQ=255)
                      gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAACCCGGAAGTCGGGTCGTCGAGAAAc                                                            >  1:583871/1‑90 (MQ=255)
                      gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACCCCCGAAGGCGGGTCGTCGGGaaaa                                                            >  2:613300/1‑89 (MQ=255)
                       tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                                           <  1:1266528/90‑1 (MQ=255)
                       tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                                           <  1:374700/90‑1 (MQ=255)
                        gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCGGGTTGTCGGGAAAAAg                                                          >  1:1811641/1‑90 (MQ=255)
                        gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGTCGGGTTGTCGAGAAAAAg                                                          >  2:1189456/1‑90 (MQ=255)
                              agtctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTc                                                    <  2:131363/88‑1 (MQ=255)
                               ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                   <  2:1867045/90‑1 (MQ=255)
                               ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                   <  1:685978/90‑1 (MQ=255)
                                    tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAACAGAAAGTCCGGGTc                                              >  2:506954/1‑90 (MQ=255)
                                       cATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGCAAGCCGGGAGGCGGGTCGGCGGGAAAAAGAAAGACCGGGTCGTg                                           >  1:291613/1‑90 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGTGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                          <  2:319005/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                          <  1:1903764/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                          <  2:233861/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTAGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                          <  1:832850/90‑1 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGTGGGGCAACGCCGGGAGTCGGGTGGTCGAGAAAAAGAAAGTCCGGGTCAttc                                          >  1:319005/1‑88 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTCGTCGAGAAAAAAAAAATCCGGGTCATGc                                          >  2:685978/1‑90 (MQ=255)
                                        aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGGGGGGTCGTCGGGAAAAAAAAAGTCCGGGTCAttc                                          >  2:1136371/1‑88 (MQ=255)
                                         tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                      <  1:413663/61‑1 (MQ=255)
                                         tGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTtc                                                                      >  2:413663/1‑61 (MQ=255)
                                          gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAGAAAAAGAAAGGCGGGGgcatgcat                                        >  1:1602637/1‑90 (MQ=255)
                                           cTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGTCGGGAAAAAGAATGTCCGGGTGATTCATc                                       >  1:233861/1‑90 (MQ=255)
                                             gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                     <  1:174716/90‑1 (MQ=255)
                                             gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAAAAAAAAGGCCGGGTCATGCCTCAg                                     >  1:902902/1‑90 (MQ=255)
                                               tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCg                                   >  1:1262057/1‑90 (MQ=255)
                                                     gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGACGGTTCGTCGGGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAgg                             >  1:1621513/1‑90 (MQ=255)
                                                       aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGGGAAAAAGAAAGGCCGGGTCATGCATCAGCGAACAGgcc                           >  2:174716/1‑90 (MQ=255)
                                                              aTTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgccagcgcca                    <  1:506954/90‑1 (MQ=255)
                                                                    ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt              <  2:20162/90‑1 (MQ=255)
                                                                         gTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTa         <  1:1189455/90‑1 (MQ=255)
                                                                               gggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCt   >  2:1810114/1‑90 (MQ=255)
                                                                                ggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTa  <  2:1811642/90‑1 (MQ=255)
                                                                                ggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTa  <  1:156833/90‑1 (MQ=255)
                                                                                  |                                                                                       
AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTA  >  NZ_CP009273/822881‑823050

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: