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

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAT  >  NZ_CP009273/822886‑823052
                                                                             |                                                                                         
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                               >  1:119628/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAAAGCCGGAAg                                                                               >  2:168512/1‑90 (MQ=255)
                cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa                                                               <  2:387818/90‑1 (MQ=255)
                     acagtttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                          <  1:892043/86‑1 (MQ=255)
                      gCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTCGGCGGGAAAAAGaaa                                                         >  1:532791/1‑89 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                     <  2:908260/90‑1 (MQ=255)
                              ctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                                 <  2:140918/90‑1 (MQ=255)
                                 ccATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAt                                              <  2:166793/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGGGAAACAGAAGGTCCGGTTCAttc                                            >  2:331609/1‑88 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGGGAAAAAGAAAGTCCGGGTCAttc                                            >  2:681699/1‑88 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGc                                            >  2:62093/1‑90 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTTGTCGAGAAAAAAAAAGTCCGGTTCATGc                                            >  2:782977/1‑90 (MQ=255)
                                     gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCATGCAt                                          >  1:363647/1‑90 (MQ=255)
                                     gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGAGGGGTCGTCGGGGAAAAAAAAGGCCGGGTGATGCAt                                          >  1:934318/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                       <  2:227212/90‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCGGTTCGTCGTGAAAAAAAAAGTCCGGGTTATGCATTAg                                       >  2:791128/1‑90 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  <  1:119388/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                  <  1:159176/90‑1 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAgg                               >  2:666908/1‑90 (MQ=255)
                                                gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGCATCAGCGAACAgg                               >  1:54531/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGGTCGTCGGGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACAgggc                             >  1:53115/1‑88 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  1:344563/90‑1 (MQ=255)
                                                              cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc                 <  1:791128/90‑1 (MQ=255)
                                                                          ggggCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGGTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGGTTAAAACCt     >  1:827980/1‑90 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAAt  >  2:406586/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 < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: