Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,953,2610TG38.1% 36.8 / 10.2 42intergenic (‑88/+111)ppiC/BW25113_RS25790peptidylprolyl isomerase PpiC/type II toxin‑antitoxin system PemK/MazF family toxin
Reads supporting (aligned to +/- strand):  ref base T (6/20);  new base G (16/0);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATT  >  NZ_CP009273/3953173‑3953349
                                                                                        |                                                                                        
tGATCGTGCCTTGAGTGAATAATTGTCAGGGGGCGATTTTAGCATGATCCACAGATGCTAACGTAATTGTCATTCAGGGAGTGGCTGGTa                                                                                         <  2:934002/90‑1 (MQ=255)
tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGTCTGGTa                                                                                         >  2:88836/1‑90 (MQ=255)
tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGTCTGGTa                                                                                         <  1:88836/90‑1 (MQ=255)
tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTa                                                                                         <  2:820915/90‑1 (MQ=255)
tGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTa                                                                                         <  2:147085/90‑1 (MQ=255)
     gTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACatat                                                                                    <  1:1015333/90‑1 (MQ=255)
       gCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACatatat                                                                                  <  1:123381/90‑1 (MQ=255)
           tGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCa                                                                              >  1:895366/1‑90 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaa                                                                             <  2:1036464/90‑1 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCaa                                                                             >  2:220849/1‑90 (MQ=255)
            gAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCaa                                                                             >  1:554721/1‑90 (MQ=255)
               tGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCaaaaa                                                                          >  2:64370/1‑89 (MQ=255)
                gAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGc                                                                         >  1:553098/1‑90 (MQ=255)
                      ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGAATATATGGCAAAAAACCggga                                                                   >  2:313302/1‑88 (MQ=255)
                      ttGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCggga                                                                   >  1:166226/1‑88 (MQ=255)
                        gTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAGCCCGGaaaa                                                                 >  1:1053606/1‑90 (MQ=255)
                          cAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATc                                                               <  2:323499/90‑1 (MQ=255)
                          cAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATTGCCAAAACCCGGaaaaac                                                               >  1:482309/1‑88 (MQ=255)
                                gcgATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCt                                                         <  2:699909/90‑1 (MQ=255)
                                       tAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAgg                                                  <  1:313302/90‑1 (MQ=255)
                                             gATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTc                                            >  1:360255/1‑90 (MQ=255)
                                                     gATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCGGAAAAACCGGGCTTCGCAGGTCCCTCAAgggggg                                    >  1:660477/1‑88 (MQ=255)
                                                      aTGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGgaga                                   <  2:287116/90‑1 (MQ=255)
                                                       tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCAAAAACCCGGAAAAAACGGGGTTCGCAGGGCCCTCCAGGGGGGaa                                  >  2:455079/1‑90 (MQ=255)
                                                       tGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGAAATATATGGCAAAAAACCGGGAAAACCGGGGTTCGCCGGGCCCTCCAGGGGAGaa                                  >  2:627630/1‑90 (MQ=255)
                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGGCCCTCAAGGGGAGaaa                                 >  1:405130/1‑90 (MQ=255)
                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGGCATATATGGCAAAAGCCCGGAAAAATCGGGGTTCGCAGGGGCCTCAAGGGGGGaaa                                 >  2:938903/1‑90 (MQ=255)
                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGAGAAATTTTGGAAAAGACCGGAAAATCCCGGCTTCGGATGGGCCTCAAGGGGGGaaa                                 >  2:380916/1‑90 (MQ=255)
                                                        gCTAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCCGAAAATCCGGGCTTCGCAGGGGCCTCAAAGGGAGaaa                                 >  1:78453/1‑90 (MQ=255)
                                                          tAACGTAATTGTCAGTCAGGGAGTGGCTGGGACATATATGGCAAAAACCCCGGAAAACCGGGGTTCGCAGGGCCCCCCAGGGGGGAAAAc                               >  2:900547/1‑90 (MQ=255)
                                                            aCGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACtt                             <  1:777987/90‑1 (MQ=255)
                                                              gTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATTCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAg                           <  1:64370/90‑1 (MQ=255)
                                                               tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAgg                          <  1:152469/90‑1 (MQ=255)
                                                               tAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAAATTAgg                          >  1:410652/1‑90 (MQ=255)
                                                                 aTTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGc                        <  1:455079/90‑1 (MQ=255)
                                                                    gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                     <  1:508511/90‑1 (MQ=255)
                                                                    gtcagtcaGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCctc                     <  1:900547/90‑1 (MQ=255)
                                                                        gtcagGGAGTGGCTGGGACATATATTGCAAAAACCCGGGAAAACCGGGCTTCGCAGGGCCCTCAAGGGGGGAAAACTTAGGGGCTCTATg                 >  2:508511/1‑90 (MQ=255)
                                                                         tcagGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATgg                <  2:410652/90‑1 (MQ=255)
                                                                          cagGGAGTGGCTGGTACATATATGGCAAAAGCACGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATggg               <  1:367555/90‑1 (MQ=255)
                                                                                      ggTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAt   >  1:321725/1‑90 (MQ=255)
                                                                                       gTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCAtt  <  1:917529/90‑1 (MQ=255)
                                                                                        |                                                                                        
TGATCGTGCCTTGAGTGAATAATTGTCAGGCGGCGATTTTAGCATGATCCGCAGATGCTAACGTAATTGTCAGTCAGGGAGTGGCTGGTACATATATGGCAAAAGCCCGGAAAATCCGGGCTTCGCAGGTCCCTCAAGGGGAGAAAACTTAGGGCCTCTATGGGCTACAAAGTCATT  >  NZ_CP009273/3953173‑3953349

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