Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9550TG61.1% 13.2 / 19.4 18T495P (ACC→CCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (11/0);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 3.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGG  >  NZ_CP009273/822870‑823023
                                                                                     |                                                                    
gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg                                                                  <  2:524352/90‑1 (MQ=255)
      tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCa                                                            <  1:772972/90‑1 (MQ=255)
                 ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAgg                                                 >  1:851783/1‑89 (MQ=255)
                      tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCgggt                                            >  2:649666/1‑88 (MQ=255)
                        aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCCGGAGGCGGTtcg                                          >  1:590380/1‑90 (MQ=255)
                          cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTtcgtc                                        >  2:364294/1‑90 (MQ=255)
                                cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa                                  <  2:50578/90‑1 (MQ=255)
                                     cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                             <  1:389456/90‑1 (MQ=255)
                                         tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGGGAAAAAAAAAGTc                         >  1:45412/1‑90 (MQ=255)
                                          gtctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                        <  2:699753/89‑1 (MQ=255)
                                            ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg                      <  2:471667/90‑1 (MQ=255)
                                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGGGAAACAGAAAGTCCGGGTCATGc               >  2:242609/1‑90 (MQ=255)
                                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGc               >  2:235176/1‑90 (MQ=255)
                                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg          <  2:183729/90‑1 (MQ=255)
                                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGAAGGTGGGTTGTCGGGAAAAAGAAAGTCCGGGTTATGCCTCaa          >  2:824447/1‑89 (MQ=255)
                                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCCTGCATCAg          >  1:183729/1‑90 (MQ=255)
                                                             tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTGCGGGTCATGGATCAGCGAAc     >  2:853240/1‑90 (MQ=255)
                                                                gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAAGCATCAGCGAACAgg  >  2:264829/1‑90 (MQ=255)
                                                                                     |                                                                    
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGG  >  NZ_CP009273/822870‑823023

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