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

GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822870‑823044
                                                                                     |                                                                                         
gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg                                                                                       <  1:516018/90‑1 (MQ=255)
 cTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAggg                                                                                      >  2:137098/1‑90 (MQ=255)
                gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAg                                                                       >  2:529971/1‑90 (MQ=255)
                gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGTCAACGCCGGAAg                                                                       >  1:430325/1‑90 (MQ=255)
                     gtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGGACGCCGGAAGGCggg                                                                  >  1:545349/1‑89 (MQ=255)
                                   gACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAg                                                    <  1:618794/90‑1 (MQ=255)
                                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  1:137098/90‑1 (MQ=255)
                                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                             <  2:430325/90‑1 (MQ=255)
                                              agaCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGtt                                         <  1:228708/88‑1 (MQ=255)
                                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                    <  2:545349/90‑1 (MQ=255)
                                                          tGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGAAACAAAAAGTTCGGGTCATGCATCAGct                             >  2:228708/1‑89 (MQ=255)
                                                                gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCGGGGCGTCGGGAAAAAGAAAGTGCGGGTGATTCATCCGCGGAAAgg                       >  1:452189/1‑90 (MQ=255)
                                                                gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAAAAgg                       >  2:118240/1‑90 (MQ=255)
                                                                 aggCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgc                      <  1:207067/88‑1 (MQ=255)
                                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGTCGGTTCGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAGCGAACAgggc                     >  2:534489/1‑88 (MQ=255)
                                                                            cACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGcc           <  1:534489/90‑1 (MQ=255)
                                                                             aCGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAAAAAAAAGTCCGGTTCATGCATCAGCGAACAGGCCAGcgcccgccc          >  2:207067/1‑89 (MQ=255)
                                                                                     tGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaa  <  2:584357/90‑1 (MQ=255)
                                                                                     |                                                                                         
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAA  >  NZ_CP009273/822870‑823044

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