Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,9630TG55.0% 7.6 / 26.7 20D492A (GAC→GCC) ybhFATP‑binding cassette domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (11/0);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 7.14e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA  >  NZ_CP009273/822886‑823051
                                                                             |                                                                                        
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAg                                                                              >  1:159143/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAg                                                                              >  2:56527/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAg                                                                              >  2:372697/1‑90 (MQ=255)
               ccGTGACGCCTTTCTCTACCATGCTGTTGATGTGCGGCCAAAATTCACGGCGGGTGGGGGGGGCAACACCGGAAGACGGGttgtcgcgga                                                               >  2:427241/1‑86 (MQ=255)
                 gTGACGCTTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGAAAGCGGGTCGTCGGGaaaa                                                             >  2:719973/1‑89 (MQ=255)
                     cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                         <  2:154764/90‑1 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                                    <  2:190471/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGc                                           <  1:611060/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGTGGGGGGCGAAGCCGGGAGGGGGGTCCTCGAGGAAAAGAAGGTCCGGGTCATGc                                           >  1:677634/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAg                                      >  1:676074/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAGGCGGGTTGTCGGGAAAAAAAAAGTCCGGGTCATGCATCAg                                      >  1:142733/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCCGAAGGCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCCTGCATccg                                      >  1:317018/1‑88 (MQ=255)
                                          tGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCCGGAGGCGGGTCGTCGAGAAAAAGAAAGTGCGGGTGATGCATCAGCg                                    >  1:566296/1‑90 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                 <  1:668888/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                 <  2:26940/90‑1 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                                 <  1:372697/90‑1 (MQ=255)
                                               tgCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAAAGGGGTCGTCGGGGAACAAAAAGTGCGGTTTATTAATCAACGGACAg                               >  1:602307/1‑90 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa        <  2:159143/90‑1 (MQ=255)
                                                                           ggTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTa   <  2:15567/90‑1 (MQ=255)
                                                                            gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa  >  1:478837/1‑90 (MQ=255)
                                                                             |                                                                                        
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA  >  NZ_CP009273/822886‑823051

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