Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,953,2610TG53.3% 5.7 / 16.0 15intergenic (‑88/+111)ppiC/BW25113_RS25790peptidylprolyl isomerase PpiC/type II toxin‑antitoxin system PemK/MazF family toxin
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (8/0);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAG  >  NZ_CP009273/3953173‑3953344
                                                                                        |                                                                                   
tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTa                                                                                    <  1:256179/90‑1 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaa                                                                        >  1:376049/1‑90 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaa                                                                        >  1:62158/1‑90 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaa                                                                        >  2:314507/1‑90 (MQ=255)
            gAGTGAATAATTGTAAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTAACGATATgggca                                                                        >  1:263386/1‑87 (MQ=255)
                  ataatTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGGAAAAAccc                                                                  >  1:227518/1‑90 (MQ=255)
                                             gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCGGGCTTCGCAGGGCCCTc                                       >  1:95313/1‑90 (MQ=255)
                                                       tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGAAAAACCGGGCTTCGCAGGTGCCTCAAGGGgagga                             >  1:386589/1‑88 (MQ=255)
                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGaaa                            <  2:62158/90‑1 (MQ=255)
                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCGGGCTTCGCAGGGCCCTCAAGGGGAGAAAAc                          >  1:190829/1‑90 (MQ=255)
                                                                       agtcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAACCCGGGCTTCGCAGGTCCCTCCAGGGGAGAAAACTTAGGGGCTCTAt             >  1:123715/1‑90 (MQ=255)
                                                                               aGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACa     <  1:163810/90‑1 (MQ=255)
                                                                                  ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAggg                    <  2:65282/72‑1 (MQ=255)
                                                                                  ggCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAg  >  2:276050/1‑90 (MQ=255)
                                                                                  ggCTGGGACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAATTAggg                    >  1:65282/1‑72 (MQ=255)
                                                                                        |                                                                                   
TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAG  >  NZ_CP009273/3953173‑3953344

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