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

GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822886‑823045
                                                                             |                                                                                  
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAg                                                                        >  1:115157/1‑90 (MQ=255)
gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGGAg                                                                        >  2:80885/1‑90 (MQ=255)
   tggtgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCg                                                                     >  2:132972/1‑90 (MQ=255)
          cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAGTCGGgtcgtc                                                              >  2:387880/1‑90 (MQ=255)
                     cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa                                                   <  1:329883/90‑1 (MQ=255)
                        cTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGt                                                <  2:282541/90‑1 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  1:387880/90‑1 (MQ=255)
                          ttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  1:475880/90‑1 (MQ=255)
                                aCCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCa                                        <  2:163031/90‑1 (MQ=255)
                                   aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGGTCGTCGAGAAACAGAATGTCCCGTTCATgg                                     >  2:88080/1‑89 (MQ=255)
                                     gCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGGAt                                   >  2:354970/1‑90 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg                                <  1:132972/90‑1 (MQ=255)
                                        gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGTCGGGTTGTCGAGAAAAAAAAAGTCCGGGTCATGCATCAg                                >  1:176912/1‑90 (MQ=255)
                                             tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc                           <  2:102511/90‑1 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGAAGGCGGGTCGTCGAGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGgcc                      >  1:317525/1‑90 (MQ=255)
                                                  aGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGgcc                      >  2:196583/1‑90 (MQ=255)
                                                                   ggggAGGGGGGCAACGCCGGAAGGCGGTTCGTCGGGAAACAGAAAGTCCGGGTCATGCATCAGCGAACAGGCCCGCGCCCGCCCCCGttt     >  2:434160/4‑90 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  1:80865/90‑1 (MQ=255)
                                                                      gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa  <  2:458241/90‑1 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              <  2:198545/39‑1 (MQ=255)
                                                                             tCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc                                              >  1:198545/1‑39 (MQ=255)
                                                                             |                                                                                  
GGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAA  >  NZ_CP009273/822886‑823045

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