Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG61.9% ‑5.7 / 18.9 21D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base G (13/0);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 4.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822875‑823037
                                                                                        |                                                                          
aTCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTc                                                                           <  1:99842/90‑1 (MQ=255)
       aaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCg                                                                    <  2:301325/90‑1 (MQ=255)
           gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                >  2:548302/1‑90 (MQ=255)
           gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                >  2:485805/1‑90 (MQ=255)
           gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                >  2:401429/1‑90 (MQ=255)
                       tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGTCGGGttgtcgg                                                    >  1:116348/1‑89 (MQ=255)
                       tCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGGGGGTTGTCga                                                    >  1:472494/1‑90 (MQ=255)
                            gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGGGaaaa                                               >  1:426149/1‑89 (MQ=255)
                             tGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACa                                              <  2:464714/90‑1 (MQ=255)
                                acagtttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                           <  2:444808/86‑1 (MQ=255)
                                    tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAGAAAGTc                                       >  2:180977/1‑90 (MQ=255)
                                              aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGGAACGCCGGAAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTTATGc                             >  2:321982/1‑90 (MQ=255)
                                                   gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                        <  2:504723/90‑1 (MQ=255)
                                                   gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAAAAAGTCCGGGTCATTCATCAg                        >  2:566478/1‑90 (MQ=255)
                                                        tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                   <  1:566478/90‑1 (MQ=255)
                                                        tgtgCACCCAAAATTCAAGGCGGGTGGGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                   <  1:517366/90‑1 (MQ=255)
                                                           gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGGAGGCGGGTCGGCGAGAAAAAGAAAGGCCCGGTCATGCATCAGCGAACAgg                >  1:254645/1‑90 (MQ=255)
                                                           gCAGCCAAAATTCACGGCGGGGGGGGGGGGGAAAGCCGGAAGGCGGGTCGTCGAGAAAAAGAAAGTCCGGGTCATGCATCAGCGAACaag                >  1:216387/1‑88 (MQ=255)
                                                           gCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGGCGGGTCGTCGGGAAAAAGAATGTCCGGGTTATGCATCAGCGAACAgg                >  2:260875/1‑90 (MQ=255)
                                                             aGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGGCGAGAAAAAGAAAGTTCGGGTCCTGCATCACCGGACAAgcc              >  1:208430/1‑90 (MQ=255)
                                                                         cggcggGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGc  <  1:485805/90‑1 (MQ=255)
                                                                                        |                                                                          
ATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGC  >  NZ_CP009273/822875‑823037

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: