Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,648,6370AC56.2% 8.5 / 14.6 16V209G (GTG→GGG) yhiDMgtC/SapB family protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (0/9);  total (3/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGACC  >  NZ_CP009273/3648556‑3648714
                                                                                 |                                                                             
atTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGcc                                                                       >  1:268766/1‑90 (MQ=255)
                gAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCCCACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATc                                                       <  2:341484/90‑1 (MQ=255)
                      cAGGCTTACGGTGAGTCTGGCTCCGCTCCCACACAGATTAGCTAATTGAAACGCCTTCCCCCCCTGCCATACCTTTTAATAATCGCAACa                                                 <  1:331622/90‑1 (MQ=255)
                      cAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACa                                                 <  2:294854/90‑1 (MQ=255)
                          cTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGtc                                             <  2:76051/90‑1 (MQ=255)
                                               cTGCCACACAGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGa                       <  1:96248/90‑1 (MQ=255)
                                                   cacacaGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGc                    <  1:413847/90‑1 (MQ=255)
                                                      acaGATTAGCTAATTGAACCCCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTgg                 <  1:213760/90‑1 (MQ=255)
                                                      acaGATTAGCTAATTGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTgg                 <  2:259185/90‑1 (MQ=255)
                                                           ttAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATc            >  2:104112/1‑90 (MQ=255)
                                                             aGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCaa          <  1:104112/90‑1 (MQ=255)
                                                                   ttGAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCAt                         <  1:436185/69‑1 (MQ=255)
                                                                   ttGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCAt                         >  2:436185/1‑69 (MQ=255)
                                                                    tGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAc   <  1:355678/90‑1 (MQ=255)
                                                                     gAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAcc  <  1:148634/90‑1 (MQ=255)
                                                                     gAAACGCCTTTCCCCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGAcc  <  2:369017/90‑1 (MQ=255)
                                                                                 |                                                                             
ATTGATTATCTGGTGTGAATTTCAGGCTTACGGTGAGTCTGGCTACGCTGCCACACAGATTAGCTAATTGAAACGCCTTTCACCCCTGCCATACCTTTTAATAATCGCAACAGGTCTTCTATCGACGTTGTGGCATGAAGCTGGATATCAATGGCGACC  >  NZ_CP009273/3648556‑3648714

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