Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273957,6210GT71.4% 63.4 / 38.2 42Q57H (CAG→CAT) rpsA30S ribosomal protein S1
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (15/15);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAA  >  NZ_CP009273/957483‑957750
                                                                                                                                          |                                                                                                                                 
gTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAg                                                                                                                                   >  1:142982/1‑139 (MQ=255)
 tCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGCCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAgg                                                                                                                                  >  1:845856/1‑139 (MQ=255)
  ccTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAtgg                                                                                                                                 >  2:438517/1‑139 (MQ=255)
  ccTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAtgg                                                                                                                                 >  1:708249/1‑139 (MQ=255)
       aaaGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAg                                                                                                                             <  1:831437/138‑1 (MQ=255)
         aGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTg                                                                                                                          <  2:567484/139‑1 (MQ=255)
          gAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTgg                                                                                                                         <  2:817280/139‑1 (MQ=255)
            aaTCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGaa                                                                                                                       <  2:122110/139‑1 (MQ=255)
             aTCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGaaa                                                                                                                      >  1:356971/1‑139 (MQ=255)
                                      ttcgtGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGt                                                                                             >  2:7420/1‑139 (MQ=255)
                                               ttgttgttgCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    <  1:332534/139‑1 (MQ=255)
                                               ttgttgttgCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    <  2:142982/139‑1 (MQ=255)
                                                tgttgttgCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGAc                                                                                   >  2:435990/1‑139 (MQ=255)
                                                  ttgttgCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACg                                                                                  <  1:7420/138‑1 (MQ=255)
                                                       gCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGa                                                                                                                              <  1:47137/89‑1 (MQ=255)
                                                       gCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGa                                                                                                                              >  2:47137/1‑89 (MQ=255)
                                                            cGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    <  1:139857/126‑1 (MQ=255)
                                                            cGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    >  2:139857/1‑126 (MQ=255)
                                                            cGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGAc                                                                       >  2:106801/1‑139 (MQ=255)
                                                                      gtagtaCTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTg                                                             <  2:505444/139‑1 (MQ=255)
                                                                              ggTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTctg                                                     >  1:226020/1‑139 (MQ=255)
                                                                               gTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    >  2:452007/1‑107 (MQ=255)
                                                                               gTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    <  1:452007/107‑1 (MQ=255)
                                                                               gTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACg                                                                                  <  1:2028/109‑1 (MQ=255)
                                                                               gTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACg                                                                                  >  2:2028/1‑109 (MQ=255)
                                                                                  gACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGt                                                 >  2:502615/1‑139 (MQ=255)
                                                                                    cGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    <  2:638995/102‑1 (MQ=255)
                                                                                    cGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGa                                                                                    >  1:638995/1‑102 (MQ=255)
                                                                                         gTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGt                                            >  1:847243/1‑137 (MQ=255)
                                                                                         gTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGt                                            >  1:618854/1‑137 (MQ=255)
                                                                                                aTCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCt                                   <  2:847243/139‑1 (MQ=255)
                                                                                                         cGCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                          <  1:438517/139‑1 (MQ=255)
                                                                                                          gCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         <  2:172543/139‑1 (MQ=255)
                                                                                                          gCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         <  1:762329/139‑1 (MQ=255)
                                                                                                          gCCATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         <  1:502615/139‑1 (MQ=255)
                                                                                                           ccATCCCGGCTGAGCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         >  2:285542/1‑138 (MQ=255)
                                                                                                                        gCAGTTCAAAAACGCCCATGGCGAGCTGGAAATCCAGGCAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGCAGCCTGGATCACg           >  2:192558/1‑139 (MQ=255)
                                                                                                                        gCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACg           >  2:573789/1‑139 (MQ=255)
                                                                                                                            ttCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCt                                   <  2:368594/111‑1 (MQ=255)
                                                                                                                            ttCAAAAACGCCCATGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCt                                   >  1:368594/1‑111 (MQ=255)
                                                                                                                             tCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         >  2:638023/1‑120 (MQ=255)
                                                                                                                             tCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACg                         <  1:638023/120‑1 (MQ=255)
                                                                                                                               aaaaaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGaaa    <  2:845856/139‑1 (MQ=255)
                                                                                                                                 aaaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                          <  2:55616/115‑1 (MQ=255)
                                                                                                                                 aaaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCAc                          >  1:55616/1‑115 (MQ=255)
                                                                                                                                 aaaCGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGaaaaa  <  1:729578/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                 
GTCCTTAAAAGAAATCGAAACCCGCCCGGGTTCTATCGTTCGTGGCGTTGTTGTTGCTATCGACAAAGACGTAGTACTGGTTGACGCTGGTCTGAAATCTGAGTCCGCCATCCCGGCTGAGCAGTTCAAAAACGCCCAGGGCGAGCTGGAAATCCAGGTAGGTGACGAAGTTGACGTTGCTCTGGACGCAGTAGAAGACGGCTTCGGTGAAACTCTGCTGTCCCGTGAGAAAGCTAAACGTCACGAAGCCTGGATCACGCTGGAAAAA  >  NZ_CP009273/957483‑957750

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑