Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG63.2% 6.0 / 41.8 38D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/14);  new base G (24/0);  total (24/14)
Fisher's exact test for biased strand distribution p-value = 1.03e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.10e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT  >  NZ_CP009273/822886‑823052
                                                                             |                                                                                          
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGGAg                                                                               >  2:299759/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGTGGGGGCTACGCCGGGAg                                                                               >  1:451737/1‑90 (MQ=255)
   tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAACCGGAAGGCg                                                                            >  2:696750/1‑90 (MQ=255)
       gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTtc                                                                        <  1:125850/90‑1 (MQ=255)
          cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGGCGGgtcgtc                                                                     >  1:508803/1‑90 (MQ=255)
            tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCga                                                                   <  1:53615/90‑1 (MQ=255)
               ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGGCCGGAAGCGGGTTGTCgagga                                                                >  1:1042420/1‑88 (MQ=255)
                 gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGaaaa                                                              >  2:53615/1‑89 (MQ=255)
                    aCGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGa                                                           <  2:1150490/90‑1 (MQ=255)
                         tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAACAGAAAGTc                                                      >  2:144076/1‑90 (MQ=255)
                         tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGGCGGGTCGTCGAGAAAAAGaaaggc                                                      >  1:990402/1‑86 (MQ=255)
                                aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCa                                               >  2:63797/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                            <  2:269101/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTCATGc                                            >  2:1109342/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTTGTCGGGAAAAAGAAAGTCCGGGTCATGc                                            >  2:781360/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGGCCGGGTCATGc                                            >  1:1163/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGGGAAAAAGAATGTCCGGGTCCTgg                                            >  1:750442/1‑89 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAATGTCCGGGTCATGc                                            >  2:83312/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGc                                            >  1:205245/1‑90 (MQ=255)
                                     gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCGGGTCATGCAt                                          >  2:228081/1‑90 (MQ=255)
                                      cTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATc                                         <  2:897853/90‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                       <  1:543319/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  <  1:781360/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  <  2:750442/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGGTCGTCGCGAAAAAGAAAGTCCGGGTTCTGCATCAGCGaaa                                  >  2:822987/1‑89 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGTTCATGCATCAGCGAACAgg                               >  1:284489/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgg                               >  1:320062/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACCCCCCAAGGCGGTTCGCTGAGAAACAAAATATCCGGGTCATGGATCAGCGCACAgg                               >  1:1140847/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGTAGGCGGGTCGGCGTAAAACAGAAAAGGCGGGTTATGCAGCAGCGGACaag                               >  1:897853/1‑88 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGAAGTCGGTTTGTCGAGAAACAGAATGTCCGGGTTCTTCATCAGCGCACAgg                               >  2:911656/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                             <  1:918416/90‑1 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTCGTCGGGAAAAAGAAAGTGCGGGGCATGCCTCAGCGGAAAGgcc                             >  1:931361/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAGAATGTCCGGGTTATGCATCAGCGAACAGgcc                             >  2:1042041/1‑90 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  1:1042987/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  1:940332/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  1:83312/90‑1 (MQ=255)
                                                               ggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCt                <  2:725932/90‑1 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa         <  1:1109342/90‑1 (MQ=255)
                                                                          ggggCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCCCTGTTTAAAAccc     >  1:686864/1‑89 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  >  2:721353/1‑58 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  <  1:721353/58‑1 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa          >  2:3784/1‑82 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa          <  1:3784/82‑1 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt  <  2:931937/90‑1 (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 < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: