Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,350,4580TG53.8% 9.8 / 17.2 13T304P (ACC→CCC) sapApeptide ABC transporter substrate‑binding protein SapA
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (7/0);  total (7/6)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTG  >  NZ_CP009273/1350369‑1350543
                                                                                           |                                                                                     
gtaataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGgt                                                                                       <  1:1314341/89‑1 (MQ=255)
     aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGttta                                                                                  >  1:1193489/1‑88 (MQ=255)
        tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTTaatg                                                                                >  2:563478/1‑88 (MQ=255)
          gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAAAGc                                                                              >  1:331035/1‑90 (MQ=255)
          gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGcc                                                                             <  1:257270/90‑1 (MQ=255)
          gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTAGCGGGGTTTAATGcc                                                                             >  2:1234429/1‑90 (MQ=255)
                   aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCg                                                                    <  1:295511/90‑1 (MQ=255)
                         tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTc                                                              <  2:180787/90‑1 (MQ=255)
                              attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTTAAATCCCGAAAAGCGGCGGTCATccc                                                         >  2:510681/1‑90 (MQ=255)
                                     gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGc                                                  <  1:563478/90‑1 (MQ=255)
                                                    tGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCAGGGCGCAGCGTTAAAAGCAAg                                   >  1:421383/1‑90 (MQ=255)
                                                       cGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcgc                                <  1:494629/90‑1 (MQ=255)
                                                                                ggTTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGGCGTCACGCAAAAAGGATaac        >  2:560904/1‑88 (MQ=255)
                                                                                       gCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggctg  <  1:952166/90‑1 (MQ=255)
                                                                                           |                                                                                     
ATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTG  >  NZ_CP009273/1350369‑1350543

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