Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273117,9940AT53.3% 5.2 / 17.0 15L15L (CTT→CTAaroParomatic amino acid transporter AroP
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base T (8/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATG  >  NZ_CP009273/117905‑118081
                                                                                         |                                                                                       
ccTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGAGATAAGCTGAATATGGCGGTTTTTa                                                                                         <  1:260677/90‑1 (MQ=255)
               aCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAAGCCGCGCTTTTgc                                                                          >  2:163991/1‑90 (MQ=255)
                  gAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTtgctgc                                                                       >  2:214295/1‑90 (MQ=255)
                    ggCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTc                                                                     >  2:219103/1‑90 (MQ=255)
                    ggCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGGTTTTTAAGCCGCGCCTTTgctggct                                                                     >  1:175312/1‑87 (MQ=255)
                       gcTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGcc                                                                  >  2:10578/1‑90 (MQ=255)
                                         ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTccct                                                >  2:22568/1‑88 (MQ=255)
                                         ggTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCcat                                                <  1:163991/90‑1 (MQ=255)
                                                                  aTAAGCTGAATATGGCGGGTTTTTAAGCCGCGCTTTAGCCGCTCCCCGGGCTGTTGGCCTTCCCTCATTAAAACTCGTGCGGGGggtggt                       >  1:247536/1‑88 (MQ=255)
                                                                      gCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGtttttt                   <  1:109021/90‑1 (MQ=255)
                                                                      gCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGtttttt                   <  2:355010/90‑1 (MQ=255)
                                                                              tGGCGGTTTTTTAGGCCGCGCTTTTGCTGCTCGCCGGGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTTATCTAc           >  1:5149/1‑90 (MQ=255)
                                                                                gCGGTTTTTTAGGCCGCGCTTTAGCTGCTCGCCGGGCTGTTGGCCTTCCATCATGAAAACTCGTGCGGGGGTTGGTTTTTTTTTCTTCGc         >  1:317164/1‑90 (MQ=255)
                                                                                 cGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCa        <  1:185666/90‑1 (MQ=255)
                                                                                 cGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAg       <  1:361445/90‑1 (MQ=255)
                                                                                       ttAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATg  <  1:187699/90‑1 (MQ=255)
                                                                                         |                                                                                       
CCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGCTCGCCGTGCTGTTGACCTTCCATCATGAAACCTCGTGCGGTGGTTGTTTTTTTGATCTACGCAGTGATG  >  NZ_CP009273/117905‑118081

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