Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273822,4140TC58.3% 5.7 / 10.7 12K99E (AAA→GAA) ybhSABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (0/7);  total (4/8)
Fisher's exact test for biased strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTG  >  NZ_CP009273/822347‑822495
                                                                    |                                                                                  
ggggtgGCGTGGGCGCCCCCCTTCTTTCCCGAAAAACCCCCCGGAATACCCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATc                                                             <  2:313697/88‑1 (MQ=255)
gcggtgGCGTGGGCGCGCTCCTTCTGTCCCGCAAAACCCCCGGAAATACCCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATc                                                             <  2:510657/90‑1 (MQ=255)
gcggtgGCGTGGGCCCGCCCCTTCTTTCCCGCAAAACCCCCCGAAATACCCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATc                                                             <  2:733691/90‑1 (MQ=255)
gcggggGCGTGGGCGCGCTCCTTCTTTCCCGAAAAACCCCCCGAAATACCCCCCAGCCCGCGATTTTCCCCCGCCTGCATTTTGGCGATc                                                             <  2:64417/90‑1 (MQ=255)
               cgcTCCATCTGTCCCGCAAAACCCCCCGAAATACCCCCCAGACCGCGAATTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTa                                              <  1:22032/90‑1 (MQ=255)
                     aTCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTg                                        >  1:595264/1‑90 (MQ=255)
                         gTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtgg                                    >  1:535131/1‑90 (MQ=255)
                           cccGAAAAACCCACCGAAAAAACCCCCAGCCCGCAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtg                                  <  2:535131/89‑1 (MQ=255)
                                 aaaaCCCCCGGAAATACCCCCCAGCCCGCAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTc                             <  2:373628/90‑1 (MQ=255)
                                      ccccGGAAAAACCCCCCGCCCCGAAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                        <  2:279214/88‑1 (MQ=255)
                                                        gACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATgg      >  1:953737/1‑90 (MQ=255)
                                                            gcgAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGtgtg  >  1:774179/1‑90 (MQ=255)
                                                                    |                                                                                  
GCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAACCACCAGACCGCGAATTTTCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTG  >  NZ_CP009273/822347‑822495

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